fix normal-mixture_k example

Author: Ben Goodrich

basic_estimators/normal-mixture_k-sim.R

@@ -1,17 +1,3 @@
-K <- 3;
-
-N1 <- 500;
-N2 <- 300;
-N3 <- 200;
-N <- N1 + N2 + N3;
-mu1 <- -2.0;
-mu2 <- 0.5;
-mu3 <- 3.5;
-sigma1 <- 1.0;
-sigma2 <- 0.5;
-sigma3 <- 2;
-
-y <- c(rnorm(N1,mu1,sigma1),rnorm(N2,mu2,sigma2),rnorm(N3,mu3,sigma3));
-
 library("rstan");
-fit <- stan(file="normal_mixture_k_prop.stan", data=c("K","N","y"), iter=1000, chains=1, init=0);
+source("normal-mixture_k.data.R")
+fit <- stan(file="normal_mixture_k.stan", data=c("K","N","y"), iter=1000, chains=1, init=0);

basic_estimators/normal-mixture_k.data.R

@@ -0,0 +1,14 @@
+K <- 3;
+
+N1 <- 500;
+N2 <- 300;
+N3 <- 200;
+N <- N1 + N2 + N3;
+mu1 <- -2.0;
+mu2 <- 0.5;
+mu3 <- 3.5;
+sigma1 <- 1.0;
+sigma2 <- 0.5;
+sigma3 <- 2;
+
+y <- c(rnorm(N1,mu1,sigma1),rnorm(N2,mu2,sigma2),rnorm(N3,mu3,sigma3));