Package PLS: fonction RMSEP

Bonjour,

Pourquoi dans le package PLS, quand on calcule la RMSE du modèle de prédiction en utilisant la fonction RMSEP(),

la valeur donnée correspond a celle obtenue par l’équation :

Code:

RMSE<-sqrt(sum((mydata$Y-Ypredit)^2)/N)

alors que l'on devrait avoir :

Code:

RMSE<-sqrt(sum((mydata$Y-Ypredit)^2)/N-2)

La différence concerne le nombre de dégrées de liberté, comme l’équation correspond a Y=aX+b, de devrait-on pas
avoir N-2 comme dénominateur a la place de simplement N ?

Merci beaucoup !

Le n-2 s'applique aux données qui ont servis à construire le modèle. Comme le RMSEP est fait pour de "nouvelles" données il me semble relativement cohérent que son calcul se fasse sur N et non N-2.

Merci beaucoup !