data(bottles.cont)
plot(bottles.cont)
dev.qual(bottles.cont)
dev.quant(bottles.cont)
harm.pow(bottles.cont)
nef <- get.Nef(bottles.cont)
fac <- factor(rep(c("beer","whisky"), each=20))
pca(nef, fac)
pca3(nef, fac)
pca(tps(nef, fac)
from <- c(-0.5, 0.25)
to <- -from
tps.iso(nef, fr=from, to=to)
tps.vf(nef, fr=from, to=to)
tps.grid(nef, fr=from, to=to)
Run the code above in your browser using DataLab