-
0 votesanswersviews
使用jagam在空间样条模型中链接函数问题
我正在尝试基于地理坐标（x和y）在一些二进制数据（y）中建模空间模式 . 最后我将这个代码构建到一个BUGS模型中，所以我在包mgcv中使用jagam来为我开发代码 . 我能够使用logit链接函数（在下面的代码中使用ilogit - mod 1）获得合理的结果（即，预测显示预期的模式） . 但是这个模型似乎不能用我喜欢的链接函数，cloglog，（mod 3）或者我写出长期的logit转换（... -
2 votesanswersviews
如何从ZIP或ZINB模型中获取贝叶斯p值的新样本
希望有人可以帮助我这个，因为我真的卡住了，没有找到我的编码错误！ 我在JAGS中使用零膨胀泊松/负二项式GLM（没有随机效应）（使用R2Jags），并且参数估计，先验，初始值和链收敛都很好 . 所有结果都完全符合例如pscl-package的估计值，包括我在模型中计算的皮尔逊残差... 我唯一无法工作的是从模型中采样一个新样本以获得贝叶斯p值来评估模型拟合 . 我之前拟合的“正常”泊松和负二项模型... -
2 votesanswersviews
用RJAGS约束贝叶斯线性回归中的参数
我正在使用RJAGS构建一个分层贝叶斯线性回归模型，我想约束三个参数值的总和，以正态分布1.3 . 那是： 该模型是： Y = B1 * X1 B2 * X2 B3 * X3 ... BN * XN 和， B1 B2 B3~dnorm（1.3,1 /（0.2）^ 2） 有可能吗？使用一行代码来分配上一行中的参数总和似乎不起作用 . 第二个最好的选择是完全约束参数（B1 B2 B3 = 1.3），但... -
2 votesanswersviews
在ubuntu 14.04中安装rjags
在Ubuntu 14.04中从RStudio安装rjags，我收到错误消息， configure: error: "JAGS module directory /usr/lib/JAGS/modules-3 does not exist." 在usr / lib / JAGS中有模块-4 . 我也无法从控制台安装rjags，因为我收到错误消息， Problem with hea... -
4 votesanswersviews
混合帕累托和普通斯坦模型不起作用
我正试图通过rstan学习斯坦（因为我熟悉R） . 我试过运行一个简单的混合Pareto和Normal模型 . 编译很好（据我所知），但它没有采样，给我错误： “尝试100次后，（ - 2，2）之间的初始化失败 . 请尝试指定初始值，减少约束值的范围或重新参数化模型 . 调用采样器时发生错误;未完成抽样“ 我只想说我已经尝试了各种方法来参数化，并尝试设置初始值，但都无济于事 . 我的Rrstan... -
0 votesanswersviews
R中的残差和包络图
我使用JAGS拟合回归模型，现在我想做一个残差的模拟包络来检查这个模型的拟合 . 假设下面的矩阵是残差矩阵，其中每行是观察，每列是模拟 . 第一列是原始数据集的残差，其他是模拟残差 resid <- matrix(runif(330,0,2),ncol =11) resid <- apply(resid,2,sort) 为了进行模拟包络，我需要找到每个观测值的最小值，最大值和平均值，...