R_L11

/*
Задание
По данным о процентном содержании гидробионтов в желудке горбуши проведите
дискриминантный анализ для определения района вылова. Таблица данных INZs.csv. Все
ли рыбы были правильно классифицированы?

link https://yadi.sk/d/0oIEomPxkd1x7w
*/

View(INZc)
str(INZc)
INZc.pca<-princomp(scale(INZc[,2:11]))
plot(INZc.pca)

INZ.p<-predict(INZc.pca)
INZ.p
plot(INZ.p[,2:11],type="n",xlab="PC1",ylab="PC2")
text(INZ.p[,2:11],labels=abbreviate(INZc[,1],1,method="both.sides"))

biplot(INZc.pca)
loadings(INZc.pca)

#Дискриминантный анализ

INZ1=INZc[seq(1,nrow(INZc),11),]
INZ2=INZc[-seq(1,nrow(INZc),11),]
INZ.lda=lda(INZ1[,2:11],INZ1[,1])
lda(INZ1[,2:11],INZ1[,1])
INZ.ldap=predict(INZ.lda,INZ2[,2:11])$class