[R] compare linear regressions

Jeff Newmiller jdnewmil at dcn.davis.ca.us
Tue Nov 8 15:58:37 CET 2011


Just because your model specification is valid does not mean your data can be analyzed using that particular model specification.
---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
--------------------------------------------------------------------------- 
Sent from my phone. Please excuse my brevity.

Holger Taschenberger <Holger.Taschenberger at mpi-bpc.mpg.de> wrote:

>Thank you. I was confused because the output of the line
>"summary(lm(formula
>=..." reads:
>
>"Coefficients: (1 not defined because of singularities)",
>
>which did not look like a normal message (which can safely be ignored)
>to me.
>
>--Holger
>
>On Tue, 8 Nov 2011 15:06:28 +0100
>Peter Konings <peter.l.e.konings at gmail.com> wrote:
>
>> On Tue, Nov 8, 2011 at 2:52 PM, Holger Taschenberger <
>> Holger.Taschenberger at mpi-bpc.mpg.de> wrote:
>> <snip>
>> 
>> > summary(lm(formula = y ~ x - 1 + x:g, dta))
>> > ##################
>> > I assume that the last line is the correct way to specify a linear
>model
>> > without intercept. But I'm not certain about that. Can someone
>please
>> > confirm?
>> 
>> 
>> Yes, that's true. See chapter 11 of the "Introduction To R" manual
>that was
>> installed with R for an overview of model specification in R.
>> 
>> HTH
>> Peter.
>
>______________________________________________
>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