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.