On Jan 18, 2011, at 8:45 PM, Lao Meng wrote:

Hi all:
Here's a question about result of loglinear analysis.
There're 2 factors:area and nation.The raw data is in the attachment.

I fit the saturated model of loglinear with the command:
glm_sat<-glm(fre~area*nation, family=poisson, data=data_Analysis)

After that,I extract the coefficients:
result_sat<-summary(glm_sat)
result_coe<-result_sat$coefficients

I find that all the coeffients are 1 or very near to 1.

I didn't get that result with that code. Did you perhaps attach some other data structure that you have failed to inform us about?

> result_coe
               Estimate Std. Error   z value      Pr(>|z|)
(Intercept)  5.40707901 0.06797962  79.53971  0.000000e+00
area        -0.12665730 0.01498742  -8.45091  2.890419e-17
nation      -0.42998381 0.01615198 -26.62112 3.867012e-156
area:nation  0.04823113 0.00317437  15.19392  3.879947e-52




How does this happen?Why all the coeffients are 1 or very near to 1?

Thanks!

My best
<area_nation.txt>______________________________________________
R-help@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.

David Winsemius, MD
West Hartford, CT

______________________________________________
R-help@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.

Reply via email to