[R] forecast.lm() and NEWDATA

PIKAL Petr petr.pikal at precheza.cz
Tue Nov 5 14:12:03 CET 2013


Hi

> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of Ryan
> Sent: Monday, November 04, 2013 4:19 PM
> To: David Winsemius
> Cc: r-help at r-project.org
> Subject: Re: [R] forecast.lm() and NEWDATA
> 
> Hi David (and everyone)
> 
> Thank you for your reply.
> 
> I see I copied down the code wrong for the regression. I was using the
> "+" seperator, and not ",".  The regression was working.
> 
> I made a mistake with the NEWDATA, where I also used "+", instead of
> ",", however I see both work but return very different results. I have
> corrected the mistake.
> 
> However, I am still receiving the same error when I try to run the
> forecast.lm.
> "
> 
> forecast.lm(REGGY, h=5)
> 
> i receive the following error
> "Error in as.data.frame(newdata) :
>   argument "newdata" is missing, with no default""
> 

Where is forecast.lm from? Anyway, it probably requires an argument for newdata so maybe

forecast.lm(REGGY, newdata=somedata, h=5)

will work.

Regards
Petr

> 
> 
> I am sorry to say, I don't understand your reference to the string
> output when you mentioned
> 
> "
> 
> You should post the output of str() on the data-objects that has the 12
> variables and if it was modified the data argument pasted to `lm()`
> when you made REGGY.
> 
> "
>   As of right now, NEWDATA includes all the independent variables from
> the REGGY Regression.
> 
> 
> Should the NEWDATA data.frame include the dependent variable?
> 
> Also, is it necessary to reformat all variables as time series? (I have
> not currently done so)
> 
> Regards
> Ryan
> 
> On 2013-11-01 07:28 PM, David Winsemius wrote:
> > On Nov 1, 2013, at 6:50 AM, Ryan wrote:
> >
> >> Good day all.
> >>
> >> I am hoping you can help me (and I did this right). I've been
> working in R for a week now, and have encountered a problem with
> forecast.lm().
> >>
> >> I have a list of 12 variables, all type = double, with 15 data
> entries.
> >> (I imported them from tab delimited text files, and then formatted
> as.numeric to change from list to double)
> >> (I understand that this leaves me rather limited in my degrees of
> freedom, but working with what I have, sadly. )
> >>
> >> I have a LM model, such that
> >> REGGY = lm(formula=Y~A,B,C,...,I,J)
> > This looks wrong. Separating independent predictors with commas would
> be highly unusual.
> >> which I am happy with.
> >>
> >> I have
> >> NEWDATA = data.frame(A+B+C+D....+I+J)
> > This also looks wrong. Separating arguments to data.frame with "+"-
> signs is surely wrong.
> >> When i try to run
> >>
> >> forecast.lm(REGGY, h=5)
> >>
> >> i receive the following error
> >> "Error in as.data.frame(newdata) :
> >>   argument "newdata" is missing, with no default"
> > If your code prior to calling forecast on the REGGY-object was really
> what you showed here, I am not surprised. You should post the output of
> str() on the data-objects that has the 12 variables and if it was
> modified the data argument pasted to `lm()` when you made REGGY.
> (Beginners should name their data arguments.)
> >
> >> When I run
> >> forecast.lm(REGGY, NEWDATA, h=5)
> >> I receive the confidence intervals of the 15 data entries I already
> possess. I understand that by including NEWDATA, the "h=5" is ignored,
> but without NEWDATA, I receive the error message.
> >>
> >> Can anyone help me please?
> >>
> >> Regards
> >> Ryan
> >>
> >> P.S The forecast is trying to predict the next 5 values for Y from
> the regression model pasted above. I'm a bit rusty with regressions,
> but I think I've covered my bases as well as I can, and from what I
> understand of the R code, I'm following the right steps.
> > Not if what you posted here was your code. I think you missed a few
> crucials points about R syntax.
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list