# Load data set
data(abr1)
# Select data set
dat <- abr1$neg
# number of observations and variables
dim(dat)
# Transform data
dat.log <- preproc(dat, method = "log")
dat.sqrt <- preproc(dat, method = "sqrt")
dat.asinh <- preproc(dat, method = "asinh")
op <- par(mfrow=c(2,2), pch=16)
matplot(t(dat),main="Original",type="l",col="blue",
ylab="Intensity")
matplot(t(dat.log),main="Log",type="l",col="green",
ylab="Intensity")
matplot(t(dat.sqrt),main="Sqrt",type="l",col="red",
ylab="Intensity")
matplot(t(dat.asinh),main="ArcSinh)",type="l",col="black",
ylab="Intensity")
par(op)
mtext("Data set", line=2.5, font=3, cex=1.5)
Run the code above in your browser using DataLab