Dear Bob,
I am migrating our recent discussion from Google groups. Thanks for pointing me here.
Indeed, I am trying to simulate the design matrix using a distribution. Currently, I am using this R code:
## Generate response
response <- rep(x = c(1, 0), times = 25)
## Generate covariates conditioned on response
generate_covariates <- function(response_variable){
## covariate_i: Bernoulli random variable drawn with probability ranging from 0 to 1
prob_given_1 <- runif(n = 1, min = 0, max = 1)
prob_given_0 <- runif(n = 1, min = 0, max = 1)
## generate class-conditional distribution
covar_given_1 <- rbinom(n = length(which(response == 1)), size = 1, prob = prob_given_1)
covar_given_0 <- rbinom(n = length(which(response == 0)), size = 1, prob = 1 - prob_given_0)
return(c(covar_given_1, covar_given_0))
}
covariate_list <- sapply(c(1:10), function(x){
covars <- generate_covariates(response_variable = response)
return(covars)
}, simplify = TRUE)
## Put together simulated response and covariates
simulated_data <- data.frame(response, covariate_list, stringsAsFactors = FALSE)
Once the simulated_data is generated, I run a logistic regression on it using Stan.
My question is can I do the simulation process above in Stan itself?
Thanks for any insights!
Best,
Jy
On May 8, 2017, at 6:57 PM, Bob Carpenter carp@alias-i.com wrote:
You are never going to be able to simulate the design matrix
without giving it a distribution—if it’s just covariates,
there’s nothing to simulate from. Usually the covariates don’t
get modeled because everything’s conditionally independent given
their value.
Which interface are you using?
- Bob
P.S. We’re shutting down this list soon. We’ve switched to:
discourse.mc-stan.org
On May 7, 2017, at 8:48 PM, Jy jyoti.shankar@gmail.com wrote:
Dear Bob,
Is this possible now in Stan - i.e. a way to simulate the design matrix that will then be provided > as input to stan? If yes, would you please point me to documentation of how to do this?
Thanks!
JyOn Saturday, 17 October 2015 17:46:12 UTC-4, Bob Carpenter wrote:
What I’m talking about here is a way to generate a single fake data
set to use as input to Stan to make sure my program is doing the
right thing. Right now, there’s no easy way (that I know of) to convert
Stan output (in any interface) into something that can be fed back in directly
as data.
- Bob