/dports/biology/gcta/gcta_1.26.0_src/ |
H A D | mlm_assoc.cpp | 241 float *Vi_X=new float[n]; in mlma_calcu_stat() local 269 cblas_sgemv(CblasRowMajor, CblasNoTrans, n, n, 1.0, Vi, n, X, 1, 0.0, Vi_X, 1); in mlma_calcu_stat() 270 Xt_Vi_X=cblas_sdot(n, X, 1, Vi_X, 1); in mlma_calcu_stat() 272 beta[i]=se[i]*cblas_sdot(n, y, 1, Vi_X, 1); in mlma_calcu_stat() 280 delete[] Vi_X; in mlma_calcu_stat() 291 float *Vi_X=new float[n*col_num]; in mlma_calcu_stat_covar() local 325 …asRowMajor, CblasNoTrans, CblasNoTrans, n, col_num, n, 1.0, Vi, n, X, col_num, 0.0, Vi_X, col_num); in mlma_calcu_stat_covar() 326 …RowMajor, CblasTrans, CblasNoTrans, col_num, col_num, n, 1.0, X, col_num, Vi_X, col_num, 0.0, Xt_V… in mlma_calcu_stat_covar() 328 … cblas_sgemv(CblasRowMajor, CblasTrans, n, col_num, 1.0, Vi_X, col_num, y, 1, 0.0, Xt_Vi_y, 1); in mlma_calcu_stat_covar() 340 delete[] Vi_X; in mlma_calcu_stat_covar()
|
H A D | est_hsq.cpp | 642 eigenMatrix Vi_X(_n, _X_c), Xt_Vi_X_i(_X_c, _X_c), Hi(_r_indx.size(), _r_indx.size()); in reml() local 645 …double lgL = reml_iteration(Vi_X, Xt_Vi_X_i, Hi, Py, varcmp, reml_priors_var_flag | reml_priors_fl… in reml() 654 if (est_fix_eff) _b = Xt_Vi_X_i * (Vi_X.transpose() * _y); in reml() 874 eigenMatrix Vi_X(_n, _X_c), Xt_Vi_X_i(_X_c, _X_c), Hi(_r_indx.size(), _r_indx.size()); in lgL_reduce_mdl() local 879 double lgL = reml_iteration(Vi_X, Xt_Vi_X_i, Hi, Py, varcmp, false, no_constrain); in lgL_reduce_mdl() 884 double gcta::reml_iteration(eigenMatrix &Vi_X, eigenMatrix &Xt_Vi_X_i, eigenMatrix &Hi, eigenVector… in reml_iteration() argument 929 logdet_Xt_Vi_X = calcu_P(_Vi, Vi_X, Xt_Vi_X_i, _P); // Calculate P in reml_iteration() 1245 double gcta::calcu_P(eigenMatrix &Vi, eigenMatrix &Vi_X, eigenMatrix &Xt_Vi_X_i, eigenMatrix &P) in calcu_P() argument 1247 Vi_X = Vi*_X; in calcu_P() 1248 Xt_Vi_X_i = _X.transpose() * Vi_X; in calcu_P() [all …]
|
H A D | bivar_reml.cpp | 405 eigenMatrix Vi_X(_n, _X_c), Xt_Vi_X_i(_X_c, _X_c), Hi(_r_indx.size(), _r_indx.size()); in lgL_fix_rg() local 412 double lgL = reml_iteration(Vi_X, Xt_Vi_X_i, Hi, Py, varcmp, false, no_constrain, true); in lgL_fix_rg()
|
H A D | gcta.h | 244 …double reml_iteration(eigenMatrix &Vi_X, eigenMatrix &Xt_Vi_X_i, eigenMatrix &Hi, eigenVector &Py,… 254 double calcu_P(eigenMatrix &Vi, eigenMatrix &Vi_X, eigenMatrix &Xt_Vi_X_i, eigenMatrix &P);
|