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: http://r.789695.n4.nabble.com/access-the-se-of-a-forecast-tp4628847p4629088.h

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.

[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] 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 thr

[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 funct

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: http://r.789695.n4.nabble.com/pr

[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 pre

[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 pre

[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 pre

[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 pre