R语言迅速计算多基因评分(PRS)

Polygenic Risk Scores in R

最朴素的理解PRS:

GWAS分析结果中,有每个SNP的beta值、se值、P值,因为GWAS分析中将SNP变为0-1-2编码,所以这些显著的SNP的beta值,就可以用于预测。

比如:GWAS分析中,显著的SNP效应值为:

SNP1: 0.3

SNP2: 0.2

SNP3: -0.1

对于target data(目标群体),检测了3个个体,3个SNP的分型分别为:

ID1 0 0 1

ID2 1 0 2

ID3 2 2 1

那么个体1的多基因评分为:00.3 + 00.2 + 1*-0.1 = -0.1

个体2的多基因评分为:0.3 + 0 + -0.1 = 0.2

个体3的多基因评分为:0.6 + 0.4 + -0.1 = 0.9

用数学公式表示:

  • beta是效应值
  • G是0-1-2的编码
  • m是m个SNP

实际项目的PRS计算

实际中的项目,考虑的因素比较多,比如:

  • 数据质控
  • 群体结构
  • LD值(clumping)
  • beta矫正值
  • 通过P值筛选最优组合

相关软件实现PRS分析

  • plink
  • biqsnpr,一个R包
  • PRSice,应用最广泛,通过C+T的策略
  • LDpred,通过贝叶斯收缩的模型
  • PRS-CS
  • JAMPred
  • Lassosum

之前写过PRS的操作流程,可以作为参考:

多基因风险预测模型1--先立Flag

多基因风险预测模型2--相关概念和软件

不会安装使用PRSice-2软件就太不讲究了

相关推荐
菜鸟小贤贤几秒前
pyhton+yaml+pytest+allure框架封装-全局变量接口关联
开发语言·python·macos·自动化·pytest
南东山人16 分钟前
python问题解决-外部模块明明安装了,却总是无法找到
开发语言·python
Adolf_199324 分钟前
Django 自定义路由转换器
后端·python·django
B站计算机毕业设计超人24 分钟前
计算机毕业设计Python+大模型美食推荐系统 美食可视化 美食数据分析大屏 美食爬虫 美团爬虫 机器学习 大数据毕业设计 Django Vue.js
大数据·爬虫·python·深度学习·机器学习·课程设计·推荐算法
ThetaarSofVenice27 分钟前
【Java从入门到放弃 之 Java程序基础】
java·开发语言·python
XMYX-030 分钟前
使用 Django 构建支持 Kubernetes API 测试连接的 POST 接口
python·kubernetes·django
使者大牙32 分钟前
【单点知识】基于PyTorch进行模型部署
人工智能·pytorch·python·深度学习
L Jiawen38 分钟前
【Python · PyTorch】循环神经网络 RNN(基础概念)
pytorch·python·rnn
int WINGsssss38 分钟前
对pytorch的底层nccl库进行插桩
人工智能·pytorch·python
qq_189370491 小时前
鸢尾花Iris训练数据和测试数据的分割和训练数据的散点图矩阵绘制
python·线性代数·机器学习·散点图矩阵