Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/wpf/13.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
如何使用rstanarm设置不同的优先级_R_Stan_Rstan_Rstanarm - Fatal编程技术网

如何使用rstanarm设置不同的优先级

如何使用rstanarm设置不同的优先级,r,stan,rstan,rstanarm,R,Stan,Rstan,Rstanarm,假设我有一个形式为y=a{i}+b{i,1}*x{1}+b{2}*x{2},其中i=1,2,…,12,我想用rstanarm来估计这个模型 可以为每个截距设置不同的优先级吗?(假设前4个具有normal(位置=0,刻度=1,自动刻度=TRUE),后4个具有normal(位置=1,刻度=2,自动刻度=TRUE),最后4个具有student(df=1,位置=0,刻度=NULL,自动刻度=TRUE))。我还想为b{I,1}和最后的b{2}~normal(location=3,scale=1,autos

假设我有一个形式为
y=a{i}+b{i,1}*x{1}+b{2}*x{2}
,其中
i=1,2,…,12
,我想用
rstanarm
来估计这个模型

可以为每个截距设置不同的优先级吗?(假设前4个具有
normal(位置=0,刻度=1,自动刻度=TRUE)
,后4个具有
normal(位置=1,刻度=2,自动刻度=TRUE)
,最后4个具有
student(df=1,位置=0,刻度=NULL,自动刻度=TRUE)
)。我还想为
b{I,1}
和最后的
b{2}~normal(location=3,scale=1,autoscale=TRUE)设置相同的优先级


rstanarm是否有可能做到这一点?

rstanarm支持的模型中最多有一个截距,但您可以通过在公式中包含-1并将虚拟变量上的系数视为系数来抑制截距。对于系数,可以执行以下操作

previor=student_t(df=c(rep(Inf,8),rep(1,4)),
位置=c(代表(0,4),代表(1,4),代表(0,4)),
比例=c(代表(1,4),代表(2,4),代表(1,4)),
自动缩放=真)

但是,似乎您打算建立某种层次模型,在这种情况下,全局参数偏差的先验值只能是多元正态的。请参阅
?装饰前

谢谢您的回答。1) 如果我理解正确,添加“-1”也不会导致多重共线性,对吗?只是解释发生了变化。2) 在我的例子中,我想说的是,如果你可以用位置和比例参数代替向量,你也可以传递一个具有不同分布的向量(如果我错了,请纠正我,但在你的回答中,你可以通过操纵df将t分布转换为正态)是,在R公式中使用
-1
+0
会删除截距,因此它可以包含一个更高级别的因子作为伪变量。否,在rstanarm包中,您不能将函数列表或字符向量传递给
Previor
参数。只能将
位置
比例
测向
等向量传递给一个函数,如
学生
。但是,由于当自由度为无穷大时,
student\u t
相当于
normal
,因此这相当于在您最初给出的示例中使用“不同”函数。使用brms R包中的
brm
功能,您可以为Stan估计的模型中的不同参数指定不同的先验族。在rstanarm中,您不能这样做。但是如果你要用二项数据来计算β先验分布,那么你可以用解析的方法来计算后验分布。嗨,本,我在这里简单地问一些问题,如果你能回答的话,我会正式地问。我有一个前后控制设计。也就是说,对对照组(c)进行预测试(c前)和后测试(c后)。同样,对治疗组(t)进行预测试(前测试)和后测试(后测试)。所以我们在两个时间点(时间因子)测试了两组(组因子)。问题:我想知道使用
rstanarm
我们可以拟合两个因素(即时间和组)的多水平回归,并且可以估计它们的相互作用?