Mode and Credible Interval Calculation for The Difference between Two Beta Distributions — sumBetadiff • phase1b

A function to summerize the characters of a betadiff distribution dbetadiff

Usage

sumBetadiff(parY, parX, level = 0.9, cutB, cutW)

Arguments

parY: two parameters of Y's beta distribution (Treatment)
parX: two parameters of X's beta distribution (Control)
level: credible interval
cutB: a cutoff value for the upper bound
cutW: a cutoff value for the lower bound

Value

sumbetadiff gives the mode, credible interval of the distribution function, along with the probabilities of above cutB and below cutW.

Examples

parX <- c(1, 52)
parY <- c(5.5, 20.5)
sumBetadiff(parX, parY, 0.9, cutB = 0.8, cutW = 0.2)
#> $mode
#> [1] -0.1704321
#> 
#> $ci
#> [1] -0.33643341 -0.07270616
#> 
#> $go
#> [1] 3.379515e-43
#> 
#> $nogo
#> [1] 0.9999999
#>