R语言stan进行基于贝叶斯推断的回归模型

可以从许多统计软件包中运行Stan。到目前为止，我一直在从R运行Stan。

我们围绕stan进行一些咨询，帮助客户解决独特的业务问题。

简单线性回归

第一步是为Stan模型编写文件。这包含一个文件linreg.stan：

 

 视频：线性回归中的贝叶斯推断与R语言预测工人工资数据案例

，时长09:58

R的Stan

 
data {
  int N;
  vector[N] x;
  vector[N] y;
}
 
model {
  y ~ normal(alpha + beta * x, sigma);
}

该文件的第一部分称为数据，它声明了将作为输入传递给Stan的标量，向量和矩阵。

接下来，我们可以通过运行以下R代码来模拟数据集，并使用Stan和我们的文件linreg.stan来拟合模型：

stan(file = 'linreg. ', data