See sampling, for which this function is
a wrapper, for additional options. See prepData to
create a suitable data list. See
print.stanfit for ways of getting tables
summarizing parameter posteriors.
# NOT RUN {dl <- prepData(phyActFlowPropensity[,c(1,2,3)])
dl$varCorrection <- 2.0
# }# NOT RUN {pcStan('unidim_adapt', data=dl)
# }# NOT RUN {# takes more than 5 seconds# }