/* * weibull.stan * fit tree diameter distribution to Weibull distriution */ data { int N; // number of trees vector[N] D; // tree diameters } parameters { real alpha; // shape parameter real sigma; // scale parameter } model { D ~ weibull(alpha, sigma); alpha ~ normal(0, 10); sigma ~ normal(0, 40); }