MSPE=anova(fitted.model)["Residuals", "Sum Sq"]/(nrow(DF)-ncol(DF))*
(1+1/nrow(DF)+(newdata-mean(DF[,2]))^2/(var(DF[,2])*(nrow(DF)-1)))
Asumiendo que los datos están en el data.frame "DF" y que el dato usado para la prediccion esta en "newdata" he usado: