R/onoff_singlethres.R

In biosignalEMG: Tools for Electromyogram Signals (EMG) Analysis

onoff_singlethres <- function(data, channel, eemg, t = 0.05, data.name) {
    if (missing(data)) 
        stop("'data' argument is not specified")
    if (!is.emg(data)) 
        stop("an object of class 'emg' is required")
    if (missing(channel)) {
        if (missing(data.name)) 
            data <- extractchannel(data) else data <- extractchannel(data, data.name = data.name)
    } else {
        if (missing(data.name)) 
            data <- extractchannel(data, channel) else data <- extractchannel(data, channel, data.name)
    }
    if (!is.numeric(t)) 
        stop("The threshold 't' must be a numeric value.")
    if (missing(eemg)) {
        eemg <- envelope(data, method = "MA", wsize = 60)
    } else {
        if (!is.emg(eemg)) 
            eemg <- envelope(data, method = "MA", wsize = 60)
    }
    detected <- as.numeric(eemg$values > t)
    detected[is.na(detected)] <- 0
    return(detected)
}