[R-sig-eco] R-sig-ecology Digest, (80/1): Binomial GLM: Is it the right way?
Abdoul Dia
dia.abdoul at gmail.com
Sat Nov 1 16:15:49 CET 2014
Dear Frate, I don't know if I do understand well your problem but
Binomial suppose that you have a binary response variable. For
proportions, you should use other models otherwise you should specifiy
proportions of success vs failure.
Hope this could help.
Cheers
Abdoul
On 14-11-01 10:29 AM, Ludovico Frate wrote:
> Dear Class, thank you for your reply.
> Maybe I have not been very clear in the exposition. I'll try to reformulate.
>
> "the likelihood of colonizing or becoming extinct" is a proportion so the
> binomial distribution should be still ok!
>
> Considering only colonizing species, or arrivals, for a given plot. This
> means that a species is not new in overall sampled flora but only for that
> plot (i.e. *local* colonizer).
>
> PLOTID;TRAITS;NUMBER OF NEW COLONIZERS;NUMBER OF SUCCESS;PROPORTION
> 1; ChF; 4; 0; 0
> 1; ChC; 4; 0; 0
> 1; ChR; 4; 0; 0
> 1; ChSC;4; 0; 0
> 1; ChS; 4; 1; 0,25
> 1; Ph; 4; 0; 0
> 1; Geo; 4; 0; 0
> 1; HeB; 4; 0; 0
> 1; HeC; 4; 0; 0
> 1; HeR; 4; 0; 0
> 1; HeS; 4; 3; 0,75
> 1; Th; 4; 0; 0
> 2; ChF; 17; 0; 0
> 2; ChC; 17; 0; 0
> 2; ChR; 17; 1; 0,058823529
> 2; ChSC;17; 0; 0
> 2; ChS; 17; 2; 0,117647059
> 2; Ph; 17; 0; 0
> 2; Geo; 17; 1; 0,058823529
> 2; HeB; 17; 0; 0
> 2; HeC; 17; 0; 0
> 2; HeR; 17; 2; 0,117647059
> 2; HeS; 17; 7; 0,411764706
> 2; Th; 17; 4; 0,235294118
> 3; ChF; 16; 0; 0
> 3; ChC; 16; 0; 0
> 3; ChR; 16; 2; 0,125
> 3; ChSC;16; 0; 0
> 3; ChS; 16; 0; 0
> 3; Ph; 16; 0; 0
> 3; Geo; 16; 2; 0,125
> 3; HeB; 16; 0; 0
> 3; HeC; 16; 0; 0
> 3; HeR; 16; 1; 0,0625
> 3; HeS; 16; 9; 0,5625
> 3; Th; 16; 2; 0,125
> . . . .
> . . . .
> . . . .
> 24 . . .
>
>
> for example, in the table, the plot 1 has 4 new species found of which 1 is
> a ChS and 3 are HeS with a proportion of 0.25 and 0.75 respectively.
>
> In a GLM model, the response variable should be PROPORTION (that I have
> improperly called PROBABILITY)
>
> An hypothetical model could be:
> M0<-glm(PROPORTION~TRAITS, family = binomial, weights = NUMBER OF NEW
> COLONIZERS, data = ...)
>
>
> Does it make sense?
>
>
> For the second question, I agree with you that "the probability of a species
> colonizing a plot not only depends on which life form that particular
> species has".
>
> Regards,
> Ludovico
>
>
>
> --
> View this message in context: http://r-sig-ecology.471788.n2.nabble.com/Re-R-sig-ecology-Digest-80-1-Binomial-GLM-Is-it-the-right-way-tp7579191p7579192.html
> Sent from the r-sig-ecology mailing list archive at Nabble.com.
>
> _______________________________________________
> R-sig-ecology mailing list
> R-sig-ecology at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
More information about the R-sig-ecology
mailing list