@@ -125,79 +125,79 @@ test_that("survParamPlot works.", {
125125
126126# logisticCoefs
127127
128- test_that(" logisticCoefs works." , {
129-
130- skip_on_cran()
131-
132- d1 <- defData(varname = " x1" , formula = 0 , variance = 1 )
133- d1 <- defData(d1 , varname = " b1" , formula = 0.5 , dist = " binary" )
134-
135- coefs <- log(runif(2 , min = .8 , max = 1.2 ))
136-
137- # ## Prevalence
138-
139- d1a <- defData(d1 , varname = " y" ,
140- formula = " t(..B) %*% c(1, x1, b1)" ,
141- dist = " binary" , link = " logit"
142- )
143-
144- tPop <- round(runif(1 , .2 , .5 ), 2 )
145- B <- logisticCoefs(defCovar = d1 , coefs = coefs , popPrev = tPop )
146-
147- dd <- genData(100000 , d1a )
148- expect_equal(dd [, mean(y )], tPop , tolerance = .025 )
149-
150- # ### Comparisons
151-
152- d1a <- defData(d1 , varname = " rx" , formula = " 1;1" , dist = " trtAssign" )
153- d1a <- defData(d1a , varname = " y" ,
154- formula = " t(..B) %*% c(1, rx, x1, b1)" ,
155- dist = " binary" , link = " logit"
156- )
157-
158- # ## Risk ratio
159-
160- rr <- runif(1 , .1 , 1 / tPop )
161- B <- logisticCoefs(d1 , coefs , popPrev = tPop , rr = rr , trtName = " rx" )
162-
163- dd <- genData(100000 , d1a )
164- expect_equal(dd [rx == 0 , mean(y )], tPop , tolerance = .025 )
165- expect_equal(dd [rx == 1 , mean(y )]/ dd [rx == 0 , mean(y )], rr , tolerance = 0.025 )
166-
167- # ## risk difference
168-
169- rd <- runif(1 , - tPop , 1 - tPop )
170- B <- logisticCoefs(d1 , coefs , popPrev = tPop , rd = rd , trtName = " rx" )
171-
172- dd <- genData(100000 , d1a )
173- expect_equal(dd [rx == 0 , mean(y )], tPop , tolerance = .025 )
174- expect_equal(dd [rx == 1 , mean(y )] - dd [rx == 0 , mean(y )], rd , tolerance = 0.025 )
175-
176- # ## AUC
177-
178- d1a <- defData(d1 , varname = " y" ,
179- formula = " t(..B) %*% c(1, x1, b1)" ,
180- dist = " binary" , link = " logit"
181- )
182-
183- auc <- runif(1 , 0.6 , 0.95 )
184- B <- logisticCoefs(d1 , coefs , popPrev = tPop , auc = auc )
185-
186- dx <- genData(500000 , d1a )
187- expect_equal(dx [, mean(y )], tPop , tolerance = .025 )
188-
189- form <- paste(" y ~" , paste(d1 [, varname ], collapse = " + " ))
190-
191- fit <- stats :: glm(stats :: as.formula(form ), data = dx )
192- dx [, py : = stats :: predict(fit )]
193-
194- Y1 <- dx [y == 1 , sample(py , 1000000 , replace = TRUE )]
195- Y0 <- dx [y == 0 , sample(py , 1000000 , replace = TRUE )]
196- aStat <- mean(Y1 > Y0 )
197-
198- expect_equal(aStat , auc , tolerance = 0.025 )
199-
200- })
128+ # test_that("logisticCoefs works.", {
129+ #
130+ # skip_on_cran()
131+ #
132+ # d1 <- defData(varname = "x1", formula = 0, variance = 1)
133+ # d1 <- defData(d1, varname = "b1", formula = 0.5, dist = "binary")
134+ #
135+ # coefs <- log(runif(2, min = .8, max = 1.2))
136+ #
137+ # ### Prevalence
138+ #
139+ # d1a <- defData(d1, varname = "y",
140+ # formula = "t(..B) %*% c(1, x1, b1)",
141+ # dist = "binary", link = "logit"
142+ # )
143+ #
144+ # tPop <- round(runif(1, .2, .5), 2)
145+ # B <- logisticCoefs(defCovar = d1, coefs = coefs, popPrev = tPop)
146+ #
147+ # dd <- genData(100000, d1a)
148+ # expect_equal(dd[, mean(y)], tPop, tolerance = .025)
149+ #
150+ # #### Comparisons
151+ #
152+ # d1a <- defData(d1, varname = "rx", formula = "1;1", dist = "trtAssign")
153+ # d1a <- defData(d1a, varname = "y",
154+ # formula = "t(..B) %*% c(1, rx, x1, b1)",
155+ # dist = "binary", link = "logit"
156+ # )
157+ #
158+ # ### Risk ratio
159+ #
160+ # rr <- runif(1, .1, 1/tPop)
161+ # B <- logisticCoefs(d1, coefs, popPrev = tPop, rr = rr, trtName = "rx")
162+ #
163+ # dd <- genData(100000, d1a)
164+ # expect_equal(dd[rx==0, mean(y)], tPop, tolerance = .025)
165+ # expect_equal(dd[rx==1, mean(y)]/dd[rx==0, mean(y)], rr, tolerance = 0.025)
166+ #
167+ # ### risk difference
168+ #
169+ # rd <- runif(1, -tPop, 1 - tPop)
170+ # B <- logisticCoefs(d1, coefs, popPrev = tPop, rd = rd, trtName = "rx")
171+ #
172+ # dd <- genData(100000, d1a)
173+ # expect_equal(dd[rx==0, mean(y)], tPop, tolerance = .025)
174+ # expect_equal(dd[rx==1, mean(y)] - dd[rx==0, mean(y)], rd, tolerance = 0.025)
175+ #
176+ # ### AUC
177+ #
178+ # d1a <- defData(d1, varname = "y",
179+ # formula = "t(..B) %*% c(1, x1, b1)",
180+ # dist = "binary", link = "logit"
181+ # )
182+ #
183+ # auc <- runif(1, 0.6, 0.95)
184+ # B <- logisticCoefs(d1, coefs, popPrev = tPop, auc = auc)
185+ #
186+ # dx <- genData(500000, d1a)
187+ # expect_equal(dx[, mean(y)], tPop, tolerance = .025)
188+ #
189+ # form <- paste("y ~", paste(d1[, varname], collapse = " + "))
190+ #
191+ # fit <- stats::glm(stats::as.formula(form), data = dx)
192+ # dx[, py := stats::predict(fit)]
193+ #
194+ # Y1 <- dx[y == 1, sample(py, 1000000, replace = TRUE)]
195+ # Y0 <- dx[y == 0, sample(py, 1000000, replace = TRUE)]
196+ # aStat <- mean(Y1 > Y0)
197+ #
198+ # expect_equal(aStat, auc, tolerance = 0.025)
199+ #
200+ # })
201201
202202test_that(" logisticCoefs throws errors." , {
203203
0 commit comments