我使用来自农产包装的HSD.test进行了后特别测试,但结果是MSD值而不是HSD。在某些情况下,使用aov拟合的方差分析结果显着性(p<0.05),但HSD.test不会输出HSD或MSD值。知道为什么吗拜托?
HSD.test(model, trt = "Groups")
> summary(aov)
Df Sum Sq Mean Sq F value Pr(>F)
Groups 2 29696 14848 36.87 5.3e-11 ***
Residuals 57 22953
问题:我想学习如何改变图基在R中的测试计算平均值和分配相应字母的样本顺序。下面是一个非常简单的例子。
我玩了虹膜数据,发现不同物种的Sepal.Length存在差异。这是盒子的情节:
我进行了一次方差分析,发现差异在统计学上是显著的。
> fit <- lm(Sepal.Length ~ Species, data = iris)
> summary(aov(fit))
Df Sum Sq Mean Sq F value Pr(>F)
Species 2 63.21 31.606 119.3 <2e
在运行Tukey测试时,我遇到了一个奇怪的错误。我希望有人能够帮助我,因为我已经尝试了很多。这是我的数据框架:
Name Score
1 A 2.29
2 B 2.19
这是我的Tukey测试代码:
#TUKEY HSD TEST
tukey = pairwise_tukeyhsd(endog=df['Score'].astype('float'),
groups=df['Name'],
我尝试使用statsmodels模块执行多个样本比较和Tukey HSD,但我一直收到这个错误消息,“当p >= .9时,ValueError:v必须大于1”。我试着在互联网上查找可能的解决方案,但没有效果。任何熟悉这个模块的人都可以帮助我破解我做错了什么来提示这个错误。我使用的是Python版本2.7x和spyder。下面是我的数据和print语句的示例。谢谢!
import numpy as np
from statsmodels.stats.multicomp import (pairwise_tukeyhsd,MultiComparison)
###--- Here are t
我是R的新手,我有一个关于如何将循环中的结果保存到文件中的问题。下面是一个示例:
# Read in the data, set up variables
sdata<-read.csv("sdata.csv", header=T)
m=ncol(sdata)
x=matrix(0,m,4)
row.names(x) <- variable.names(sdata)
colnames(x) <- c( "Ground", "111d", "125d", "Ground")
for (i in
我正在尝试执行Tukey的HSD测试,看看我的数据中的几个组的值平均值是否存在显著差异。例如,在这里,我试图看看变量'acad_se_communicate_needs‘按组'Class’是否存在平均差异。但是,我在结果中遇到了NaN值。这是怎么回事,我该怎么解决呢?
我已经使用状态模型函数来完成这个任务。我避免了需要为每个组将数据分割成不同数据的方法,因为我必须对多个变量执行此分析。而且,这些方法对我来说真的很难理解。
from statsmodels.stats.multicomp import pairwise_tukeyhsd
from statsmodels.stat
我正在寻找一种方法来从Python中的数据帧执行ANOVA和HSD测试。我试着在论坛和教程上阅读一些例子,但我没有实现将其应用到我的工作中。
这是一个简单的Pandas数据帧:
Date Density Hour Repetition Glucose
A HD AM 1 6.7
A HD AM 2 6.8
A HD PM 2 9.6
A HD PM 3 11.9
B HD AM 1 23
B
我关于在Tukey HSD测试后显示同构子集的问题在过去已经被问过(请参阅,在讨论的最后),但从未得到回答。我现在也遇到了同样的问题,我迫切需要一个解决方案,因为这个函数似乎是R唯一可以显示这样的组的函数。如果不是,请告诉我。简而言之,我得到的是:
Groups, Treatments and means a 140095-001 36.79 b 150004-001 32.1 b 136936-021 31.97 bc 137219-004 31.39 bc 136673-017