Re: [R] Odp: Help me with prediction in linear model
Hi r-help-boun...@r-project.org napsal dne 24.07.2010 11:48:11: Thanks Murphy and pikal, I need another help,for fitting first fourier transformation ,i used following thing .Please advise on this beer_monthl has 400+ records EXample: head(beer_monthly) beer 1 93.2 2 96.0 3 95.2 4 77.1 5 70.9 6 64.8 time-seq(1956,1995.2,length=length(beer_monthly)) Even this is not working time-seq(1956,1995.2,length=length(beer_monthly)) Error in seq.default(1956, 1995.2, length = length(beer_monthly)) : object 'beer_monthly' not found sin.t-sin(2*pi*time) cos.t-cos(2*pi*time) beer_fit_fourier=lm(beer_monthly[,1]~poly(time,2)+sin.t+cos.t) #this is not working gives me result without error beer_fit_fourier=lm(beer_monthly[,1]~time+time2+sin.t+cos.t) #it is working #prediction is not workinng tpred_four - data.frame(time = seq(1995, 1998, length = 20)) predict(beer_fit_fourier, newdata = tpred_four) You'd rather read some documentation about lm and predict. You do not have variables sin.t and cos.t in your newdata. Please do read docs, posting guide and provide reproducible example. Is there any way to fit first fourier frequency , Please assist. Thanks in advance -- View this message in context: http://r.789695.n4.nabble.com/Help-me-with- prediction-in-linear-model-tp2297313p2300991.html Sent from the R help mailing list archive at Nabble.com. __ R-help@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. __ R-help@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.
Re: [R] Odp: Help me with prediction in linear model
Thanks Murphy and pikal, I need another help,for fitting first fourier transformation ,i used following thing .Please advise on this beer_monthl has 400+ records EXample: head(beer_monthly) beer 1 93.2 2 96.0 3 95.2 4 77.1 5 70.9 6 64.8 time-seq(1956,1995.2,length=length(beer_monthly)) sin.t-sin(2*pi*time) cos.t-cos(2*pi*time) beer_fit_fourier=lm(beer_monthly[,1]~poly(time,2)+sin.t+cos.t) #this is not working beer_fit_fourier=lm(beer_monthly[,1]~time+time2+sin.t+cos.t) #it is working #prediction is not workinng tpred_four - data.frame(time = seq(1995, 1998, length = 20)) predict(beer_fit_fourier, newdata = tpred_four) Is there any way to fit first fourier frequency , Please assist. Thanks in advance -- View this message in context: http://r.789695.n4.nabble.com/Help-me-with-prediction-in-linear-model-tp2297313p2300991.html Sent from the R help mailing list archive at Nabble.com. __ R-help@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.
[R] Odp: Help me with prediction in linear model
Hi r-help-boun...@r-project.org napsal dne 21.07.2010 17:46:02: Hi R-community, I have the code as follows,i Fitted model as follows lbeer-log(beer_monthly) t-seq(1956,1995.2,length=length(beer_monthly)) #beer_monthly contains 400+ entries not needed t2=t^2 beer_fit_parabola=lm(lbeer~t+t2) beer_fit_parabola=lm(lbeer~t+I(t^2)) you can use poly(t,2) but coefficients are not directly applicable in equation. Below is not working for me. Please help me in preparing the new data set for the below prediction predict(beer_fit_parabola,newdata=data.frame(t=seq(1995,1998,length=20),t2=seq (1995,1998,length=20)) predict(beer_fit_parabola,newdata=data.frame(t=seq(1995,1998,length=20)) shall be enough Regards Petr #it is listing all 400+ entries ,but not 20 ahead prediction. Thanks In advance for your help -- View this message in context: http://r.789695.n4.nabble.com/Help-me-with- prediction-in-linear-model-tp2297313p2297313.html Sent from the R help mailing list archive at Nabble.com. __ R-help@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. __ R-help@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.