#-- Author: Omar S. G—mez G.
#-- Fecha: Mar.2011
#-- Versiion: 1.0

PrintMeta <- function(file, meta){
	#postscript(file, width = 4, height = 5,horizontal = FALSE, onefile = FALSE, paper = "special",
	postscript(file, width = 8.5, height = 3,horizontal = FALSE, onefile = FALSE, paper = "special",
       			family = "Helvetica")

	forest.meta(meta, comb.fixed=FALSE, comb.random=TRUE, fontsize=8, plotwidth=unit(2.8,"inches"), 
				print.Q=TRUE, print.pval.Q=TRUE, just="left",xlim="s",leftcols=c("studlab", "n.e",
				"mean.e","sd.e","n.c","mean.c","sd.c"), lab.e.attach.to.col="sd.e", 
				lab.c.attach.to.col="sd.c",digits=4)

	#plot.meta(meta, comb.fixed=FALSE, comb.random=TRUE, text.random="", xlab="")

	dev.off()
 	return(TRUE)
}

PrintCC <- function(file, mean, studies, title, efect){
	postscript(file, width = 5.5, height = 5, horizontal = FALSE, onefile = FALSE, paper = "special", 
				family = "Helvetica", bg="transparent", pointsize=10)
	qcc(mean, type="xbar.one", std.dev="MR", labels=studies, title=title, ylab=efect, xlab="", 
		axes.las=2, nsigmas=3)
	dev.off()
 	return(TRUE)
}


#-----------
#Operacionalizaciones de Causa
C_Op_A <- "P.C.E."
C_Op_B <- "Cob. Dec."

#Operacionalizaciones de Efecto
E_Op <- "Porcentaje de faltas observadas"

ds 			<- read.table("validacion/p01/res_comp/exp_res_part01.txt", header=T)
ds$J		<- 1-(3/(4*(ds$n.e+ds$n.c-2)-1))
ds$Swithin	<- sqrt( ((ds$n.e-1) * ds$sd.e^2 + (ds$n.c-1) * ds$sd.c^2) / (ds$n.e+ds$n.c-2) )
ds$d		<- (ds$mean.e - ds$mean.c) / ds$Swithin
ds$g		<- ds$J * ds$d
ds$Vd		<- ( (ds$n.e + ds$n.c)/(ds$n.e * ds$n.c) ) + ds$d^2/(2*(ds$n.e + ds$n.c))
ds$Vg		<- ds$J^2*ds$Vd
ds$SEg		<- sqrt(ds$Vg)

meta01 <- metagen(ds$g, ds$SEg, sm="SMD", data=ds, studlab=paste(studlab), label.e=C_Op_A, 
				  label.c=C_Op_B, complab=paste(C_Op_A," vs. ",C_Op_B), outclab=E_Op )

meta01$n.e<-addvar(meta01,ds,"n.e")
meta01$mean.e<-addvar(meta01,ds,"mean.e")
meta01$sd.e<-addvar(meta01,ds,"sd.e")
meta01$n.c<-addvar(meta01,ds,"n.c")
meta01$mean.c<-addvar(meta01,ds,"mean.c")
meta01$sd.c<-addvar(meta01,ds,"sd.c")


PrintMeta("validacion/p01/res_comp/_meta_part01.eps", meta01)

PrintCC("validacion/p01/res_comp/_cc_a_part01.eps", ds$mean.e, ds$studlab, paste("Control Chart for ",C_Op_A), E_Op)
PrintCC("validacion/p01/res_comp/_cc_b_part01.eps", ds$mean.c, ds$studlab, paste("Control Chart for ",C_Op_B), E_Op)