powered by
getVarbayesResp(data, model)
# get resp for only a subsample, as this operation is rather long. temp <- getVarbayesResp(pendat[1:10,], vbpen[[2]])
Run the code above in your browser using DataLab