回答情况：
提问时间：
问题标签：
我在R中遇到了来自gamlss的一些令人困惑的行为，文档列表predict.gamlss和?predict.gamlss返回RStudio中的函数文档。但是，该函数在键入predict.gamlss并试图运行时不会自动完成，因此将返回Error: 'predict.gamlss' is not an exported object from 'namespace:gamlss'。gamlss版本为5.1.6。
浏览 1修改于2020-05-05得票数 1
回答已采纳
在R中使用GAMLSS包时，有许多不同的方法可以将一个分布适合于一组数据。我的数据是一个单一的值向量，我在这些值上拟合一个分布。我的问题是:使用()和gamlss()之间的主要区别是什么，因为它们对于参数值和不同的蠕虫图给出了相似但不同的答案。此外，使用函数with ()可用于安装了gamlss()的对象，而不适用于安装了fitDist()的对象。是否有任何方法来生成与fitDist()函数匹配的参数的置信区间？
浏览 5提问于2022-01-31得票数 1
data_1 <- read.csv('data_1.csv', sep = ',')mZAGA <- gamlss(cent ~ treat * cat * season, family = ZAGA(), data = data_1)
我想比较一下模型中主要影响和交互的重要性，这似乎不是gamlss中直接提供的特性。本文描述了如何将gamlss对象转换为nlme::lme对象(，第17页)，这似乎是一种很好的方法，因为统计数据包更广泛
浏览 4修改于2021-12-13得票数 -1
实际上，我正在研究GAMLSS模型(用于位置、规模和形状的通用加法模型)。我的问题是:正确地说，这些模型是GLM和线性回归的推广？如果使用带指数族分布的GAMLSS和σ上的同一性链接函数，则得到GLM？作为特例，如果我使用正态分布和两个Idenity链接函数，我得到一个线性回归？
感谢每个人。
浏览 9提问于2022-02-05得票数 0
我一直在研究gamlss软件包来拟合半参数模型，并在ga()函数中发现了一些奇怪的东西。即使模型被指定为具有伽马分布，并使用REML拟合，模型的输出也是高斯的，用GCV拟合。例子：library(gamlss)data(rent)gn3 <- gamlss(R~ga(~s(Fl)+
浏览 1修改于2022-04-13得票数 1
所使用的R代码如下：函数pvc()是变系数函数。当使用包“gamlss”(2020年2月3日)上可用的下列数据集和代码时，也会出现此错误：先考虑一个因素现在考虑到这个因素
r2 <- gamlss(R~pvc(F1，by=
浏览 1提问于2020-05-26得票数 1
我正在使用包gamlss并指定了家庭BEINF。我想知道如何提取$mu.coefficients的p值。0.004517611 26 0.003602046 48 0.002993016 4
> fit.3 = gamlss(y ~ factor(index), family = BEINF, data = mydata)
**********
浏览 1修改于2012-05-20得票数 0
我正在使用gamlss()函数运行一个通用的加法混合模型。我在我的数据上使用了fitDist()，它建议我使用零膨胀的泊松。我的响应变量是'deg‘，是计数数据，但有很多零。fitDEG <- fitDist(deg, data=node_dat, k = 2, type = "counts", try.gamlss = TRUE)
Family:mDEG_zip <- gamlss(formula = deg ~ pb(SE_score) + TL + species
浏览 51修改于2021-06-11得票数 0
我正在尝试使用rpy2在python中重现这个R code： library(gamlss)library(gamlss.add)
fit <- fitDist(x, k = 2, type = "realpl
浏览 9修改于2020-07-15得票数 0
回答已采纳
我想使用gamlss软件包来拟合一个模型，从该包中更多可用的发行版中受益。然而，我很难正确地说明我的随机效应，或者至少我认为这是一个错误，因为如果我将lmer模型的输出与高斯分布进行比较，而gamlss模型与高斯分布的输出不同。如果将无随机效应的lm模型与高斯分布和无随机效应的gamlss模型进行比较，则输出是相似的。df <- na.omit(df)df$Year <- as.fact
浏览 4修改于2021-10-20得票数 1
我的y变量(n=30,000)分布非常重(正和负)，fitDist GAMLSS函数为此选择ST4家族。我试图用一个解释变量x (pb平滑)来评估基于GAMLSS的回归，但是y上的尾太重了，以至于在50个周期后，甚至在改装(time consuming+++)之后，收敛都不会达到。因此，我使用y变换( orderNorm变换，bestNormalize包)进行归一化，使其易于快速收敛，然后从GAMLSS对象中预测拟合值。然而，这些匹配的"orderNormalized“值是一个GAMLSS对象，因此不能使用bes
浏览 3修改于2019-12-10得票数 3
回答已采纳
我已经使用R中的gamlss包运行了带有随机截获的beta回归模型。调用如下所示：
uik.max.model <- gamlss(formula = max.opp.vote ~ n_commute + vote_commute + uik_dummy +
浏览 71提问于2018-01-26得票数 2
回答已采纳
values are higher than 0 x2 = table$T
mod.GA2 <- gamlss(y~x1+x2,sigma.fo=~1,family=GA)
gamlss中的错误(公式=表$
浏览 1提问于2018-08-21得票数 4
最近，我尝试使用gamlss的dredge函数对一个带有dredge的glm进行平均化，得到了与另外两个用户和相同的错误。下面是一个可重复的例子和错误：library(MuMIn)
structure(list(offset_var = c(1.28599085563929model <-gamlss(response_var~pred1*pred2 + pred1*pred3 + pred1*pred4 + pred1*pred5 + pred6 + offset(最新的包文档仍然
浏览 7修改于2022-06-05得票数 1
回答已采纳
对于gamlss函数，我得到了以下错误ERROR : NA's in the working vector or weights for parameter mu：fit <- gamlss::gamlss
浏览 24提问于2020-03-27得票数 1
我想从gamlss R函数预测二进制类概率/类标签，如何利用预测函数得到它们？我有以下示例代码X1 <- rnorm(500)Y <- ifelse(X1>0.2& X2=="A",1,0)
training <- sample(1:n, 400)
浏览 0修改于2018-02-19得票数 3
回答已采纳
我正在努力将gamlss结果收集到一个数据框架中。这将继续示例library(tidyverse)library(gamlss)
mrofit = map(data, ~lm(mpg ~ hp + wt + disp, data = .)),使用gamlss的断例
# GAMLSS model.frame workaround for
浏览 0提问于2018-03-06得票数 0
mtcars %>% model.frame()library(broom) nest(-am) %>%.) +
labs(x = "Miles Per Gallon", y = "Predicted V
浏览 0修改于2018-03-02得票数 3
回答已采纳
我正在使用lms函数在gamlss包中进行厘米回归。运行lms()函数后，得到一个对象，但我不知道如何显示厘米曲线。
如何提取参数，这样我就可以自己绘制出厘米曲线。
浏览 10修改于2013-03-01得票数 0
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号