linear_reg.stan

gistfile1.txt

data {
  int<lower=0> N;          // n data
  vector[N] x;             // regressor
  vector[N] y;             // regressand
}
parameters {
  real alpha;               // intercept
  real beta;                // coefficient
  real<lower=0> sigma;      // variance
}
model {
  y ~ normal(alpha + beta * x, sigma);
}