[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