Searched refs:Yk_arr (Results 1 – 4 of 4) sorted by relevance
/dports/math/R-cran-MCMCpack/MCMCpack/src/ |
H A D | cMCMChregress.cc | 136 Matrix<double> *Yk_arr = new Matrix<double>[NGROUP]; in cMCMChregress() local 142 Yk_arr[k] = Matrix<double>(nobsk[k],1); in cMCMChregress() 150 Yk_arr[k](m,0)=Y_vect[posk_arr[k][m]]; in cMCMChregress() 167 tXYk_arr[k] = t(Xk_arr[k])*Yk_arr[k]; in cMCMChregress() 168 tWYk_arr[k] = t(Wk_arr[k])*Yk_arr[k]; 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() 363 delete[] Yk_arr; in cMCMChregress()
|
H A D | cHMMpanelFE.cc | 195 Matrix<double> *Yk_arr = new Matrix<double>[nsubj]; in HMMpanelFE_impl() local 199 Yk_arr[k] = Matrix<double>(nobsk[k], 1); in HMMpanelFE_impl() 204 Yk_arr[k](l,0) = Y[posk_arr[k][l]]; in HMMpanelFE_impl() 219 newY[s] = Yk_arr[s] - delta0; in HMMpanelFE_impl() 220 Yres[s] = Yk_arr[s]; in HMMpanelFE_impl() 233 Yres[s] = Yk_arr[s] - Xk_arr[s]*beta; in HMMpanelFE_impl() 240 newY[s] = Yk_arr[s] - Zk_arr*delta_holder[s][0]; in HMMpanelFE_impl() 279 Matrix<> Yj = Yk_arr[s]((delta_count - nstate[s][j]), 0, (delta_count - 1), 0); in HMMpanelFE_impl() 365 delete [] Yk_arr; in HMMpanelFE_impl()
|
H A D | HMMmultivariateGaussian.cc | 102 Matrix<double> *Yk_arr = new Matrix<double>[nsubj]; in MultivariateGaussian_impl() local 106 Yk_arr[k] = Matrix<double>(nobsk[k], 1); in MultivariateGaussian_impl() 111 Yk_arr[k](l,0) = Y[posk_arr[k][l]]; in MultivariateGaussian_impl() 119 tXYk_arr[k] = t(Xk_arr[k])*Yk_arr[k]; in MultivariateGaussian_impl() 158 Matrix<> e = t(Yk_arr[s]-Xk_arr[s]*beta)*(Yk_arr[s] - Xk_arr[s]*beta); in MultivariateGaussian_impl() 230 Matrix<> e = t(Yk_arr[s]-Xk_arr[s]*beta_st)*(Yk_arr[s] - Xk_arr[s]*beta_st); in MultivariateGaussian_impl() 248 loglike += lndmvn(Yk_arr[s], Mu, Sigma); in MultivariateGaussian_impl() 282 delete [] Yk_arr; in MultivariateGaussian_impl() 347 Yk_arr[k] = Matrix<double>(nobsk[k], 1); in HMMmultivariateGaussian_impl() 352 Yk_arr[k](l,0) = Y[posk_arr[k][l]]; in HMMmultivariateGaussian_impl() [all …]
|
H A D | cHMMpanelRE.cc | 137 Matrix<double> *Yk_arr = new Matrix<double>[nsubj]; in GaussianPanelRE_impl() local 143 Yk_arr[k] = Matrix<double>(nobsk[k], 1); in GaussianPanelRE_impl() 152 Yk_arr[k](l,0) = Y[posk_arr[k][l]]; in GaussianPanelRE_impl() 166 tXYk_arr[k] = t(Xk_arr[k])*Yk_arr[k]; in GaussianPanelRE_impl() 170 tWYk_arr[k] = t(Wk_arr[k])*Yk_arr[k]; in GaussianPanelRE_impl() 216 (Yk_arr[s] - Xk_arr[s]*beta - Wk_arr[s]*bi(_,s)); in GaussianPanelRE_impl() 383 loglike += lndmvn(Yk_arr[s], Mu, Sigma); in GaussianPanelRE_impl() 415 delete [] Yk_arr; in GaussianPanelRE_impl() 508 Yk_arr[k] = Matrix<double>(nobsk[k], 1); in HMMpanelRE_impl() 516 Yk_arr[k](l,0) = Y[posk_arr[k][l]]; in HMMpanelRE_impl() [all …]
|