重测序关系矩阵构建方式汇总

样本间亲缘关系矩阵(kinship matrix)和同源性矩阵(IBS matrix)构建的方式

1. 可以使用plink的--make-rel计算个体之间的亲缘关系(强调个体之间的遗传相似性)

bash 复制代码
/opt/software/plink --bfile vcf_bfile--make-rel --out relatedness_matrix # 得到亲缘关系距离矩阵:
# relatedness_matrix.rel

2. kinship

bash 复制代码
# 利用tassel计算
run_pipeline.pl -Xmx1536m-Xms512m -SortGenotypeFilePlugin -inputFile 你的vcf文件 -outputFile outvcf -fileType VCF
run_pipeline.pl-Xmx1536m -Xms512m -importGuess outvcf -KinshipPlugin -methodCentered_IBS -endPlugin -export tassel_kinship.txt -exportType SqrMatrix
# 利用gcta计算
使用 --make-grm-alg 1 或 --make-grm 0
gcta --make-grm --make-grm-alg 1 --out snp.gcta --bfile vcf_bfile snp --autosome-num 90

3. IBS

bash 复制代码
/opt/software/plink --bfile  vcf_bfile --make-bed --out IBS_matrix --maf 0.05 --recode --double-id --allow-extra-chr --chr-set 90  --distance square ibs
要计算遗传距离,使用1-ibs

群体关系矩阵如何构建?

转换方法:平均IBS(个体对间均值)

计算所有个体两两之间的IBS均值,反映群体内遗传相似性。

python 复制代码
#算法示例
import numpy as np

def pairwise_ibs(genotype_matrix):
    n_individuals = genotype_matrix.shape[0]
    ibs_matrix = np.zeros((n_individuals, n_individuals))
    
    for i in range(n_individuals):
        for j in range(i+1, n_individuals):
            ibs_sum = 0
            for snp in range(genotype_matrix.shape[1]):
                ibs_sum += calculate_ibs(
                    genotype_matrix[i, snp, 0], genotype_matrix[i, snp, 1],
                    genotype_matrix[j, snp, 0], genotype_matrix[j, snp, 1]
                )
            ibs_matrix[i, j] = ibs_sum / genotype_matrix.shape[1]
            ibs_matrix[j, i] = ibs_matrix[i, j]  # 对称矩阵
    return ibs_matrix

# 示例基因型矩阵(个体数×SNP数×2等位基因)
genotype_data = np.array([
    [['A', 'A'], ['G', 'G']],  # 个体1
    [['A', 'G'], ['G', 'G']],  # 个体2
    [['T', 'T'], ['A', 'G']]   # 个体3
])
ibs_matrix = pairwise_ibs(genotype_data)
print("群体IBS矩阵:\n", ibs_matrix)
相关推荐
Eloudy4 小时前
复矩阵与共轭转置矩阵乘积及其平方根矩阵
人工智能·算法·矩阵
m0_631354454 小时前
VTK开发day2:切片矩阵
人工智能·算法·矩阵
go54631584654 小时前
在本地环境中运行 ‘dom-distiller‘ GitHub 库的完整指南
人工智能·深度学习·神经网络·算法·矩阵·github
PerfumerKarma10 小时前
【WebGPU学习杂记】数学基础拾遗(2)变换矩阵中的齐次坐标推导与几何理解
学习·线性代数·矩阵
牵牛老人1 天前
OpenCV学习探秘之二 :数字图像的矩阵原理,OpenCV图像类与常用函数接口说明,及其常见操作核心技术详解
opencv·学习·矩阵
lqjun08271 天前
相机内外参矩阵:从3D世界坐标到2D像素坐标变换
数码相机·3d·矩阵
恣艺1 天前
LeetCode 1074:元素和为目标值的子矩阵数量
算法·leetcode·矩阵
Alfred king1 天前
面试150 搜索二维矩阵
线性代数·矩阵·二分查找·数组
Alfred king2 天前
面试150 建立四叉树
矩阵··数组·分治
石去皿2 天前
QKV 为什么是三个矩阵?注意力为何要除以 √d?多头注意力到底有啥用?
人工智能·线性代数·机器学习·矩阵