1

我绘制了一些采用这种格式的数据

Time    SD  Average Situation   n   se

我使用了这段代码

library(ggplot2)
ggplot(data, aes(x=Time, y=Average, colour=Situation)) + 
    geom_errorbar(aes(ymin=Average-se, ymax=Average+se), width=.1) +
    geom_line() +
    geom_point()

我想在 Exp 和 Control 之间进行 T.test，并且我想使用星号在图中显示重要性*。我很难找到执行此操作的代码。我知道如何进行正常的 T.test，但我真的希望那些带有星号的图表能够显示重要性。

4

1 回答 1

0

编辑我已经根据您的评论重写了这个，以给出一个对我有用的例子：

假设您想将 t 检验中的 p 值添加到您的图表中......

set.seed(1)
### data following OPs format
data <- data.frame(Time=seq(10),
                   Average=abs(rnorm(10)),
                   Situation=rep(letters[1:2], 5)
                   )
### get p value
pval <- t.test(data$Average[data$Situation=="a"],
               data$Average[data$Situation=="b"])$p.value
### strip to 3 significant digits
pval <- signif(pval, 3)
### add asterisk if <0.05
pval <- ifelse(pval<0.05, paste0(pval, " *"), pval)
pval <- paste0("t-test \n p=", pval)

### plot as per OP
ggplot(data, aes(x=Time, y=Average, colour=Situation)) + 
    geom_line() +
    geom_point()+
### annotate it near lower left corner
    annotate("text",
             x=0.25*max(data$Time),
             y=0.25*max(data$Average),
             label=pval)

给予

在此处输入图像描述 查看?annotate更多选项。如果您想要 p 值的一系列符号，那么您可以通过调用switch而不是ifelse.

一个更大的问题可能是 at-test是否适合与时间相关的结果......

于 2013-10-05T16:40:48.437 回答