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]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
