[R] understanding the output
Ronaldo Reis Jr.
chrysopa at insecta.ufv.br
Thu Jun 20 18:40:24 CEST 2002
Em Quinta 20 Junho 2002 09:11, Peter Dalgaard BSA escreveu:
> "Ronaldo Reis Jr." <chrysopa at insecta.ufv.br> writes:
> > Estimate Std. Error t value Pr(>|t|)
> > (Intercept) 4.8660 0.9808 4.961 0.000101 ***
> > gen2g 3.8504 1.3870 2.776 0.012459 *
> > ---
>
> ...
>
> > > tapply(peso,gen,mean)
> >
> > 1g 2g
> > 4.8660 8.7164
> >
> > The values is the same in estimate.
> > But the SE, when I try using tapply the SE is different os the SE of
> > estimates. look
> >
> > > tapply(peso,gen,sd)/sqrt(tapply(peso,gen,length))
> >
> > 1g 2g
> > 0.8665567 1.0829789
> >
> > Why this difference? What SE I can use with mean???
>
> (I assume you are only looking at the intercept, i.e. the mean of
> group "1g"). glm() assumes a constant SD and uses the pooled SD.
>
> Try for illustration:
>
> sds <- tapply(peso,gen,sd)
> lens <- tapply(peso,gen,len)
> dfs <- lens - 1
> sd.pool <- sqrt(sum(sds^2*dfs)/sum(dfs))
> sd.pool/sqrt(lens)
>
> It's a matter of taste which SE you'd put on a plot. Some people feel
> that you "shouldn't mix in the other groups" in a mainly descriptive
> plot, but as soon as any modelling is involved (e.g., age-adjusted
> means) you pretty much have to extract SE's from the model.
Its Ok,
I choice to use the model SE. So, my table of mean +- SE maybe:
trat mean SE
1g 4.8660 0.9808
2g 8.7164 1.3870
or
trat mean SE
1g 4.8660 0.9808
2g 8.7164 1.3870+0.9808
or
trat mean SE
1g 4.8660 0.9808
2g 8.7164 0.9808
Or neither??
I'm a bit confused
And for poisson errors? The mean to use is antilog of estimate, and the
associated error, is the glm calculated or need to be transformed?
Sorry by these newbie questions :))
Thanks
Ronaldo
--
[It is] best to confuse only one issue at a time.
-- K&R
--
| //|\\ [*****************************][*******************]
|| ( õ õ ) [Ronaldo Reis Júnior ][PentiumIII-600 ]
| V [ESALQ/USP-Entomologia, CP-09 ][HD: 30 + 10 Gb ]
|| / l \ [13418-900 Piracicaba - SP ][RAM: 128 Mb ]
| /(lin)\ [Fone: 19-429-4199 r.229 ][Video: SiS620-8Mb ]
||/(linux)\ [chrysopa at insecta.ufv.br ][Modem: Pctel-onboar]
|/ (linux) \[ICQ#: 5692561 ][SO: CL 7.0 (2.2.19)]
|| ( x ) [*****************************][*******************]
||| _/ \_Powered by Conectiva Linux 7.0 D+:) | Lxuser#: 205366
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list