Message-ID: <546404572.1122955.1586898034318@mail.yahoo.com>
Date: 2020-04-14T21:00:34Z
From: varin sacha
Subject: Span for loess regression
Dear R-experts,
I am trying to find the best span for my loess regression. Here below a reproducible example. I don't get the result. Am I missing something ?
Many thanks for your help.
############################################################
a<-c(2,3,4,3,2,6,5,7,4,5,12,13,21,6,4,5,6,7)
b<-c(12,13,32,14,23,21,76,34,12,32,43,23,12,32,43,12,32,11)
model <- loess(b ~ a)
bestLoess <- function(model, spans = c(.05, .95)) {
?f <- function(span) {
???????? mod <- update(model, span = span)
???????? loessGCV(mod)[["gcv"]]
???? }
???? result <- optimize(f, spans)
???? result
?}
#######################################################################