R中的car包实现了带非正的Box-Cox (BCN)。这不是标准的Box-Cox变换，尽管它是在此基础上构建的。 emmeans包可以处理转换后的数据。如果定义了适当的函数，它可以进行反向变换。这是emmeans使用的接口。 make.tran返回一个列表，其中包含变换、其逆、雅可比矩阵和用于测试域的函数。 如何将BCN转换添加到内置转换中？只是简单地创建一个包含适当内容的新列表(就像make.tran指定的那样)并在调用emmeans</e
浏览 46提问于2019-10-12得票数 2
回答已采纳
我在我的数据上尝试了log，sqrt root和arcsine变换，但都不起作用。我尝试使用boxcox，但得到的错误响应变量必须为正 md<-lm(Score1~Location1+Site1+Trial1+Stage1)
plot(md, which
浏览 159提问于2021-09-21得票数 0
R中是否存在等效于幂变换的python函数？这是R：powerTransform(X, family = "bcnPower")中的 from scipy.stats import skew,boxcox_normmaxsklearn.preprocessing import power_transform我试图使用(vec))
print(bo
浏览 0提问于2019-08-27得票数 2
1回答
我可以使用glmer.nb模型和emmeans，使用type = "response"反向转换估计值和置信区间。model = glmer.nb (response ~ p1 + p2 + (1|block))我认为glmmTMB
浏览 19提问于2019-08-02得票数 1
我现在知道，对于我的方差分析测试，我需要使用III型平方和，这涉及到使用lm()而不是aov()进行拟合。mod<-lm(Snavg~StudentEthnicity*StudentGender)
Anova(mod, type="III&q
浏览 1提问于2011-10-12得票数 12
1回答
= new Car(width / 2, 500, -Math.PI / 2 + Math.PI / 2); var flagger =if (car.y < car.height / 2) { }
if (car.x < car.height
浏览 59提问于2019-12-10得票数 5
回答已采纳
1回答
emmeans的输出类型为"response“，提供了我正在比较的两组的平均值和置信区间。然而，我想要的是组之间差异的平均值和CI (即估计值)。emmeans只提供比率(使用type=“响应”)或对数比率(使用type=“链接”)，我不确定如何将其更改为绝对值。如果您在不使用对数(x+1)转换的情况下运行模型，则emmeans提供估计的差异和围绕此差异的CI，而不是比率。当我的响应变量被对数(x+1)转换时，我如何也这样做
浏览 10提问于2019-09-08得票数 0
1回答
为了满足lmer中线性模型的假设，我必须转换变量响应(例如变量1)，使用这里建议的重尾数据的方法，并演示如下：接下来，我更新引用网格，以便按照这里的建议对转换进行解释( )但是，在使用以下命令后，to方法不会被转换回原来的比例：
fitemm <- as.data.frame(<em
浏览 12提问于2022-06-12得票数 1
回答已采纳
我一直在使用lme4中的lmer函数，它给出了独立变量的系数，这样我就可以建立预测模型，但它没有给出每个因素的F统计量。我可以使用lsmeans成对地运行tukey来检查特定的比较，但我想首先证明使用F-statistics运行这些比较是合理的。
浏览 3提问于2017-12-16得票数 0
1回答
在这种模式下：mod <- lm(conc~source + percent, data=pigs)
percentemmean SE df lower.CL upper.CL我正在使用emmeans包，它给出了percent但是，如果我想预测percent值8、10和13怎么办？如何在emmea
浏览 13提问于2022-10-31得票数 1
回答已采纳
我想用离散预测器和不平衡数据来计算glm模型中的边际均值。利用函数emmeans的emmeans包得到边缘均值，给出了设置weights="cell"和weights="proportional"的不同结果。包装文档称，"proportional"使用的权重与平均的因子组合的频率(在原始数据中)成比例，而"cells"则根据被平均的细胞频率使用权重。但我不明白这到底是什么意思？！model <- glm
浏览 5提问于2021-03-22得票数 1
回答已采纳
我希望将带有自定义因子的变量从ANOVA模型赋给emmeans()语句。在这里，我使用来自R的橙子数据集来使代码可复制。这是我的模型，我通常是如何计算因子存储的方法的：oranges$store<-as.factor(oranges$store)means<-emmeans(model, pairwise ~ store, adjust="tuk
浏览 5提问于2021-06-29得票数 2
回答已采纳
1回答
我试图使用R来运行后临时比较之后，一个重要的交互混合方法Anova。我想做的后特别类似SPSS EMMEANS=TABLES(组*时间)比较(组) ADJ(BONFERRONI)，使用估计的边际平均，但不假定方差相等。aov_car(depvar ~ group*timept + Error(id/(timept)), data=myData)
如果我使
浏览 0提问于2018-06-13得票数 0
回答已采纳
但会产生以下错误： UseMethod("cld")中出错:没有适用于"c('emm_list'，'list')“类的对象的”cld“方法” arcdata 代码如下 install.packages("car
浏览 27提问于2019-05-05得票数 2
1回答
我有一个超过100个线性模型的列表，我想为每个模型取估计的平均值和标准误差。mtcars %<>% cyl = as.factor(cyl)使用emmeans，我可以很容易地得到每个模型的估计均值和标准误差。mod <-
浏览 0提问于2018-09-18得票数 3
回答已采纳
1回答
我使用简单的平方变换转换了一个变量(例如leaf_area)，然后对包含交互作用的以下模型进行了拟合：fit.emm <- emmeans(fit, ~ genotype*soil_type + date_measurement, type = "response")
尽管如此，如下面的例子所示，30均值是土壤MT和测量日期2
浏览 6提问于2022-07-11得票数 2
1回答
我在R中使用了更大的函数，我也想使用“hoc”函数来进行后期的分析和绘图。然而，函数不能处理类“大”的对象。有一种方法可以从类“”构造类"glm“的对象library(biglm)data(trees)ff <-log(Volume)~f
a <- glm(ff,data=tree
浏览 1提问于2018-09-09得票数 1
回答已采纳
我在这里从car数据集中介绍了这个示例，并省略了错误项，因为这个示例不需要它：df=Cars93[c(1:15),]
df=Cars93这些是我看过的来源：和。不过，最后一个引用并不报告交互。
浏览 4提问于2020-07-06得票数 1
回答已采纳
我想使用带有lsmeans对象的contrast()来执行一些特定的计划比较，但是我找不到一个方法来执行我想要的比较。我想要比较一个因素是否在两个因素水平中的任何一个水平上产生影响。
浏览 17提问于2019-10-03得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
腾讯会议
云直播
对象存储活动推荐
腾讯云开发者
