Very simple question on creating an intercept-only model: How do I do it? I would have thought my code below would work, but it is returning an error.
> InterceptModel <- stan_glm( + Accept_Reject ~ 1, + data = pubdata, + family = binomial(link = "logit"), + prior_intercept = NULL, + cores=3, + QR = FALSE, + chains = 3, iter = 50000, + diagnostic_file=file.path(tempdir(), "df.csv"))
Error in read.table(file = file, header = header, sep = sep, quote = quote, :
no lines available in input
Error in dat[, -1, drop = FALSE] : incorrect number of dimensions