Nothing
#' Remove aliquots
#'
#' This function removes a list of aliquots from a \code{\linkS4class{TLum.Analysis}} object.
#'
#' @param object
#' \code{\linkS4class{TLum.Analysis}} (\bold{required}): object containing the initial TL curves.
#' @param list
#' \link{numeric} (\bold{required}): list containing the position of the aliquots to remove.
#'
#' @return
#' This function provides a \code{\linkS4class{TLum.Analysis}} without the aliquots specified in the list.
#'
#'
#' @author David Strebler, University of Cologne (Germany).
#'
#' @export mod_remove.aliquot
mod_remove.aliquot <- function(
object,
list
){
# ------------------------------------------------------------------------------
# Integrity Check
# ------------------------------------------------------------------------------
if (missing(object)){
stop("[mod_update.dType] Error: Input 'object' is missing.")
}else if (!is(object,"TLum.Analysis")){
stop("[mod_update.dType] Error: Input 'object' is not of type 'TLum.Analysis'.")
}
if(missing(list)){
stop("[mod_align.peaks] Error: Input 'list' is missing.")
}
else if(!is.numeric(list)){
stop("[mod_align.peaks] Error: Input 'list' is not of type 'numeric'.")
}
# ------------------------------------------------------------------------------
nRecord <- length(object@records)
#temp.id <- 0
temp.records <- list()
for(i in 1:nRecord){
temp.curve <- object@records[[i]]
temp.position <- temp.curve@metadata$POSITION
if(!(temp.position %in% list)) {
#temp.id <- temp.id+1
#temp.curve@metadata$ID <- temp.id
temp.records <- c(temp.records, temp.curve)
}
}
#----------------------------------------------------------------------------------------------
# Generate TLum.Analysis
#----------------------------------------------------------------------------------------------
new.records <- temp.records
new.protocol <- object@protocol
new.history <- c(object@history,
as.character(match.call()[[1]]))
new.plotData <- list()
new.plotHistory <- object@plotHistory
new.plotHistory[[length(new.plotHistory)+1]] <- new.plotData
new.analysis <- set_TLum.Analysis(records = new.records,
protocol = new.protocol,
history = new.history,
plotHistory = new.plotHistory)
#--------------------------------------------------------------------------------------------------------
#Plot results
#--------------------------------------------------------------------------------------------------------
#----------------------------------------------------------------------------------------------
#Return results
#----------------------------------------------------------------------------------------------
return(new.analysis)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.