setwd("C:/Users/Lisa/Dropbox/CSUMB stats class")

library("rjags")
library("R2jags")
library("coda")

x <- read.csv(file = "Norm_tst300.csv", header = FALSE)
N <- 300
tmp <- list('x' = x,
            'N' = N)
jags <- jags.model('Normalmodel.bug',
                   data = tmp,
                   n.chains = 4,
                   n.adapt = 100)

update(jags, 1000)

samples <- coda.samples(jags,
             c('mu', 'tau'),
             1000)

plot(samples)
autocorr.plot(samples)

jags2 <- jags.model('Normalmodel2.bug',
                   data = tmp,
                   n.chains = 4,
                   n.adapt = 100)

update(jags2, 1000)

samples2 <- coda.samples(jags2,
                        c('mu', 's2'),
                        1000)

plot(samples2)
autocorr.plot(samples2)

mean(x)
sd(x)