[R] A simple question about summary.glm

michal33 michal_sa at hotmail.com
Tue Jul 31 23:48:22 CEST 2007


Thanks Uwe for your question,

Yesterday I found out what I wanted to know: 
In the summary table the values refer to the different of each treatment to
the first treatment alphabetical (in my case it was the control group and
therefore served as Tukey test, which is what I wanted in the first place).



Uwe Ligges wrote:
> 
> 
> 
> michal33 wrote:
>> Hello,
>> 
>> I am new to R and have tried to search similar questions but could not
>> find
>> exactly what I am looking for, but I apologize if the question was
>> already
>> asked.
>> 
>> I have 10 different treatments and want to know whether they affect the
>> sex
>> ratios of insect emergence. After running the glms I got this table:
>> 
>>       Df Deviance Resid. Df Resid. Dev      F   Pr(>F)   
>> NULL                    133     9250.3                   
>> sex    1    481.5       132     8768.9 7.7212 0.006314 **
>> trt    9   1099.1       123     7669.7 1.9585 0.049780 * 
>> 
>> But now I would like to know WHICH of the treatments was significant. I
>> tried to use Tukey test but for some reason it does not work. 
>> My question is:
>> I used the following function: 
>>> summary(file.name, corr=F)
>> and got the following table:
>> Deviance Residuals: 
>>     Min       1Q   Median       3Q      Max  
>> -14.118   -4.808   -1.466    2.033   33.882  
>> Coefficients:
>>               Estimate Std. Error   t value Pr(>|t|)    
>> (Intercept)  8.696e+00  1.893e+00     4.594 1.06e-05 ***
>> sexm        -3.791e+00  1.364e+00    -2.779  0.00631 ** 
>> trtccc      -1.050e+00  4.325e+00    -0.243  0.80859    
>> trtcga3      2.450e+00  4.325e+00     0.566  0.57211    
>> trtcga4     -2.300e+00  4.325e+00    -0.532  0.59584    
>> trtg         1.550e+00  2.497e+00     0.621  0.53593    
>> trtga4      -5.550e+00  4.325e+00    -1.283  0.20183    
>> trtp         5.422e+00  2.566e+00     2.113  0.03658 *  
>> trtpg       -1.850e+00  2.497e+00    -0.741  0.46019    
>> trtw        -3.634e-17  2.497e+00 -1.46e-17  1.00000    
>> trtwg       -3.750e+00  2.497e+00    -1.502  0.13573    
>> 
>> What do the stars  mean?
> 
> 
> Well, you omitted the last lines from the output which include:
> 
> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> 
> 
> i.e. *** for values < 0.001; * for values in [0.01, 0.05) etc.
> 
> Uwe Ligges
> 
> 
> 
>   Is it the same as Tukey test that tells me which
>> treatment is different from which? i.e. is trtp (with *) significantly
>> different to the control (which, by the way do not appear in this list
>> and I
>> do not know why)? 
>> 
>> Thanks
>> Michal
>> 
>>
> 
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
> 
> 

-- 
View this message in context: http://www.nabble.com/A-simple-question-about-summary.glm-tf4167757.html#a11935815
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list