[R] understanding the output

Ronaldo Reis Jr. chrysopa at insecta.ufv.br
Thu Jun 20 13:55:00 CEST 2002


Hi all,
I try to understand the output of the summary(model).
I make a simple model.

> summary(m1anova)

Call:
glm(formula = peso ~ gen)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-4.3114  -2.3788  -0.9167   2.1581   5.7856  

Coefficients:
            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 *  
---
Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 

(Dispersion parameter for gaussian family taken to be 9.618818)

    Null deviance: 247.27  on 19  degrees of freedom
Residual deviance: 173.14  on 18  degrees of freedom
AIC: 105.92

Number of Fisher Scoring iterations: 2

I need the mean and se to make an plotbar (I know to make a plotbar)
The mean is easy to get.

> 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???

When I use others errors that dont gaussian (poisson, binomial, etc), the mean 
is easy transformed (exp(x) for poisson) but how to use the SE? I can to use 
always the SE on the summary(model)? Or need transformation?

Thanks
Ronaldo
 
-- 
Hello, GORRY-O!!  I'm a GENIUS from HARVARD!!
--
|   //|\\   [*****************************][*******************]
|| ( õ õ )  [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