|
#母分散が未知＝t分布
# （１）簡単にする方法
# 参考 赤間・山口[2006],p109-110
#平均が50のデータ・正規分布の乱数100個を発生
x1 = rnorm(100,50,1);x1
t.test(x1)
x2 = rnorm(100,50,10);x2
t.test(x2)
#平均が50のデータ・正規分布の乱数10000個を発生
x3 = rnorm(10000,50,1);x3
t.test(x3)
x4 = rnorm(10000,50,10);x4
t.test(x4)
# （２）定義にしたがい導出
# 基本データ
m1 = mean(x1);v1 =var(m1);n1=length(x1);df1=n1-1
# 95%信頼区間
t1 = qt(1-0.05/2,df1);t1 #統計量
u1 = (1/(n1-1))*sum((x1-m1)^2)
T1 = m1-t1*sqrt(u1/n1)
T2 = m1+t1*sqrt(u1/n1)
T1;T2
|
