Eu posso estar enganado então me corrijam em caso afirmativo. O
leave-ONE-out é mais simples do que imagina. Grosseiramente falando, crie
um laço for para deixar um caso de fora e ajuste o modelo. Depois estude
uma medida de ajuste. No caso de lm, não requer de fato fazer a exaustiva
tarefa de ajustar n modelos (deixando um caso de fora), pois se chega as
medidas leave-one-out por projeções matriciais, etc. Mas num caso mais
geral é algo como:
da ## seu data.frame
L <- vector(mode="list", length=nrow(da))
for(i in 1:nrow(da)){
L[[i]] <- sua_funcao_R(..., data=da[-i,])
}
## Medidas de deviance leave-one-out.
sapply(L, deviance)
## log-veross leave-one-out.
sapply(L, logLik)
Código não testado.
À disposição.
Walmes.
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
