Eu estou enfrentando um problema com a suavização de curvas ajustadas pelas funções lowess
ou loess
no R, no que diz respeito a escolha do parâmetro de suavização f
no lowess
ou span
no loess
Eu procurei algumas dicas na web e encontrei a resolução de um problema similar para a função smooth.spline()
da seguinte maneira
tuneSpline = function(x,y,span.vals=seq(0.1,1,by=0.05),fold=10){
require(bootstrap)
fun.fit <- function(x,y,span) {smooth.spline(x = x,y = y,spar = span)}
fun.predict <- function(fit,x0) {predict(fit,x0)$y}
mae <- sapply(span.vals, function(span){
y.cv <- bootstrap::crossval(
x,y,fun.fit,fun.predict,span=span,ngroup = fold
)$cv.fit
fltr <- which(!is.na(y.cv))
mean(abs(y[fltr]-y.cv[fltr]))
})
span.vals[which.min(mae)]
}
attach(cars)
tuneSpline(speed,dist,fold = length(dist))
# 0.75
tentei modificar esta rotina da seguinte maneira:
tuneSpline = function(x,y,span.vals=seq(0.1,1,by=0.05),fold=10){
require(bootstrap)
fun.fit <- function(x,y,f) {lowess(x = x,y = y, f= f)}
fun.predict <- function(fit,x0) {predict(fit,x0)$y}
mae <- sapply(span.vals, function(f){
y.cv <- bootstrap::crossval(
x, y, fun.fit, fun.predict, f=f, ngroup = fold
)$cv.fit
fltr <- which(!is.na(y.cv))
mean(abs(y[fltr]-y.cv[fltr]))
})
span.vals[which.min(mae)]
}
No entanto quando rodo código no mesmo conjunto de dados, dá o seguinte erro:
attach(cars)
tuneSpline(speed,dist,fold = length(dist))
# Error in UseMethod("predict") :
# método não aplicável para 'predict' aplicado a um objeto de classe "list"
Agradeceria se pudessem ajudar, obrigado
Ps: não postei dados próprios pq não tenho ainda
obrigado desde já
lowess
retorna uma lista dex
ey
, enquantosmooth.spline
retorna um objeto de classesmooth.spline
. A funçãocrossvall
sabe o que fazer com o resultado da segunda, mas não da primeira. Tem certeza que esta solução serve para o seu problema? Talvez o ideal seja postar o seu problema original e não a sua tentativa de adaptar outra solução a ele.loess
dentro de uma função, e então utilizaroptimize
para o parâmetro em questão. Se você tem os dados de entrada, o parâmetro a ser otimizado (dentro de uma faixa) e o critério para considerar ótimo, encontrar o valor ideal não deve ser complicado.