使用ggboxplot定位选定的p值
ggboxplot是一个R语言中的数据可视化包，用于绘制箱线图。箱线图是一种常用的统计图形，用于展示一组数据的分布情况，包括中位数、四分位数、异常值等。
使用ggboxplot可以帮助我们定位选定的p值，即在一组数据中找到特定的p值所对应的位置。具体步骤如下：
- 导入ggplot2和ggpubr包：在R语言中，首先需要导入ggplot2和ggpubr包，这两个包提供了绘制箱线图的函数。
library(ggplot2)
library(ggpubr)- 准备数据：将需要绘制箱线图的数据准备好，可以是一个向量、一个数据框或一个矩阵。
data <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10)- 绘制箱线图：使用ggboxplot函数绘制箱线图，并设置参数来定位选定的p值。
ggboxplot(data, ylab = "Value", notch = TRUE, add = "jitter", pvalue = TRUE, pval.method = "t.test", pval.adjust.method = "bonferroni")参数解释：
- data：需要绘制箱线图的数据。
- ylab：y轴标签。
- notch：是否绘制缺口。
- add：是否添加抖动点。
- pvalue：是否显示p值。
- pval.method：计算p值的方法，这里使用t.test表示使用t检验计算p值。
- pval.adjust.method：p值校正方法，这里使用bonferroni表示使用Bonferroni校正。
- 结果解读：根据绘制的箱线图，可以观察到中位数、四分位数、异常值等信息。如果设置了显示p值，还可以根据p值来判断数据的显著性差异。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云服务器（CVM）：提供弹性计算能力，满足各类业务需求。产品介绍链接
- 腾讯云数据库（TencentDB）：提供高性能、可扩展的数据库服务。产品介绍链接
- 腾讯云人工智能（AI）：提供丰富的人工智能服务，包括图像识别、语音识别、自然语言处理等。产品介绍链接
- 腾讯云物联网（IoT）：提供全面的物联网解决方案，帮助连接和管理物联网设备。产品介绍链接
- 腾讯云移动开发（Mobile）：提供移动应用开发和运营的一站式解决方案。产品介绍链接
- 腾讯云存储（COS）：提供安全、可靠的云存储服务，适用于各类数据存储需求。产品介绍链接
- 腾讯云区块链（Blockchain）：提供高性能、可扩展的区块链服务，支持多种场景应用。产品介绍链接
- 腾讯云元宇宙（Metaverse）：提供虚拟现实、增强现实等技术支持，构建沉浸式体验。产品介绍链接
相关·内容
我需要在ggboxplot为纸张生成的盒图上报告确切的p值，但是对于非常小的p值，stat_compare_means只给出< XX。在下面的例子中，我怎样才能得到精确的p值而不是最小的值呢？
library(ggpubr)
data(ToothGrowth)
data_to_change <- which(ToothGrowth$dose==0.5)
ToothGrowth$len[data_to_change] <- ToothGrowth$len[data_to_change]-10
ggboxplot(ToothGrowth, x = "dose"
浏览 11提问于2022-11-14得票数 1
我试着用显示p值的方式制作盒图。
my_comparisons <- list( c("0.5", "1"), c("1", "2"), c("0.5", "2") )
ggboxplot(ToothGrowth, x = "dose", y = "len")+
stat_compare_means(comparisons = my_comparisons, method = "wilcox.test")
结果是
但是如何在计算出
浏览 9提问于2021-12-07得票数 3
使用ToothGrowth数据集(内置到R中)，我使用了以下代码。
library(ggplot2)
library(tidyverse)
library(ggpubr)
p <- ggboxplot(ToothGrowth, x = "supp", y = "len",
color = "supp", palette = "jco",
add = "jitter",
facet.by = "dose", short.panel
浏览 0提问于2019-05-09得票数 6
回答已采纳
我想要将p值添加到ggg图中。他们应该有三个3小数位，如果p< 0.001，那么它应该显示("p < 0.001")。"p“应该用斜体表示。不知怎么的，一个人必须用sprintf来做这件事，但究竟是怎样做到的呢？
library(ggplot2)
library(tidyverse)
library(ggpubr)
ggboxplot(ToothGrowth, x = "dose", y = "len")+
stat_compare_means(label = "p.format")
浏览 7提问于2022-09-12得票数 2
回答已采纳
2回答
让我们想象一下我有这样的命令
library(ggpubr)
ggboxplot(ToothGrowth, x = "supp", y = "len",
color = "supp", palette = "jco",
add = "jitter")
p + stat_compare_means(method = "t.test")
我想把y值画成一个计算值：2^len，并在len原始值上留下比较方法。
有什么帮助吗？
浏览 3提问于2020-10-16得票数 1
我想限制y轴的末端，而不限制y-axis上该图中显示的数据的范围，从0到32。具体地说，对于较高的值，不应显示轴，但应包括显著性比较，并且不应显示红色部分。
library(ggpubr)
data("ToothGrowth")
head(ToothGrowth)
# Visualize: Specify the comparisons you want
my_comparisons <- list( c("0.5", "1"), c("1", "2"), c("0.5", "2&#
浏览 0提问于2021-09-24得票数 3
3回答
我试图迭代地注释许多情节。我希望在所有图上都有一个注释，但它们都有不同的y轴标度。如何在不使用y-值的情况下注释一个图？必须有一个独立于数据的方法，但我想我缺少一个关键的词汇，所以我的搜索是徒劳的。
通常情况下，我会将annotate("text", x = 5, y = 10)添加到ggplot对象中，但这是行不通的，因为y的范围在不同的组之间有很大的差异。
请记住，大约有30个地块可供跟踪，下面是一个使用mtcar和ggpubr的玩具示例。
library(ggpubr)
ggboxplot(mtcars, x = "gear", y = "wt&#
浏览 11提问于2022-06-13得票数 2
我想使用包ggplot中的R函数stat_compare_means，将Kruskal测试的p值绘制到我的ggpubr。
但是，如果只运行函数，则绘制的值与值不同：
kruskal.test(value ~ type, data = Profile_melt)
我绘制p值的代码是：
ggplot(Profile_melt, aes(type, value)) +
geom_boxplot(aes(fill = factor(type), alpha = 0.5),
outlier.shape = NA, show.legend = FALSE) +
浏览 8提问于2019-12-26得票数 3
回答已采纳
我有一个动画，可以从1到50个关键点。在第50个键之后，我希望键旋转180度
我可以使用以下命令在1个选定的关键点中旋转1个控制器
//mel
rotate -y 50 -ws -p 5 6 7;
但是非常耗时，我如何才能对所有选定的控制器和关键点执行此操作？
因此，有一种在关键点上执行mel命令的直接方法。如何继续做这样的事情呢？
浏览 3提问于2016-02-26得票数 1
我正在尝试用R ggboxplot创建一个带有'star‘符号的分组箱形图。然而，我在添加误差条时遇到了困难，特别是在上下水平线的情况下。这是我的代码和输出。
library(ggpubr)
library("ggsignif")
s=15
p <- ggboxplot(data, x = "MRI", y = "featurevalue",
color = "black",fill="filename",palette = "aaas",
浏览 27提问于2020-02-10得票数 0
1回答
有时，在没有胡须的情况下显示一个盒子情节是合适的。在geom_boxplot (ggplot2)中，我们可以通过coef=0来实现这一点。是否有一种方法可以在ggboxplot (gpubrv0.5.0，编写时的当前版本)中实现这一点？我注意到ggboxplot与geom_boxplot有很多共同之处，比如在每种情况下使用outlier.shape=NA来抑制异常值的能力。似乎应该有一个简单的方法来抑制胡须。
浏览 4提问于2022-11-23得票数 1
回答已采纳
这种标记均值的尝试会导致每个组得到一个标签的中间值。
ggboxplot(DataW, x = "Factor-r-Level", y = c("Var1","Var2","Var2"),
merge = TRUE, palette = "jco", add="mean") +
stat_summary(fun.y=mean, colour="darkred", geom="point", hape=18,
size=3,show_guide = FALS
浏览 2提问于2020-09-26得票数 1
回答已采纳
我被困在如何绘制分组箱形图的平均p值上。
这是我数据的一部分：
ID<-c('E5b','R6',"S22","E5b","R6","S22","E5b","R6","S22","E5b","R6","S22","E5b","R6","S22","E5b","R6",
"S22","E
浏览 544提问于2021-02-14得票数 0
1回答
我需要帮助来调整柱状图中的列宽度的大小。我使用下面的R脚本获得一个带有p值的条形图，但是条形图(列)的宽度太宽了。谢谢 library(ggpubr)
library(rstatix)
df <- ToothGrowth
df$dose <- as.factor(df$dose)
head(df, 6)
# Statistical test
stat.test <- df %>%
t_test(len ~ supp) %>%
add_significance()
stat.test
# Box plots with p-values
bxp
浏览 24提问于2021-08-25得票数 0
回答已采纳
下面是我目前拥有的创建下图的代码。相反，我正在尝试制作一个图表，其中蓝色的箱形图(双关)在x轴上的标题为“双关语”，旁边是黄色的箱形图“mv”。 ggboxplot(
d, x = "Condition.Name", y = c("pun", "mv"),
merge = TRUE, palette = "jco",
title = NULL,
xlab = '',
ylab = ''
) ?
浏览 347提问于2020-10-09得票数 1
回答已采纳
1回答
我从教程中获得了以下代码，用于使用ggpubr中的函数，但输出图的统计测试结果下面没有括号(或线)。 library(ggpubr)
data("ToothGrowth")
p <- ggboxplot(ToothGrowth, x = "dose", y = "len",
color = "supp", palette = "jco",
add = "jitter")
p + stat_compare_means(aes(group = supp))
浏览 17提问于2021-11-12得票数 0
1回答
将抖动添加到ggpaired
、、、、
包提供了函数ggpaired来绘制成对数据。它不像ggboxplot那样有add='jitter'选项。
有没有可能通过不同的方式达到类似的效果？
文档中的示例：
require(ggplot)
require(ggpubr)
before <-c(200.1, 190.9, 192.7, 213, 241.4, 196.9, 172.2, 185.5, 205.2, 193.7)
after <-c(392.9, 393.2, 345.1, 393, 434, 427.9, 422, 383.9, 392.3, 352.2)
d <- data.frame(
浏览 19提问于2017-12-31得票数 1
回答已采纳
如何在下图中更改stat_compare_means的字体大小？即，更改"Kruskal-Wallis，p= 1.5e-09“和其他p值字体大小？我想使用比默认字体更小的字体...
按照数据示例...
library(ggpubr)
data("ToothGrowth")
compare_means(len ~ dose, data = ToothGrowth)
# Visualize: Specify the comparisons I want
my_comparisons <- list( c("0.5", "1"), c
浏览 8提问于2018-02-01得票数 5
回答已采纳
我正在尝试复制一个盒子图使用R g图。然而，我有困难添加错误栏，特别是与上和下水平线。在我的图表中，抖动点的分布也不均匀。这是我的代码和输出
survey <- data.frame(sample=rep(c("sample1","sample1", "sample1", "sample1", "sample1", "sample1", "sample1", "sample1", "sample1", "sample1",
浏览 3提问于2020-01-02得票数 1
回答已采纳
最初，我问这个如何压制ggboxplot制作的盒子里的胡须。(无法使用设置geom_boxplot选项的预期方法。)出现了一个很好的解决方案，它符合原来的问题。然而，需要解决的更广泛的问题是如何在盒图上抑制晶须，但仍然保留对ggpubr中漂亮的添加的访问，例如能够自动计算统计测试结果并将它们放在方格图上。
为了得到我想要的东西，我修改了@Julian_Hn的解决方案。既然我已经提出了更广泛的问题，那么有两个更有知识的人可能会帮上忙：
有办法使解决方案更有效吗？，如何添加更改x-值范围的能力？(我尝试了使用ggpar和coord_cartesian的各种方法，但没有效果。我可能对如何使用像eff
浏览 1提问于2022-11-25得票数 1
回答已采纳
1回答
我使用的是ggboxplot图和"facet.by“参数。默认设置为水平面，但我希望我的情节是垂直的(即。它们都有相同的x轴：“年龄”)。我已经扫描了这样和其他网站，但没有看到这个问题被问到。
我还尝试使用相同的"facet.by“分类变量过滤我的数据，然后使用牛仔的plot_grid和gridExtra的grid.arrange函数。牛仔图和gridExtra的问题是，它们都删除了我的scale_y_beak()设置。我尝试使用free_y，但现在它为我的绘图创建了三个不同的y轴，我只想将一个y轴应用于所有三个绘图。
我的问题是:是否有一种方法来使用ggboxplot图中内置
浏览 11提问于2022-05-08得票数 0
我正在生成不同组的多面图，我想将不同组的平均值与对照样本进行比较。我有一个广泛的值范围，并希望以编程方式将“重要性标签”放置在每个箱线图的顶部。
这里有一些玩具数据。
library(ggpubr)
ggboxplot(mtcars, x = "am", y = "mpg") +
facet_wrap(~vs) +
stat_compare_means(aes(label = ..p.signif..),
method = "t.test", ref.group = "0")
给出了下面
浏览 110提问于2021-02-18得票数 1
我想从HTML输出中的引用中删除一个带有按钮的选定子项。
关键字是聊天消息或帖子。
现在我想要的是，用户可以用按钮删除所选的消息/帖子。
我可以全部删除，但不能完全删除选定的子项。
结果和内容来自“for in”循环。
删除选定的(子)与角度和AngularFire是没有问题的，但我希望这没有角度。
// get ref
var dbRef = firebase.database().ref().child("user/chat");
// child added
dbRef.on("child_added", function(snapshot) {
var
浏览 23提问于2017-02-22得票数 0
回答已采纳
1回答
在使用ggboxplot时，有没有办法自定义图例中的标签？我有一个带有图例的条形图，其中的标签取自数据框列。我想要更改这些标签，而不必更改列的标签。类似于ggplot中的scale_fill_discrete(labels=c(......)。
浏览 2提问于2017-08-10得票数 0
1回答
我用鸡蛋做了一个多重图：：egg排5位数。我想知道怎样才能使b，d，c，e在垂直位置上的地块面积对齐？
PS:图b和c必须在y轴上以不同的小数点表示。。
library(egg)
library(ggplot2)
data("ToothGrowth")
data("mtcars")
P1 <- ggplot(mtcars, aes(x = wt, y = mpg, color=cyl))+
geom_point() # Add correlation coefficient
浏览 6提问于2020-02-04得票数 2
回答已采纳
1回答
绘制组内和组间的差异
、、、
我目前正在使用ggpubr来绘制我的值。由于这些方面的原因，我成功地根据某些因素创建了子图。最后，使用函数stat_compare_means() i只绘制了组之间的差异()。
生成我的数字的最小代码是
library(ggpubr)
data("ToothGrowth")
my_comparisons <- list( c("0.5", "1"), c("1", "2"), c("0.5", "2") )
p <- ggboxplot(ToothGrowth, x
浏览 5提问于2020-01-20得票数 1
1回答
我编写了一个for循环来生成一个图表矩阵，但是我编写了如下代码
par(mfrow = c(4, 4))
for (i in c(6:17)) {
print(ggpubr::ggboxplot(logdat,
x = "Diagnostic",
y = names(logdat)[i] ,
color = "Diagnostic",
add
浏览 1提问于2018-10-02得票数 0
回答已采纳
1回答
我使用mtcars数据集作为示例来运行以下代码。 library(ggplot2)
library(ggpubr)
ggboxplot(mtcars, x = "cyl", y = "drat", fill = "cyl",
facet.by = "am", width = 0.5, outlier.shape = NA,
bxp.errorbar = TRUE, bxp.errorbar.width = 0.2) +
stat_compare_means(aes(label =
浏览 21提问于2019-06-05得票数 0
回答已采纳
1回答
我正在尝试报告下面的图，一些多重成对比较。为了做到这一点，我编写了以下代码 p <- ggboxplot(data, x = "COND", y = "P3(400-450).Pz",
color = "SES", palette = "npg",
facet.by = "SES", short.panel.labs = FALSE)
# Use only p.format as label. Remove method name.
p + sta
浏览 40提问于2021-11-10得票数 0
回答已采纳
我希望使用ggpubr中的表达式功能，使我的APA测试统计信息能够在使用ggpubr中的geom_bracket进行组间比较的文本中斜体。
我的问题是，我需要删除生成的前0，但这似乎不是在字符串中呈现的(因此，第二个示例显示gsub()不能直接取出它。
library(ggpubr)
library(ggplot2)
# Example close to what I currently have in code with leading 0 in-front of p value
ggboxplot(df, x = "dose", y = "len") +
浏览 13提问于2022-06-20得票数 1
回答已采纳
1回答
我使用了常规的“盒图”函数来检查我的数据的极值值。对于项目中的演示文稿，我使用ggpubr包(它建立在ggplot2上)创建了相同的方块图。
据我所知，胡须应代表相同的区域，在这两个污点。令我惊讶的是，在ggpubr框图中出现了一个极值，而r- not图中没有显示它。
代码R-框图：
boxplot(data$vtrust_post1, data$vtrust_post2)
代码Ggpubr:
vtrust_a_b_long %>%
ggboxplot(x = "drive", y = "count", bxp.errorbar = TRUE)
对于
浏览 0提问于2020-11-25得票数 1
回答已采纳
0回答
我的问题是如何在if语句后插入多个操作。例如：
vartest <- var.test(var1 ~ group, alternative='two.sided', conf.level=.95, data=data1)
ttest <- t.test(var1~group, alternative='two.sided', conf.level=.95, var.equal=FALSE, data=data1)
if (vartest$p.value>0.05) {
if (ttest$p.value<=0.05) {
ca
浏览 2提问于2018-07-10得票数 0
回答已采纳
我有数字选择器，我用字符串数组填充它。
final String[] power = {"0.00", "0.25", "0.50", "0.75","1.00"};
npPowerSecond.setMinValue(0);
npPowerSecond.setMaxValue(power.length - 1);
npPowerSecond.setDisplayedValues(power);
如上所示，我显示值。当我打开数字选择器时，它将显示"0.00“作为选定的值。现在我选择"0.50“。现在
浏览 6提问于2017-04-05得票数 1
回答已采纳
我想将一些缺陷修复到我正在尝试用.pdf创建的RMarkdown文档中
正如你在这里所看到的，我有一系列的表格，它们超出了纸张的页边距限制。在这里
轴上的标签被覆盖。我不知道是否有必要，但对于表，我刚刚遇到了RMardown文档
```{r model, echo = FALSE, results = 'asis'}
Tables_models，~cat(knitr：：knit_print(.x))
While graphs come from this for cycle
```javascript
for (i in 5:ncol(data)) {
P
浏览 2提问于2021-11-10得票数 0
回答已采纳
我在错误信息上遇到了一些问题，我无法解决。
我使用以下代码在R中运行Kruskal Wallis，然后使用Dunn进行成对比较：
res.kruskal <- mydata %>% kruskal_test(values ~ group)
res.kruskal
stat.test <- mydata %>% dunn_test(values ~ group, p.adjust.method = "hochberg")
stat.test
stat.test <- stat.test %>% add_xy_position(x = "
浏览 4提问于2020-05-20得票数 2
回答已采纳
我有以下数据框架：
library(tidyverse)
dat <- structure(list(feat = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L), .Label = c("A", "KyteDoo", "FoldIndex"
), class = "factor"), dose = c("0.3mM", "1mM", "3mM", "10mM",
浏览 8提问于2022-01-05得票数 0
回答已采纳
我正在尝试创建一个显示多个比较的箱线图。我制作了一个玩具数据集，它给我带来的错误似乎与我在较大的数据集中得到的错误相同。 library(tidyverse)
library(ggsignif)
library(ggpubr)
dat <- data.frame(measurement = c("750","850","900", "300","200","400", "20", "30", "50"),
浏览 137提问于2020-04-10得票数 1
回答已采纳
解决了，现在我想让这些点更明显。我选择pch=21有一个充满颜色的黑色圆圈。但是，盒子变成了相同的配色方案。我怎么才能把盒子的颜色改回来呢？我希望在图A中有和图B相同的方框颜色。
library(tidyverse)
library(ggpubr)
mtcars$cyl=factor(mtcars$cyl)
p1=mtcars %>% ggplot(aes(x=cyl, y=mpg, fill=cyl))+
geom_boxplot(show.legend = F, aes(fill=cyl))+
geom_point(position=position_jitterdod
浏览 0提问于2018-09-25得票数 4
回答已采纳
我试图使用ggboxplpot将p值添加到盒图中，但是当我有多个y =值时，stat_compare_means()似乎无法工作。
下面是来自palmerpenguin数据集的示例代码
library(palmerpenguins)
library(tidyverse)
library(ggplot2)
library(ggpubr)
#Load data
data(package = 'palmerpenguins')
#Remove NA data
df_clean <- na.omit(penguins)
#Group dataset according
浏览 3提问于2022-01-30得票数 1
回答已采纳
是否可以在下面的ggpboxplot中更改p值的字体系列？字体系列应为"Times New Roman“。 使用theme_bw(base_family = "Times New Roman")不会更改p值。 library(ggpubr)
library(ggplot2)
data("ToothGrowth")
# Box plot
p <- ggboxplot(ToothGrowth, x = "supp", y = "len",
color = "supp", palette
浏览 240提问于2020-04-09得票数 3
回答已采纳
1回答
我正在尝试复制Github的定价页面，特别是团队计划中最受欢迎的边界。
如果您查看源代码，您将看到<p class="...">Most Popular</p>有position: absolute; top: 6px;。现在，据我所知，absolute定位从屏幕左上角开始采用x和y值。
为什么<p class="...">Most popular</p>标记与这种样式完全匹配呢？
浏览 3提问于2021-09-08得票数 0
1回答
我有一个假设的数据样本，如下所示： df <- read.table(header = TRUE, text =
"time goal book pen temp1 temp2 weight hight
18 13 18 15 13 15 18 13
18 13 20 16 16 15 20 19
20 12 20 14 18 20 17 15
15 19 18 16 18 20 15 15
16 17 14 12 17 20 20 15
17 12 18 16 17
浏览 27提问于2020-11-08得票数 0
回答已采纳
函数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
回答已采纳
如何检查是否在那里选定了什么？我至少需要‘事实’：选择与否(并且很好地知道如何阅读真实的设置)。
我查看了NotesViewColumn，没有在那里找到任何东西，可能是没有文档的特性：
浏览 1提问于2015-04-27得票数 0
回答已采纳
2回答
我想从列表框中选择记录并保存页面。保存后，所选记录将被填充到列表的顶部。未选中的记录将显示在“选定记录”列表下面。
如何标识列表框中的非选定记录以在下一次更新中选择记录。
如何从要存储在变量中的列表中标识选定的记录数据。
下面是非选定记录的元素
<li ng-class="{'list-group-item-info': agent.$selected}" ng-repeat="agent in newAgents" style="padding: 5px 15px;" class="list-gro
浏览 2提问于2015-12-09得票数 2
回答已采纳
按照示例，我可以向每个小平面网格添加一些文本，但在一个面板上，我需要在每组箱图上放置两个单独的标签，类似于stat_compare_means生成的标签的间距。我当前的逻辑在每个面板上只能容纳一个标签。如何防止它们重叠或为每个标签设置单独的坐标？
ID<-c('E5b','R6',"S22","E5b","R6","S22","E5b","R6","S22","E5b","R6","S22",&
浏览 8提问于2021-02-16得票数 1
1回答
我已经修改了下面从引用到我的数据集中的代码：
Day<-c(1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2)
Group<-c("A","A","A","B","B","B","C","C","C","A","A","A","A","B","B","B","C","C",
浏览 1提问于2021-01-05得票数 0
回答已采纳
我拟合了一个线性混合模型(带有lmer函数)。“比例”作为我的反应变量，“阶段”作为一个固定的因素(有5个层次)和被试作为一个随机因素。
LMM.fit<-lmer(Proportion~Stage+(1|Subject),data=tab)
Chisq Df Pr(>Chisq)
Stage 290.07 4 < 2.2e-16 ***
采用hoc函数进行了临时两两比较，以确定不同阶段的差异发生在何处：
emmeans(LMM.fit, pairwise ~ Stage)$contrasts
contrast esti
浏览 3提问于2021-09-28得票数 0
回答已采纳
我一直在使用下面的R脚本，，但是X轴的宽度太宽了。有人能帮我调整一下X轴的宽度吗？谢谢
library(ggpubr)
library(rstatix)
df <- ToothGrowth
df$dose <- as.factor(df$dose)
head(df, 6)
# Statistical test
stat.test <- df %>%
t_test(len ~ supp) %>%
add_significance()
stat.test
bxp <- ggboxplot(df, x = "supp", y
浏览 2提问于2021-08-25得票数 0
回答已采纳
大家好，我正在比较不同分布的权重。我的数据格式如下： Weight_G1 Gaussian_One
0.714484 Non_Treated
0.6992871 Non_Treated
0.7144726 Non_Treated
0.7062965 Non_Treated
0.7331027 Non_Treated
0.6927765 Non_Treated
0.5939587 Treated1
0.6230649 Treated1
0.5935826 Treated1
0.59120414 Treated1
0
浏览 22提问于2021-04-30得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
