全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2324 11
2022-05-31

微信图片_1.png
我本来想做一个上图所示的2*4列联表进行卡方检验
根据fertilization和TPOAb情况分成4组
fertilization是一个二分类变量,不是0就是1.
TPOAb虽然是具体数值,但是以5.61为界可以分为两类
所以我输入以下代码
微信图片_2.png
但是返回的结果却是一个data.frame。
微信图片_3.png
请问这是怎么回事呀,该如何解决呢?
二维码

扫码加我 拉你入群

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

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

全部回复
2022-6-1 08:19:32
代码中,ABCD是什么啊?
如果是分组的名称,可以把A改成“A”,其他类推。也就是说分组的名称需要加引号。
二维码

扫码加我 拉你入群

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

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

2022-6-9 00:22:22
谢谢您,上面的问题已经解决了!我还有个问题想请教您:
微信图片_4.png
如QQ图图所示,该数据存在几个极大极小值,并且这几个极值可能会影响我后面的分析,所以我想把它们删除,请问我该如何才能删除这些极值呢?非常感谢!
二维码

扫码加我 拉你入群

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

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

2022-6-10 16:21:59
几种方法:
1、这好像是qqnorm()做的图吧,从图中可以大概估计,假设大于40的是极值。
2、或者你用分位数函数:
upperLim<- quantile(y, 0.975)
看一下0.975(这是一个主观标准,也可以改为0.95或其他)分位对应的y值。
3、也可以用evir包的函数:
upperLim <- findthresh(y, 3) #确定y有3个极值时的阈值。

然后:
new<-data[which(data$y<upperLim), ]就可以删除存在极值的样本了。
二维码

扫码加我 拉你入群

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

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

2022-6-12 10:32:39
llb_321 发表于 2022-6-10 16:21
几种方法:
1、这好像是qqnorm()做的图吧,从图中可以大概估计,假设大于40的是极值。
2、或者你用分位数 ...
是的,感谢大神!!!
二维码

扫码加我 拉你入群

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

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

2022-6-14 13:36:40
大神您好,我又来打扰您了
我想做交互作用,需要安装epiR包,但是提示需要加载Rtools(如下图)
微信图片_5.png
然后我就去下载Rtools,但是貌似也没有安装成功(如下图)
微信图片_6.png
请问这是怎么回事,或者有什么其他的方法可以做相加交互作用分析吗~感谢您~

二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群