# Chapter 4 Methods

We describe our methods in this chapter. Let’s run some code

a=c(1,2,3)
print(a)
## [1] 1 2 3

schools_dat=list(J=8,y = c(28,  8, -3,  7, -1,  1, 18, 12),
sigma = c(15, 10, 16, 11,  9, 11, 10, 18))

fit <- stan(file = './schools.stan', data = schools_dat)