search for: evd

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], model = model, cmax = cmax, r =...

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...

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 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. T...

As I tried to draw random sample from Hüsler-Reiss density with the following code, library(evd) margins <- cbind(0, 1, seq(-0.5,0.5,0.1)) x <- rbvevd(101, dep = 1.7, model = "hr", mar1 = margins) [,1] [,2] [1,] -0.56662298 0.8448505 [2,] 0.55824918 1.0217279 [3,] 0.24741124 0.6684668 [4,] -0.31547985 0.1781680 [5,] 0.69466134 0.2299363 [6,...

...I hope you can help me: I?ve got a dataframe (df) within which I am looking for Peak Over Threshold values as well as the length of the events. An event starts when walevel equals 5.8 and it should end when walevel equals the lower threshold value (5.35). I tried ?clusters (?)? from ?evd package?, and varied r (see example) but it did not work for all events (again see example). walevel <- c(5.75, 5.75, 5.75, 5.76, 5.80, 5.82, 5.85, 5.86, 5.91, 5.91, 5.88, 5.92, 5.99, 6.02, 6.00, 6.08, 6.11, 6.10, 6.10, 6.11, 6.14, 6.12, 6.15, 6.17, 6.15, 6.08, 6.01, 5.95, 5.89, 5.8...

Hello, I would like to sample from a Gumbell (minimum) distribution. I have installed package {evd} but the Gumbell functions there appear to refer to the maximum case. Unfortunately, setting the scale parameter negative does not appear to work. Is there a separate package for the Gumbell minimum? -- _____________________________ Rich Gwozdz Fire and Mountain Ecology Lab College of Forest...

Hi, I use the function fpot of packages evd. If I call the fit that I obtain "fit", I want to modify the name of the axis and the main title that is produced by plot(fit1). To do this, I want to create a new function where I would have the names modified. The problem is that I can't find the function that produce these plots......

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 1...

...changes axis label size). What I would like to change is the label that says: "Residuals vs. Fitted" and "Q-Q normal". So what other cex commands are there? THX for any help. Bernhard Müller ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Eidgenössisches Volkswirtschaftsdepartement EVD Forschungsanstalt Agroscope Reckenholz-Tänikon ART Ökologische Landbausysteme Reckenholzstrasse 191, CH-8046 Zürich Tel. ++41 (0)44 / 377 74 47 FAX ++41 (0)44 / 377 72 01 mailto:bernhard.mueller@art.admin.ch http://www.reckenholz.ch [[alternative HTML version deleted]]

hello anyone faced this error before? dtrace: processing aborted: Abort due to systemic unresponsiveness I used iosnoop -evD > iosnoop.log thanks. __________________________________________________________ Sent from Yahoo! Mail - a smarter inbox http://uk.mail.yahoo.com

...notypes, so the pairwise comparison is not that usefull and clear to look at. How would I do that? Is there an other test that does give grouping as output? Thnaks for your help. Bernhard Müller, Master Student ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Eidgenössisches Volkswirtschaftsdepartement EVD Forschungsanstalt Agroscope Reckenholz-Tänikon ART Ökologische Landbausysteme Reckenholzstrasse 191, CH-8046 Zürich Tel. ++41 (0)44 / 377 74 47 FAX ++41 (0)44 / 377 72 01 mailto:bernhard.mueller@art.admin.ch http://www.reckenholz.ch [[alternative HTML version deleted]]

I am finding real trouble in installing evd library in R for linux -- View this message in context: http://www.nabble.com/library-instal-tp19453453p19453453.html Sent from the R help mailing list archive at Nabble.com.

Dear R Helpers Which package is available for estimatine the parameters of three parameter FRECHET distribution. Also, how to generate the random numbers for Frechet using these three estimated parameters. Thanking in advance Maithili

...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 extRemes program (in R), and I've already uploaded my data, and found the probability using the fitted GEV distribution (I used the pgev function from the evd package). I want to compare the results from the GEV distribution to those of the GPD distribution. I've already fitted the GPD distribution to my data (at a threshold of 400 cents, where there are 52 entries a year (for most years, at least)) and have obtained the supporting graphs. M...

Hi R lovers! Are there any means to define and plot a function given the equation that specifies the function? For example I'd like to plot and work with the Gumbel Distribution density defined by Lambda(x)=exp(-exp(-x)) My question may appear very simple but I haven't got an idea yet about how to do that. I could plot something with x a vector/set of value but I don't know how to

...statistics, distributed computing, and bioinformatics. Contents: Time Series in R 1.5.0 Naive Time Series Forecasting Methods Rmpi: Parallel Statistical Computing in R The grid Graphics Package Lattice Programmer's Niche geoRglm: A Package for Generalised Linear Spatial Models Querying PubMed evd: Extreme Value Distributions ipred: Improved Predictors Changes in R Changes on CRAN Upcoming Events which makes a total of 44 pages. For the editorial board, Kurt Hornik -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-announce mailing list -- Read http://www.ci...

...statistics, distributed computing, and bioinformatics. Contents: Time Series in R 1.5.0 Naive Time Series Forecasting Methods Rmpi: Parallel Statistical Computing in R The grid Graphics Package Lattice Programmer's Niche geoRglm: A Package for Generalised Linear Spatial Models Querying PubMed evd: Extreme Value Distributions ipred: Improved Predictors Changes in R Changes on CRAN Upcoming Events which makes a total of 44 pages. For the editorial board, Kurt Hornik -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-announce mailing list -- Read http://www.ci...

Dear useRs, With the following set of command i managed to create the copula density ================================================ >library(MASS) >library(evd) >X<-matrix(sample(1:3000), ncol=2) >U=cbind(rank(X[,1])/(nrow(X)+1),rank(X[,2])/(nrow(X)+1)) > mat1=kde2d(U[,1],U[,2],n=150) image(mat1$x,mat1$y,mat1$z,col= rev(heat.colors(1500)),xlab="",ylab="") contour(mat1$x,mat1$y,mat1$z,add= TRUE,levels = pretty(c(0,4), 11)...

...9;, 'DAAG', 'date',   'deldir', 'descr', 'deSolve', 'devtools', 'digest', 'diptest', 'doBy', 'DoE.wrapper', 'e1071', 'effects',   'ENmisc', 'epiR', 'eRm', 'evaluate', 'evd', 'FactoMineR', 'fArma', 'fAssets', 'fBasics', 'fdrtool',   'fExoticOptions', 'fExtremes', 'fGarch', 'fields', 'flexclust', 'fMultivar', 'fNonlinear', 'fOptions',   'forecast', 'f...