# 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")