Searched refs:prev_varcmp (Results 1 – 4 of 4) sorted by relevance
/dports/biology/gcta/gcta_1.26.0_src/ |
H A D | reml_within_family.cpp | 59 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 D | est_hsq.cpp | 899 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 D | bivar_reml.cpp | 291 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 D | gcta.h | 246 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);
|