r - Rの対数確率プロット？

Question

x軸が確率で、y軸が対数スケールであるRで、このような対数確率プロットを作成する方法を知っている人はいますか。heR.Miscパッケージを読んでダウンロードしましたが、使い方がわかりません。！

Answer 1

#create log probablity plot
#MPM 131201
#Make some dummy data
set.seed(21)
Dt<-as.data.frame(rlnorm(625, log(10), log(2.5)))
names(Dt)<-"Au_ppm"

#Create probablity scale lines and 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 some value grid lines then convert to 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 up lattice packages - latticeExtra for nice log scale
require(lattice)
require(latticeExtra)

#Use qqmath to make the plot (note lattice does not work for weighted data - shame about 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)
        }

    )