/**********************************************************************
 * model: 
 *    y ~ normal(mu, sigma) 
 *    mu ~ normal(1.9, 0.95) 
 *    sigma ~ scaled inverse chi-square(1, 0.1) 
 **********************************************************************/

data {
    int<lower=0> n;           // sample size 
    real y[n];                // data - midge wing length
}
parameters { 
    real mu;                  // mean
    real<lower=0> sigma;      // standard deviation 
}
model {
    mu ~ normal(1.9, 0.95);   // 2nd argument is std.dev
    sigma^2 ~ scaled_inv_chi_square(1, 0.1); 
    y ~ normal(mu, sigma);   
}

/**********************************************************************
 * THE END 
 **********************************************************************/