[R] Problems with Cochrane-Orcutt procedures
John Fox
jfox at mcmaster.ca
Tue Jun 17 18:51:44 CEST 2008
Dear Tolga,
I'm afraid that I don't see an error. (I expect in any event that the
Cochrane-Orcott and Prais estimators are now only of historical interest.)
Regards,
John
------------------------------
John Fox, Professor
Department of Sociology
McMaster University
Hamilton, Ontario, Canada
web: socserv.mcmaster.ca/jfox
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
On
> Behalf Of tolga.i.uzuner at jpmorgan.com
> Sent: June-17-08 11:13 AM
> To: r-help at r-project.org
> Subject: [R] Problems with Cochrane-Orcutt procedures
>
> Hi John,
>
> Hi Folks/Prof. Fox,
>
> I found some code John Fox had written sometime back on the
Cochrane-Orcutt
> and Prais procedures here:
> https://stat.ethz.ch/pipermail/r-help/2002-January/017774.html
>
> I thought I would try it out and get the following errors below. Was
> wondering if anyone had any immediate opinions why this might be ?
>
> The linear model is the object regrCMSlm .
>
> Thanks,
> Tolga
>
> > regrCMSlm
>
> Call:
> lm(formula = regrCMS[, 1] ~ regrCMS[, 2])
>
> Coefficients:
> (Intercept) regrCMS[, 2]
> 25.7067 -0.3409
>
> > summary(regrCMSlm)
>
> Call:
> lm(formula = regrCMS[, 1] ~ regrCMS[, 2])
>
> Residuals:
> 09/20/07 11/28/07 02/01/08 04/09/08 06/16/08
> 10.0593 0.3588 -1.1459 0.1340 -9.8520
>
> Coefficients:
> Estimate Std. Error t value Pr(>|t|)
> (Intercept) 25.70673 0.85300 30.14 <2e-16 ***
> regrCMS[, 2] -0.34092 0.02205 -15.46 <2e-16 ***
> ---
> Signif. codes: 0 b
More information about the R-help
mailing list