similar to: Package ismev, gpd.fit, and interpretation for statistics of extreme values

Dear all, I'm trying to produce a log-likelihood profile for a GPD estimate with the ismev package. When following the examples, everything is working fine, i.e. I get a nice parabolic log-likelihood curve. When using my own data (intensity of hurricanes at landfall), the log-likelihood curve looks like a step function. Am I doing something wrong here? Is my data not suited for the GPD fit?

Hello, I want to compute the P-value from the joint cumulative distribution of an n-dimensional order statistic in R, using the formula found on http://cmgm.stanford.edu/%7Ekimlab/multiplespecies/Supplement/methods_network.html My data consists of three different techniques (G2D, POCUS and RANDOM), and each has associated with it a number of rankings (integer between 0 and 1000), like for

Hi, I'm exploring the ismev package and in particular the pp.fit function. Documentation says that this function performs Maximum-likelihood fitting for the point process model, including generalized linear modelling of each parameter. Actually I'm afraid not to understand what it means with "Point process model". Any suggestion on possible references for further infos? thanks

Up to now, I have recognized problems with "gpd(..)", the function from the package "evir" I think that all these functions that estimate the parameters xi, beta for the GPD by given threshold mu use the function "optim(..)" ( gpd, fitgpd, ...) "Error" example: data1 <- rgpd(1000, xi= -1.5, mu=1000, beta=100) so the created poinnts take place in about

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

Greetings I am attempting to load exTremes and this is what I get... *> utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- trying URL ' http://cran.uk.r-project.org/bin/windows/contrib/2.9/extRemes_1.60.zip' Content type 'application/zip' length 676182 bytes (660 Kb) opened URL downloaded 660 Kb package 'extRemes' successfully

[This email is either empty or too large to be displayed at this time]

Hi all, I'd like to export the results of GPD or GEV analysis generated with the Extremes Toolbox for plotting in Grapher or Excel (for manipulation by our publications group). Is it possible to access the data used to generate the plots in extRemes, or do I need to code the analysis from scratch and write to file? Thanks, Ben Stay connected to the people that matter most with a

Hi, I'd like to export the results of GPD or GEV analysis generated with the Extremes Toolbox for plotting in Grapher or Excel (for manipulation by our publications group). Is it possible to access the data used to generate the plots in extRemes, or do I need to code the analysis from scratch and write to a file? Stay connected to the people that matter most with a smarter inbox.

Hello, First of all, let me apologize that my statistics background is modest at best. I am doing some extreme value analysis on model output (WRF) which have the following dimensions: speed(time,lat,lon) I am trying to fit the GPD (gpd.fit) to each point (time,lat,lon) to get a return level plot with values at each grid point. (Map with return level by location.) Here is some code I

Hi, The lib call (generated through LLVM framework) for 64 bit operation returns 64 bit value. But only the lower 32 bits are consumed. Even the users of higher part are consuming the lower 32 bit value only. Consider the case in Sparc or ARM long long lla, llb; int ia; void fun() { lla = llb << ia; } For such targets in the DAG after LegalizeTypes the higher 32

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],

Dear R-help! I need to draw contour lines in a plot of wave heights (Hs) versus peak periods (Tp) showing the joint probabilities of 1-year wave heights~peak periods, 10-year wave heights~peak periods and 100-year wave heights~peak periods. I've used the contourplot() function in the plot I've added in this mail. You can use the added dataset "Rhelpdata.txt" to reproduce a

Hola, De esta manera sí funciona: #----------------------------------------------- #................... vec1a <- as.vector(1:3) n1a <- 3 n2a <- 3 #P1a <- matrix(rnorm(30),nrow=3) #P2a <- matrix(rnorm(30),nrow=3) Ka <- 10 n.val <- 30 # En forma de funcion producto <- function(n1a,n2a,Ka,vec1a,n.val){ # A eliminar una vez compruebes que la vectorización funciona

Hi Yan, Before I try to answer your question: please send R questions to the r- help mailing list (cc'd here). You will likely get your answer more quickly (since I'm not the only one looking at it) and it might help other people who get stuck on a similar situation by finding the answer on a public forum (just like you did). On Sep 8, 2009, at 5:44 AM, Yan Hui wrote: > Dear

Alright, going ot ask this question in advance of doing bone headed tests. I have one application that does not understand the idea of /usr/local/ and insists on populating /usr. Of course it is closed source and commercial and it is the one they want to use. Below is the list of junk it shoves in /usr/share. Between the Type Reference guide and another email , I believe I can move this tree

rpcclient localhost -d 4 -U'root%Passwd' -c 'adddriver "Windows NT x86"

Hello, I am a new user of R. I am trying to use the packages fBasics and fExtremes when i am running the examples I get few error. Could someone tell me what is happenig? Thank you beforehand. from Fbasics packages: xmpfBasics() Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file '/usr/lib/R/library/fBasics/demoIndex'

Hola a tod en s, tengo que hacer una operación con matrices que lleva un doble bucle. He intentado vectorizarlo pero sin mucho éxito con la función "Vectorize". ¿sabríais de alguna manera de evitar los bucles o de que funcionase Vectorize? Adjunto un ejemplo #................... vec1a <- as.vector(1:3) n1a <- 3 n2a <- 3 P1a <- matrix(rnorm(30),nrow=3) P2a <-

Hi, I am trying to migrate print queues from a windows server to a samba share. 1. I ran printmig.exe on the windows server. 2. Backed up the printers on the windows servers to a cab file. 3. Tries to restore the printers to the sambashare by specifying the target as //sambaShare. But I get the following error: 2006:01:25 15:35:38 Access Granted to: \\sambaShare 2006:01:25 15:35:38