t分布のグラフ

例1.　自由度10のt分布

DF<-function(x)dt(x,10)
plot(DF,-4,4)

例2.　t二乗分布の塗りつぶし指定

　以下は塗りつぶしのための関数tpaint()の定義です。

tpaint <- function(x1, x2,color)
{
xvals <- seq(x1, x2, length=50)
dvals <- c(DF(xvals), 0, 0)
xvals <- c(xvals, x2, x1)
polygon(xvals, dvals, col=color)
}

　もし-4～-2、2～4の範囲を塗りつぶすのであれば次のようにします。

tpaint(-4,-2,color="pink")
tpaint(2,4,color="pink")