data-olive: Olive oil samples from Italy

In cepp: Context Driven Exploratory Projection Pursuit

Description

This data is from a paper by Forina, Armanino, Lanteri, Tiscornia (1983) Classification of Olive Oils from their Fatty Acid Composition, in Martens and Russwurm (ed) Food Research and Data Anlysis. We thank Prof. Michele Forina, University of Genova, Italy for making this dataset available.

• region Three super-classes of Italy: North, South and the island of Sardinia

• area Nine collection areas: three from North, four from South and 2 from Sardinia

• palmitic, palmitoleic, stearic, oleic, linoleic, linolenic, arachidic, eicosenoic fatty acids percent x 100

Usage

data(olive)

Format

A 572 x 10 numeric array

Examples

data(olive)
head(olive)
##Permutation
OlivesT <- as.matrix(olive[, -c(1:2)])
OlivesF <- OlivesT
#You should set.seed here so as to "fix" the benchmark
OlivesF[, 'palmitic']	<- OlivesF[sample(572,572), 'palmitic']
OlivesF[, 'palmitoleic'] <- OlivesF[sample(572,572), 'palmitoleic']
OlivesF[, 'stearic'] 	<- OlivesF[sample(572,572), 'stearic']
OlivesF[, 'oleic'] 	<- OlivesF[sample(572,572), 'oleic']
OlivesF[, 'linoleic'] 	<- OlivesF[sample(572,572), 'linoleic']
OlivesF[, 'linolenic']	<- OlivesF[sample(572,572), 'linolenic']
OlivesF[, 'arachidic']	<- OlivesF[sample(572,572), 'arachidic']
OlivesF[, 'eicosenoic']	<- OlivesF[sample(572,572), 'eicosenoic']
##
oil1 <- pp(r=2, n=50, oth=OlivesF, data=OlivesT, k=2)
##In practice try at least >10 starting values
F1   <- basis_random(8)
##Increase iterations to >2000 for useful results
o1 <- optim(par=F1, fn=oil1, gr=basis_nearby(), method='SANN',
            control=list(fnscale=-1, maxit=50, trace=6))

Example output

Loading required package: trust
  region          area palmitic palmitoleic stearic oleic linoleic linolenic
1      1 North-Apulia      1075          75     226  7823      672        36
2      1 North-Apulia      1088          73     224  7709      781        31
3      1 North-Apulia       911          54     246  8113      549        31
4      1 North-Apulia       966          57     240  7952      619        50
5      1 North-Apulia      1051          67     259  7771      672        50
6      1 North-Apulia       911          49     268  7924      678        51
  arachidic eicosenoic
1        60         29
2        61         29
3        63         29
4        78         35
5        80         46
6        70         44
Warning message:
In crossprod(x[, j], x[, i]) * x[, i] :
  Recycling array of length 1 in array-vector arithmetic is deprecated.
  Use c() or as.vector() instead.

sann objective function values
initial       value -0.027156
final         value -17.594911
sann stopped after 49 iterations
There were 50 or more warnings (use warnings() to see the first 50)

cepp documentation built on May 2, 2019, 3:44 p.m.