parody.R
In parody: Parametric And Resistant Outlier DYtection

## ----setuppp,echo=FALSE-------------------------------------------------------
suppressPackageStartupMessages({
library(parody)
library(BiocStyle)
})

## ----do1----------------------------------------------------------------------
library(parody)

## ----dun21--------------------------------------------------------------------
lead <- c(83, 70, 62, 55, 56, 57, 57, 58, 59, 50, 51, 
        52, 52, 52, 54, 54, 45, 46, 48,
        48, 49, 40, 40, 41, 42, 42, 44, 
        44, 35, 37, 38, 38, 34, 13, 14)

## ----dobox1-------------------------------------------------------------------
boxplot(lead)

## ----dobox2-------------------------------------------------------------------
calout.detect(lead,alpha=.05,method="boxplot", scaling=
function(n,alpha)1.5)

## ----dobox3-------------------------------------------------------------------
calout.detect(lead,alpha=.05,method="boxplot",ftype="ideal")

## ----dor----------------------------------------------------------------------
calout.detect(lead,alpha=.05,method="GESD",k=5)

## ----domemd-------------------------------------------------------------------
calout.detect(lead,alpha=.05,method="medmad",scaling=hamp.scale.3)

## ----dosh---------------------------------------------------------------------
calout.detect(lead,alpha=.05,method="shorth")

## ----dot----------------------------------------------------------------------
data(tcost)
ostr = mv.calout.detect(tcost)
ostr

## ----doparis------------------------------------------------------------------
thecol = rep("black", nrow(tcost))
thecol[ostr$ind] = "red"
pairs(tcost, col=thecol, pch=19)

## ----dobipi-------------------------------------------------------------------
pc = prcomp(tcost)
biplot(pc)

## ----mvmv---------------------------------------------------------------------
ftcost = tcost[-c(9,21),]
fpc = prcomp(ftcost)
biplot(fpc)