1標本t検定のスクリプト
只今練習中~
2標本xとyについて
sumxy <- function(x,y,sxy)
{
xs<-summary(x)
ys<-summary(y)
xt<-t.test(x)
yt<-t.test(y)
print(xs);print(ys);print(xt);print(yt)
}
sumxy (x,y,sxy)
より
x<-c( );y<-c( )にデータを入れると、
各標本の要約・t値・p値・95%CIのみを出力させるプログラム
sumxy <- function(x,y,sxy)
{
cat1<-c("Min.","1st Qu","Median","Mean", "3rd Qu.","Max.")
xs<-summary(x)
ys<-summary(y)
cat(cat1,"\n","x",xs,"\n","y",ys,"\n\n",file = "TEST001") #各標本のまとめ
tx<-t.test(x)
cat2<-c("tvalue","pvalue","95%CI")
txt<- tx$statistic
txp<- tx$p.value
txc1<- tx$conf.int[1]
txc2<- tx$conf.int[2]
x1 <- sprintf("%.3f\t%.3f\t%.3f\t%.3f",txt,txp,txc1,txc2)
ty<-t.test(y)
tyt<- ty$statistic
typ<- ty$p.value
tyc1<- ty$conf.int[1]
tyc2<- ty$conf.int[2]
y1 <- sprintf("%.3f\t%.3f\t%.3f\t%.3f",tyt,typ,tyc1,tyc2)
cat(cat2,"\n","x",x1,"\n","y",y1,"\n",file = "TEST001", append =T)
}
sumxy (x,y,sxy)