Unit testing for postprob() in Mixture scenario

[audreyyeoCH]

the following did not apply and can benefit from troubleshooting and unit testing:

a = postprob(x = 16,
n = 23,
p = 0.60,
parE = c(0.6, 0.4), weights = 1)
b = postprob(x = 16,
n = 23,
p = 0.60,
parE = c(2, 4), weights = 1)

does 0.5*(a + b ) = postprob(x = 16,
n = 23,
p = 0.60,
par = rbind(c(0.6, 0.4),
c(2, 4)),
weights = c(0.5, 0.5))

Currently this is not the case.

• unit test for the above
• trouble shoot line 165 in pbetaMix() to make sure matrix multiplication is correct
• For the above unit testing, can do weights (1,0) = for case a, and weights (0,1) = for case b.
• There should be a check or assert for weights not equalling to 1, and a way to calculate that and print comment to user that the function has calculate the weights as user input was not summing to 1.