# R program to display the observed cane infections.

thedata = read.table("scylvFC.txt", header=T)

K = 30

thedata[1:5,]

fact = 10 / K  # to make total display be 10 seconds

plot(thedata$x,thedata$y, xlim=c(0,50), ylim=c(0,50), pch="-")
Sys.sleep(6*fact)
points(thedata$x[thedata$S6==1],thedata$y[thedata$S6==1], col="red", pch="*")
Sys.sleep(4*fact)
points(thedata$x[thedata$S10==1],thedata$y[thedata$S10==1], col="red", pch="*")
Sys.sleep(4*fact)
points(thedata$x[thedata$S14==1],thedata$y[thedata$S14==1], col="red", pch="*")
Sys.sleep(5*fact)
points(thedata$x[thedata$S19==1],thedata$y[thedata$S19==1], col="red", pch="*")
Sys.sleep(4*fact)
points(thedata$x[thedata$S23==1],thedata$y[thedata$S23==1], col="red", pch="*")
Sys.sleep(7*fact)
points(thedata$x[thedata$S30==1],thedata$y[thedata$S30==1], col="red", pch="*")