if (FALSE) {
data(tecator)
x=tecator$absorp.fdata[1:129]
y=tecator$y$Fat[1:129]
res1=fregre.pc(x,y,1:5)
# time consuming
res.infl1=influence(res1)
res2=fregre.basis(x,y)
res.infl2=influence(res2)
res<-res1
res.infl<-res.infl1
mat=cbind(y,res$fitted.values,res.infl$DCP,res.infl$DCE,res.infl$DP)
colnames(mat)=c("Resp.","Pred.","DCP","DCE","DP")
pairs(mat)
}
Run the code above in your browser using DataLab