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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.