Laura, puedes acceder directamente al modelo de esta forma <modelo>$finalModel
y ejecutar predict
tal como lo vienes haciendo normalmente. Veamos un ejemplo:
df <- read.table(text = 'Cobranza Facturacion
100 200
200 400
300 600
100 220
100 210', header=TRUE)
cobranza.lm <- lm(Cobranza ~ Facturacion, data=df)
Esto es crear un modelo lineal básico para predecir una cobranza en función a su facturación, ahora si quiero predecir este valor a partir de un dato de facturación de 150$
, podemos hacer esto:
newdata <-data.frame(Facturacion = c(150))
predict(cobranza.lm, newdata, interval="predict")
Esto nos retornará, el valor y el intervalo pedido:
fit lwr upr
1 69.61691 52.76363 86.47019
Si usamos el paquete caret
, se puede hacer lo mismo que pides, pero de forma ligeramente distinta:
library("caret")
cobranza.caret.lm <- train(Cobranza ~ Facturacion, data=df, method='lm')
predict(cobranza.caret.lm$finalModel, newdata, interval = "predict")
Como puedes ver, construimos el mismo modelo anterior, al que podemos acceder desde cobranza.caret.lm$finalModel
, y usarlo directamente en predict
tal como hicimos antes. Deberíamos obtener los mismos intervalos:
fit lwr upr
1 69.61691 52.76363 86.47019
Nota Importante: el parámetro interval
es propio de la Clase del modelo definido, particularmente lo implementan los modelos lm
y nls
, para cualquier otro modelo que no lo implemente, el parámetro será ignorado.