R/plotCrit.R

In geozoning: Zoning Methods for Spatial Data

#########################################################
#' plotCrit
#'
#' @details  reads loopQ1-5 results, filters results by keeping th best criteria ) and plots them together with corresponding costs.
#' @param m1 dataset with loopQ1 results
#' @param m2 dataset with loopQ2 results
#' @param m3 dataset with loopQ3 results
#' @param m4 dataset with loopQ4 results
#' @param m5 dataset with loopQ5 results
#' @param ONE single plot
#' @param title plot title
#'
#' @return a vector of probabilities corresponding to best results
#' @importFrom graphics plot points lines axis
#' @export
#'
#' @examples
#' # not run
#' 
plotCrit=function(m1=NULL,m2=NULL,m3=NULL,m4=NULL,m5=NULL,ONE=FALSE,title="Gaussian field simulation")
################################################################################################################
{

# 
# plot results if (max(crit)-crit) <=maxd

if (!ONE) par(mfrow=c(3,2))

maxd=0.8
# no need to give all mk arguments
q=1:5
mask=rep(TRUE,length(q))
for (k in q)
{
mk=get(paste("m",k,sep=""))
if(is.null(mk)) mask[k]=FALSE
}
# keep only non null mk
q=q[mask]
nl=length(q)+1
maxy=0
miny=Inf
best=c()

for (k in q)
{
mk=get(paste("m",k,sep=""))
mbkname=paste("mb",k,sep="")
mask=(-mk[,"crit"]+ mk[1,"crit"])<maxd & mk[,"nq"]==k

mbk=mk[mask,]
if(nrow(mbk)>5) mbk=mbk[1:6,]
assign(mbkname,mbk)
maxy=max(maxy,max(mbk[1,"crit"]))
miny=min(miny,min(mbk[,"crit"]))
}
miny=min(miny,4)

firstq=q[1]
mb1=get(paste("mb",firstq,sep=""))
plot(2:nl,rep(mb1[1,"crit"],nl-1),type="n",xlim=c(1.5,nl),ylim=c(miny,ceiling(maxy)),ylab="Criteria",xlab="Number of labels",xaxt="n",cex=2,main=title)
axis(1,at=2:nl,cex=2)
step=0.05

for (k in q)
{
mbkname=paste("mb",k,sep="")
mbk=get(mbkname)
n=nrow(mbk)
kx=k+1+seq(-n/2,n/2,length=n)*step
points(kx,mbk[,"crit"],pch=20,cex=1.5)
if (k!=q[length(q)]) lines(c(k+1.5,k+1.5),c(miny*0.8,ceiling(maxy)*1.2),lty=2)
mdb=mbk[,"cost"]
#points(kx,mdb,col="blue")
crit=mbk[,"crit"]-1

numic=(ncol(mbk)-1-k+1):(ncol(mbk)-1)
best=c(best,mbk[1,numic])
}

#legend(x=5.5,y=3,leg=c("Crit","Cost"),col=c("black","blue"),pch=c(20))

return(unlist(best))
}