search for: yearbind

...). I've pasted the rest of the code below for reference. I'd appreciate any suggestions as to what the issue here is! Thanks, Varina transect=read.csv("Transect data.csv") attach(transect) replicates=c(1,2) yearfactor=as.factor(replicates) yearframe=data.frame(yearfactor) yearbind=cbind(transect$TP10, transect$TP11) yearmodel=lm(yearbind~1) rmanova=anova(yearmodel, idata=yearframe, idesign=~yearfactor) summary(rmanova) [[alternative HTML version deleted]]