循环 chisq.test

1. x1<-c(10,20,30,30)
   
2. x2<-c(1,10,2,10)
   
3. x3<-c(2,10,5,40)
   
4. aaa<-data.frame(x1,x2,x3)
   
5. aaa<-t(aaa)
   
6. colnames(aaa)<-c("a","b","c","d")
   
7. aaa

复制代码

数据框aaa中每一行是四格表的a b c d，循环进行卡方检验，
如果n>=40 and T>=5      进行卡方检验
如果n>=40 and 1<=T<5  进行校正卡方检验
如果n<40   and T<1         不进行检验

结果要求在数据框aaa后面加上两列，分别是卡方值和P值，谢谢！

You can take a look at this video.

https://zh.coursera.org/learn/spss-ruanjian/lecture/WARgh/5-2-du-li-yang-ben-si-ge-biao-zi-liao-de-qia-fang-jian-yan