LDA and QDA

简单的R实现：

library(MASS)
Iris <- data.frame(rbind(iris3[,,1], iris3[,,2], iris3[,,3]), Sp = rep(c("s", "c", "v"), rep(50, 3)))
y_hat <- predict(lda(Sp~.,Iris, prior = c(1, 1, 1)/3),Iris)$class
sum(y_hat == Iris$Sp)/150

一般numpy中对于ndarray的处理对python内置list也是可以作用的。（构造函数实现了转型）
下面对sklearn中线性判别及二次判别进行了调用，并与DIY二次判别进行比较，
这里使用了python字符串执行及实例化函数exec及eval提高了DIY函数的可重用性：

from __future__ import division
import numpy as np 
from sklearn.datasets import load_iris 
from sklearn.discriminant_analysis import LinearDiscriminantAnalysis, QuadraticDiscriminantAnalysis
from functools import partial

iris = load_iris()
data, target = iris.data, iris.target

X, y = data, target
clf = LinearDiscriminantAnal