Nothing
#' @title Consumer's Risk for Multi-state RDT with Multiple Failure Modes (MFM)
#'
#' @description Define the consumer risk function which gets the probability of passing the test when the lower level reliability requirements are not satisfied under different failure modes (for Multi-state RDT, Multiple Failure Modes).
#'
#' @param n RDT sample size
#' @param cvec Maximum allowable failures for each separate period
#' @param pivec Failure probability for each seperate period
#' @param Rvec Lower level reliability requirements for each cumulative period from the begining of the test.
#' @return Probability for consumer's risk
#' @examples
#' pi1 <- pi_MCSim_beta(M = 1000, seed = 10, a = 1, b = 1)
#' pi2 <- pi_MCSim_beta(M = 1000, seed = 10, a = 2, b = 18)
#' MFM_consumerrisk(n = 10, cvec = c(1, 1), pivec = cbind(pi1, pi2), Rvec = c(0.8, 0.7))
#' @export
#' @family MSRDT for MFM functions
#' @seealso \code{\link{MFM_core}} for getting the core probability of passting the test;
#' \code{\link{MFM_Indicator}} for getting the binary indicator;
#' \code{\link{MFM_optimal_n}} for getting the optimal test sample size;
MFM_consumerrisk <- function(n, cvec, pivec, Rvec){
sum1 <- rep(NA, length(1:dim(pivec)[1]))
sum2 <- rep(NA, length(1:dim(pivec)[1]))
for(i in 1:dim(pivec)[1]){
sum2[i] <- MFM_core(n, cvec, pivec[i,])
sum1[i] <- sum2[i] * MFM_Indicator(pivec[i,], Rvec)
}
return(1 - (sum(sum1) / sum(sum2)))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.