Determine a cutoff value between the upper and lower mode of a bimodal dataset

This function finds the cutoff value between the modes of bimodally distributed data at the indicated probability value. As a matter of course, this function runs the cleanData and datamodel functions for you, negating the need to run those separately.

Usage

findcutoff(x = cutoffvalue:::exampledata, proba = 0.5, i = model$indexLower)

Arguments

x

Your dataset specified as "DatasetName$ColumnName" or converted to a numeric list with a name (e.g., "yourrawdata <- as.numeric(yourrawdata$columnname)"). Regardless of how you import or specify it, data should be a single column of log-transformed data.

proba

The probability value you want the cutoff value to represent

i

Determines the location, i.e., index of the (first) minimum of the model$mydata$mu

Value

Determines and reports the cutoff value for the indicated probability value.

Examples

cutoff <- findcutoff(cutoffvalue:::exampledata)
#> number of iterations= 18 

#> Cutoff Value: 0.1137328