[R] GLM with a multinomial response

Petr Smilauer petrsm at jcu.cz
Sat Nov 2 10:33:12 CET 2002

Dear all,
  does anyone know how to setup a generalized linear model for data looking
the following example. This would be an experiment on food preferences for,
snails, and I would measure for each specimen its choice made on food type
plus one or more quantitative (or factorial) explanatory variables, here
by the snail size.
SpeciesSelected    SnailSize
    Spc1                      2.7
    Spc2                      3.5
    Spc1                      2.9
    Spc3                      4.6


I believe that even in the situation where my explanatory variables were
I shall not use the loglinear models (i.e. modelling counts and moving the
descriptor to the right side, studying its interaction terms with the
predictors), because
the total count has an upper bound determined by the experimental design and I
model the probability of selecting particular plant species given snail size.
Doing separate models with binomial distribution for each food type is not
correct, either,
because the constrain of probabilities summing up to 1.0 is violated.

Does anyone have any thoughts on this issue? This shall not be so rare type
of analysis, given the various customer preference studies. Of course, I can
classification trees (or neural-network models), but I would like to use more
description of dependence.

Thank you for any help you might provide!

Best wishes

Petr Smilauer
Faculty of Biological Sciences
Ceske Budejovice, Czech Republic
Course Multivariate Analysis of Ecological Data <http://regent.bf.jcu.cz>

r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch

More information about the R-help mailing list