
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@fmrp.usp.br
On Mar 23, 2016, at 3:50 PM, Wagner Tassinari <wtassinari@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 ?