search for: shorlin

...ist.,data=d) kruskal.test(s~d$Dist.,data=d) df<-df.residual(a1) t<-qt(.975,df) k<-data.frame(s=seq(10.5,max(s,na.rm=T),1)) p<-predict(a1,k,se=T) est<-exp(p$fit) low<-exp(p$fit-t*p$se.fit) hight<-exp(p$fit-t*p$se.fit) plot(d$Dist.,s,ylim=c(min(s),max(s)),xlab="Distance from Shorline (m)", ylab="Length of D. aethiops (mm)", pch=7) lines(k$s,est,lty=1,col="red") lines(k$s,low,lty=2,col="blue") lines(k$s,high,lty=2,col="blue") [[alternative HTML version deleted]]