Searched refs:Xt_Vi_X_i (Results 1 – 3 of 3) sorted by relevance
/dports/biology/gcta/gcta_1.26.0_src/ |
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() 738 cout << setiosflags(ios::fixed) << _b[i] << "\t" << sqrt(Xt_Vi_X_i(i, i)) << endl; in reml() 795 …for (i = 0; i < _X_c; i++) o_reml << setprecision(6) << _b[i] << "\t" << sqrt(Xt_Vi_X_i(i, i)) << … 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() 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 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);
|