search for: strubb

...example) lme(Bird_abundance = study_area + count +year+ X1 + X2 + X3,random = ~count|study_area). However, I have found no analogue design for a GAM. For now, I have averaged my bird abundances but I wondered whether a more subtle and elegant strategy exists...? Many thanks, Diederik Diederik Strubbe Evolutionary Ecology Group Department of Biology, University of Antwerp Universiteitsplein 1 B-2610 Antwerp, Belgium http://webhost.ua.ac.be/deco tel : 32 3 820 23 85 [[alternative HTML version deleted]]

...s : Warning messages: 1: In Ops.factor(pts$x, xlim[1]) : - not meaningful for factors 2: In Ops.factor(pts$y, ylim[1]) : - not meaningful for factors 3: In Ops.factor(pts$z, zlim[1]) : - not meaningful for factors Any idea about what is wrong and what to do about it? Thanks! Diederi Diederik Strubbe Evolutionary Ecology Group Department of Biology, University of Antwerp Universiteitsplein 1 B-2610 Antwerp, Belgium http://webhost.ua.ac.be/deco tel : 32 3 820 23 85 [[alternative HTML version deleted]]

...; mean <- mean(data$duck) > var <- var(data$duck) > proportion_non_zero [1] 0.1153846 > mean [1] 1.906593 > var [1] 37.35587 (I have no idea how to simulate a zero-inflated overdispersed Poisson variable, but the data used can be found at http://www.ua.ac.be/main.aspx?c=diederik.strubbe&n=23519). First, I create a (strong) pattern in the dataset by: data$LFAP200 <- data$LFAP200 + (data$duck*data$duck) I try to analyze these data by fitting several possible distributions (Poisson PO, zero-inflated Poisson ZIP, negative binomial type I and type II NBI NBII and zero-inflat...

...03 0 0.418671 0.087187 0.504914 2 E 4 2003 6 0.418671 0.087187 0.504914 2 E 5 2003 4 0.418671 0.087187 0.504914 3 E 1 2003 0 0.418671 0.087187 0.504914 3 E 2 2003 0 0.418671 0.087187 0.504914 3 E 3 2003 9 0.418671 0.087187 0.504914 3 E 4 2003 0 0.418671 0.087187 0.504914 3 E 5 2003 0 Diederik Strubbe Evolutionary Ecology Group Department of Biology, University of Antwerp Universiteitsplein 1 B-2610 Antwerp, Belgium http://webhost.ua.ac.be/deco tel : 32 3 820 23 85 [[alternative HTML version deleted]]

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems