Normal Distribution (Code)

########################################################################

# Normal Distribution. Mean: unknown. Variance: unknown.

########################################################################

model

{

for(i in 1:n)

{

time[i] ~ dnorm(mu, inv.sigma2) # likelihood

}

mu ~ dnorm(0, 0.000001) # approximates the improper prior

inv.sigma2 ~ dgamma(0.0001, 0.0001) # approximates the improper prior

sigma2 <- 1/inv.sigma2

sigma <- sqrt(sigma2)

}

# Initial values

list( mu = 0, inv.sigma2 = 1)

# Data 1

list(n = 20)

# Data 2

time[]

182

201

221

234

237

251

261

266

267

273

286

291

292

296

296

296

326

352

359

365

END