Searched refs:beta_run (Results 1 – 3 of 3) sorted by relevance
/dports/math/R-cran-MCMCpack/MCMCpack/src/ |
H A D | cMCMChregress.cc | 186 Matrix<double> beta_run(NP,1,false); // Unicolumn matrix of fixed effects in cMCMChregress() local 188 beta_run(p)=beta_vect[p*NSAMP]; in cMCMChregress() 236 beta_run=myrng.rmvnorm(big_V*small_v,big_V); in cMCMChregress() 249 Matrix<double> small_vk=(t(Wk_arr[k])*(Yk_arr[k]-Xk_arr[k]*beta_run))/V_run; in cMCMChregress() 273 e+=crossprod(Yk_arr[k]-Xk_arr[k]*beta_run-Wk_arr[k]*bk_run[k]); in cMCMChregress() 293 Matrix<double> Y_hat=Xk_arr[k](m,_)*beta_run+Wk_arr[k](m,_)*bk_run[k]; in cMCMChregress() 308 beta_vect[p*NSAMP+(isamp-1)]=beta_run(p); in cMCMChregress() 327 Matrix<double> Y_hat=Xk_arr[k](m,_)*beta_run+Wk_arr[k](m,_)*bk_run[k]; in cMCMChregress()
|
H A D | cMCMChlogit.cc | 191 Matrix<double> beta_run(NP,1,false); // Unicolumn matrix of fixed effects in cMCMChlogit() local 193 beta_run(p)=beta_vect[p*NSAMP]; in cMCMChlogit() 253 Matrix<double> logit_theta_hat=Xk_arr[k](m,_)*beta_run+Wk_arr[k](m,_)*bk_run[k]; in cMCMChlogit() 299 beta_run=myrng.rmvnorm(big_V*small_v,big_V); in cMCMChlogit() 312 Matrix<double> small_vk=(t(Wk_arr[k])*(logit_thetak_run[k]-Xk_arr[k]*beta_run))/V_run; in cMCMChlogit() 336 e+=crossprod(logit_thetak_run[k]-Xk_arr[k]*beta_run-Wk_arr[k]*bk_run[k]); in cMCMChlogit() 377 beta_vect[p*NSAMP+(isamp-1)]=beta_run(p); in cMCMChlogit() 395 Matrix<double> logit_theta_hat=Xk_arr[k](m,_)*beta_run+Wk_arr[k](m,_)*bk_run[k]; in cMCMChlogit()
|
H A D | cMCMChpoisson.cc | 180 Matrix<double> beta_run(NP,1,false); // Unicolumn matrix of fixed effects in cMCMChpoisson() local 182 beta_run(p)=beta_vect[p*NSAMP]; in cMCMChpoisson() 242 Matrix<double> log_lambda_hat=Xk_arr[k](m,_)*beta_run+Wk_arr[k](m,_)*bk_run[k]; in cMCMChpoisson() 288 beta_run=myrng.rmvnorm(big_V*small_v,big_V); in cMCMChpoisson() 301 Matrix<double> small_vk=(t(Wk_arr[k])*(log_lambdak_run[k]-Xk_arr[k]*beta_run))/V_run; in cMCMChpoisson() 325 e+=crossprod(log_lambdak_run[k]-Xk_arr[k]*beta_run-Wk_arr[k]*bk_run[k]); in cMCMChpoisson() 366 beta_vect[p*NSAMP+(isamp-1)]=beta_run(p); in cMCMChpoisson() 383 Matrix<double> log_lambda_hat=Xk_arr[k](m,_)*beta_run+Wk_arr[k](m,_)*bk_run[k]; in cMCMChpoisson()
|