事象がおきたかどうかについて、ある値が高いか低いかでわけるとき R ではどうするか
データを用意
読み込み
d2 <- read.table(pipe("pbpaste"), header=TRUE)
subsetを使って切り分け、そのグループをt検定(あるいはt.test(LHratio~ACS,d2) だけでもいいっぽい)
survive <- subset(d2,ACS==0)$LHratio acs <- subset(d2,ACS==1)$LHratio t.test(survive,acs)
結果、どうやら有意差があるっぽい
Welch Two Sample t-test data: survive and acs t = -2.4542, df = 75.68, p-value = 0.01642 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -0.84383119 -0.08776028 sample estimates: mean of x mean of y 2.116852 2.582648
2群の箱表示(これも boxplot(LHratio~ACS,d2)とできる)
boxplot(survive,acs,names=c("ACS(-)","ACS(+)"),boxwex=0.5)