Bonjour,
Je réponds ci-dessous à certaines questions/problèmes laissés en suspend à la fin du TP.
- Pour les données tecator, il faut prendre 20 noeuds pour les splines:
res=funopare.kernel.cv(Y,CURVES,CURVES,semimetric='deriv',q=2,nknot=20,range=c(0,1) )
plot(Y,res$Estimated.values)
abline(0,1)
Cela permet de voir que l'on obtient de meilleurs résultats avec la dérivée d'ordre 2
- D'autre part, pour appliquer d'autres pseudo métriques (pca, ...), il peut être nécessaire de taper CURVES=as.matrix(CURVES) afin de s'assurer que CURVES soit bien au format matriciel.
Cordialement,
Laurent DELSOL