Hi all, I am trying to interparete the result of the following output from lm; fit1 =lm(Feed _Intake ~ weight + season + weight*season) Season has three classes(x,y,z) Reults are Estimate (Intercept) 21.51559 weight 2.13051 factor(season)y 10.59739 factor(season)z 1.30421 weight:factor(season)y 10.1 weight:factor(season)z 21.70288 My question are what is the estimate of season x? Could it be possible to change the output in the following way? factor(season)x factor(season)y weight:factor(season)x weight:factor(season)y Thanks in adavance
Try this: DF$season <- relevel(DF$season, 'y') fit1 <- lm(Feed_Intake ~ weight + season + weight*season, data = DF) On Mon, Jan 18, 2010 at 2:00 PM, Ashta <sewashm at gmail.com> wrote:> Hi all, > I am trying to interparete ?the result of the following output from ?lm; > > > fit1 =lm(Feed _Intake ~ weight + season + weight*season) > Season has three classes(x,y,z) > > Reults are > > Estimate (Intercept) ? ? ? ? ? ? ? 21.51559 > weight ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 2.13051 > factor(season)y ? ? ? ? ? ? ? ? ? ? ?10.59739 > factor(season)z ? ? ? ? ? ? ? ? ? ? ? ?1.30421 > weight:factor(season)y ? ? ? ? ?10.1 > weight:factor(season)z ? ? ? ? ?21.70288 > > My question are ?what is the estimate of season x? > > Could it be possible to change the output in the following way? > > factor(season)x > factor(season)y > weight:factor(season)x > weight:factor(season)y > > Thanks in adavance > > ______________________________________________ > 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. >-- Henrique Dallazuanna Curitiba-Paran?-Brasil 25? 25' 40" S 49? 16' 22" O
Season X is taken as the reference category. So the output
"factor(season)y 10.59739" means the feed_intake is higher by 10.59
units in Season Y _compared to_ Season X.
Change your levels in season. E.g.
season <- factor(season, levels=c("Z", "X",
"Y")
which means that Z will be taken as the reference category. Also read
help(contrasts).
Regards, Adai
Ashta wrote:> Hi all,
> I am trying to interparete the result of the following output from lm;
>
>
> fit1 =lm(Feed _Intake ~ weight + season + weight*season)
> Season has three classes(x,y,z)
>
> Reults are
>
> Estimate (Intercept) 21.51559
> weight 2.13051
> factor(season)y 10.59739
> factor(season)z 1.30421
> weight:factor(season)y 10.1
> weight:factor(season)z 21.70288
>
> My question are what is the estimate of season x?
>
> Could it be possible to change the output in the following way?
>
> factor(season)x
> factor(season)y
> weight:factor(season)x
> weight:factor(season)y
>
> Thanks in adavance
>
> ______________________________________________
> 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.