clean

Author: audreyyeoCH

R/dbetabinom.R

@@ -96,7 +96,7 @@ h_getBetamixPost <- function(x, n, par, weights) {
       lbeta(a = par[, 1], b = par[, 2])
   )
   # We compute the updated weights of the posterior
-  postWeights <- exp(log(weights) + log(tmp)) - sum(weights * tmp)
+  postWeights <- weights * tmp / sum(weights * tmp)
   assert_numeric(postWeights)
   list(
     par = postPar,

examples/ocPredprobDist.R

@@ -1,4 +1,3 @@
-\donttest{
 # Here we illustrate an example for Decision 1 with the following assumptions :
 # Efficacy Looks and Futility looks are identical at sample size of 10, 20 and 30.
 # True response rate or truep of the treatment group = 40%