没有LLS.fun函数

1. library(locfit);
   
2. data(ethanol);
   
3. attach(ethanol);
   
4. summary(ethanol)
   
5. Y=NOx;
   
6. X=E;
   
7. Y=Y[order(X)]
   
8. X=sort(X)
   
9. plot(X,Y);
   
10. plot(X,Y); #画原始数据的散点图
    
11. X2=X^2; X3=X^3; X4=X^4;
    
12. fit1<-lm(Y~X)
    
13. lines(X,fit1$fitted.values,lty=1,col=1,lwd=3)
    
14. fit2<-lm(Y~X+X2)
    
15. lines(X,fit2$fitted.values,lty=2,col=2,lwd=3)
    
16. fit3<-lm(Y~X+X2+X3)
    
17. lines(X,fit3$fitted.values,lty=3,col=3,lwd=3)
    
18. fit4<-lm(Y~X+X2+X3+X4)
    
19. lines(X,fit4$fitted.values,lty=4,col=4,lwd=3)
    
20. legend(0.53,4,c("线性","二次","三次","四次"),lty=c(1,2,3,4),col=c(1,2,3,4))
    
21. 
    
22. 
    
23. 
    
24. # Rule-of-Thumb Global Bandwidth Selector;
    
25. X2=X^2; X3=X^3; X4=X^4;
    
26. fit1 <- lm(Y~X+X2+X3+X4);
    
27. coefE=c(fit1$coeff);
    
28. resids=fit1$residuals;
    
29. sigmaE=sqrt(var(resids));
    
30. CK=1.719 # Epanechnikov kernel used.
    
31. temp=cbind(2,3*2*X,4*3*X^2)%*%as.vector(coefE[-(1:2)]);
    
32. den=sum(temp^2);
    
33. h.ROT=CK*(sigmaE^2/den)^(1/(2*1+3));
    
34. print(h.ROT)  #选择窗宽
    
35. 
    
36. #final fitting;
    
37. xfine=seq(0,4,length=50);
    
38. ypred1 <- rep(0,length(xfine));
    
39. for(i in 1:length(xfine)){
    
40.     ypred1[i] <- LLS.fun(xfine[i],X,Y,h=h.ROT);   
    
41. }
    
42. plot(X,Y);
    
43. lines(xfine,ypred1,lty=1,col=3);#拟合数据

复制代码

关键是，提示找不到LLS.fun函数，怎么办？？