Re: [R] lmom and lmomRFA Upper and Lower Bounds Simulation Questions

' is likely to be better. For inference about 'zz_gev', 'CI_gev$Qhat' is likely to be better. J. R. M. Hosking jrmhosk...@gmail.com __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] Best Distribution

Gumbel fit evdistq(qualn3, pelln3(samlmu(Prec), bound=0), col='blue') # lognormal legend(topleft,c(Gumbel,lognormal),lty=1,col=c(green,blue)) J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] R 14.0, ggplot2: could not find function initRefFields

is an implementation of the Grammar of Graphics. It seems likely that whatever you were hoping to achieve with ggplot2 can also be done natively in SPSS. J. R. M. Hosking plugin requires R 14.0 (14.2 or 15 does not work). When I install R it gives an warning: library(ggplot2) Warning message: package

Re: [R] GEV distribution fitted by L-moment graph

on this? Thanks. [[alternative HTML version deleted]] The help file for function evplot() in package lmom has an example that does exactly what you ask. J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] Adding gamma and 3-parameter log normal distributions to L-moments ratio diagram lmrd()

of the 3-parameter lognormal distribution). Is this what you want? If not, I don't know what else you want to know about the specification of the parameter (lmom). J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] Gamma distribution parameter estimation

, by fitting a 3-parameter gamma distribution (Pearson type III in the terminology of package lmom), gives a visually much better fit: evdistq(quape3, pelpe3(samlmu(x)), col='blue') J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] covariance matrix TL-moments

(2007). Some theory and practical uses of trimmed L-moments. Journal of Statistical Planning and Inference, 137, 3024-3039. J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Strange space characters in character strings

spaces. In the original HTML table they are defined as non breaking spaces i.e. nbsp; So my question is WHAT ARE THEY? Is there a way to show the binary (hex) values of these characters? charToRaw(...) will show them gsub([[:space:]], , ...) may remove them J. R. M. Hosking Here is my

Re: [R] lmomRFA-package: regsimq()

will see what I can do. J. R. M. Hosking 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

Re: [R] lmomRFA package: error bounds/confidence intervals

to be a way of getting as close as possible to what confidence intervals would be (or ought to be, given consideration (b) above). J. R. M. Hosking Thank you in advance, Tonja __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] is there a function to find the quantile of the mean of a vector?

((x.sort - x.mean)^2) / x.length Another user pointed out a function suited for this purpose, findInterval() p.mean - findInterval(x.mean, xsort) / x.length Thanks for your help, David No need to sort: x - runif(1000) p.mean - mean(x = mean(x)) J. R. M. Hosking On Thu, Jun 17, 2010

Re: [R] lmom: plotting log Pearson Type III

of the LPE3, then use it ... quaLPE3 - function(f, para) exp(quape3(f, para)) evdistq(quaLPE3, parLPE3, col='magenta') (and thank you for providing such clear sample code) J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] tapply for function taking of 1 argument?

as ? is to tapply ) Thanks for your help. coef(lm(data ~ -1 + as.factor(groups), weights=weights)) Not the fastest, but IMO more comprehensible than the constructions involving anonymous functions. J. R. M. Hosking __ R-help@r-project.org mailing list

Re: [R] Problems with fitdistr

supplied by optim Any help or suggestions are most welcomed Use function pelwei() in package lmom. J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] estimate inverse gaussian in R

in package lmom. J. R. M. Hosking __ 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.

Re: [R] lmomco package and confidence limits?

distribution is the one fitted to the actual data sim - regsimq(rfit$qfunc, nrec = rdat$n, f = 1 - 1 / c(2,5,10,25,50,100)) # Compute error bounds for quantiles of the site's # frequency distribution sitequantbounds(sim, rfit) J. R. M. Hosking

Re: [R] Log logistic Distribution - Parameter estimation

parameter Log Logistic distribution to this data? Thanking in advance Maithili Log logistic is, after reparametrization, a subset of the generalized logistic distribution used in package lmom, so function pelglo in that package, together with a bit of algebra, should get you there. J. R. M

Re: [R] lmom - Estimating Normal Distribution Parameters using lmom package

inappropriate, they can be very far apart. Your data, which have sample skewness 22 and sample L-skewness 0.93, fall into this category. J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] RDCOMClient: how to close Excel process?

- COMCreate(Excel.Application) wk - xl$Workbooks() sh-wk$Open(normalizePath(sample_file.xls))$Sheets()$Count() wk$Close() xl$Quit() rm(sh) rm(wk) rm(xl) gc() J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] threshold distribution

0.80008 1.11947 1.09484 0.81494 0.68696 + 0.82364 0.84390 0.71402 0.80293 1.02873 + )) Read 39 items y - (x-mean(x))/sd(x) library(lmom) pelln3(samlmu(y)) zeta mu sigma -1.5362134 0.2554631 0.5896735 J. R. M. Hosking __ R-help

Re: [R] Estimating paramters of 3 Parameter Gamma Distribution

-parameter gamma. For the relationship between the two distributions, see the help for function cdfpe3 in package lmom. J. R. M. Hosking Please guide Regards Maithili __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] Fw: Estimating Parameters of Weibull and Pareto distribution using LMOM package

) xialphak 1.0072.993343-0.838561 pelwei offers similar options for the Weibull distribution. J. R. M. Hosking Please help me. Regards and thanking in advance Maithili __ R-help@r-project.org

Re: [R] Fw: Fitting GUMBEL Distribution - CDF function and P P Plot

) exp(-exp(-(x - para[1])/para[2])) J. R. M. Hosking Maithili Shiva wrote: Dera R Helpers, I am re-posting my query. Please guide me. Maithili --- On Fri, 3/13/09, Maithili Shiva maithili_sh...@yahoo.com wrote: I am trying to fit the Gumbel distribution to a data. I am using lmom

Re: [R] Parameter Estimation - Generalized Extreme Value Distribution

Maithili Shiva wrote: Dear R helpers, How do you estimate the (Location, Scale, Shape) parameters of Generalized Extreme Value distribution using R? ... Package lmom, function pelgev. J. R. M. Hosking __ R-help@r-project.org mailing list https

Re: [R] Problems with NA's

errors. My advice is to check the matching of parentheses in the line if (y [3] != NA){(print (no)} J. R. M. Hosking __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R