大数跨境
0
0

R语言计算Logistic的efect和OR值以及置信区间

R语言计算Logistic的efect和OR值以及置信区间 育种数据分析之放飞自我
2023-11-26
0
导读:各位小伙伴,大家好,我是邓飞,今天介绍一下,如何使用R语言进行logistic分析,并且计算OR值和置信区间

各位小伙伴,大家好,我是邓飞,今天介绍一下,如何使用R语言进行logistic分析,并且计算OR值和置信区间。

OR值的定义:

OR 值(odds ratio)是一种统计量,用于度量两个事件的概率发生的相对大小。在二分类 GWAS(基因组关联研究)分析中,OR 值用于衡量某个基因变异(或基因型)与某个疾病(或特征)之间的关联程度。
OR 值大于 1 表示基因变异和疾病间存在正相关,即该变异可能增加患病的概率;OR 值小于 1 则表示基因变异和疾病之间存在负相关,即该变异可能减少患病的概率。一般情况下,OR 值越大表示基因变异和疾病间的关联程度越强。
在二分类 GWAS 分析中,通过计算每个基因变异的OR值,可以评估其与疾病之间的关联程度,从而推断基因变异对疾病风险的贡献。这可以帮助研究人员了解潜在的遗传基础以及疾病风险的变化。

plink的Logisitic模型的GWAS分析计算结果如下:

R语言的解决方案:

m1 = glm(phe.V3 ~ rs3131972_A,family = "binomial",data=dd )summary(m1)
# 计算OR值exp(coef(mod))## 置信区间exp(confint(mod))
# 一步到位:OR值和置信区间library(questionr)odds.ratio(mod)

结果:


手动计算OR值:


一步到位的OR值和置信区间:


搞定!

1,快来领取 | 飞哥的GWAS分析教程


2,飞哥汇总 | 入门数据分析资源推荐


3,数量遗传学,分享几本书的电子版


4,R语言学习看最新版的电子书不香嘛?


5,书籍及配套代码领取--统计遗传分析导论


6,飞哥的学习圈子


【声明】内容源于网络
0
0
育种数据分析之放飞自我
本公众号主要介绍动植物育种数据分析中的相关问题, 算法及程序代码.
内容 912
粉丝 0
育种数据分析之放飞自我 本公众号主要介绍动植物育种数据分析中的相关问题, 算法及程序代码.
总阅读370
粉丝0
内容912