# leggi i dati di input (lista 1D)
# e fissa valori iniziali di mu e sigma
d <- read.table("sample.dat")
media <- mean(d)
std   <- sqrt(mean(d^2)-media^2)
tau <- 1/std^2
init <- sprintf("list(mu=%f,tau=%f);", media,tau)
write (init, file="mu0_tau0.in")