Differences

This shows you the differences between two versions of the page.

--- wald_test [2023/12/07 13:37] – hkimscil
+++ wald_test [2023/12/07 23:51] (current) – [Wald test in logistic regression] hkimscil
@@ Line 69: / Line 69: @@
 ilogit     <- function(x)      exp(x)/(1+exp(x))
-iter <- 100000
+iter <- 10000
-n <- 100
+n <- 350
 p.cancer <- 0.08
 p.mutant <- 0.39
-orgroup <- rep (NA, iter)
+logor <- rep (NA, iter)
+pp0 <- rep (NA, iter)
+pp1 <- rep (NA, iter)
+op0 <- rep (NA, iter)
+op1 <- rep (NA, iter)
+or <- rep (NA, iter)
 for(i in 1:iter){
   c <- runif(n, 0, 1)
@@ Line 83: / Line 89: @@
   da <- data.frame(gene, canc)
   tab <- table(da)
-  pp0 <- tab[1,1] / (tab[1,1] + tab[1,2])
+  pp0[i] <- tab[1,1] / (tab[1,1] + tab[1,2])
-  pp1 <- tab[2,1] / (tab[2,1] + tab[2,2])
+  pp1[i] <- tab[2,1] / (tab[2,1] + tab[2,2])
-  op0 <- odds(pp0)
+  op0[i] <- odds(pp0[i])
-  op1 <- odds(pp1)
+  op1[i] <- odds(pp1[i])
-  ortemp <- odds.ratio(pp0, pp1)
+  or[i] <- odds.ratio(pp0[i], pp1[i])
-  orgroup[i] <- log(ortemp)
+  # stats <- c(pp0, pp1, op0, op1, ortemp)
+  logor[i] <- log(or[i])
 }
-hist(orgroup)
+hist(logor,breaks = 50)
 </code>