TukeyHSD在R中出现错误
我正在处理混合设计的ANOVA，并想运行TukeyHSD进行事后检验。 我不断收到错误提示： "Error in UseMethod("TukeyHSD") : no applicable method for 'TukeyHSD' applied to an object of cla...
嵌套ANOVA独特因子水平
我正在运行一个嵌套ANOVA分析，并具有以下设置：2个区域，一个是参考区域，一个是暴露区域（列名为CI = Control/Impact）。两个时间期间（影响前和影响后，列名为BA），在之前的时期有1年，在之后的时期有3年。这些年份是嵌套的。 我的问题是：如果我使用原始年份（在玩具数据集中的...
带有区组设计和重复测量的ANOVA
我试图对同一生长季内在两个场地进行的田间试验进行一些统计分析。 在两个场地（Site，水平：HF | NW）上，实验设计为RCBD，每个Site内有4个块（Block，水平：1 | 2 | 3 | 4）。共有4种处理方法-3种不同形式的氮肥和一种对照（无氮肥）（Treatment，水平：AN...
R中的aov()函数错误项：Error(id)和Error(id/timevar)规范有什么区别？
aov(depvar~timevar+Error(id)) 和 aov(depvar~timevar+Error(id/timevar)) 公式规范有何区别？这两种变体会产生略微不同的结果。 同样的问题曾在此处提问过：https://stats.stackexchange.com/quest...
ANOVA: 自由度几乎全部相等为1
我有一个数据集，起始如下： > d.weight R N P C D.weight 1 1 0 0 GO 45.3 2 2 0 0 GO 34.0 3 3 0 0 GO 19.1 4 4 0 0 ...
如何使用plm包比较R中的2个模型？
我正在使用R中的plm包运行一个固定效应模型，我想知道如何比较哪个模型更适合。 例如，这是我构建的两个模型的代码： library(plm) eurofix <- plm(rlogmod ~ db+gdp+logvix+gb+i+logtdo+fx+ld+euro+core, ...
如何在R中使用anova()函数比较两个glm模型时提取p值
因此，我正在尝试比较两个模型fit1和fit2。 最初，我只是做了anova（fit1，fit2），这产生了我理解的输出（包括p值）。 然而，当我将我的模型从基于lm（）的模型切换到基于glm（）的模型时，anova（fit1，fit2）现在产生了剩余自由度，剩余偏差和Df偏差，我很难解释...
Scikit Learn：如何检查系数的显著性
我尝试使用SKLearn对一个包含大约600个虚拟变量和少量区间变量（数据集中有300K行）的大型数据集进行逻辑回归，得到的混淆矩阵看起来很可疑。我想检查返回系数和ANOVA的显著性，但我不知道如何访问它。这真的可能吗？对于包含大量虚拟变量的数据，最佳策略是什么？非常感谢！
在R中自定义对比：对比系数矩阵或对比矩阵/编码方案？如何实现？
自定义对比在分析中非常广泛使用，例如："这个三级因子的第1级和第3级的DV值是否显著不同？" 直观地，这种对比是以单元均值为基础表达的： c(1,0,-1) 一种或多种这些对比项，以列的形式绑定，组成一个对比系数矩阵，例如： mat = matrix(ncol = 2, byrow ...