#' @export
zelig2factor.bayes <- function (
formula,
factors = 2,
burnin = 1000, mcmc = 20000,
verbose=0,
...,
data
) {
loadDependencies("MCMCpack", "coda")
if (missing(verbose))
verbose <- round((mcmc + burnin)/10)
if (factors < 2)
stop("Number of factors needs to be at least 2")
x <- as.matrix(model.response(model.frame(formula, data=data, na.action=NULL)))
list(
.function = "MCMCfactanal",
.hook = "McmcHookFactor",
formula = formula,
x = x,
burnin = burnin,
mcmc = mcmc,
verbose= verbose,
data = data,
factors = factors,
...
)
}
#' @S3method param factor.bayes
param.factor.bayes <- function (...) {
}
#' @S3method param factor.bayes
qi.factor.bayes <- function(obj, x=NULL, x1=NULL, y=NULL, num=1000, param=NULL) {
stop('There is no qi function for the "factor.bayes" model')
list(
"Expected Values: E(Y|X)" = NA
)
}
#' @S3method describe factor.bayes
describe.factor.bayes <- function(...) {
list(
authors = c("Ben Goodrich", "Ying Lu"),
text = "Bayesian Factor Analysis",
year = 2013
)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.