# The MIT License
# Copyright (c) 2007 The GGobi Foundation
# http://www.ggobi.org/book/code-license.txt


# Read data
d.olive<-read.csv("olive.csv",row.names=1)
d.olive.sub <- subset(d.olive,region!=1,
   select=c(region,palmitic:eicosenoic))

# Load library
library(classifly)

# Run classifly
classifly(d.olive.sub,region~linoleic+oleic+arachidic, lda)
classifly(d.olive.sub,factor(region)~linoleic+oleic+arachidic, svm, 
  probability=TRUE, kernel="linear")

# More examples
classifly(d.olive, region ~ eicosenoic+linoleic , lda)
classifly(d.olive, region ~ eicosenoic+linoleic+arachidic, lda)

# 
d.olive.sub <- subset(d.olive,region==1,
   select=c(area,palmitic:eicosenoic))
d.olive.sub[d.olive.sub[,1]!=3,1]<-1
d.olive.sub[d.olive.sub[,1]==3,1]<-2
classifly(d.olive.sub,factor(area)~palmitic+palmitoleic+stearic+
  linoleic+oleic+linolenic+arachidic+eicosenoic, svm, 
  probability=TRUE, kernel="linear")
classifly(d.olive.sub,factor(area)~palmitic+palmitoleic+stearic+
  linoleic+oleic+linolenic+arachidic+eicosenoic, svm, 
  probability=TRUE, kernel="radial")