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

样本间亲缘关系矩阵(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)
相关推荐
我爱C编程10 小时前
基于拓扑结构检测的LDPC稀疏校验矩阵高阶环检测算法matlab仿真
算法·matlab·矩阵·ldpc·环检测
CVer儿12 小时前
svd分解求旋转平移矩阵
线性代数·算法·矩阵
张晓~1833994812117 小时前
数字人分身+矩阵系统聚合+碰一碰发视频: 源码搭建-支持OEM
线性代数·矩阵·音视频
YuTaoShao2 天前
【LeetCode 热题 100】73. 矩阵置零——(解法二)空间复杂度 O(1)
java·算法·leetcode·矩阵
szekl3 天前
HDMI 2.0 4×2矩阵切换器412HN——多信号输入输出的高清解决方案
linux·矩阵·计算机外设·电脑·ekl
引量AI13 天前
TikTok 矩阵如何快速涨粉
大数据·人工智能·矩阵·tiktok矩阵·海外社媒
Ven%14 天前
矩阵阶数(线性代数) vs. 张量维度(深度学习):线性代数与深度学习的基石辨析,再也不会被矩阵阶数给混淆了
人工智能·pytorch·深度学习·线性代数·矩阵·tensor·张量
云云32114 天前
亚矩云手机赋能Vinted矩阵运营:破解二手电商多账号与本地化困局
网络·人工智能·智能手机·矩阵·自动化
云云32115 天前
Subway Surfers Blast × 亚矩阵云手机:手游矩阵运营的终极变现方案
大数据·人工智能·线性代数·智能手机·矩阵·架构
点云侠15 天前
PCL 点云旋转的轴角表示法
人工智能·线性代数·算法·计算机视觉·矩阵