构建H逆矩阵Julia代码展示
?
Julia是新一代的计算科学语言, 速度非常快, 比R语言、Python要快很多, 内存占用也比较少. 比C语言、Fortran要简单很多,是一门未来的明星语言....代码思路:
依赖的包: 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...meandiagA22=mean(diagA22) print("Means_off_diag\t","Means_diag:\n","G\t",meanoffdiagG,"\t",meandiagG, "\nA22...Adjust G, and the value of alpha and beta is:",alpha,beta,"\n\n\n") G=alpha .+ beta .* G # # 在将a