Home
last modified time | relevance | path

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

/dports/math/R-cran-MCMCpack/MCMCpack/src/
H A DcMCMChregress.cc136 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 DcHMMpanelFE.cc195 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 DHMMmultivariateGaussian.cc102 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 DcHMMpanelRE.cc137 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 …]