Come ottenere un R-quadrato per una vestibilità loess?


15

Come calcolare la statistica R-quadrato ( ) in R per e / o output della funzione? Ad esempio per questi dati:r2loesspredict

cars.lo <- loess(dist ~ speed, cars)
cars.lp <- predict(cars.lo, data.frame(speed = seq(5, 30, 1)), se = TRUE)

cars.lpha due array fitper il modello e se.fitper l'errore standard.


Nella regressione lineare è uguale alla correlazione al quadrato tra i valori osservati e i valori adattati - che ne dici? R2
Macro

Risposte:


10

Il mio primo pensiero è stato di calcolare una misura pseudo R2 come segue:

ss.dist <- sum(scale(cars$dist, scale=FALSE)^2)
ss.resid <- sum(resid(cars.lo)^2)
1-ss.resid/ss.dist

Qui, otteniamo un valore di 0.6814984 ( ), vicino a quello che si otterrebbe da un GAM : cor(cars$dist, predict(cars.lo))^2

library(mgcv)
summary(gam(dist ~ speed, data=cars))

loessMultiple R-squaredairqualityairloessR2

inserisci qui la descrizione dell'immagine

R2

spanR2


8
R2R2

r2

3
Non necessariamente "più preciso". In effetti, usare Loess per ottenere la precisione in un modello predittivo sarebbe insensato. Penso che riferirsi a Loess come un "modello" trasmetta un possibile fraintendimento su come funziona e su come deve essere utilizzato: è davvero uno strumento grafico ed esplorativo per aiutare a vedere modelli e tendenze. Poiché è davvero solo una finestra mobile, agisce come un modello di quartiere spaziale abbastanza complicato in cui il valore adattato in un punto dipende da quali punti vicini esistono nel set di dati e dai valori lì.
whuber

2
r2
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.