log probability plot in R? -
does know how create log probability plot 1 in r x-axis probability , y-axis in log-scale. read , downloaded package her.misc package don't know how use it. !
#create log probablity plot #mpm 131201 #make dummy data set.seed(21) dt<-as.data.frame(rlnorm(625, log(10), log(2.5))) names(dt)<-"au_ppm" #create probablity scale lines , associated labels - prbgrd <- qnorm(c(0.001,0.01, 0.05, 0.10,0.20,0.30,0.40, 0.50, 0.60, 0.70,0.80,0.90,0.95,0.99,0.999)) prbgrdl<-c("0.1","1","5","10","20","30","40","50","60","70","80","90","95","99","99.9") #create value grid lines convert logs valgrd<-c(seq(0.001,0.01,0.001),seq(0.01,0.1,0.01),seq(0.1,1,0.1),seq(1,10,1),seq(10,100,10)) valgrd<-log10(valgrd) #load lattice packages - latticeextra nice log scale require(lattice) require(latticeextra) #use qqmath make plot (note lattice not work weighted data - shame that) qqmath(~ au_ppm, data= dt, distribution = function(p) qnorm(p), main = "normal probablity / log (base 10) plot", pch=20, cex=0.5, xlab="normal distribution scale (%)", scales=list(y=list(log=10,alternating=1),x = list(at = prbgrd, labels = prbgrdl, cex = 0.8)), yscale.components=yscale.components.log10ticks, panel=function(x,...){ panel.abline(v=prbgrd ,col="grey",lty=3) panel.abline(h=valgrd,col="grey",lty=3) panel.qqmath(x,distribution=qnorm) } )
Comments
Post a Comment