## Jonkers and Kučera (2015) present data on the abundance (shell influx) of the foraminiferan Turborotalita quinqueloba


day_in_year<-c(0.9023162,1.9475052,1.9111508,17.947441,17.61407,19.799334,33.986637, 33.615097,39.078434,49.624847,50.02184,58.024895,65.629875,65.61569,76.569275,81.67198,82.0228,95.510284,97.67119,98.09835,114.81702,113.7471,113.76201,133.79655,129.08539,145.89676,145.54774,152.80917,161.92287,171.76254,177.9617,191.43936,194.37025,209.6631,209.6791,228.60011,225.7012,247.55856,251.5441,254.8538,270.53964,270.5327,286.1364,286.54538,289.83145,302.56824,302.5624,308.39548,318.93423,318.5678,327.66803,334.57318,334.956,346.61084,350.22012,350.6073)
log_number<-c(1.434,0.454,-0.294,-0.123,0.513,0.506,0.109,-0.041,-0.056, -0.415, 0.258,0.034,-0.887,-1.179,-0.655,-0.596,-0.872,-0.678,-1.187, 0.108,-0.671,-0.199,0.107,0.099,0.601,1.73,2.044,1.55,1.69,1.782,1.916,1.909,2.253,2.118,2.447,2.012,2.327,2.35,1.907,2.55,3.01,2.86,1.63,2.55,2.707,2.445,2.325,2.422,1.906,1.86,1.726,1.397,1.779,1.741,1.0525,1.524)
plot(day_in_year,log_number)
lines(lowess(day_in_year, log_number, f = .2, iter=5, delta=1),col="red",lwd=2)
m<-lowess(day_in_year, log_number)
lines(m,col="green",lwd=2)
model<-loess(log_number ~ day_in_year)
xv<-seq(0,366,0.1) 
yv<-predict(model,xv) 
lines(xv,yv,col="blue",lwd=2)