similar to: fgev + negative location parameter for positive data?

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

hi, rgev function gives me random deviates and I have a data set which I am fitting to an EVD,IS there a way I can plot both observed and ideal evd on the same plot thankyou Rangesh

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 everyone, is there a function to fit a frechet distribution? The only thing I found is gev.fit from ismev which fits a generalized extreme value distribution (if shape>1 => Frechet) . Is there a function to only fit a frechet? Thank you Benjamin

hello everybody i want to use the maximum likelihood method to estimate FRECHET parameters of my sample data. Should it work with fitdistr in the package MASS? I only find how to do it for GEV, Gumbel, and almost all other distributions, but FRECHET? I would be very happy if somebody can tell me how to do fit the FRECHET distribution! Thanks Nadja Riedwyl

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

Hi guys, I'm trying to use lmomco package. first I did the manual calculation on what is the estimates scale and location parameter given L-CV=0.2, L1=1000 L-moments and k (shape parameter) =- 0.1. so what i get is: location: 821.0445 scale: 260.7590 shape: -0.1000 #I assign this as GEV vectors using vec2par GEVpara2<-vec2par(c( 821.0445 , 260.7590 ,-0.1),'gev') #then I

I am estimating a large model by groups. How do you save the results and?returns the associated quantiles? For this example I need a data frame n?? ?xi??????? mu????????beta 1?? 0.1033614? 2.5389580 0.9092611 2? ?0.3401922? 0.5192882 1.5290615 3?? 0.5130798? 0.5668308 1.2105666 I also want to apply gevrlevelPlot() for each "n" or group. ? #Example n <- c(1, 1, 1, 1, 1, 1, 2, 2, 2,

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   > head(dat,20) A B C D E F G H I J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5 45.6 4 50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7 38.3 5 39.3 30.6 46.9 23.8 25.8

i need HELPPP!! how do i calculate the RMSE value for two GEV models?first GEV is where the three parameters are constant.2nd GEV model a 4 parameter model with the location parameter is allowed to vary linearly with respect to time while holding the other parameters at constant. is there any programming code for this? i really really need help. please reply to me as soon as possible. thanks in

hello, i'am is new in R software.i have try to make a function but it can't give what it should.i dunno what have to do next. Can somebody help me to solve it.i'll very appreciate... ##GEV simulation(Non-stationary) dsim<-function(n, alpha, beta,sca,sha){ t <- 1:n location <- alpha + beta*t inv.df<-function(x) location + -(sca/sha)+(sca/sha)*(-log(x))^(-sha) u<-runif(n)

Hello, and please advise regarding any errors/omissions on my part. However, the documentation for R's tcplot function in the evd package appears to contain an error. I am using evd version 2.2-4 in R 2.14.0 with Windows XP. > data(portpirie) > mrlplot(portpirie) ## No Error > tlim <- c(3.6, 4.2) > tcplot(portpirie, tlim) ## Line from documentation Error in fpot(data, u[1],

Hi, I am using the clusters function in the evd package and was wondering if anyone knew how to set the row names as date/time objects. My data has 1 column of date and time and another of wave height (H). My two options are: 1 - When in import the data using read.table make sure that the row.names is set = 1 giving row names to the H data. 2 - Import H from the data file and set the row names

Hi All, I am a newcomer to R. Could anyone explain me how to define link functions for either mu/sigma to allow for quadratic trends in the same, when fitting non-stationary GEV distributions? Thanks -- View this message in context: http://www.nabble.com/QUADRATIC-TREND-FOR-LINK-FUNCTIONS-ON-NON-STATIONARY-GEV-tp23360751p23360751.html Sent from the R help mailing list archive at Nabble.com.

Hi, I'm new to R, and am trying to embed R into another application. I'm calling gev.fit() from the ismev package, and it is crashing somewhere inside it. gdb is not catching it, and valgrind is not showing any memory corruption issues. I suspect it's memory corruption, because it doesn't crash in exactly the same spot each time. I'm running R 2.12.2 on a 64 bit linux (Ubuntu

Hi, I am using the clusters function in the evd package in order to determine storm events from a wave time series. So far I have the code working as I want it for wave height on its own but I would now like to include the period as well. The input data is in the form of: H t H t H t so every height measurement has a corresponding period. The storms are defined when the wave height exceeds a

Hello, I am using the lmomco package (lmom.ub and pargev) to compute the GEV parameters (location, scale, and shape), which are used to estimate return values. I was wondering how/if I can calculate upper and lower confidence (CI_u, CI_l) intervals for each return frequency using the GEV parameters to fill-in the table below? Xi (location) = 35.396 Alpha (scale) = 1.726 Kappa (shape) =

Dear All- I am looking at the bivariate point process model for extreme values as described in Section 8.3.2 of Coles book. I am stuck at a very low level in being unable to reproduce the initial step of transforming the data to unit Frechet margins. In particular, I am unable to reproduce the wave-surge pairs distribution after transformation to Frechet scale (Figure 8.6 of Coles

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,

Dear R helpers, I need to use KS and AD test for Generalized Pareto and Generalized extreme value. E.g. if I need to use KS for Weibull, I have teh syntax ks.test(x.wei,"pweibull", shape=2,scale=1) Similarly, for AD I use ad.test(x, distr.fun, ...) My problem is fir given data, I have estimated the parameters of GPD and GEV using lmom. But I am not able to find out the distribution