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