PFloat Verify::mult_weights(Eigen::ArrayXd w)
{
PFloat ret = PFloat(w[0]);
for(int i = 1; i < w.size(); i++)
{
PFloat tmp = PFloat(w(i));
ret = ret*tmp;
}
return ret;
}
GLMMBelief::GLMMBelief(const std::vector<int>& items,
const Eigen::MatrixXd& X,
const Eigen::MatrixXd& Zt,
const Eigen::SparseMatrix<double>& Lambdat,
const Eigen::VectorXi& Lind,
const Eigen::ArrayXd& response,
const Eigen::ArrayXd& weights):
ContinuousBeliefBase(items),
numObservations_(response.size()), numFixed_(X.cols()),
X_(X), Zt_(Zt), Lambdat_(Lambdat), Lind_(Lind),
response_(response), weights_(weights)
{
initializeParameterDependents();
}