R/Combined_cutoff_percent.R
In manateelab/EMAeval-R-Package: Identification of Careless Responses and Responders in EMA Data
Defines functions
Combined_cutoff_percent
Documented in
Combined_cutoff_percent
#' Identification of Careless Response Percentages Using Standard Deviation and Time per Item
#'
#' \code{Combined_cutoff_percent} returns dataframe of percent of assessment per ID flagged by two cutoff values: Standard Deviation per assessment and Time per Item.
#'
#' This function creates the dataframe that includes the ID and the percentage of assessments that met the cutoff criteria for \emph{Standard Deviation per assessment}, \emph{Time per Item}, and \emph{Percent of Items at Mode }.
#' The percentage is calculated as number of assessments that met both criteria by the total number of assessments for the individual ID.
#' If an assessment has a Standard Deviation \strong{less than or equal to} the cutoff value, a Time per Item \strong{less than or equal to} the cutoff value, or the Percent of Items scores at the Mode is \strong{greater than or equal to} the cutoff value, it will be flagged and placed in the dataframe.
#' This function uses the Boolean logic \emph{"OR"} in order to identify assessments that meet any of the listed criteria.
#' This function allows users to customize cutoff conditions for each separate cutoff to \strong{"<", ">", "<=", ">="}. The default values for the cutoff conditions will be set as described before.
#' This function allows users to customize the Boolean logic to \strong{"AND", "OR"}. The default values for the Boolean logic will be set as \emph{"OR"} between all three criteria.
#'
#' @param data dataframe to be analyzed.
#' @param SD.cutoff numeric cutoff value for Standard Deviation, 'default' value set to 5.
#' @param SD.condition character string of condition that is desired for comparison of data to SD cutoff value, i.e <, >=, etc. 'Default' logic set to "<=".
#' @param TPI.cutoff numeric cutoff value for Time per Item, 'default' value set to 1.
#' @param TPI.condition character string of condition that is desired for comparison of data to TPI cutoff value, i.e <, >=, etc. 'Default' logic set to "<=".
#' @param Perc.Mode.cutoff numeric cutoff value for the Percent of Item scores at Mode, 'default' value set to 0.7.
#' @param Perc.Mode.condition character string of condition that is desired for comparison of data to Percent of Items at Mode cutoff value, i.e <, >=, etc. 'Default' logic set to ">=".
#' @param Combined.logic character string of logic that is desired for classifying the use of the SD, TPI, and Percent of Items at Mode cutoff values, 'default' logic set to 'OR'.
#' @param ttc.colnames vector of column names of "Start Time" and "End Time" to calculate time to complete, also can be "Completion Time" if already calculated.
#' @param number.items integer, number of items per assessment.
#' @param mandatory.response logical value based on whether response to items in each assessment were mandatory to complete.
#' @param item.colnames vector of column names of all items/questions to be used to calculate item score Standard Deviation responses.
#' @param ID.colname character string of column name for ID of assessment.
#' @return The item \code{"ttc.colnames"} must be the names of columns, corresponding to Start Time and End Time, ordered \code{ttc.colnames = c("StartTime", "EndTime")} start time first, followed by end time. If the data includes assessment duration, then list the column name that corresponds with assessment completion: \code{ttc.colnames = "SurveyDuration"}
#' @return The item \code{"mandatory.response"} should be \strong{\code{TRUE}} if participants were required to complete all items per assessment. This item should be \strong{\code{FALSE}} if not all items were required to be completed per assessment.
#' @return If there is variability in the items that are asked, mark the item \code{"mandatory.response"} as \strong{\code{FALSE}} and ensure that \code{"item.colnames"} includes all items.
#' @seealso \code{\link{Combined_cutoff}} for a similar function, presenting a dataframe that includes the ID and index of assessments that met the cutoff criteria specified rather than a dataframe of percentages.
#' @seealso See the following functions for more information on Careless Response Identification in EMA: \code{\link{flagging_df}}, \code{\link{flagging_plots}}, \code{\link{longstringr}}, \code{\link{SD_cutoff}}, \code{\link{TPI_cutoff}}, and \code{\link{Perc_Mode_cutoff}}.
#' @references Jaso, B.A., Kraus, N.I., Heller, A.S. (2020) \emph{Identification of careless responding in ecological momentary assessment: from post-hoc analyses to real-time data monitoring.}
#'
#' @export
#'
Combined_cutoff_percent <- function(data, SD.cutoff, SD.condition, TPI.cutoff, TPI.condition, Perc.Mode.cutoff, Perc.Mode.condition, Combined.logic, ttc.colnames, number.items, mandatory.response, item.colnames, ID.colname){
if (missing(TPI.cutoff)){
warning("TPI.cutoff missing.\nTPI cutoff value automatically set to 1 second. Specify TPI cutoff value to change from 1 second with variable 'TPI.cutoff =...' \n")
TPI.cutoff <- 1
}
if (missing(SD.cutoff)){
warning("SD.cutoff missing.\nSD cutoff value automatically set to 5. Specify SD cutoff value to change from 5 with variable 'SD.cutoff =...' \n")
SD.cutoff <- 5
}
if (missing(Perc.Mode.cutoff)){
warning("Perc.Mode.cutoff missing.\nPercent of Item scores at Mode cutoff value automatically set to 0.7. Specify Percent of Item scores at Mode cutoff value to change from 0.7 with variable 'Perc.Mode.cutoff =...' \n")
Perc.Mode.cutoff <- 0.7
}
if (missing(SD.condition)){
warning("SD.condition missing.\nSD condition (<, >, <=, >=) automatically set to '<='. Specify SD condition to change from '<=' with variable 'SD.condition =...' \n")
SD.condition <- "<="
}
if (SD.condition == "="){
warning("SD.condition cannot be set to '='.\nSD condition (<, >, <=, >=) automatically set to '<='. Specify SD condition to change from '<=' with variable 'SD.condition =...' \n")
SD.condition <- "<="
}
if (missing(TPI.condition)){
warning("TPI.condition missing.\nTPI condition (<, >, <=, >=) automatically set to '<='. Specify TPI condition to change from '<=' with variable 'TPI.condition =...' \n")
TPI.condition <- "<="
}
if (TPI.condition == "="){
warning("TPI.condition cannot be set to '='.\nTPI condition (<, >, <=, >=) automatically set to '<='. Specify TPI condition to change from '<=' with variable 'TPI.condition =...' \n")
TPI.condition <- "<="
}
if (missing(Perc.Mode.condition)){
warning("Perc.Mode.condition missing.\n Percent of Items at Mode condition (<, >, <=, >=) automatically set to '>='. Specify Percent of Items at Mode condition to change from '>=' with variable 'Perc.Mode.condition =...' \n")
Perc.Mode.condition <- ">="
}
if (Perc.Mode.condition == "="){
warning("Perc.Mode.condition cannot be set to '='.\n Percent of Items at Mode condition (<, >, <=, >=) automatically set to '>='. Specify Percent of Items at Mode condition to change from '>=' with variable 'Perc.Mode.condition =...' \n")
Perc.Mode.condition <- ">="
}
if (missing(Combined.logic)){
warning("Combined.logic missing.\nThe logic for combining SD, TPI, and Percent of Items at Mode cutoffs has automatically been set to 'or'. Specify the logic to change from 'or' with variable 'Combined.logic =...' \n")
Combined.logic <- "or"
}
if (missing(ttc.colnames)){
ttc.colnames <- stop("ttc.colnames missing. Specify Start-Date & End-Date / Completion Time with ttc.colnames\n")
}
if (isTRUE(mandatory.response) & missing(number.items)){
number.items <- stop("number.items missing. Specify number of items with the variable 'number.items =...'\n")
}
if (missing(mandatory.response)){
mandatory.response <- stop("mandatory.response missing.\nSpecify if response to all items was mandatory, resulting in with no NAs for item responses.
\nIf responses are mandatory, mandatory.response = TRUE \nIf responses are NOT mandatory, mandatory.response = FALSE")
}
if (isFALSE(mandatory.response) & missing(item.colnames)){
item.colnames <- stop("item.colnames missing.\nSpecify column names of items with the variable 'item.colnames =...'\n")
}
if (missing(item.colnames)){
item.colnames <- stop("item.colnames missing.\nSpecify number of items with the variable 'item.colnames =...'\n")
}
if (missing(ID.colname)){
ID.colname <- stop("ID.colname missing.\nSpecify column name for the identifying information (Participant ID) with the variable 'ID.colname =...'\n")
}
newDF <- c()
for(idx in 1:nrow(data)){
# number of items if response to items not mandatory
if(isFALSE(mandatory.response)){
item_length <- length(which(!is.na(data[idx, item.colnames])))
} else {
item_length <- number.items
}
# calc TTC for the assessment
data_TTC <- c()
if (length(ttc.colnames[!is.na(ttc.colnames)]) > 1){
starttime <- lubridate::ymd_hms(data[idx, which(colnames(data) %in% ttc.colnames[1])])
endtime <- lubridate::ymd_hms(data[idx, which(colnames(data) %in% ttc.colnames[2])])
data_TTC <- as.numeric(endtime-starttime, units = "secs")
} else if (length(ttc.colnames[!is.na(ttc.colnames)]) == 1){
data_TTC <- (data[idx, ttc.colnames[!is.na(ttc.colnames)]])
} else {
data_TTC <- NA
}
# TPI
data_TPI <- (as.numeric(data_TTC) / as.numeric(item_length))
# SD
data_SD <- c()
if (length(item.colnames[!is.na(item.colnames)]) == 0){
data_SD <- NA
} else {
data_SD <- suppressWarnings(sd(as.numeric(data[idx,item.colnames]), na.rm = TRUE))
}
# % items at MODE
data_mode <- c()
if (length(item.colnames[!is.na(item.colnames)]) == 0){
data_mode <- NA
} else {
mode <- suppressWarnings(as.vector(apply(data[idx, which(colnames(data) %in% item.colnames[!is.na(item.colnames)])], 1, DescTools::Mode, na.rm = TRUE)))
Number_items <- length(which(!is.na(data[idx, item.colnames])))
Number_items_Mode <- length(which(data[idx, item.colnames] == mode[1]))
Mode_Frequency_Percent <- Number_items_Mode/Number_items
data_mode <- Mode_Frequency_Percent
}
# check to see if TPI is <= cutoff and paste into dataframe if so
data_point <- c()
flag_point <- c()
#pipe in conditionals specified
Combined.logic <- tolower(Combined.logic)
SD_TPI_Mode_cond <- paste("SD:", SD.condition, " TPI:", TPI.condition, " Perc.Mode:", Perc.Mode.condition, sep = "")
switch(Combined.logic,
"and" = switch(SD_TPI_Mode_cond,
# Different comb with SD TPI when Perc. Mode <=
"SD:<= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# Different comb iwth SD TPI when Perc. Mode >=
"SD:<= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# Different comb with SD TPI when Perc. Mode <
"SD:<= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# Different comb with SD TPI when Perc. Mode >
"SD:<= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next}
),
"or" = switch(SD_TPI_Mode_cond,
# comb SD and TPI with Perc Mode <=
"SD:<= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# comb SD and TPI with Perc Mode >=
"SD:<= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# comb SD and TPI with Perc Mode <
"SD:<= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# comb SD and TPI with Perc Mode >
"SD:<= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next}
)
)
newDF <- rbind(newDF, flag_point)
}
newDF <- as.data.frame(newDF)
colnames(newDF)[1:2] <- c("ID", "Index_of_Flagged_Assessment")
newDF2 <- cbind(ID = rep(NA,length(unique(newDF$ID))), Percent_Flagged = rep(NA,length(unique(newDF$ID))))
uniqueID1 <- unique(newDF$ID)
for(idx in 1:length(uniqueID1)){
newDF2[idx,1] <- uniqueID1[idx]
newDF2[idx,2] <- ((length(which(newDF$ID == uniqueID1[idx])) / length(which(data[,ID.colname] == uniqueID1[idx])) ) * 100)
}
newDF2 <- as.data.frame(newDF2)
newDF2 <- na.omit(newDF2)
newDF2
}
manateelab/EMAeval-R-Package documentation built on Oct. 13, 2021, 6:48 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions Combined_cutoff_percent
Documented in Combined_cutoff_percent
#' Identification of Careless Response Percentages Using Standard Deviation and Time per Item
#'
#' \code{Combined_cutoff_percent} returns dataframe of percent of assessment per ID flagged by two cutoff values: Standard Deviation per assessment and Time per Item.
#'
#' This function creates the dataframe that includes the ID and the percentage of assessments that met the cutoff criteria for \emph{Standard Deviation per assessment}, \emph{Time per Item}, and \emph{Percent of Items at Mode }.
#' The percentage is calculated as number of assessments that met both criteria by the total number of assessments for the individual ID.
#' If an assessment has a Standard Deviation \strong{less than or equal to} the cutoff value, a Time per Item \strong{less than or equal to} the cutoff value, or the Percent of Items scores at the Mode is \strong{greater than or equal to} the cutoff value, it will be flagged and placed in the dataframe.
#' This function uses the Boolean logic \emph{"OR"} in order to identify assessments that meet any of the listed criteria.
#' This function allows users to customize cutoff conditions for each separate cutoff to \strong{"<", ">", "<=", ">="}. The default values for the cutoff conditions will be set as described before.
#' This function allows users to customize the Boolean logic to \strong{"AND", "OR"}. The default values for the Boolean logic will be set as \emph{"OR"} between all three criteria.
#'
#' @param data dataframe to be analyzed.
#' @param SD.cutoff numeric cutoff value for Standard Deviation, 'default' value set to 5.
#' @param SD.condition character string of condition that is desired for comparison of data to SD cutoff value, i.e <, >=, etc. 'Default' logic set to "<=".
#' @param TPI.cutoff numeric cutoff value for Time per Item, 'default' value set to 1.
#' @param TPI.condition character string of condition that is desired for comparison of data to TPI cutoff value, i.e <, >=, etc. 'Default' logic set to "<=".
#' @param Perc.Mode.cutoff numeric cutoff value for the Percent of Item scores at Mode, 'default' value set to 0.7.
#' @param Perc.Mode.condition character string of condition that is desired for comparison of data to Percent of Items at Mode cutoff value, i.e <, >=, etc. 'Default' logic set to ">=".
#' @param Combined.logic character string of logic that is desired for classifying the use of the SD, TPI, and Percent of Items at Mode cutoff values, 'default' logic set to 'OR'.
#' @param ttc.colnames vector of column names of "Start Time" and "End Time" to calculate time to complete, also can be "Completion Time" if already calculated.
#' @param number.items integer, number of items per assessment.
#' @param mandatory.response logical value based on whether response to items in each assessment were mandatory to complete.
#' @param item.colnames vector of column names of all items/questions to be used to calculate item score Standard Deviation responses.
#' @param ID.colname character string of column name for ID of assessment.
#' @return The item \code{"ttc.colnames"} must be the names of columns, corresponding to Start Time and End Time, ordered \code{ttc.colnames = c("StartTime", "EndTime")} start time first, followed by end time. If the data includes assessment duration, then list the column name that corresponds with assessment completion: \code{ttc.colnames = "SurveyDuration"}
#' @return The item \code{"mandatory.response"} should be \strong{\code{TRUE}} if participants were required to complete all items per assessment. This item should be \strong{\code{FALSE}} if not all items were required to be completed per assessment.
#' @return If there is variability in the items that are asked, mark the item \code{"mandatory.response"} as \strong{\code{FALSE}} and ensure that \code{"item.colnames"} includes all items.
#' @seealso \code{\link{Combined_cutoff}} for a similar function, presenting a dataframe that includes the ID and index of assessments that met the cutoff criteria specified rather than a dataframe of percentages.
#' @seealso See the following functions for more information on Careless Response Identification in EMA: \code{\link{flagging_df}}, \code{\link{flagging_plots}}, \code{\link{longstringr}}, \code{\link{SD_cutoff}}, \code{\link{TPI_cutoff}}, and \code{\link{Perc_Mode_cutoff}}.
#' @references Jaso, B.A., Kraus, N.I., Heller, A.S. (2020) \emph{Identification of careless responding in ecological momentary assessment: from post-hoc analyses to real-time data monitoring.}
#'
#' @export
#'
Combined_cutoff_percent <- function(data, SD.cutoff, SD.condition, TPI.cutoff, TPI.condition, Perc.Mode.cutoff, Perc.Mode.condition, Combined.logic, ttc.colnames, number.items, mandatory.response, item.colnames, ID.colname){
if (missing(TPI.cutoff)){
warning("TPI.cutoff missing.\nTPI cutoff value automatically set to 1 second. Specify TPI cutoff value to change from 1 second with variable 'TPI.cutoff =...' \n")
TPI.cutoff <- 1
}
if (missing(SD.cutoff)){
warning("SD.cutoff missing.\nSD cutoff value automatically set to 5. Specify SD cutoff value to change from 5 with variable 'SD.cutoff =...' \n")
SD.cutoff <- 5
}
if (missing(Perc.Mode.cutoff)){
warning("Perc.Mode.cutoff missing.\nPercent of Item scores at Mode cutoff value automatically set to 0.7. Specify Percent of Item scores at Mode cutoff value to change from 0.7 with variable 'Perc.Mode.cutoff =...' \n")
Perc.Mode.cutoff <- 0.7
}
if (missing(SD.condition)){
warning("SD.condition missing.\nSD condition (<, >, <=, >=) automatically set to '<='. Specify SD condition to change from '<=' with variable 'SD.condition =...' \n")
SD.condition <- "<="
}
if (SD.condition == "="){
warning("SD.condition cannot be set to '='.\nSD condition (<, >, <=, >=) automatically set to '<='. Specify SD condition to change from '<=' with variable 'SD.condition =...' \n")
SD.condition <- "<="
}
if (missing(TPI.condition)){
warning("TPI.condition missing.\nTPI condition (<, >, <=, >=) automatically set to '<='. Specify TPI condition to change from '<=' with variable 'TPI.condition =...' \n")
TPI.condition <- "<="
}
if (TPI.condition == "="){
warning("TPI.condition cannot be set to '='.\nTPI condition (<, >, <=, >=) automatically set to '<='. Specify TPI condition to change from '<=' with variable 'TPI.condition =...' \n")
TPI.condition <- "<="
}
if (missing(Perc.Mode.condition)){
warning("Perc.Mode.condition missing.\n Percent of Items at Mode condition (<, >, <=, >=) automatically set to '>='. Specify Percent of Items at Mode condition to change from '>=' with variable 'Perc.Mode.condition =...' \n")
Perc.Mode.condition <- ">="
}
if (Perc.Mode.condition == "="){
warning("Perc.Mode.condition cannot be set to '='.\n Percent of Items at Mode condition (<, >, <=, >=) automatically set to '>='. Specify Percent of Items at Mode condition to change from '>=' with variable 'Perc.Mode.condition =...' \n")
Perc.Mode.condition <- ">="
}
if (missing(Combined.logic)){
warning("Combined.logic missing.\nThe logic for combining SD, TPI, and Percent of Items at Mode cutoffs has automatically been set to 'or'. Specify the logic to change from 'or' with variable 'Combined.logic =...' \n")
Combined.logic <- "or"
}
if (missing(ttc.colnames)){
ttc.colnames <- stop("ttc.colnames missing. Specify Start-Date & End-Date / Completion Time with ttc.colnames\n")
}
if (isTRUE(mandatory.response) & missing(number.items)){
number.items <- stop("number.items missing. Specify number of items with the variable 'number.items =...'\n")
}
if (missing(mandatory.response)){
mandatory.response <- stop("mandatory.response missing.\nSpecify if response to all items was mandatory, resulting in with no NAs for item responses.
\nIf responses are mandatory, mandatory.response = TRUE \nIf responses are NOT mandatory, mandatory.response = FALSE")
}
if (isFALSE(mandatory.response) & missing(item.colnames)){
item.colnames <- stop("item.colnames missing.\nSpecify column names of items with the variable 'item.colnames =...'\n")
}
if (missing(item.colnames)){
item.colnames <- stop("item.colnames missing.\nSpecify number of items with the variable 'item.colnames =...'\n")
}
if (missing(ID.colname)){
ID.colname <- stop("ID.colname missing.\nSpecify column name for the identifying information (Participant ID) with the variable 'ID.colname =...'\n")
}
newDF <- c()
for(idx in 1:nrow(data)){
# number of items if response to items not mandatory
if(isFALSE(mandatory.response)){
item_length <- length(which(!is.na(data[idx, item.colnames])))
} else {
item_length <- number.items
}
# calc TTC for the assessment
data_TTC <- c()
if (length(ttc.colnames[!is.na(ttc.colnames)]) > 1){
starttime <- lubridate::ymd_hms(data[idx, which(colnames(data) %in% ttc.colnames[1])])
endtime <- lubridate::ymd_hms(data[idx, which(colnames(data) %in% ttc.colnames[2])])
data_TTC <- as.numeric(endtime-starttime, units = "secs")
} else if (length(ttc.colnames[!is.na(ttc.colnames)]) == 1){
data_TTC <- (data[idx, ttc.colnames[!is.na(ttc.colnames)]])
} else {
data_TTC <- NA
}
# TPI
data_TPI <- (as.numeric(data_TTC) / as.numeric(item_length))
# SD
data_SD <- c()
if (length(item.colnames[!is.na(item.colnames)]) == 0){
data_SD <- NA
} else {
data_SD <- suppressWarnings(sd(as.numeric(data[idx,item.colnames]), na.rm = TRUE))
}
# % items at MODE
data_mode <- c()
if (length(item.colnames[!is.na(item.colnames)]) == 0){
data_mode <- NA
} else {
mode <- suppressWarnings(as.vector(apply(data[idx, which(colnames(data) %in% item.colnames[!is.na(item.colnames)])], 1, DescTools::Mode, na.rm = TRUE)))
Number_items <- length(which(!is.na(data[idx, item.colnames])))
Number_items_Mode <- length(which(data[idx, item.colnames] == mode[1]))
Mode_Frequency_Percent <- Number_items_Mode/Number_items
data_mode <- Mode_Frequency_Percent
}
# check to see if TPI is <= cutoff and paste into dataframe if so
data_point <- c()
flag_point <- c()
#pipe in conditionals specified
Combined.logic <- tolower(Combined.logic)
SD_TPI_Mode_cond <- paste("SD:", SD.condition, " TPI:", TPI.condition, " Perc.Mode:", Perc.Mode.condition, sep = "")
switch(Combined.logic,
"and" = switch(SD_TPI_Mode_cond,
# Different comb with SD TPI when Perc. Mode <=
"SD:<= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# Different comb iwth SD TPI when Perc. Mode >=
"SD:<= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# Different comb with SD TPI when Perc. Mode <
"SD:<= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# Different comb with SD TPI when Perc. Mode >
"SD:<= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD <= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD >= SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD < SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) & (!is.na(data_SD) & data_SD > SD.cutoff) & (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next}
),
"or" = switch(SD_TPI_Mode_cond,
# comb SD and TPI with Perc Mode <=
"SD:<= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode <= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# comb SD and TPI with Perc Mode >=
"SD:<= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>=" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode >= Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# comb SD and TPI with Perc Mode <
"SD:<= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:<" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode < Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next},
# comb SD and TPI with Perc Mode >
"SD:<= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:<= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI <= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:>= Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI >= TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:< Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI < TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:<= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD <= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:>= TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD >= SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:< TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD < SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
}else {next},
"SD:> TPI:> Perc.Mode:>" = if((!is.na(data_TPI) & data_TPI > TPI.cutoff) | (!is.na(data_SD) & data_SD > SD.cutoff) | (!is.na(data_mode) & data_mode > Perc.Mode.cutoff)){
data_point <- idx
ID <- data[idx,ID.colname]
flag_point <- cbind(ID, data_point)
} else {next}
)
)
newDF <- rbind(newDF, flag_point)
}
newDF <- as.data.frame(newDF)
colnames(newDF)[1:2] <- c("ID", "Index_of_Flagged_Assessment")
newDF2 <- cbind(ID = rep(NA,length(unique(newDF$ID))), Percent_Flagged = rep(NA,length(unique(newDF$ID))))
uniqueID1 <- unique(newDF$ID)
for(idx in 1:length(uniqueID1)){
newDF2[idx,1] <- uniqueID1[idx]
newDF2[idx,2] <- ((length(which(newDF$ID == uniqueID1[idx])) / length(which(data[,ID.colname] == uniqueID1[idx])) ) * 100)
}
newDF2 <- as.data.frame(newDF2)
newDF2 <- na.omit(newDF2)
newDF2
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.