dmontecino

## gist:b804853e4b36a57990a7108a35201cf5
### Categorical missing data in Stan ###
library(rstan)


N <- 10000
N_missing <- 100
K <- 3 # number of categories
x <- sample( seq(from=1,to=K), size=N, replace=TRUE ) # unordered cateogrical covariate

# Simulating bivariate respose as function of a categorcial variable
	### Categorical missing data in Stan ###
	library(rstan)


	N <- 10000
	N_missing <- 100
	K <- 3 # number of categories
	x <- sample( seq(from=1,to=K), size=N, replace=TRUE ) # unordered cateogrical covariate

	# Simulating bivariate respose as function of a categorcial variable