Codes/Figures                        Données                         Compléments                     Erratum





Comp9. Température mensuelle moyenne à Nottingham



# 9.2.2 Régression sur fonctions trigonométriques

require("zoo")
yt=function(t) 48.9426-9.3016*cos(2*pi*t/12)-6.9901*sin(
2*pi*t/12)+1.5393*sin(2*pi*t/12)+0.495*sin(2*pi*t/12)

date.1 = as.Date(nottem)

yt.zoo=zoo(cbind(nottem,yt(1:length(date.1))),date.1)

plot(yt.zoo)

Fig9-1b