Home
last modified time | relevance | path

Searched refs:Xt_Vi_X_i (Results 1 – 3 of 3) sorted by relevance

/dports/biology/gcta/gcta_1.26.0_src/
H A Dest_hsq.cpp642 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 Dbivar_reml.cpp405 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 Dgcta.h244 …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);