我正在尝试执行在网站上找到的一些代码:http://r-statistics.co/Loess-Regression-With-R.html,使用给定的示例函数来找到最佳的平滑范围。
由于某些原因,我总是使用示例中给出的确切代码获得未找到的错误对象sse
。
任何人都可以帮助/解释为什么这吗? 谢谢
# define function that returns the SSE
calcSSE <- function(x){
loessMod <- try(loess(uempmed ~ index, data=economics, span=x), silent=T)
res <- try(loessMod$residuals, silent=T)
if(class(res)!="try-error"){
if((sum(res, na.rm=T) > 0)){
sse <- sum(res^2)
}
}else{
sse <- 99999
}
return(sse)
}
# Run optim to find span that gives min SSE, starting at 0.5
optim(par=c(0.5), calcSSE, method="SANN")