#======================================================================= plot.profmle <- function(mleobj) { if( (class(mleobj))[1] != "profile.mle" & class(mleobj)[1] != "profile.mle2" ) stop( "Object should have class \'profile.mle\' or \'profile.mle2\'") mleprof <- mleobj@profile npar <- length(mleprof) for(i in 1:npar) { tmp <- mleprof[i][[1]] varname <- names(mleprof[i]) y <- tmp[,1]^2/2 x <- (tmp[,2][,i]) plot(x, y, type="l", xlab=varname, ylab="Log-Verossimilhança Negativa", col="red", ) abline(h=2, col="blue", lty=2) } }