Function for replacing outliers for single observations by NA.
removeSingleOut(TP, singleOut, trait = attr(x = singleOut, which ="trait"))
Arguments
TP: An object of class TP.
singleOut: A data.frame with at least the columns plotId and timePoint with values corresponding to those in TP. If a column outlier is present, as in the output of detectSingleOut, only plotId x timePoint combinations for which outlier = 1 will be set to NA. If no column outlier is present, all observations in singleOut will be set to NA.
trait: The trait that should be set to NA. Can be ignored when using the output of detectSingleOut as input.
Returns
An object of class TP, the input with the outlier replaced by NA.
Examples
## Create a TP object containing the data from the Phenovator.PhenovatorDat1 <- PhenovatorDat1[!PhenovatorDat1$pos %in% c("c24r41","c7r18","c7r49"),]phenoTP <- createTimePoints(dat = PhenovatorDat1, experimentName ="Phenovator", genotype ="Genotype", timePoint ="timepoints", repId ="Replicate", plotId ="pos", rowNum ="y", colNum ="x", addCheck =TRUE, checkGenotypes = c("check1","check2","check3","check4"))## First select a subset of plants, for example here 9 plants.plantSel <- phenoTP[[1]]$plotId[1:9]# Then run on the subsetresuVatorHTP <- detectSingleOut(TP = phenoTP, trait ="EffpsII", plotIds = plantSel, confIntSize =3, nnLocfit =0.1)## Replace the studied trait by NA for the plants marked as outliers.phenoTPOut <- removeSingleOut(phenoTP, resuVatorHTP)
See Also
Other functions for detecting outliers for single observations: detectSingleOut(), detectSingleOutMaize(), plot.singleOut()