R/qat_plot_lim_rule_sigma_2d.R

In qat: Quality Assurance Toolkit

qat_plot_lim_rule_sigma_2d <-
function(flagvector, filename, measurement_vector=NULL, sigma_factor=NULL, meanofvector=NaN, sdofvector=NULL, measurement_name="",directoryname="", plotstyle=NULL) {
## functionality: plot lim-static-rule
## author: André Düsterhus
## date: 03.08.2011
## version: A0.1
## input: flagvector from qat_analyse_lim_rule_static, measurement_vector, min_value, max_value, directoryname, filename, plotstylelist
## output: plot
	if (is.null(plotstyle)) {
		# if no plotstyle available, use standard plotstyle
		plotstyle<-qat_style_plot()
	}
#	library("fields")
	# set up savepath of the plot
	path <- paste(directoryname,filename,".png", sep="")
	png(filename=path,width=800,height=600, pointsize=12, bg=plotstyle$basecolor)
		par(font.lab=2,oma=c(0,0,2,0), mar=c(5.1,5.1,5.1,8), font=2)
		image(1:dim(flagvector)[1], 1:dim(flagvector)[2], flagvector, col=c(plotstyle$plotcolormain,plotstyle$plotcolorbackground,plotstyle$plotcolorminor), xlab="", ylab="", col.lab=plotstyle$fontcolor, col.main=plotstyle$fontcolor, col.sub=plotstyle$fontcolor,fg=plotstyle$frontcolor, col.axis=plotstyle$fontcolor, font.axis=2, zlim=c(-1,1))
		if(!is.na(sigma_factor)&&!is.na(meanofvector)&&!is.na(sdofvector)) {
			bordertext2 <- paste("Minimum errors: ",length(which(flagvector==-1)), " Maximum errors: ",length(which(flagvector==1)))
		} else {
			bordertext2<-""
		}
		if(!is.na(sigma_factor)) {
			bordertext<- paste("Sigma factor: ", sigma_factor, " ")
		} else {
			bordertext <- ""
		}
		if(!is.na(meanofvector)) {
			bordertext <- paste(bordertext,"Mean of vector: ", meanofvector, " ")
		} else {
			bordertext<- paste(bordertext,"")
		}
		if(!is.na(sdofvector) != 0) {
			bordertext <- paste(bordertext,"SD of vector: ", sdofvector)
		} else {
			bordertext<- paste(bordertext,"")
		}
		if (bordertext != "") {
			mtext(bordertext, side=3, line=2, font=2, col=plotstyle$fontcolor)
			mtext(bordertext2, side=3, line=1, font=2, col=plotstyle$fontcolor)

		}
		if(measurement_name != "") {
			bordertext3<-paste("Data: ",measurement_name, sep="")
			mtext(bordertext3, side=3, line=3, font=2, col=plotstyle$fontcolor)
		}
		title(main=list("LIM-Rule sigma", col=plotstyle$fontcolor), outer=TRUE)
		image.plot(1:dim(flagvector)[1], 1:dim(flagvector)[2], flagvector, col=c(plotstyle$plotcolormain,plotstyle$plotcolorbackground,plotstyle$plotcolorminor), xlab="", ylab="", col.lab=plotstyle$fontcolor, col.main=plotstyle$fontcolor, col.sub=plotstyle$fontcolor,fg=plotstyle$frontcolor, col.axis=plotstyle$fontcolor, font.axis=2, legend.only=T, axis.args=list(at=c(-1,0,1), labels=c("min", "ok", "max"), font=2), zlim=c(-1,1))
	dev.off()		
}