cat("Test n")
matrix_with <-
matrix(c(12, (335-12), 4, (341-4)),
     nrow = 2,
     dimnames = list(Treated = c("Positive", "Negative"),
                     Control = c("Positive", "Negative")))
print(matrix_with)
cat("Proportions postive: TX followed by controln")
print(12/335)
print(4/341)
print(fisher.test(matrix_with, alternative = "two.sided"))
#------------
    Fisher's Exact Test for Count Data

data:  matrix_with
p-value = 0.04477
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
  0.9348216 13.4298907
sample estimates:
odds ratio
  3.125192