run-model8-1.R

library(rstan)

d <- read.csv('input/data-salary-2.txt')
N <- nrow(d)
K <- 4
data <- list(N=N, X=d$X, Y=d$Y, KID=d$KID)
fit1 <- stan(file='model/model8-1.stan', data=data, seed=1234)

save.image('output/result-model8-1.RData')