发帖
每发1个新帖
可以获得0.5个丁当奖励
回帖

【求助】haplo.stats 计算进行单倍体分析，能计算OR吗

只看楼主
页码直达：
直达末页

楼主

longjx12345

入门站友

1楼

这个帖子发布于7年零246天前，其中的信息可能已发生改变或有所发展。

问题已解决悬赏丁当:2

haplo.stats 计算进行单倍体分析，能计算OR和95%可信区间吗？我见别人的文章的能给出，如果能怎么计算，如果不能别人是怎么算出来的？文章理都没交代清楚

邀请讨论

搜索

2013-03-30 08:28 浏览 : 2837 回复 : 7

投票
收藏 1
打赏
引用
分享
- 微信扫一扫
- 新浪微博
- 丁香客
- 复制网址
举报
- 广告宣传推广
- 政治敏感、违法虚假信息
- 恶意灌水、重复发帖
- 违规侵权、站友争执
- 附件异常、链接失效
- 其他

• 全科医生转岗培训

wudaliming

丁香园准中级站友

2楼

你用haplo.stats做过单倍型分析没？如果做过的话就会发现OR和95%置信区间和P值是一起计算出来的。haplo.stats是在R语言环境下操作的，楼主会操作R软件不？

2013-03-30 12:12

投票
收藏
打赏
引用
分享
- 微信扫一扫
- 新浪微博
- 丁香客
- 复制网址
举报
- 广告宣传推广
- 政治敏感、违法虚假信息
- 恶意灌水、重复发帖
- 违规侵权、站友争执
- 附件异常、链接失效
- 其他

• 从老家公立医院辞职后2年又想回去正常吗？

楼主

longjx12345

入门站友

3楼

2013-03-31 09:39

投票
收藏
打赏
引用
分享
- 微信扫一扫
- 新浪微博
- 丁香客
- 复制网址
举报
- 广告宣传推广
- 政治敏感、违法虚假信息
- 恶意灌水、重复发帖
- 违规侵权、站友争执
- 附件异常、链接失效
- 其他

• 公立还是私立？编制还是合同？大城市还是小乡镇？这些问题一次性解决！

楼主

longjx12345

入门站友

4楼

我的只能算出分值和P值，用的是haplo.score，好像 haplo.cc是可以给出所有结果，但是不能校正性别年龄 haplo.cc命令为

cc.hla <- haplo.cc(y=Y01.vec,geno=Geno,miss.val=NA, locus.label=NA, control = haplo.glm.control(haplo.freq.min = 0.02))

cc.hla
haplo.score的命令为，可以校正性别年龄

Y01.hap.reg.FD <- haplo.score(y=Y01.vec,geno=Geno,trait.type="gaussian",miss.val=NA, locus.label=NA, x.adj=X.adj,skip.haplo)

Y01.hap.reg.FD

我想修改haplo.cc命令校正性别年龄，加了adj=X.adj, X.adj <- data.mat[,c(2:3)]已经赋值，但是出错了
cc.hla <- haplo.cc(y=Y01.vec,geno=Geno,miss.val=NA, locus.label=NA,x.adj=X.adj, control = haplo.glm.control(haplo.freq.min = 0.02))
错误于haplo.cc(y = Y01.vec, geno = Geno, miss.val = NA, locus.label = NA, :
变元((x.adj = X.adj)) 没有用
求高手帮我把haplo.cc命令改成可以校正性别年龄的，谢谢

2013-03-31 09:40

投票
收藏
打赏
引用
分享
- 微信扫一扫
- 新浪微博
- 丁香客
- 复制网址
举报
- 广告宣传推广
- 政治敏感、违法虚假信息
- 恶意灌水、重复发帖
- 违规侵权、站友争执
- 附件异常、链接失效
- 其他

longjx12345 编辑于 2013-03-31 10:10

• 耗材再次砍，骨科，心内，普外等集体采购！

遗传发育

【求助】haplo.stats 计算进行单倍体分析，能计算OR吗

提示

丁香园旗下网站

关于丁香园

官方链接

遗传发育

分享到：

微信扫一扫

【求助】haplo.stats 计算进行单倍体分析，能计算OR吗

微信扫一扫

微信扫一扫

微信扫一扫

微信扫一扫

提示

丁香园旗下网站

关于丁香园

官方链接