更正Python中Means - Tukey HSD的多重比较
在Python中，Mean - Tukey HSD是指一种多重比较方法，用于在一组数据中比较多个组之间的平均值差异。它是Tukey的Honest Significant Difference (HSD)方法的一种变体，用于确定哪些组之间存在显著差异。
具体步骤如下：
- 首先，收集一组数据，并将其分为多个组。
- 计算每个组的平均值和方差。
- 使用Tukey HSD方法来计算每对组之间的差异是否显著。该方法使用方差分析来确定平均值之间的差异是否超过了预设的显著性水平。
- 如果两个组的平均值之差超过了显著性水平，则认为它们之间存在显著差异。
Mean - Tukey HSD方法的优势包括：
- 它能够处理多个组之间的比较，并提供准确的差异结果。
- 它可以控制整体类型I错误率（即错误地拒绝了无效的假设）。
- 它考虑了样本量和组间差异，因此比其他方法更精确。
应用场景： Mean - Tukey HSD方法适用于那些需要比较多个组之间平均值差异的场景，例如实验结果的分析、市场调研数据的比较等。
推荐的腾讯云相关产品和产品介绍链接地址： 腾讯云提供了多个与数据分析和云计算相关的产品，如腾讯云数据分析平台（https://cloud.tencent.com/product/da），腾讯云云服务器（https://cloud.tencent.com/product/cvm）等。您可以根据实际需求选择适合的产品来支持您的数据分析和云计算工作。
相关·内容
我正在处理一个不平衡的设计/样本，并且最初学习了aov()。我现在知道，对于我的方差分析测试，我需要使用III型平方和，这涉及到使用lm()而不是aov()进行拟合。
问题是使用lm()进行后期测试(特别是Tukey的HSD)。我所做的所有研究都表明，在multcomp包中使用simint是可行的，但现在它已经更新了，该命令似乎不可用。它似乎也依赖于通过aov()进行计算。
基本上，我找到的所有针对R的Tukey HSD测试都假定您使用aov()而不是lm()进行比较。为了得到不平衡设计所需的类型III平方和，我必须使用：
mod<-lm(Snavg~StudentEthnicity*S
浏览 1提问于2011-10-12得票数 12
我有一个包含4个条件(A，B，C，D)的日期集。我观察到运行单向方差分析的结果是，在4种情况下，我的因变量(反应时间，RT)呈线性增长。 我想运行一个后自组织测试，看看在Tukey HSD后自组织测试中，RT从A到B，从B到C，以及从C到D的增加是否显著。 为了在Python中运行测试，我使用了以下代码： #Multiple Comparison of Means - Tukey HSD
from statsmodels.stats.multicomp import pairwise_tukeyhsd
print(pairwise_tukeyhsd(df["RT"], df[
浏览 104提问于2020-11-10得票数 1
回答已采纳
我一直在环顾四周，我对图基调整的方式感到很困惑。根本没有提到Tukey，但是在中有这样的说法：“对于大多数对比度()结果来说，调整通常是另一回事，这取决于创建了哪种类型的对比。例如，成对比较默认地调整= " Tukey "，即Tukey HSD方法。”
据我所知，Tukey HSD本质上是一系列配对t测试，并对I型错误进行调整。但是emmeans函数是计算估计的边际平均(EMMs)，我假设它不是成对的t检验；然后将Tukey调整应用于EMMs输出，将不等同于Tukey HSD后特别测试。
第二个相关的问题是"tukey.emmc“函数是什么？
我想我的第二个问题是，tu
浏览 5提问于2020-11-25得票数 2
2回答
我的问题是，我如何在类型2的Anova (混合效果模型)上运行帖子？到目前为止，我使用的是"lme4“包中的glmer()，"car”包中的Anova()，并尝试运行"agricolae“包中的HSD测试。
在搜索了一段时间后，这是我能找到的最好的结果，但是，在这样做的时候，我收到了一条错误消息。有没有人知道如何解决这个问题，或者我做错了什么？或者是一种不同的方式？
library(lme4)
totaldiversity.model <- glmer(totaldiversity ~ focalspecies + (1|site), family = "
浏览 107提问于2018-06-20得票数 1
回答已采纳
1回答
我使用来自农产包装的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
浏览 9提问于2022-10-15得票数 0
我有一组数据，我已经用一个广义线性模型进行了分析，该模型在3向交互中有三个分类因素(factorA，factorB，factorC)，第四个连续因素(factorD)简单地添加到模型中。我试图从模型中获得一组Tukey字母组(即紧凑的字母显示)，但还没有找到一种方法来成功地包含交互。我对包含factorD不感兴趣，只对交互中的三个感兴趣。
我已经得到了Tukey调整后的成对比较：
lsmeans(my.glm, factorA*factorB*factorC)
但我不能想出如何从它产生一个紧凑的字母显示。这可以用multcomp包来完成，但我只能找到使用该包的主效果来做这件事的方法，而不是交互
浏览 4提问于2013-08-01得票数 1
我在一个函数中使用rstatix包中的tukey_hsd()后，遇到了add_xy_position()的问题。下面是我的函数： make_tukey_test <- function (data,variable,grouping_variable){
data %>%
tukey_hsd(variable ~ grouping_variable)
} 当我使用以下代码调用该函数时，它工作得很好，并且测试结果保存在data.frame中： test <- make_tukey_test(data = dat, variable = dat$num_var,
浏览 161提问于2021-10-26得票数 0
回答已采纳
当我查看我的数据集的方差分析时，我得到的P值为0.018，这在alpha=0.05中很重要。然而，当我在方差分析上做Tukey的HSD测试时，所有的东西都被归入A组，有人能解释一下这里发生了什么吗？ 以下是我的数据： dataHodag <- read.table(text="
Plot Rate Rep HollowHeartPercent
1 2155 1 1 0.0
2 2252 1 2 0.2
3 2343 1 3 0.2
4 2433
浏览 29提问于2019-10-15得票数 1
回答已采纳
3回答
问题：我想学习如何改变图基在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
浏览 0提问于2015-04-24得票数 4
我想从函数内部的Tukey测试中生成标签。在那里，它会抛出一个错误，而如果我调用代码，而不是函数中的代码，它就会运行得非常好。
其他线程表示，解决方案是使用"as.character“转换变量，但当我尝试时，它没有工作。
遗憾的是，我找不到错误，如果你能帮我的话，那就太好了。错误和跟踪在最后一行。
最小工作实例：
require(plyr)
require(multcomp)
require(multcompView)
require(datasets)
data(iris)
iris
Dataset <- iris
############################
浏览 4提问于2019-10-30得票数 1
回答已采纳
1回答
我正在尝试复制一个用于单向混合效应方差分析的Tukey检验的教科书示例(来自Statistics，William L.Hays p 581-583)，但我使用lme & glht获得的p值没有意义
该研究重复测量了4个处理和10个受试者
数据
subject=c(1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6,
6, 7, 7, 7, 7, 8, 8, 8, 8, 9, 9, 9, 9, 10, 10, 10, 10)
treatment=c("a1", "a2&
浏览 21提问于2018-02-20得票数 2
我想知道我的代码出了什么问题？
我的代码如下：
a=aov(tem.range~Butterfly.colonies, data=my.data)
summary(a)
tHSD <- TukeyHSD(a, ordered = FALSE, conf.level = 0.95)
tHSD
generate_label_df <- function(HSD, flev){
# Extract labels and factor levels from Tukey post-hoc
Tukey.levels <- HSD[[flev]][,4]
Tukey.labe
浏览 0提问于2017-06-07得票数 0
我正在尝试使用Agricolae包中的HSD.test对我的数据进行图基测试。但是，我的结果总是单词NULL。 我的数据： structure(list(Sample = c(1011L, 1012L, 1014L, 1024L, 1025L,
1026L), Collection = c(1L, 1L, 1L, 1L, 1L, 1L), Irrigation = structure(c(3L,
3L, 3L, 5L, 5L, 5L), .Label = c("Rate1", "Rate2", "Rate3", "R
浏览 24提问于2019-06-21得票数 0
1回答
我已经尝试运行一个多变量数据集的Tukey HSD。但是，当我对单个变量运行相同的测试时，结果完全相反。
在运行多个变量时，我观察到在ANOVA输出中出现了以下错误：
87个影响中有8个是不可估计的估计影响可能是不平衡的
在运行单变量时，我观察到在ANOVA输出中出现了以下错误：
估计的影响可能是不平衡的
这与我收到的完全相反的Tukey HSD输出有任何关系吗？另外，我如何继续解决这个问题呢？
我使用了aov()，在我的数据集中有将近500000个数据点。
具体来说，下面的代码给了我一个不同的结果：
code1：
lm_test1 <- lm(y ~ x1+ x2, d
浏览 1提问于2016-04-18得票数 0
好的，我想我已经很接近了，但是当我试图在最后构建我的箱形图时，我得到了一个错误。我的目标是在每个盒子图的上方放置字母，表示时间点之间的统计关系。我在这个网站上看到了两个关于这一点的讨论，可以从他们的代码中复制结果，但不能将其应用于我的数据集。
套餐
library(ggplot2)
library(multcompView)
library(plyr)
以下是我的数据：
dput(WaterConDryMass)
structure(list(ChillTime = structure(c(1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 3L, 3L, 3L, 3L, 4L,
浏览 24提问于2017-06-23得票数 2
回答已采纳
我对R语言比较陌生，并且使用它来分析我的数据。我正在使用一个名为'agricolae‘的包和它的HSD.test组件来输出基于处理它的anova的tukeyHSD结果。我简化了下面的代码，但本质上讲，tukey代码是在一个for循环中运行的，该循环遍历了一系列因素/交互。它工作得很好，但是，我只希望循环在相应因子的anova结果显着的地方运行。(即<= 0.05)。
f1 <- as.formula(paste(m,'~ Treatment + Genotype + GENOTYPExTREATMENT', sep = ''))
anova.
浏览 3提问于2020-07-31得票数 2
在运行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'],
浏览 1提问于2018-10-03得票数 0
1回答
我想对数据帧列表执行Tukey Post Hoc测试。作为结果，我希望有信件表明哪些群体彼此之间存在显着差异。agricolae包的HSD.test()可以做到这一点，但是我不知道如何将它同时应用到多个变量上。这将为我节省大量时间，因为我的数据集包含许多变量。
这是我数据的一部分：
category <- c(rep("young", 3), rep("Middle", 4), rep("old", 5))
fat <- c(1857.87, 1953.90, 1440.70, 1553.81, 1785.91, 1893.82, 1
浏览 0提问于2018-03-08得票数 2
在使用rstatix进行了Tukey测试之后，我想用Tukey组向数据帧中添加另一列。
library(tidyverse)
library(rstatix)
ggplot(data = iris,
aes(x = Species,
y = Petal.Length)) +
geom_boxplot()
Tukey <- iris %>%
tukey_hsd(Petal.Length ~ Species) %>%
add_significance() %>%
# ????? <------- Help her
浏览 15提问于2021-12-26得票数 1
我想将Tukey.HSD后期测试的结果添加到ggplot2箱线图中。包含一个我想要的手动示例(即，绘图上的字母是手动添加的；共享一个字母的组是无法区分的，p>whatever)。
有没有一个基于AOV和Tukey HSD post-hoc分析的自动添加字母到箱线图的功能？
我认为编写这样的函数并不太难。它看起来像这样：
set.seed(0)
lev <- gl(3, 10)
y <- c(rnorm(10), rnorm(10) + 0.1, rnorm(10) + 3)
d <- data.frame(lev=lev, y=y)
p_base <-
浏览 1提问于2013-09-13得票数 11
回答已采纳
我尝试使用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
浏览 2提问于2018-11-04得票数 5
1回答
我有一个不平衡的实验，在三个地点(L，M，H)，我们测量四种不同植被类型(a，b，c，d)的参数(met)。所有的植被类型都出现在所有三个地点。植被类型在L和M上重复4次，在H上重复8次。
因此，简单的anova和TukeyHSD是行不通的。包农业(HSD.test)和DTK (DTK.test)只工作于单向设计，然后有多个复合...mcp函数中的Tukey测试是计算Tukey-Kramer对比度，还是给出常规的Tukey对比度？我假设是第一种情况，因为该软件包旨在测试不平衡设计的多个比较，但我不确定，因为两种方法产生的p值实际上是相同的。那么什么测试才是合适的呢？
此外，有没有更合适的方法来
浏览 2提问于2012-09-12得票数 7
回答已采纳
我试图对我的数据执行图基的HSD测试或LSD测试。我有两个因素，收集(2个处理)和灌溉(5个处理)，并想做测试蔗糖的反应，从每个组合，所以10个总处理。
数据：
structure(list(Collection = structure(c(1L, 1L, 1L, 1L, 1L, 2L
), .Label = c("1", "2"), class = "factor"), Irrigation = structure(c(1L,
2L, 3L, 4L, 5L, 1L), .Label = c("Rate1", "Rate
浏览 1提问于2019-06-19得票数 0
1回答
我想知道是否有可能在Python中找到Studentized统计量(q*)作为Tukey的HSD计算的输入，类似于插值这样的表()或从连续分布中提取。
我已经在这里找到了一些指导()，但在如何输入df (自由度)或k(样本组的#)方面迷失了方向。
我正在寻找类似于临界F或临界t统计量的东西，这些统计量可以通过
Scipy.stats.f.isf(α，df-中间，df-内部)
或
Scipy.stats.t.isf(α，df)。
浏览 7提问于2015-04-16得票数 4
1回答
我试图评估几个信号的频域。为此，我使用了中给出的PSD实现。作为比较，我使用了the中提供的signal.periodogram函数：
from scipy.signal import tukey
import scipy as sp
f, Pxx_den = sp.signal.periodogram(a_gtrend_orig,12,window=tukey( len(a_gtrend_orig) ))
然而，当我在自我实现的PSD旁边画这幅图时，它们看起来却有很大的不同：
由于使用相同的窗口函数，而周期图函数也应该使用FFT，这种差异从何而来？
浏览 0提问于2018-10-09得票数 2
回答已采纳
我做了5x2的交叉验证实验，之后，我做了Tukey-HSD成对比较，比较了5种技术的10个准确性，如下所示。
上图是在执行以下R命令后生成的：
data <- read.table("experimento-geral.txt", head=TRUE, sep=",", dec=".");
data$metodo <- factor(data$metodo);
summary(data$acuracia)
aov.data <- aov(data$acuracia ~ data$metodo, data=data)
a
浏览 4提问于2014-11-01得票数 0
我得到了这些原始数据在Spss中使用，我很困惑，因为我已经习惯了R。一项实验监测了厌食症女孩在各种治疗后体重增加的情况。女孩被分配到三组中的一组。第1组不接受治疗，第2组接受认知行为治疗。第3组采用家庭治疗。研究人员想知道，与对照组相比，这两个治疗组是否会增加体重。
这是数据
第1组<- c(-9.3，-5.4,12.3，-2，-10.2，-12.2,11.6，-7.1,6.2,9.2,8.3,3.3,11.3，-10.6，-4.6，-6.7,2.8,3.7,15.9，-10.2) 2<-c(-1.7，-3.5,14.9,3.5,17.1，-7.6,1.6,11.7,6.1，-4
浏览 3提问于2013-10-15得票数 2
有没有办法在与Kruskall wallis测试进行成对比较后获得显著性代码？对于重要性代码，我指的是分配给人群的字母代码，以指示差异显着的地方。
对于传统的anova，可以使用agricolae库中的HSD.test执行这样的测试，但是对于非参数化的anova，我还找不到任何东西。
一个小玩具示例：
dv <- c(runif(100, 5.0, 10))
iv <- as.factor( c(rep("I", 10), rep("II", 10), rep("III", 10), rep("IV",
浏览 0提问于2013-06-05得票数 4
回答已采纳
1回答
我是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
浏览 1提问于2014-08-15得票数 0
1回答
我正在把代表图基的HSD结果的字母放在我制作的盒子里。我在这里找到了一个关于如何做它的描述：，但是在运行函数之后，我会得到一个错误。
这是我的数据：
dput(WaterConDryMass)
structure(list(ChillTime = structure(c(1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L), .Label = c("Pre-chill",
"6", "13", "24", "Pos
浏览 3提问于2017-06-21得票数 2
回答已采纳
我有4个样本，代表他们的均值和标准偏差。如果我认为所有样本的样本大小为100，我如何使用R中的R库进行单向方差分析？
样品均为正态分布且相互独立。
这些是他们的价值观：
S1: mean1=5.5，sd1=0.6
S2: mean2=5.6，sd2=0.5
S3: mean3=5.9，sd3=0.5
S4: mean4=6.1，sd4=0.6
我如何执行后专案对明智的比较(图基-HSD)使用R以及？
浏览 2提问于2014-06-25得票数 1
我正在尝试执行Tukey的HSD测试，看看我的数据中的几个组的值平均值是否存在显著差异。例如，在这里，我试图看看变量'acad_se_communicate_needs‘按组'Class’是否存在平均差异。但是，我在结果中遇到了NaN值。这是怎么回事，我该怎么解决呢？
我已经使用状态模型函数来完成这个任务。我避免了需要为每个组将数据分割成不同数据的方法，因为我必须对多个变量执行此分析。而且，这些方法对我来说真的很难理解。
from statsmodels.stats.multicomp import pairwise_tukeyhsd
from statsmodels.stat
浏览 0提问于2019-07-29得票数 0
回答已采纳
2回答
R限制法很慢
、、、
下面是代码示例
cq = c(15.5193589369923,15.6082097993496,15.3048276641115,15.887944383963,15.3813224284544,14.9723432922121,14.8742353464212,15.0448020475332,15.1584221729435,15.3692219904727,15.2369219681739,15.0804950645883,15.0836397511495,14.8821059462034,14.6827696388115,14.5701385743889,14.85062481036
浏览 1提问于2016-07-29得票数 3
回答已采纳
函数stat_pvalue_manual()将使用ggboxplot将p值添加到盒图中。然而，打印的p值有时有很多位长.我想把小数点的数目限制在3，我怎么做呢？
从下面的图中，您将看到，versicolor和virginica的p值为5小数位，我如何调整下面的代码来报告小数点后3位(即0.009位)？
library(ggplot2)
library(tidyverse)
library(dplyr)
library(rstatix)
library(ggpubr)
test <- iris
test$Species <- as.factor(test$Species)
test
浏览 10提问于2022-03-13得票数 0
回答已采纳
我正在寻找一种保存结果的方法，将Tukeyhsd的结果保存到一只熊猫数据中。见下文：
import matplotlib.pyplot as plt
import statsmodels.formula.api as smf
import statsmodels.stats.multicomp as multi
mcDate = multi.MultiComparison(df['Glucose'], df['Date'])
Results = mcDate.tukeyhsd()
print(Results)
Multiple Compar
浏览 1提问于2016-11-09得票数 6
回答已采纳
我希望将带有自定义因子的变量从ANOVA模型赋给emmeans()语句。在这里，我使用来自R的橙子数据集来使代码可复制。这是我的模型，我通常是如何计算因子存储的方法的：
library(emmeans)
oranges$store<-as.factor(oranges$store)
model <- lm (sales1 ~ 1 + price1 + store ,data=oranges)
means<-emmeans(model, pairwise ~ store, adjust="tukey")
现在，我想给一个变量( lsmeans )赋值，定义计
浏览 5提问于2021-06-29得票数 2
回答已采纳
2回答
R:图基跟随单向方差分析
、、、
我在R中进行了一次单向的方差分析，但是当我尝试做一个图基临时测试来了解哪些治疗方法不同时，我会不断地收到错误信息。(我希望对结果进行排序(a，ab，b，bcd...etc)。
数据详细信息：
data = "abh2“
X=6个处理：“处理”
y=水分读数“湿润”(每次处理的n=63，total=378)
我做了一次单向方差分析：
anov <- anova(lm(moist~treatment, data=abh2))
.#结果表明，我可以移动到post (p<0.05)：
Analysis of Variance Table
Response: moist
浏览 0提问于2017-02-02得票数 1
回答已采纳
我试着用图基的后置分析从方差分析中添加错误条到我使用ggplot2的R线图中。谢谢提前，我很高兴张贴更多的信息，将帮助我解决这个问题！
到目前为止，我的代码如下：
#### ANOVA ####
aov_CW <- aov(CW ~ Subject, data = VisualAcuity)
summary(aov_CW)
model.tables(aov_CW, "means")
#### Tukey HSD post-hoc Test ####
TukeyHSD(aov_CW, conf.level = 0.95)
#### Line Plot ####
Vis
浏览 4提问于2020-04-02得票数 0
回答已采纳
2回答
我试图绘制Tukey测试的结果，但我很难根据P值将数据分组。我试图复制的。我一直在使用SciPy单向方差分析测试和Tukey测试状态模型，但无法以同样的方式完成这些测试。
任何帮助都是非常感谢的。
我还找到了这个
浏览 5提问于2017-05-15得票数 1
回答已采纳
1回答
我使用来自包的lmer来运行线性混合效应模型。我有3年的温度数据未处理(5)和处理的地块(10)。模式：
modela<-lmer(ave~yr*tr+(1|pl), REML=FALSE, data=mydata)
模型检查残差的正态性；qqnorm绘制我的数据：
'data.frame': 6966 obs. of 7 variables:
$ yr : Factor w/ 3 levels "yr1","yr2","yr3": 1 1 1 1 1 1 1 1 1 1 ...
$ pl : Factor w/
浏览 9提问于2017-04-19得票数 0
我正在寻找一种方法来从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
浏览 1提问于2014-09-11得票数 1
我关于在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
浏览 0提问于2016-11-16得票数 0
1回答
我和图基有问题..。我真的不明白为什么..。
我尝试过几个代码，我也遇到了使用as.factor解决的因子变量的问题。但是with NaN产生了我不知道如何管理
我的数据：
> zz2
id Mean group
1 1 1.17 1
2 2 1.57 2
3 3 1.16 2
4 4 1.14 1
5 5 1.81 4
6 6 1.33 1
7 7 1.27 1
8 8 1.57 1
9 9 1.56 3
10 10 1.36 1
11 11 1.09 1
浏览 3提问于2020-03-24得票数 0
回答已采纳
当调用ggpubr的stat_pvalue_manual函数时：
ggplot + stat_pvalue_manual(statistics, label = "p.adj.signif",
coord.flip = TRUE,
#hide.ns = TRUE,
tip.length = 0,
bracket.size = 0)+coord_flip()
绘制了以下图表：
但是，我只对有意义的数据感兴趣，因此我通过调用hide.ns
浏览 6提问于2022-01-24得票数 0
回答已采纳
1回答
我使用lme4软件包运行了线性混合模型。然后，要查看所有的对比，我必须使用lsmeans包运行lsmeans函数。这是我运行的函数：
library('lme4')
library('lsmeans')
lsmeans(lmer1, pairwise ~ vowel * experience * lang_sound, adjust="tukey")
但是，我看不到整个输出，因为有太多的对。有人能告诉我应该怎么做才能得到这个函数的输出吗？
我试过“sink()”，但它不起作用。可能是因为我运行的lsmeans命令没有名称。我在Windows上使用
浏览 1提问于2017-03-30得票数 2
回答已采纳
我不知道我是否在题目中使用了正确的字眼，所以请让我解释一下。比如说，我有一个包含三个因素的数据集:物种、生命和暴露，还有一个响应变量: mean_v，总共有108个层次。当我应用3路方差，然后是图基的后特别测试，它产生了非常多的统计比较。
life <- rep(c("5d", "15d", "45d"), 2, each = 18)
species <- rep(c("SP1", "SP2"), each = 54)
exposure <- rep(c("c1", "c2
浏览 1提问于2018-04-16得票数 0
回答已采纳
我想把代表图基的HSD的结果的字母放在我用ggplot制作的盒子图上。我知道其他的帖子(，)，但在我的例子中，我与facet_grid一起工作，我不知道如何继续。由于我正在创建很多的盒子情节，我希望有一些自动化的东西。任何帮助我开始工作的建议都将不胜感激。非常感谢
library(multcompView)
library(ggplot2)
#sample dataframe
df <- data.frame(var_facet=c('facet1', 'facet1', 'facet1', 'facet1', 'f
浏览 0提问于2018-10-29得票数 2
回答已采纳
我得到了下面的代码，并被要求提取原始的p值，而不是Tukey调整后的值(因为我们将在稍后阶段使用Homes Bonferroni调整多个比较)，但我不确定如何将"Tukey“替换为(我不熟悉使用R).....
res=glht(x, linfct=mcp(Letter="Tukey")
out=summary(res)
out
浏览 0提问于2018-10-29得票数 1
我试图为一个数据集运行TukeyHSD (40个观察)。结果是'mbvs‘，这两个因素是性别(0或1)和治疗(0或1)？
我试图找出男性治疗组和对照组(sex=0)之间的平均结果有什么不同？
有人能帮我查一下这个密码吗？谢谢。
浏览 2提问于2021-12-08得票数 0
回答已采纳
我有麻烦回到"CL“值给出的电子邮件包。
我的第一个目标是进行一项权力研究:我正在研究植物，我想看看我能为下一季改进我的实验设计，以使我的基因型之间有更大的差异。我的基因型是重复的，这被用作阻塞因子。我想玩的重复次数(假设手段和剩余SE将保持不变)。最后，我想在预算问题和统计权力之间找到一个折衷方案。
我在考虑置信区间(CI)，考虑到如果两个基因型的顺铂覆盖至少一个共同值，两种基因型之间没有显著差异。所以我想看看我在玩重复次数时能减少多少顺式干扰。
对于使用Tukey方法进行多个比较的CIs计算，我的猜测是：
μi±(σ/(2*√(Ni)*q
以μi为第一基因型的平均值
第一种基因型
浏览 0提问于2018-04-20得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
云直播活动推荐
腾讯云开发者
