Aggregate logit

[yoshida-gisc]

Aggregate logit

[melff]

I think if total is the total number of idividuals from which the proportion is computed then the aggregate model could be fitted as:

mclogit(cbind(total*prop.true, suburb)~distance+cost,data=Transport) 

or

mclogit(cbind(prop.true, suburb)~distance+cost,data=Transport,
        weights=total) 

[melff]

If you have only the proportions, then the results will be the same as in the case where the number of responses is the same in all suburbs. However, the inference will be much more conservative. This is not a problem of this particular implementation, but of any logit model, because the (inverse) number of the responses is used to compute the standard errors.

[yoshida-gisc]

Thank you very much for your prompt reply! I gotcha.

Just for making sure. In the data(Transport) case, total number of resp in suburb 1, 2, and 3 is 212, 212, 219, respectively. But, if we can not know the total number (that is, if we know only ratio), the code mclogit(cbind(prop.true, suburb)~distance+cost,data=Transport) will not be good and not correspond to the aggregate model, right?

So, mclogit(cbind(prop.true, suburb)~distance+cost,data=Transport) will be misleading as aggregate logit model.

[yoshida-gisc]

@melff Yes, I see this is not implementation problem. Thank you so much for your kind help and providing useful package!