TASSEL的MLM模型构建的kinship矩阵相关知识

今天，看一下TASSEL的MLM模型构建的kinship矩阵是如何计算的？

1. 导入基因型数据

数据及代码下载，请关注公众号：育种数据分析之放飞自我，进入知识星球进行相关下载和学习

这里导入vcf格式的数据：

2. 导入表型数据和协变量

表型数据：

协变量文件：

3. 构建kinship矩阵

3.1 Centered_IBS

这种方法，应该就是VanRaden的方法，中心化的IBS亲缘关系矩阵。

构建的kinship矩阵：

R语言比较：

两者一致。

3.2 Normalized_IBS

这种方法应该就是Yang的方法：

结果：

R语言对比：

两者结果完全一致。

3.3 Dominance_Centered_IBS

这个，应该是显性中心化的IBS

结果：

R语言结果对比：

结果完全一致。

3.4 Dominance_Normalized_IBS

这个应该是显性矩阵标准化的IBS矩阵。

结果：

暂时，未找到R中对应矩阵计算的方法。

4. kinship矩阵如何判断异常

1，可以将kinship做热图聚类，查看分布，类似：

2，一般对角线<0.8，或者>1.2的个体，可以判断是离群样本，可以通过PCA看一下其分布

3，非对角线为负值，一般可以认为其为0，如果负值很大，就要看一下是否异常

TASSEL的MLM模型构建的kinship矩阵相关知识相关推荐

关于构建网络安全知识库方向相关知识的学习和思考
说明:最近在看<面向知识服务的知识库结构理论与方法>蒋勋,将自己研究可能用到的知识进行梳理和摘录,并加入部分自己的主观想法,由于17年的书,有些思想或描述已过时,但可借鉴的思想还是有不少的 ...
矩阵相关知识回顾--协方差的意义
协方差协方差矩阵的解释:A geometric interpretation of the covariance matrix 这里通过探索线性变换与所得数据协方差之间的关系,提供协方差矩阵的直观 ...
差异表达基因热图怎么看_四个基因的能量代谢模型构建与验证，轻松发4分sci！...
今天分享一篇2019年6月发表于J Cell Physiol(IF:4.522)的文章,标题是Identification of a four‐gene metabolic signature pre ...
R语言构建logistic回归模型并评估模型：计算混淆矩阵、并基于混淆矩阵计算Accuray、Precision、Recall(sensitivity)、F1、Specificity指标
R语言构建logistic回归模型并评估模型:计算混淆矩阵.并基于混淆矩阵计算Accuray.Precision.Recall(sensitivity).F1.Specificity指标目录
GWAS分析中SNP解释百分比PVE | 第四篇，MLM模型中如何手动计算PVE？
系列部分: GWAS分析中SNP解释百分比PVE | 第一篇,SNP解释百分比之和为何大于1? GWAS分析中SNP解释百分比PVE | 第二篇,GLM模型中如何计算PVE? GWAS分析中SNP解释 ...
必须要GPT-3吗？不，BERT的MLM模型也能小样本学习
©PaperWeekly 原创 · 作者|苏剑林单位|追一科技研究方向|NLP.神经网络大家都知道现在 GPT-3 风头正盛,然而,到处都是 GPT-3.GPT-3 地推,读者是否记得 GPT- ...
修改Transformer结构，设计一个更快更好的MLM模型
©PaperWeekly 原创 · 作者|苏剑林单位|追一科技研究方向|NLP.神经网络大家都知道,MLM (Masked Language Model) 是 BERT.RoBERTa 的预训练 ...
二十八、电力窃漏电案例模型构建
1. 综合案例模型构建构建窃漏电用户识别模型: 构建LM神经网络模型: 构建CART决策树模型: 模型评价 2. 构建窃漏电用户识别模型 2.1 构建专家样本专家样本准备完成后,需要划分测试样本和 ...
某农商行用户画像项目——模型构建部分
文章目录三.模型构建 3.1 贷款逾期数据分析 3.1.1 数据加载与预览 3.1.2 探索性数据分析 3.1.2.1 描述性分析 3.1.2.1.1 类别型变量 3.1.2.1.2 数值型变量 3 ...