[R] logistic regression with repeated measures for species-habitat associations

2010-07-20 Thread Maureen Ryan
Hi, I have three years of species presence/absence data for a set of ~100 ponds, with a list of associated habitat characteristics (park, vegetation, hydroperiod, etc.). The datasets differ slightly by year because some ponds were dry in some years and not in others. My goal is to look at ha

[R] distribution for random effects?

2010-03-02 Thread Maureen Ryan
Hi R users, I am using the following model to analyze data from a factorial experiment (randomized complete block design with no replication within blocks): model <- glm(survival ~ density * vegetation + (1|block), data=sal2005, family="binomial") Does R use a binomial distribution in this for

[R] R, joint scaling test, quantitative genetic analysis & sensitivity to model violations

2009-02-25 Thread Maureen Ryan
Hi all,This is really a stats question as much as an R question. I'm trying to do a joint scaling test (JST - see below) on some very oddly-distributed data and was wondering if anyone can suggest a good way of dealing with model violations and/or using R to evaluate how sensitive the model is

[R] lmer, estimation of p-values and mcmcsamp

2009-02-23 Thread Maureen Ryan
sage: Error in UseMethod("HPDinterval"): no applicable method for "HPDinterval" I imagine at least some of this has to do with the fact that I don't fully understand what exactly the mcmcpvalue function is doing. Any help on either problem would be very much appreciated. Thank