similar to: (no subject)

Hi, I am performing GEV analysis on temperature/precipitation data and want to use the ks.boot function but I am unsure of how to implement it using the pgev distribution. For example: ks.test(data,pgev,shape,location,scale) ks.boot(data,(distribution parameters?),alternative="pgev",nboots=1000) Any advice? Apologies in advance if I have used the wrong email address. Regards, Louise

Can anyone show me how to add a log Pearson type III plot using the evdistq() command to an extreme value plot using the lmom package? Attached sample code below... Thanks in advance, Dave library(lmom) # annual maximum daily streamflows Mackenzie River mackenzieRiver = c(26600, 30300, 34000, 32000, 29200, 28300, 28600, 26400, 28300, 28800, 29000, 22100, 32900, 31800, 21600, 32100, 27000,

I saw a Help e-mail related to MLE. Does R have a probability weighted method (PWM) estimator function? I can't seem to find anything on PWM, unless my eyes are playing trick on me. [[alternative HTML version deleted]]

I have a small request re. R FAQ Frequently Asked Questions on R Version 1.5-10, 2002-06-13 ISBN 3-901167-51-X Kurt Hornik In section 6.2 Should I run R from within Emacs? The faq says "Yes, definitely." This led me to install emacs, and to install ess. However, I use a windows environment (mswindows2000) and I had some difficulty installing, linking and launching these

Hi everyone, In the package POT, there is a function that computes the L-moments of a given sample (samlmu). However, to compute those L-moments, one needs to obtain the total number of combinations between two numbers, which, by the way, requires the use of a factorial. See, for example, Hosking (1990 , p. 113). How does the function "samlmu" in the package POT avoids overflows? I

Hi, I am trying the following: library(ismev) library(evd) fit <- gev.fit(x,show=FALSE) ks.test(x,pgev,fit$mle[1],fit$mle[2],fit$mle[3]) but I am getting: Warning message: cannot compute correct p-values with ties in: ks.test(x, pgev, fit$mle[1], fit$mle[2], fit$mle[3]) where x is: [1] 239 38 1 43 22 1 5 9 15 6 1 9 156 25 3 100 6 [18] 5 100

Hi, I'm a second year Master's student in Applied Statistics. I am doing a project using average weekly U.S. regular gasoline prices (in cents, per gallon) from an Excel file (from the years 1990- May 2010). I want to find the probability that the average weekly U.S. regular gasoline prices (in the long term) goes over 400 cents a gallon (or $4.00 a gallon). I am using the

My apologies in advance as this seem like an easy thing to fix. I am quite new to using R. I am doing my thesis at the moment and I am trying to use the extremes package. However, when I try to access it I get this message: The following object(s) are masked from ?package:stats?: > library( extRemes) Loading required package: tcltk Loading Tcl/Tk interface ... done Loading required package:

How to adapt this piece of code but for: - gamma distribution - 3 parameter log normal More specifically, where can I find the specification of the parameter (lmom) for pelgam() and pelln3()? Lmom package info just gives: pelgam(lmom), lelln3(lmom), where lmom is a numeric vector containing the L-moments of the distribution or of a data sample. # Draw an L-moment ratio diagram and add a line for

Dear R helpers,   I have a data of 2102 observations (consisting of 0's also), to which I am trying to fit Normal distribution using "lmom" pacakage. If I use Excel, its easy to estimate the parameters of Normal distribution as simple mean and standard devaition. The results I get if I use teh excel are as   Parameters of Normal distribution :-   Mean = 22986.44 and standard

Dera R Helpers, I am re-posting my query. Please guide me. Maithili --- On Fri, 3/13/09, Maithili Shiva <maithili_shiva at yahoo.com> wrote: I am trying to fit the Gumbel distribution to a data. I am using lmom package. I am getting problem in Cumulative Distribution Function of Gumbel distribution as I am getting it as a series of 0's and 1's thereby affecting the P P

Dear R helpers I am trying to fit the Gumbel distribution to a data. I am using lmom package. I am getting problem in Cumulative Distribution Function of Gumbel distribution as I am getting it as a series of 0's and 1's thereby affecting the P P Plot. My R code is as follows. library(quantreg) library(RODBC) library(MASS) library(actuar) library(lmom) x <-

I am having major problems with SAMBA samba-3.0alpha20 in trying to connect to Windows 2000 AD. I have attached info if that helps. Any help you can give me would be greatly appreciated. Thanks in advance Clive Elsum I can get samba-3.0alpha20 working if I include reference to our NT PDC in the smb.conf file and do a net rpc join command. This joins our NT PDC domain which has a trust

Hi R-help I am trying to mask the ocean from an image plot I have made. Here is some example code: library(mapdata) image(x=110:155, y =-40:-10, z = outer(1:45, 1:30, "+"), xlab = "lon", ylab = "lat") outline <- map("worldHires", plot=FALSE) # returns a list of x/y coords xrange <- range(outline$x, na.rm=TRUE) # get bounding box yrange

Hi All Please see the Rprofile file which i have modified as follows and after that when I start R then I see that R says to me "TRUE" for all the packages implying that all loaded at once. But when i try to use commands as simple as help("lm"), it doesnt work nor any of the menu "Packages" is not working. Although the regression using lm ( Y ~ X ) is working

Hi, Interpretation problem ! so what i did is by using the: >fit1 <- fitdist(vectNorm,"beta") Warning messages: 1: In dbeta(x, shape1, shape2, log) : NaNs produced 2: In dbeta(x, shape1, shape2, log) : NaNs produced 3: In dbeta(x, shape1, shape2, log) : NaNs produced 4: In dbeta(x, shape1, shape2, log) : NaNs produced 5: In dbeta(x, shape1, shape2, log) : NaNs produced 6: In

Respected Sir I have to select a probability distribution using L-moment Ratio Diagram. I am not understanding how to plot sample TAU3 and TUA4 on L-moment Ratio Diagram. Please Guide me Best Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. Email: amnakhan493@gmail.com amna_989@hotmail.com amna_989@yahoo.com [[alternative HTML version deleted]]

Hey, I have a set of income data which I'd like to fit to a gamma distribution. How can I estimate the two parameters of the gamma distribution for a vector, e.g. c(2039L, 2088L, 5966L, 2353L, 1966L, 2312L, 3305L, 2013L, 3376L, 3363L, 3567L, 4798L, 2032L, 1699L, 3001L, 2329L, 3944L, 2568L, 1699L, 4545L) I sense this will be a very easy one-liner, but my searching didn't come up with

Dear R helpers, How do you estimate the (Location, Scale, Shape) parameters of Generalized Extreme Value distribution using R? I have tried VGAM but just not able to write the R script. Please advise. With regards Maithili

Hi all. I have problem with performance ocfs2. I just instaled ocfs2 1.4.1 in 2 nodes cluster. I use ocfs2 in mail server. This system has large amount of small files about 50kB. Ocfs2 is formated: >> mkfs.ocfs2 -T mail -N 2 /dev/sdb1 mkfs.ocfs2 1.4.1 Cluster stack: classic o2cb Overwriting existing ocfs2 partition. Proceed (y/N): y Filesystem Type of mail Filesystem label=