Como obter os intervalos de confiança para o ajuste LOWESS usando R?
Não encontrei nenhuma resposta satisfatória para os intervalos de confiança (ICs) da linha de regressão LOWESS do pacote 'stats' de R:
plot(cars, main = "lowess(cars)")
lines(lowess(cars), col = 2)
Mas não sei como desenhar um IC de 95% em torno dele ?? No entanto, eu sei que poderia obter a variação estimada de
V = s^2*sum(w^2)
onde, s2 = variação estimada do erro ew = pesos aplicados ao X. Portanto, os ICs de 95% devem ser
Y plus/minus 2*sqrt(V(Y))
Sei que há uma maneira de obter os ICs com menos flexibilidade, mas prefiro o LOWESS porque é robusto. Obrigado por suas sugestões.