setwd("") data<-read.table("Developmental_time_ancova.csv", header=T, sep=";", quote="") head(data) plot(data$Mean.developmental.time~data$Total.leg.length.3) anova(lm(data$Mean.developmental.time+data$body.length~data$Total.leg.length.3))