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. !enter image description here

#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

Popular posts from this blog

c# - SVN Error : "svnadmin: E205000: Too many arguments" -

c# - Copy ObservableCollection to another ObservableCollection -

All overlapping substrings matching a java regex -