Search Archives

Thank you for your response. In fact, it seems not to be adapted to the Poisson model... : library(evd) uvdata <- rgpd(100, loc = 0, scale = 1.1, shape = 0.2) M1 <- fpot(uvdata, 1, model="pp") library(SpatialExtremes) profile(M1...

Thank you very much for all your tips! It works perfectly know and I begin to understand how R is working. I will now see how to construct this confident interval. Have a nice evening! -- View this message in context...

Thank you very much for your answer. I didn't know this. However, I tried what you suggest and it doesn't work either : - by typing summary(fit)$s20, I get exactly what I had with summary(fit) and - by...

Hi, I try to use the function profile() of the SpatialExtremes' package to obtain the profile likelihood of parameters for an extreme values fit based on Poisson process : fit<-fpot(data, threshold, model="pp", npp=365). But when I call...

Hi, I try to use the function profile() of the SpatialExtremes' package to obtain the profile likelihood of parameters for an extreme values fit based on Poisson process : fit<-fpot(data, threshold, model="pp", npp=365). But when I call...

Hi, I try to use the function profile() of the SpatialExtremes' package to obtain the profile likelihood of parameters for an extreme values fit based on Poisson process : fit<-fpot(data, threshold, model="pp", npp=365). But when I call...

Hi, I try to use the function profile() of the SpatialExtremes' package to obtain the profile likelihood of parameters for an extreme values fit based on Poisson process : fit<-fpot(data, threshold, model="pp", npp=365). But when I call...

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

Hi everybody, I am currently trying to forecast some double seasonal time series by using the function dshw. I want to access the standard errors to build the confident interval for my forecast. I am using to following code : fit...

Thanks you very much! The way plot(fit, main="main title", xlab="X-axis lable", ylab="y-axis label") seems to work quite well, I didn't notice that I could do this. I have in fact one more problem...