Caro Wagner, boa tarde!

Não tenho experiência com esse pacote, mas tentei aplicar o procedimento que utilizo pra avaliar o preenchimento em outras situações. O dataset é restrito, mas deve servir pra testar...

### <code r>

require(mtsdi)

data(miss)

notNA   <- which(!is.na(miss), arr.ind=T)  ### valores presentes

set.seed(333); sel.pos <- notNA[sample(nrow(notNA), 20),] ### reservando 20 valores

sel.obs <- miss[sel.pos]                   ### valores reservados

miss2 <- miss

miss2[sel.pos] <- NA                       ### "exclui" valores reservados

f <- ~c31+c32+c33+c34+c35

i <- mnimput(f,miss2,eps=1e-3,ts=TRUE, method="spline",sp.control=list(df=c(7,7,7,7,7)))

summary(i)

imput <- data.frame(predict(i))

sel.pre <- imput[sel.pos]                  ### predição referente aos 20 valores reservados

cbind(sel.pre, sel.obs)                    ### comparação

library(hydroGOF)

ggof(sel.pre,sel.obs)

### </code>

Éder Comunello <comunello.eder@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
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.