I am trying to fit a simple stan_lm on some simulated data and I’m getting an error that says “please report bug”. Any ideas about what is going wrong?
X = matrix(rnorm(500), nrow = 100)
Y = 3+X%*%c(1,2,3,0.5,0.5)
data = as.data.frame(cbind(Y,X))
model <- stan_lm(V1~V2 + V3+ V4+ V5+V6, data = data, prior = R2(location = 0.2))