[R的multiplot - R的食譜錯誤的審美觀必須是一個長度......「

Question

我嘗試使用R食譜multiplot功能，使兩個地塊彼此相鄰，儘管我收到以下錯誤：

Aesthetics must either be length one, or the same length as the dataProblems:v.x

有人可以幫忙嗎？我的代碼如下：

rm(list=ls()) 
library(ggplot2) 
library(grid) 

# Multiple plot function 
multiplot <- function(..., plotlist=NULL, file, cols=1, layout=NULL) { 
    require(grid) 

    # Make a list from the ... arguments and plotlist 
    plots <- c(list(...), plotlist) 

    numPlots = length(plots) 

    # If layout is NULL, then use 'cols' to determine layout 
    if (is.null(layout)) { 
    # Make the panel 
    # ncol: Number of columns of plots 
    # nrow: Number of rows needed, calculated from # of cols 
    layout <- matrix(seq(1, cols * ceiling(numPlots/cols)), 
        ncol = cols, nrow = ceiling(numPlots/cols)) 
    } 

    if (numPlots==1) { 
    print(plots[[1]]) 

    } else { 
    # Set up the page 
    grid.newpage() 
    pushViewport(viewport(layout = grid.layout(nrow(layout), ncol(layout)))) 

    # Make each plot, in the correct location 
    for (i in 1:numPlots) { 
     # Get the i,j matrix positions of the regions that contain this subplot 
     matchidx <- as.data.frame(which(layout == i, arr.ind = TRUE)) 

     print(plots[[i]], vp = viewport(layout.pos.row = matchidx$row, 
              layout.pos.col = matchidx$col)) 
    } 
    } 
} 

# Define beta binomial function 
dbb <- function(x, N, u, v) { 
    beta(x+u, N-x+v)/beta(u,v)*(gamma(N+1)/(gamma(x+1)*gamma(N-x+1))) 
} 

#Prior parameters 
a<-1 
b<-1 

#Data 
N1<-100 
X<-10 

#New sample 
N2 <- 100 

# Generate prior, likelihood, posterior and posterior predictive distribution 
v.theta <- seq(0,1,length=1000) 
v.x <-seq(0,100,length=100) 
v.prior <- dbeta(v.theta,a,b) 
v.likelihood <- choose(N1,X)*(v.theta^X)*(1-v.theta)^(N1-X); 
v.posterior <- dbeta(v.theta,a+X,b+N1-X) 
v.posteriorPred <- dbb(v.x,N2,a+X,b+N1-X) 

# Plot as a subplot 
p1<-qplot(v.x,v.prior,geom = c("point", "line")) 
p2<-qplot(v.x,v.likelihood,geom = c("point", "line")) 

multiplot(p1, p2) 

Answer 1

問題在於矢量v.x的長度與v.prior的長度不同。這現在已經改變了！