Re: [R] Fixed zeros in tables

2006-11-27 Thread David Duffy
A.R. Criswell [EMAIL PROTECTED] wrote: Hello Andrew Robinson and R-List Thanks, Andrew, but this does not work. Puttting zero weights on structural zeros, one's elsewhere in glm() does not deliver the appropriate expected cell counts loglm() provides as one would expect. If you run the

Re: [R] Fixed zeros in tables

2006-11-26 Thread A.R. Criswell
Hello Andrew Robinson and R-List Thanks, Andrew, but this does not work. Puttting zero weights on structural zeros, one's elsewhere in glm() does not deliver the appropriate expected cell counts loglm() provides as one would expect. If you run the code provided below, you'll see loglm() delivers

Re: [R] Fixed zeros in tables

2006-11-26 Thread Andrew Robinson
Hello Andrew, I'm not sure how this is a problem. You can multiply the fitted values by the weights again, if you wish, or ignore the structural zeros altogether. The other predicted values all seem to me to be the same. I think that is a feature. Sometimes the goal for glm with structural

[R] Fixed zeros in tables

2006-11-25 Thread A.R. Criswell
Hello All R Users, Function loglm() in library MASS can be cajoled to accomodate structural zeros in a cross-classification table. An example from Fienberg demonstrates how this can be done. My question is: Can the function glm() perform the same task? Can glm() estimate a log-linear model with

Re: [R] Fixed zeros in tables

2006-11-25 Thread Andrew Robinson
Hi Andrew, try the weights argument - apply zero weight to the structural zeros, and 1 elsewhere. Cheers Andrew On Sun, Nov 26, 2006 at 11:39:39AM +0700, A.R. Criswell wrote: Hello All R Users, Function loglm() in library MASS can be cajoled to accomodate structural zeros in a