R 显示wilcox测试结果以及绘图PDF

R 显示wilcox测试结果以及绘图PDF,r,R,我编写了一个函数,该函数将为每类数据生成一个绘图,我还创建了一个函数,该函数将对所有类以及每个类单独执行wilcox测试 我想做的是获取这个测试的结果,并将其显示在每个类的每个特定PDF中。这能以一种好的方式完成吗 这是我用来执行wilcox测试的函数代码: library(ggplot2) t = read.table("data.txt", stringsAsFactors=FALSE) uni = unique(c(t$V9)) testresult<- function(colu

我编写了一个函数,该函数将为每类数据生成一个绘图,我还创建了一个函数,该函数将对所有类以及每个类单独执行wilcox测试

我想做的是获取这个测试的结果,并将其显示在每个类的每个特定PDF中。这能以一种好的方式完成吗

这是我用来执行wilcox测试的函数代码:

library(ggplot2)
t = read.table("data.txt", stringsAsFactors=FALSE)
uni = unique(c(t$V9))

testresult<- function(column, datacol) {    
  if(column=="PP") {
    for (xp in uni) { 
      for(yp in uni) {
        testx<- subset(t, V9==xp)
        testy<- subset(t, V9==xp)
        zz = wilcox.test(testx[[datacol]],jitter(testy[[datacol]]))
        print(p1)
        p1 = c(xp)
        p2 = c(yp)
        pv = c(zz$p.value)
      }
    }
  } 
}
for (x in uni) { 
  test<-subset(t, V9==x)
  dat <- melt(test, measure=c("V1", "V2" , "V3" , "V4"))

  p <- ggplot(dat, aes(Class,value , color=variable)) + 
       geom_boxplot() + 
       geom_point(position="jitter" , alpha=0.3) +  
       facet_grid(variable~., scales="free_y")
       ggsave(paste("ggplot2/",paste(x,".pdf")), p, width=25, height=15)
}
库(ggplot2)
t=read.table(“data.txt”,stringsAsFactors=FALSE)
uni=唯一(c(t$V9))

testresult您可能希望使用Swave生成混合文本和绘图的漂亮PDF报告。

您希望在每个绘图中的何处显示它?作为标题,在绘图区域本身内注释?这里有一些选择

#No sample data provided, so we use this instead.
dat <- data.frame(x = rnorm(100), y = rnorm(100))

ggplot(dat , aes(x,y)) + 
geom_point() +
opts(title = "hi world") +  #This line makes a title
geom_text(aes(max(x),max(y), label=paste("The mean of x = ", round(mean(x),3), sep = " ")), hjust = 1, size = 3)     #This line annotates within the plot region
#未提供示例数据,因此我们使用此示例数据。

谢谢你的回复,但是我怎么用这个??我是R的新手,你能给我举个例子吗?关于Swave的话题太广泛了,答案中无法涵盖。打开你最喜欢的网页搜索,键入类似于“Swave demo”的内容。应该让你开始。非常感谢,但同时我可能需要一个更简单的方法,正如我在Chase Comments上的回复中提到的,你已经提出了几个与同一问题相关的问题,并且提供了一个完全没有问题的可复制示例。这里有很多人愿意帮助你,但你目前的问题结构几乎不可能提供任何合理的答案。R提供了非常好的函数
dput()
,这使得其他人很容易复制和运行您的代码。目前,没有人可以对您的代码进行任何评估,因为我们不知道
data.txt
的结构。如果您不想继续学习Swave,请查看。感谢您的帮助,我想在另一页的PDF文件中将结果显示为表格,我希望它能为每一类数据自动生成。正如你在我的代码中看到的,我为每个(xp,yp)都有一个p值,所以我想要一个矩阵,在数据中每个类都有这些列:xp-yp p-value含义您好,这可以用一种简单的方式来完成,比如在矩阵中填充它们,然后将它们添加到绘图中??还是我应该用Swave的东西??谢谢,这是我正在处理的txt文件的一部分,我唯一想做的就是输出测试结果,正如我所说的,在一个由4列组成的表中,每个列(xp,yp),thnx提前