similar to: how to interpret the results of a simint call

Hi all, I wonder if you could help me with (what is surely) a simple R question? I've written a simple R script (attached) to perform multiple ANOVAs on the columns in a table loaded in from a a file (also attached). I loop through the list of columns for which I want to perform the ANOVAs, storing the column name in a variable. When I try to perform the ANOVA using the variable name as a

Hi everybody, I'm having trouble getting an eval-parse combination to work with lme in a for loop. Consider the code below. The call to lme outside the loop works. The call to aov inside the loop works, but the call to to lme inside the loop does not. The latter fails with Error in model.frame.default(formula = ~meanCol + Group + c1 + c3 + Subject, : variable lengths differ (found for

Edward, you also can use the package aod on CRAN, see the help page of the function negbin. Best Matthieu An example: > library(aod) > data(dja) > negbin(y ~ group + offset(log(trisk)), ~group, dja, fixpar = list(4, 0)) Negative-binomial model ----------------------- negbin(formula = y ~ group + offset(log(trisk)), random = ~group, data = dja, fixpar = list(4, 0))

Dear R-listers, I am trying to compute simultaneous confidence intervals with simint from the package multcomp. 230 measures (abundance) have been taken in 23 sites (factor) of a data.frame (donnees: a file can be sent on request, saved with save(donnees,file="donnees")). I would like to get all pairwise comparisons with : mc<- simint(ren~ID,type="Tukey",data=donnees) I

simint and TukeyHSD work for aov objects. Can someone point me to similar functions for lme objects? Dieter Menne

Dear R users; Im working with lme and Id like to have an idea of how can I get CI for the predictions made with the model. Im not a stats guy but, if Im not wrong, the CIs should be different if Im predicting a new data point or a new group. Ive been searching through the web and in help-lists with no luck. I know this topic had been asked before but without replies. Can anyone give an idea of