R/s3.R

Defines functions multiplot print.mmf summary.mmf plot.mmf

multiplot <- function(..., plotlist=NULL, file, cols=1, layout=NULL) {
  # 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))
    }
  }
}

print.mmf=function(x,...){
  summary.mmf(x)
  plot.mmf(x)
}

summary.mmf=function(object,...){
  name = names(object$thetahat)
  if(is.null(name)){
    for(i in 1:length(object[[1]])){
      print(paste0("The point estimate for the ", i, "th parameter is ", round(object[[1]][i],3),
                   ", the standard error is ", round(object[[2]][i],3)))
    }
  }else{
    for(i in 1:length(object[[1]])){
      print(paste0("The point estimate for ", name[i], " is ", round(object[[1]][i],3),
                   ", the standard error is ", round(object[[2]][i],3)))
    }
  }
  
}

plot.mmf=function(x,...){
  multiplot(x[[3]],x[[4]],cols = 2)
}
qiwei-li/mmfit documentation built on May 26, 2019, 12:33 p.m.