Dear Sir I am getting errors in using following R code for L-moment ratio
diagram
Help in this regard

> "plotlmrdia" <-
+ function(lmr,
+          nopoints=FALSE,
+         nolines=FALSE,
+         nolimits=FALSE,
+         nogev=FALSE,
+         noglo=FALSE,
+         nogpa=FALSE,
+         nope3=FALSE,
+         nogno=FALSE,
+         noexp=FALSE,
+         nonor=FALSE,
+         nogum=FALSE,
+         nouni=FALSE,
+          xlab="L-SKEWNESS",
+          ylab="L-KURTOSIS") {
+    plot(lmr$limits, xlab = xlab, ylab = ylab, type = "n",font.lab=2)
+    if(nolimits == FALSE) {
+ lines(lmr$limits,lwd=2.5,col=8)
+    }
+    if(nolines == FALSE) {
+      if(nogev == FALSE) lines(lmr$gev, col=1,lty=2,lwd=2)
+      if(noglo == FALSE) lines(lmr$glo, col=6,lwd=2)
+      if(nogno == FALSE) lines(lmr$gno, col=4, lty=2,lwd=2)
+      if(nogpa == FALSE) lines(lmr$gpa, col=2,lty=2,lwd=2)
+      if(nope3 == FALSE) lines(lmr$pe3, col=3,lwd=2)
+    }
+    if(nopoints == FALSE) {
+      if(noexp == FALSE) points(lmr$exp,pch=16,col=1)
+      if(nonor == FALSE) points(lmr$nor,pch=15,col=1)
+      if(nogum == FALSE) points(lmr$gum,pch=17,col=1)
+      if(nouni == FALSE) points(lmr$uniform,pch=18,cex=1.5,col=1)}}
>
> png("plot.png",width=480,height=400)
> par(mar=c(5,4,4,8))
>
> plotlmrdia(lmr, nopoints=FALSE, nolines=FALSE, nolimits=FALSE,
+ nogev=FALSE, noglo=FALSE, nogpa=FALSE, nope3=FALSE,
+ nogno=FALSE, noexp=FALSE, nonor=FALSE, nogum=FALSE,
+ nouni=FALSE)
Error in plot(lmr$limits, xlab = xlab, ylab = ylab, type = "n", font.lab =
2) :
  object "lmr" not found
>  grid()
Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) :
  plot.new has not been called yet
> par(xpd=TRUE)
> legend(1.15,1.05,c("LmL*","GEV","GLO","GNO","GPA","PE3","EXP","GUM","UNI","N"),
lty=c(1,2,1,2,2,1,NA,NA,NA,NA),col=c(8,1,6,4,2,3))
Error in rect(left, top, r, b, angle = angle, density = density, ...) :
  plot.new has not been called yet
> points(rep(1.3,4),c(0.32,0.22,0.12,0.02),pch=c(16,17,18,15),col=1)
Error in plot.xy(xy.coords(x, y), type = type, ...) :
  plot.new has not been called yet
> text(1.6,-0.2,"*L-moment Limits")
Error in text.default(1.6, -0.2, "*L-moment Limits") :
  plot.new has not been called yet
> par(xpd=FALSE,mar=c(5,4,4,2))
> dev.off()
windows
      2
THANK YOU

-- 
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.

	[[alternative HTML version deleted]]

