构建H逆矩阵需要考虑的参数
?...矩阵的尺度
原因:
G矩阵构建是由测序个体的基因频率构建的,它的频率可能与A矩阵的基因频率(可以追溯到基础群体base population)不一样, 这导致G矩阵和A22矩阵存在尺度上的差异, 因此需要矫正...代码思路:
依赖的包: NamedArray, 主要用于矩阵提取
参数介绍:
id_full: 为所有的ID, A矩阵的id
id_geno: 为测序的ID, G矩阵的id
Amatrix: A矩阵
Gmatrix...提取A11, A12, A21, A22
根据对角线和非对角线方程组, 计算a和b
将相关参数加进去, 构建H逆矩阵
function hmatrix_julia_adjust(id_full,id_geno...(maximum(offdiagA22),digits=3),"\t","\n\n")
# 根据方程计算alpha和beta两个参数值 meanoffdiagG=mean(offdiagG