Home
last modified time | relevance | path

Searched refs:prev_varcmp (Results 1 – 4 of 4) sorted by relevance

/dports/biology/gcta/gcta_1.26.0_src/
H A Dreml_within_family.cpp59 bool gcta::calcu_Vi_within_family(eigenMatrix &Vi, eigenVector &prev_varcmp, double &logdet, int &i… in calcu_Vi_within_family() argument
67 Vi.diagonal()=eigenVector::Constant(_n, 1.0/prev_varcmp[0]); in calcu_Vi_within_family()
68 logdet=_n*log(prev_varcmp[0]); in calcu_Vi_within_family()
71 for(i=0; i<_r_indx.size(); i++) Vi+=(_Asp[_r_indx[i]])*prev_varcmp[i]; in calcu_Vi_within_family()
H A Dest_hsq.cpp899 if (reml_bivar_fix_rg) update_A(prev_varcmp); in reml_iteration()
901 prev_varcmp = varcomp_init; in reml_iteration()
932 else if (_reml_mtd == 2) em_reml(_P, Py, prev_varcmp, varcmp); in reml_iteration()
964 varcmp = prev_varcmp; in reml_iteration()
972 varcmp = prev_varcmp; in reml_iteration()
988 prev_prev_varcmp = prev_varcmp; in reml_iteration()
989 prev_varcmp = varcmp; in reml_iteration()
1086 logdet = _n * log(prev_varcmp[0]); in calcu_Vi()
1353 if (dlogL > 1.0) varcmp = prev_varcmp + 0.316 * delta; in ai_reml()
1354 else varcmp = prev_varcmp + delta; in ai_reml()
[all …]
H A Dbivar_reml.cpp291 bool gcta::calcu_Vi_bivar(eigenMatrix &Vi, eigenVector &prev_varcmp, double &logdet, int &iter) { in calcu_Vi_bivar() argument
298 for (i = 0; i < _r_indx.size(); i++) Vi += (_Asp[_r_indx[i]]) * prev_varcmp[i]; in calcu_Vi_bivar()
382 double gcta::lgL_fix_rg(eigenVector &prev_varcmp, bool no_constrain) { in lgL_fix_rg() argument
409 varcmp[i] = fabs(prev_varcmp[_r_indx[i]]); in lgL_fix_rg()
419 void gcta::update_A(eigenVector &prev_varcmp) { in update_A() argument
423 g1 = prev_varcmp[_bivar_pos[0][i]]; in update_A()
424 g2 = prev_varcmp[_bivar_pos[1][i]]; in update_A()
H A Dgcta.h246 bool calcu_Vi(eigenMatrix &Vi, eigenVector &prev_varcmp, double &logdet, int &iter);
258 void em_reml(eigenMatrix &P, eigenVector &Py, eigenVector &prev_varcmp, eigenVector &varcmp);
259 …void ai_reml(eigenMatrix &P, eigenMatrix &Hi, eigenVector &Py, eigenVector &prev_varcmp, eigenVect…
268 … bool calcu_Vi_within_family(eigenMatrix &Vi, eigenVector &prev_varcmp, double &logdet, int &iter);
272 void update_A(eigenVector &prev_varcmp);
274 double lgL_fix_rg(eigenVector &prev_varcmp, bool no_constrain);
275 bool calcu_Vi_bivar(eigenMatrix &Vi, eigenVector &prev_varcmp, double &logdet, int &iter);