R/plot.Selectmodel.R

In FitAR: Subset AR Model Fitting

"plot.Selectmodel" <-
function(x, ...){
ans<-x #for clarity we rename object
par(mar = c(5,4,4,5)+0.1)
on.exit(par(mar=c(5, 4, 4, 2) + 0.1)) #default
if (is.list(ans)){    
    X<-Y<-a<-numeric(0)
    for (i in 1:length(ans)){
        p<-ans[[i]]$p
        X<-c(X,p)
        aic<-ans[[i]][[2]]
        Y<-c(Y,rep(aic,length(p)))
        a<-c(a,aic)
        }
    ylab<-names(ans[[1]])[2]
    plot(X,Y,xlab="lags", ylab=ylab,pch=15, cex=2)
    abline(h=a, lty=2)
    title(main=paste(attr(ans,"model"),"model selection"))
    ytic<-a
    r<-round(exp((min(Y)-ytic)/2),3)
    axis(4,at=ytic, labels=r)
    mtext(side=4, line=2.7, text="Relative Plausibility")
    }
if (is.matrix(ans)){
    X<-ans[,1]
    Y<-ans[,2]
    ylab<-dimnames(ans)[[2]][2]
    plot(X,Y,ylab=ylab,pch=15,cex=2, xlab="p")
    title(main="AR(p) model selection")
    ytic<-Y
    r<-round(exp((min(Y)-ytic)/2),3)
    axis(4,at=Y,labels=r)
    mtext(side=4, line=2.7, text="Relative Plausibility")
    }
invisible()
}