79得票7回答
从aov中提取p值

我想从R的ANOVA中提取生成的p值。 这是我运行的代码:test <- aov(asq[,9] ~ asq[,187]) summary(test) 输出: Df Sum Sq Mean Sq F value Pr(>F) asq[, 1...

38得票1回答
什么时候应该使用aov()函数,什么时候应该使用anova()函数?

我参考了很多在线文献,但它们增加了我的困惑。大部分讨论都过于技术化,涉及到不平衡设计、I、II或III因子方差分析等等术语。 我只知道 aov() 内部使用 lm(),对于具有因子的数据很有用。而 anova() 可以用于同一数据集上的不同模型。我的理解正确吗?

38得票1回答
使用statsmodels或scipy在Python中使用pandas数据框执行ANOVA?

我想使用Pandas DataFrame来分解一个变量的方差。 例如,如果我有一个名为“Degrees”的列,它被按日期、城市和夜间与白天索引,我想找出该系列中多少变化的一部分来自横向城市变化,有多少来自时间序列变化,以及有多少来自夜间与白天。 在Stata中,我会使用固定效应并查看R ...

19得票1回答
Scikit Learn:如何检查系数的显著性

我尝试使用SKLearn对一个包含大约600个虚拟变量和少量区间变量(数据集中有300K行)的大型数据集进行逻辑回归,得到的混淆矩阵看起来很可疑。我想检查返回系数和ANOVA的显著性,但我不知道如何访问它。这真的可能吗?对于包含大量虚拟变量的数据,最佳策略是什么?非常感谢!

17得票2回答
R中的aov()函数错误项:Error(id)和Error(id/timevar)规范有什么区别?

aov(depvar~timevar+Error(id)) 和 aov(depvar~timevar+Error(id/timevar)) 公式规范有何区别?这两种变体会产生略微不同的结果。 同样的问题曾在此处提问过:https://stats.stackexchange.com/quest...

16得票6回答
R错误提示说:“模型未全部拟合到相同大小的数据集上”。

我已经创建了两个广义线性模型,如下所示:glm1 <-glm(Y ~ X1 + X2 + X3, family=binomial(link=logit)) glm2 <-glm(Y ~ X1 + X2, family=binomial(link=logit)) 然后我使用anov...

15得票4回答
如何使用Anova命令(car包)进行Tukey HSD测试

我在处理一个不平衡的设计/样本,并最初学习了 aov()。现在我知道,对于我的ANOVA测试,我需要使用第III型平方和,这需要使用 lm() 进行拟合,而不是使用 aov()。 问题是如何使用 lm() 来获取事后检验(具体来说是 Tukey's HSD)。我所做的所有研究都表明,在 mu...

15得票4回答
在一个集合中检测不同元素的高效算法

想象你有一个包含五个元素(A-E)的集合,每个元素都有一些数值属性(例如“心率”的几次观测结果): A = {100, 110, 120, 130} B = {110, 100, 110, 120, 90} C = { 90, 110, 120, 100} D = {120, 100, 12...

15得票3回答
双因素方差分析的同方差性检验

我一直在使用var.test和bartlett.test来检查基本的ANOVA假设,其中包括同方差性(方差齐性)。对于单因素方差分析,该过程非常简单:bartlett.test(x ~ g) # where x is numeric, and g is a factor var.test(x...

15得票2回答
R中的重复测量/受试者内组ANOVA

我正试图使用R运行重复测量的ANOVA分析。我查看了各种网站上的各种示例,但它们似乎从来没有谈到我遇到的错误。我认为我可能误解了一些重要的东西。 我要运行的ANOVA是基于使用人类参与者的实验数据。它有一个因变量和三个自变量。所有自变量的所有水平都在所有参与者上运行,使其成为三因素重复测量/受...