[R] access the se of a forecast

2012-05-12 Thread plocq
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-dshw(eem,period1=7,period2=48,h=48) then by using summary(fit), I

Re: [R] access the se of a forecast

2012-05-12 Thread plocq
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 typing sfit-summary(fit) and then sfit$s20, I get NULL... these facts seem very strange to me.

Re: [R] access the se of a forecast

2012-05-12 Thread plocq
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:

Re: [R] modify the name of axis of an R function

2011-12-15 Thread plocq
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 with it : the fact is that I have three plots that are called by the function. I can specify my three

[R] modify the name of axis of an R function

2011-12-14 Thread plocq
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

Re: [R] profile likelihood

2011-12-10 Thread plocq
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) -- View this message in context:

[R] profile likelihood

2011-12-08 Thread plocq
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 profile(fit), I obtain the following error (even if I precise

[R] profile likelihood

2011-12-08 Thread plocq
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 profile(fit), I obtain the following error (even if I precise

[R] profile likelihood

2011-12-08 Thread plocq
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 profile(fit), I obtain the following error (even if I precise

[R] profile likelihood

2011-12-08 Thread plocq
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 profile(fit), I obtain the following error (even if I precise