R語言 t分布(不同自由度）

了解r語言幾個函數：dt，pt，qt，rt分別與dnorm，rnorm，pnorm，qnorm和rnorm對應 > * dt() 的返回值是正態分布概率密度函數(density)
> * pt()返回值是正態分布的分布函數(probability)
> * 函數qt()的返回值是給定概率p后的下百分位數(quantitle)
> * rt()的返回值是n個正態分布隨機數構成的向量

x <- seq(-4, 4, length=200)
df <- c(3, 8, 16, 61)
require(plyr)

## Loading required package: plyr

get.pt <- function(x, df) {
    prob <- dt(x, df)
    dd <- data.frame(x=x, df=factor(df), prob=prob)
    return(dd)
}
pt.df <- mdply(data.frame(x= rep(x, length(df)), df=rep(df, each=length(x))), get.pt)
require(ggplot2)

## Loading required package: ggplot2

ggplot(pt.df, aes(x, prob))+geom_line(aes(group=df, color=df), lwd=1)+geom_line(data=data.frame(x=x, prob=dnorm(x)), alpha=0.3, lwd=3, color="gray")