[R-br] Extrair P-value do resultado de uma imputação com MICE
Paulo Nogueira Starzynski
paulons em gmail.com
Sexta Abril 21 09:50:29 BRT 2017
Bom dia Elias.
Consegui extrair os valores pra depois consultar a distribuição F.
Veja abaixo.
imp <- mice(nhanes2)
fit <- with(imp, lm(chl~age+bmi))
summary(fit)
tam1 <- length(fit)
tam2 <- length(fit[[4]])
s <- summary(fit[[tam1]][[tam2]])
p.value <- 1 - pf(q = s$fstatistic["value"], df1 = s$fstatistic["numdf"],
df2 = s$fstatistic["dendf"])
p.value
Abraços,
Paulo
2017-04-20 20:30 GMT-03:00 Elias Carvalho via R-br <r-br em listas.c3sl.ufpr.br
>:
> Boa noite pessoal
>
> Estou testando o pacote mice para fazer imputação múltipla de dados
> missing, conforme exemplo:
>
> library(mice)
>
> imp <- mice(nhanes2)
> fit <- with(imp, lm(chl~age+bmi))
>
> summary(fit)
>
> O resultado se resume em:
>
> ## summary of imputation 1 :
>
> Call:
> lm(formula = chl ~ age + bmi)
>
> Residuals:
> Min 1Q Median 3Q Max
> -46.37 -19.23 -11.24 12.66 83.69
>
> Coefficients:
> Estimate Std. Error t value Pr(>|t|)
> (Intercept) 32.134 53.202 0.604 0.55232
> age2 49.058 18.153 2.702 0.01334 *
> age3 62.837 19.886 3.160 0.00472 **
> bmi 4.907 1.795 2.733 0.01246 *
> ---
> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
> Residual standard error: 34.89 on 21 degrees of freedom
> Multiple R-squared: 0.3835, Adjusted R-squared: 0.2954
> F-statistic: 4.354 on 3 and 21 DF, p-value: 0.01557
>
>
> Eu gostaria de extrair o p-value: 0.01557 dessa estrutura, mas não
> consegui com os métodos convencionais, pois a estrutura parece diferente da
> estrutura gerada por um lm ou glm.
>
> Alguém poderia me ajudar
>
> --
>
>
> *In Jesu et Maria*
> *Obrigado*
> *Prof. Elias Carvalho*
>
> *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he
> who has been able to understand the cause of things"*
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170421/e6eed917/attachment.html>
Mais detalhes sobre a lista de discussão R-br