![見出し画像](https://assets.st-note.com/production/uploads/images/112340193/rectangle_large_type_2_2c02c43905e6fa98902f696c2277b0f6.png?width=1200)
Photo by
c_midori
t検定で何をやっているか?を確認するRコード
![](https://assets.st-note.com/img/1690964495573-HygM3kkKkJ.png)
count<-0
t<-numeric(length=100)
print(-qt(0.025,50))
for (i in 1:100){
t[i]<-rt(n=1, 50)
if (t[i]> -qt(0.025,50)){
count<-count+1
}
}
print(count)
count<-0t<-numeric(length=100)
print(-qt(0.025,50))
for (i in 1:100){
t[i]<-rt(n=1, 50) if (t[i]> -qt(0.025,50)){
count<-count+1 }
}
print(count)
実行してみると、自由度50のt分布について、上側5%qt(0.025,50)の絶対値より大きい値が自由度50のt分布の乱数ひとつrt(n=1, 50)からは、100回発生させても4回しか生じませんでした。これは5%(100回中5回)より小さい値でした。
![](https://assets.st-note.com/img/1690964576314-Oz5aHjntew.png)
この記事が気に入ったらサポートをしてみませんか?