forestplot：伝説なく、各4つのcolums

Question

の凡例を与えたいが、私の森のプロットは、私はCIと各PREV1、PREV2、prev3とprev4のための凡例を置くとき、しかし、コード

rm(list=ls()) 
library(forestplot) 

test_data2<-data.frame(prev1=c(0, 70, NA), 
        prev0=c(0, 20, NA), 
        prev3=c(NA,NA , 30), 
        prev4=c(NA,NA , 30), 
        low1=c(0, 30, NA), 
        low0=c(0, 10, NA), 
        low3=c(NA, NA, 10), 
        low4=c(NA, NA, 10), 
        high1=c(30, 90, NA), 
        high0=c(20, 50, NA), 
        high3=c(NA,NA , 60), 
        high4=c(NA,NA , 60)) 

test_data2 <- rbind(rep(NA, 3), test_data2) 

col_no2<-grep("prev", colnames(test_data2)) 

row_names2<-list(
list("Author", "You", "Asiimwe", "Martin"), 
list("Year", "1988", "1984", "1964") 
) 


prev<-with(test_data2, cbind(prev1, prev0, prev3, prev4)) 
low<-with(test_data2, cbind(low1, low0, low3, low4)) 
high <- with(test_data2, cbind(high1, high0, high3, high4)) 

forestplot(row_names2, prev, low, high, 
     #zero = c(0, 0.01), 
     fn.ci_norm = list(list(fpDrawNormalCI, fpDrawCircleCI,  fpDrawDiamondCI, fpDrawPointCI), 
         list(fpDrawNormalCI, fpDrawCircleCI, fpDrawDiamondCI, fpDrawPointCI), 
         list(fpDrawNormalCI, fpDrawCircleCI, fpDrawDiamondCI, fpDrawPointCI), 
         list(fpDrawNormalCI, fpDrawCircleCI, fpDrawDiamondCI, fpDrawPointCI) 
         ), 

     xticks = c(0,10, 20, 30, 40, 50, 60, 70, 80, 90, 100), 
     txt_gp = fpTxtGp(label = list(gpar(fontfamily = ""), 
            gpar(fontfamily = "", 
              col = "black")), 
         ticks = gpar(fontfamily = "", cex=0.8), 
         xlab = gpar(fontfamily = "", cex = 0.8)), 
     boxsize=0.15, 
     line.margin = .1, # We need to add this to avoid crowding 
     lty.ci = c(1, 2, 3, 4), 
     xlab="Prevalence (%)", 
     new_page = TRUE) 

の下に使用しても動作しますエラーが表示され、動作しません。

私が使用したコードは以下の通りです。

rm(list=ls()) 
library(forestplot) 

test_data2<-data.frame(prev1=c(0, 70, NA), 
        prev0=c(0, 20, NA), 
        prev3=c(NA,NA , 30), 
        prev4=c(NA,NA , 30), 
        low1=c(0, 30, NA), 
        low0=c(0, 10, NA), 
        low3=c(NA, NA, 10), 
        low4=c(NA, NA, 10), 
        high1=c(30, 90, NA), 
        high0=c(20, 50, NA), 
        high3=c(NA,NA , 60), 
        high4=c(NA,NA , 60)) 

test_data2 <- rbind(rep(NA, 3), test_data2) 

col_no2<-grep("prev", colnames(test_data2)) 

row_names2<-list(
list("Author", "You", "Asiimwe", "Martin"), 
list("Year", "1988", "1984", "1964") 
) 


prev<-with(test_data2, cbind(prev1, prev0, prev3, prev4)) 
low<-with(test_data2, cbind(low1, low0, low3, low4)) 
high <- with(test_data2, cbind(high1, high0, high3, high4)) 

forestplot(row_names2, prev, low, high, 
     #zero = c(0, 0.01), 
     fn.ci_norm = list(list(fpDrawNormalCI, fpDrawCircleCI,  fpDrawDiamondCI, fpDrawPointCI), 
         list(fpDrawNormalCI, fpDrawCircleCI, fpDrawDiamondCI, fpDrawPointCI), 
         list(fpDrawNormalCI, fpDrawCircleCI, fpDrawDiamondCI, fpDrawPointCI), 
         list(fpDrawNormalCI, fpDrawCircleCI, fpDrawDiamondCI, fpDrawPointCI) 
         ), 

     xticks = c(0,10, 20, 30, 40, 50, 60, 70, 80, 90, 100), 
     txt_gp = fpTxtGp(label = list(gpar(fontfamily = ""), 
            gpar(fontfamily = "", 
              col = "black")), 
         ticks = gpar(fontfamily = "", cex=0.8), 
         xlab = gpar(fontfamily = "", cex = 0.8)), 
     boxsize=0.15, 
     line.margin = .1, # We need to add this to avoid crowding 
     lty.ci = c(1, 2, 3, 4), 
     xlab="Prevalence (%)", 
     new_page = TRUE, 
     legend=c("HIV-infected", "HIV-uninfected", "All participats", "Love-you"), 
     legend_args = fpLegend(pos = list("topright"), 
           r = unit(.1, "snpc") 
          )) 

誰でもこのことについてご意見がありますか？ ありがとう！

Answer 1

これはバグであることが判明しました。すばやいデバッグを可能にした徹底的なレポートのおかげです。 1.7.1 devバージョンで修正する必要があります。devtools：

devtools::install_github("gforge/forestplot", ref="develop")