[R] question about result of loglinear analysis

[David Winsemius]

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

David Winsemius, MD
West Hartford, CT