2回答
我已经使用R中的gamlss包运行了带有随机截获的beta回归模型。调用如下所示： family = BE(),我想从模型对象中提取每个分组因子(mf)的<e
浏览 71提问于2018-01-26得票数 2
回答已采纳
1回答
ggplot2中的蠕虫图残差图
、、、、
我试图在使用来自gamlss包的gamlss函数的模型上绘制蠕虫绘图残差。兴趣图如下所示：
最初，下面是引用来自wormplot_gg包的childsds函数的计算例程，但是，使用上面描述的函数表示的结果不像上面所示的示例，该示例正在应用于R中包含的数据集。library(ggplot2)library(childsds)
Dados <- Orange
Model <- gamlss(circumference~age, f
浏览 12提问于2021-06-06得票数 1
回答已采纳
1回答
我有一个连续的反应，三个范畴变量，没有随机效应。我用零调整的伽马分布拟合了一个伽马物体。链接到数据，因为它很大。, family = ZAGA(), data = data_1)
我想比较一下模型中主要影响和交互的重要性，这似乎不是gamlss中直接提供的特性。本文描述了如何将gamlss对象转换为nlme::lme对象(，第17页)，这似乎是一种很好的方法，因为统计数据包更广泛地接受lme模型。一旦混合GAMLSS被拟合并存储在一个R对象中<
浏览 4提问于2021-12-13得票数 -1
1回答
R中gamlss的随机效应规范
、、、、
我想使用gamlss软件包来拟合一个模型，从该包中更多可用的发行版中受益。然而，我很难正确地说明我的随机效应，或者至少我认为这是一个错误，因为如果我将lmer模型的输出与高斯分布进行比较，而gamlss模型与高斯分布的输出不同。如果将无随机效应的lm模型与高斯分布和无随机效应的gamlss模型进行比较，则输出是相似的。summary(fourC)
plot(four
浏览 4提问于2021-10-20得票数 1
我的y变量(n=30,000)分布非常重(正和负)，fitDist GAMLSS函数为此选择ST4家族。我试图用一个解释变量x (pb平滑)来评估基于GAMLSS的回归，但是y上的尾太重了，以至于在50个周期后，甚至在改装(time consuming+++)之后，收敛都不会达到。因此，我使用y变换( orderNorm变换，bestNormalize包)进行归一化，使其易于快速收敛，然后从GAMLSS对象中预测拟合值。然而，这些匹配的"orderNormalized“值是一个
浏览 3提问于2019-12-09得票数 3
回答已采纳
我正在使用lms函数在gamlss包中进行厘米回归。运行lms()函数后，得到一个对象，但我不知道如何显示厘米曲线。
如何提取参数，这样我就可以自己绘制出厘米曲线。
浏览 10提问于2013-03-01得票数 0
1回答
然后提取每个速度的斜率并截取它们，并将它们分组到更高的分类级别，然后通过Kullback-Leibler散度和我在引导Y值时得到的分布比较斜率的分布(截取)。也不知道如何回答我的分类水平解释了多少方差的问题，就像嵌套随机混合效应模型一样。library(gamlss)
model<-gamlss(Y~temp*fac+re(random=~1|ph
浏览 2提问于2021-12-21得票数 0
回答已采纳
4回答
我想从R中的一个定制分布中提取1000个样本。我有以下自定义发行版mu <- 1 tau <- 3 rate <- 1plot(points ~ x)我如何通过蒙特卡罗模拟从这个分布随机抽样
浏览 0提问于2018-03-20得票数 4
回答已采纳
1回答
., .id = 'shuffle')
ddd 此部分仅提取固定效果。(summary(.x)$Solutions, rownames = "model_term") 但请注意，我没有运行没有任何固定效果的模型，因此输出为空。如何使用相同或相似的代码提取随机效果？我想我可以将“解决方案”更改为其他东西，以便从我运行的没有任何固定效果的模型中提取随机效果。
浏览 37提问于2020-10-28得票数 1
我正在尝试从随机效果表中提取单个元素，该表包含在由混合效果模型的摘要调用创建的对象中。具体地说，我想提取每个2级随机效果。0.74498367 0.03515508 0.76672652 1.91896578
Number of Groups: 8
现在我可以从上面的随机效果表中提取残差modelL
浏览 1提问于2016-05-02得票数 0
0回答
我正在寻找一个类似于在nlme，lme4和brms中使用的ranef()的命令，它将允许我在我的MCMCglmm模型中提取单独的随机效果。在我的数据集中，我有40个提供者，我想为每个提供者提取随机效果，并将它们绘制在毛虫图中。任何建议都是很棒的。谢谢。
浏览 11提问于2017-12-02得票数 5
我可以使用summary(fm1)从nlme摘要中提取固定效果。但是正在苦苦挣扎着如何获得Random effects:部分。
浏览 2提问于2012-01-28得票数 8
回答已采纳
7回答
我有一个具有固定和随机效果的mer对象。如何提取随机效果的方差估计？这是我的问题的一个简化版本。study <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy)这会产生一个很长的输出--在本例中不会太长。我还检查了lme4包中的任何提取器函数，但都无济于事。请帮帮我！
浏览 1提问于2011-12-16得票数 46
回答已采纳
我已经成功地拟合了一个线性混合效果模型，我正在寻找提取单个组的随机效果分量。我知道可以使用以下命令提取随机效果的完整列表然后print(random.effects(模型))给出一个包含两列的组名称和随机效果列表，即使数据本身似乎只有一列我的问题是，是否可以通过组名“查找”与特定组相关联的随机效果，或者，如果不能，我如何才能找到与random.effec
浏览 2提问于2010-09-10得票数 2
回答已采纳
在我创建的函数中，使用R中的tryCatch()函数有问题。我想做的是：
将这些预测存储在数据框架中。我的目标是在我的模拟函数中编写一个tryCatch语句，并且在发生错误的情况下简单地第二次运行模拟/预测代码。我是在代码中编写的，我认为错误是(使用未找到mod_sim对象的预测函数)。它很可能在那里，因为cat刚好在
浏览 3提问于2020-11-01得票数 0
回答已采纳
2回答
线性混合效果模型传统上是以以下方式表示的。Ri = Xi×β+ Zi×bi +εi，其中β表示估计的固定效果，Z表示随机效果。因此，X是经典的设计矩阵。使用R，我希望能够在使用nlme包中的lme拟合模型后提取这两个矩阵。例如，数据集" rails "，也可以在nlme包中找到，包含6个随机选择的铁路轨道上的超声波旅行时间的三个单独测量。我可以拟合一个简单的模型，具有截距固定效果和每个轨道的随机效果</em
浏览 0提问于2013-11-01得票数 6
我正在尝试从我的线性混合效果模型中提取每个组的斜率和截距。该模型是在lme4库中使用lmer构建的，我可以使用jtools库中的interact_plot查看每组的结果。我知道我可以使用summary()或summ()来查看固定效果的估计值和随机效果的方差，但我看不到随机效果的估计值。因此，我无法准确计算模型的斜率和截距。> interact_plot(cond_ranin_waterxsilver, pred = LnVolu
浏览 35提问于2019-02-05得票数 1
回答已采纳
我正在使用lmer来拟合一个具有几个固定效应(包括特定于对象的变量，如年龄、短期记忆广度等)的多级多项式回归模型。和两组随机效果(主题和主题:条件)。现在，我想预测具有特定属性(年龄、短期记忆广度等)的假设对象的数据。("predict") : 我知道我可以使用暴力方法从我的模型中提取</em
浏览 2提问于2012-02-28得票数 2
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号