楼主: mumianke
2117 3

请教:循环问题。 [推广有奖]

  • 0关注
  • 0粉丝

高中生

52%

还不是VIP/贵宾

-

威望
0
论坛币
1080 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
1114 点
帖子
26
精华
0
在线时间
25 小时
注册时间
2007-3-20
最后登录
2018-7-1

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币

对每一行做卡方检验。
在样本较小时(<5),需要用Fisher精确检验来作独立性检验。
如果在data.frame中有的个数小于5,要用if 语句进行筛选。
test.txt如下:

aa   bb   cc   dd
45   24   55   76
85   52   15   48
12   41   88   59
45   14   55   86
78   25   22   75
60   3   32   11
....

最后一行是添加上验证fisher.test的
请问:如何用条件语句。
test<-read.table(file='test.txt')
result <- numeric()
for(i in 1 : nrow(test)){  
if (test$aa>4 && test$bb>4 && test$cc>4 && test$dd>4)
result<-chisq.test(matrix(unlist(test[i,]),ncol=2,byrow=T),correct=F)$p.value

else
result<-fisher.test(matrix(unlist(test[i,]),ncol=2,byrow=T))$p.value

}

结果不对。请帮帮忙,谢谢。

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:Result RESUL 卡方检验 test Est result 样本

回帖推荐

philzang 发表于2楼  查看完整内容

I add some trivial modifications on the code :result <- numeric()for(i in 1 : nrow(test)){   if (test$aa>4 && test$bb>4 && test$cc>4 && test$dd>4)result<-chisq.test(matrix(unlist(test),ncol=2,byrow=T),correct=F)$p.valueelseresult<-fisher.test(matrix(unlist(test),ncol=2,byrow=T))$p.value}  # And result is > result[1] 1.785716e-03 5.076277e-07 3.377813e-06 ...

本帖被以下文库推荐

沙发
philzang 发表于 2008-6-2 22:56:00 |只看作者 |坛友微信交流群

I add some trivial modifications on the code :

result <- numeric()

for(i in 1 : nrow(test)){  
if (test$aa>4 && test$bb>4 && test$cc>4 && test$dd>4)
result<-chisq.test(matrix(unlist(test[i,]),ncol=2,byrow=T),correct=F)$p.value

else
result<-fisher.test(matrix(unlist(test[i,]),ncol=2,byrow=T))$p.value

# And result is

> result
[1] 1.785716e-03 5.076277e-07 3.377813e-06 1.534920e-06 6.445038e-14
[6] 2.820388e-03

已有 1 人评分论坛币 学术水平 收起 理由
crystal8832 + 10 + 1 热心帮助其他会员

总评分: 论坛币 + 10  学术水平 + 1   查看全部评分

使用道具

藤椅
mumianke 发表于 2008-6-3 09:55:00 |只看作者 |坛友微信交流群
多谢,看来我得多学习了。

使用道具

板凳
wzy1713 发表于 2008-6-3 10:46:00 |只看作者 |坛友微信交流群
&nbsp; 强人太多了

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-11-6 03:51