理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

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)