http://www.commres.net/ 2026-04-17T10:46:54+00:00 http://www.commres.net/ http://www.commres.net/_media/wiki/logo.png text/html 2024-12-12T01:27:27+00:00 Anonymous (anonymous@undisclosed.example.com) http://www.commres.net/logistic_regression/examples_r?rev=1733966847&do=diff R code # log cacluation # j = 1 j <- log10(10) j # 10^j = (10) 10^j k <- log2(10) k 2^k set.seed(101) ########## # see youtube # https://youtu.be/8nm0G-1uJzA n.mut <- 23+117 n.norm <- 6+210 p.cancer.mut <- 23/(23+117) p.cancer.norm <- 6/(6+210) set.seed(1011) c <- runif(n.mut, 0, 1) # 0 = not cancer, 1 = cancer among mutant gene mutant <- ifelse(c>=p.cancer.mut, 0, 1) c <- runif(n.norm, 0, 1) # 0 = not cancer, 1 = cancer among normal gene normal <- ifelse(c>=p.cancer.norm, 0, 1) # 0 = m…