[R-br] Teste de diferça entre proporções via loop
Daniel Tiezzi
dtiezzi em usp.br
Quarta Março 23 13:46:17 BRT 2016
Boa tarde Wagner
Supondo que seus dados esteja uma tabela data_test (em anexo)
acho que você pode resolver assim
data<- read.table(“data_test.txt”, header=T, sep= “\t”)
data_df <- as.data.frame(data)
res<- list()
for (i in 1:nrow(data)) {
y<- prop.test(x = c(data$X1[i], data$X2[i]), n = c(data$N[i], data$N[i]), correct = FALSE)
x= y$p.value
res[[paste0(i)]] <- x
}
result <- t(as.data.frame(res))
data_df1 <- cbind(data_df, result)
data_df1
Espero que funcione
daniel
Daniel Tiezzi, MD, PhD
Professor Associado
Departamento de Ginecologia e Obstetrícia
Setor de Mastologia e Oncologia Ginecológica
Faculdade de Medicina de Ribeirão Preto - USP
Tel.: 16 3602-2488
e-mail: dtiezzi em fmrp.usp.br
> On Mar 23, 2016, at 3:50 PM, Wagner Tassinari <wtassinari em gmail.com> wrote:
>
> AREA X1 X2 N p_valor
> 1 1982 9416 11398 ?
> 2 6546 4084 10630 ?
> 3 3418 667 4085 ?
> 4 1503 4408 5911 ?
> 5 966 339 305 ?
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160323/a2adfaea/attachment.html>
-------------- Próxima Parte ----------
Um texto embutido e sem conjunto de caracteres especificado foi limpo...
Nome: data_test.txt
Url: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160323/a2adfaea/attachment.txt>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160323/a2adfaea/attachment-0001.html>
Mais detalhes sobre a lista de discussão R-br