To account for likely differences between
families in naturalization rates, we fitted a
generalized linear mixed model, using
PROC NLMIXED in SAS10, with the
naturalization rate per genus (that is, the
number of naturalized species in a genus as
a proportion of the total number of introduced
species in a genus) as the response
variable, a variable coding genera as containing
at least one native species or not as a
fixed-effect predictor, and a variable coding
for family as a random effect.
I would like to emulate or do something better than this PROC NLMIXED method
using R, the methods were used in a 2002 Duncan and Williams nature
publication. Any suggestions about how to proceed?
Chris Buddenhagen
[[alternative HTML version deleted]]
______________________________________________
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.
