Home
last modified time | relevance | path

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

/dports/biology/gcta/gcta_1.26.0_src/
H A Dmlm_assoc.cpp241 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 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()
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 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);