Показать результаты теста Уилкокса вместе с графиком PDF

Question

Я написал функцию, которая будет генерировать график для каждого класса данных, и я также создал функцию, которая будет выполнять тест Уилкокса для всех классов и только для каждого класса.

Что я хочуДля этого нужно взять результат этого теста и отобразить его в каждом конкретном PDF для каждого класса.Можно ли сделать это хорошим способом?

Это код функции, которую я использую для выполнения теста Уилкокса:

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)
}

Я пытался заполнить его внутри спискавнутри фрейма данных, но это не сработало.

Answer 1

Где вы хотите отобразить его на каждом графике? Как заголовок, аннотированный внутри самого участка сюжета? Вот несколько вариантов.

#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

Answer 2

Возможно, вы захотите использовать Sweave для создания хороших отчетов в формате PDF, смешивая текст и графики.