比较各组之间的基尼指数和希尔数
还 发布在 CrossValidated 但我想我可能会得到这里有更多的牵引力。 我正在寻找一个 R 软件包,可以测试两组之间的 Gini 指数和 Hill 数字的差异。 R …
glm 模型预测的置信区间
我正在尝试从 glm 获取预测的 95% 置信区间。下面是我的模型 library(caret) library(ISLR) data(Smarket) glm.fit <- glm(Direction ~ Lag1 + Lag2, …
灵敏度和特异性最大总和的最佳截止点
我想计算我的 glm 模型的敏感性特异性总和最大化阈值(约登指数): model = glm(present ~ Summer_precipitation + Summer_Temperature + Frost_days…
如何让 My.stepwise.glm 在控制台之外返回模型?
我在 RCommunity 上问了这个问题,但还没有人咬过......所以我在这里! 我当前的项目涉及预测在未来气候变化情景下某些树木是否能够生存。为了更好地…
将 glm 的所有输出呈现在一个漂亮的表格中以进行比较
我已经使用了 @StupidWold 此处 开发一个glm 并将结果存储在models 中。输出似乎是正确的。但是我想知道如何一次呈现所有输出,而不是单独调用每个输…
ggeffect() 不返回所有预测
使用安全带(包含在 R 中)数据,我想要预测“前”的边际效应。该变量有 170 个值。 ggeffect() 仅返回值 400、500、600、...、1300 处的 front 预测…
在 R 中运行 glm 时出现错误/警告 - 系数:(1 由于奇点而未定义)
当我在 model_df 上运行此 glm 时(下面的 dput): model <- glm(outcome_variable_data ~ ., data = model_df, family = 'binomial') lrm <-summary…
在 for 循环中创建多个 GLM,跳过循环中系数在 R 中不起作用的模型
我创建了一个列表,其中包含因变量的所有可能组合,并且我正在尝试使用所有这些组合创建多个 glm。 combinations = list() models = list() for(i in …
逻辑回归 - glm 中的 cbind 命令
我正在 R 中进行逻辑回归。有人可以澄清运行这两条线有什么区别吗? 1. glm(Response ~ Temperature, data=temp, family = binomial(link="logit")) 2…
是否可以拟合仅具有响应变量的线性模型?
如果这样做,我会得到两个系数(截距和年份), data <- data.frame(accidents=c(3,1,5,0,2,3,4), year=1:7) glm(accidents ~ year, family=poisson(l…
如何通过观察提取 lmer 固定效应?
我有一个 lme 对象,它是根据一些重复测量的营养摄入数据(每个 RespondentID 两个 24 小时摄入期)构建的: Male.lme2 <- lmer(BoxCoxXY ~ -1 + AgeF…
nlme 错误“组的公式无效”尽管指定了随机效应
我对此进行了一些搜索,但我发现的邮件列表帖子与未在 nlme 中指定随机效果的人相关,而我已经这样做了。我还拥有 Pinheiro 和 Bates 撰写的《S 和 S-…