discriminant analysis

• 判別分析

関

discriminant analyses

WordNet

1. an investigation of the component parts of a whole and their relations in making up the whole
2. the abstract separation of a whole into its constituent parts in order to study the parts and their relations (同)analytic thinking
3. a branch of mathematics involving calculus and the theory of limits; sequences and series and integration and differentiation
4. a form of literary criticism in which the structure of a piece of writing is analyzed
5. the use of closed-class words instead of inflections: e.g., `the father of the bride instead of `the brides father

PrepTutorEJDIC

1. (内容・状況などの)『分析』,分解;(詳細な)検討 / (化学・物理で)分析;《米》(心理学で)[精神]分析;(数学で)解析

Wikipedia preview

出典(authority):フリー百科事典『ウィキペディア（Wikipedia）』「2016/08/28 02:25:17」(JST)

wiki en

[Wiki en表示]

Not to be confused with latent Dirichlet allocation.

Ronald Fisher

Linear discriminant analysis (LDA) is a generalization of Fisher's linear discriminant, a method used in statistics, pattern recognition and machine learning to find a linear combination of features that characterizes or separates two or more classes of objects or events. The resulting combination may be used as a linear classifier, or, more commonly, for dimensionality reduction before later classification.

LDA is closely related to analysis of variance (ANOVA) and regression analysis, which also attempt to express one dependent variable as a linear combination of other features or measurements. However, ANOVA uses categorical independent variables and a continuous dependent variable, whereas discriminant analysis has continuous independent variables and a categorical dependent variable (i.e. the class label). Logistic regression and probit regression are more similar to LDA than ANOVA is, as they also explain a categorical variable by the values of continuous independent variables. These other methods are preferable in applications where it is not reasonable to assume that the independent variables are normally distributed, which is a fundamental assumption of the LDA method.

LDA is also closely related to principal component analysis (PCA) and factor analysis in that they both look for linear combinations of variables which best explain the data. LDA explicitly attempts to model the difference between the classes of data. PCA on the other hand does not take into account any difference in class, and factor analysis builds the feature combinations based on differences rather than similarities. Discriminant analysis is also different from factor analysis in that it is not an interdependence technique: a distinction between independent variables and dependent variables (also called criterion variables) must be made.

LDA works when the measurements made on independent variables for each observation are continuous quantities. When dealing with categorical independent variables, the equivalent technique is discriminant correspondence analysis.

LDA for two classes

Consider a set of observations ${\vec {x}}$ (also called features, attributes, variables or measurements) for each sample of an object or event with known class y. This set of samples is called the training set. The classification problem is then to find a good predictor for the class y of any sample of the same distribution (not necessarily from the training set) given only an observation ${\vec {x}}$.

LDA approaches the problem by assuming that the conditional probability density functions $p({\vec {x}}|y=0)$ and $p({\vec {x}}|y=1)$ are both normally distributed with mean and covariance parameters $\left({\vec {\mu }}_{0},\Sigma _{0}\right)$ and $\left({\vec {\mu }}_{1},\Sigma _{1}\right)$, respectively. Under this assumption, the Bayes optimal solution is to predict points as being from the second class if the log of the likelihood ratios is below some threshold T, so that:

$({\vec {x}}-{\vec {\mu }}_{0})^{T}\Sigma _{0}^{-1}({\vec {x}}-{\vec {\mu }}_{0})+\ln |\Sigma _{0}|-({\vec {x}}-{\vec {\mu }}_{1})^{T}\Sigma _{1}^{-1}({\vec {x}}-{\vec {\mu }}_{1})-\ln |\Sigma _{1}|\ >\ T$

Without any further assumptions, the resulting classifier is referred to as QDA (quadratic discriminant analysis).

LDA instead makes the additional simplifying homoscedasticity assumption (i.e. that the class covariances are identical, so $\Sigma _{0}=\Sigma _{1}=\Sigma$) and that the covariances have full rank. In this case, several terms cancel:

${\vec {x}}^{T}\Sigma _{0}^{-1}{\vec {x}}={\vec {x}}^{T}\Sigma _{1}^{-1}{\vec {x}}$

${\vec {x}}^{T}{\Sigma _{i}}^{-1}{\vec {\mu _{i}}}={\vec {\mu _{i}}}^{T}{\Sigma _{i}}^{-1}{\vec {x}}$ because $\Sigma _{i}$ is Hermitian

and the above decision criterion becomes a threshold on the dot product

${\vec {w}}\cdot {\vec {x}}>c$

for some threshold constant c, where

${\vec {w}}=\Sigma ^{-1}({\vec {\mu }}_{1}-{\vec {\mu }}_{0})$

$c={\frac {1}{2}}(T-{\vec {\mu _{0}}}^{T}\Sigma _{0}^{-1}{\vec {\mu _{0}}}+{\vec {\mu _{1}}}^{T}\Sigma _{1}^{-1}{\vec {\mu _{1}}})$

This means that the criterion of an input ${\vec {x}}$ being in a class y is purely a function of this linear combination of the known observations.

It is often useful to see this conclusion in geometrical terms: the criterion of an input ${\vec {x}}$ being in a class y is purely a function of projection of multidimensional-space point ${\vec {x}}$ onto vector ${\vec {w}}$ (thus, we only consider its direction). In other words, the observation belongs to y if corresponding ${\vec {x}}$ is located on a certain side of a hyperplane perpendicular to ${\vec {w}}$. The location of the plane is defined by the threshold c.

Canonical discriminant analysis for k classes

Canonical discriminant analysis (CDA) finds axes (k − 1 canonical coordinates, k being the number of classes) that best separate the categories. These linear functions are uncorrelated and define, in effect, an optimal k − 1 space through the n-dimensional cloud of data that best separates (the projections in that space of) the k groups. See “Multiclass LDA” for details below.

Fisher's linear discriminant

The terms Fisher's linear discriminant and LDA are often used interchangeably, although Fisher's original article actually describes a slightly different discriminant, which does not make some of the assumptions of LDA such as normally distributed classes or equal class covariances.

Suppose two classes of observations have means ${\vec {\mu }}_{0},{\vec {\mu }}_{1}$ and covariances $\Sigma _{0},\Sigma _{1}$. Then the linear combination of features ${\vec {w}}\cdot {\vec {x}}$ will have means ${\vec {w}}\cdot {\vec {\mu }}_{i}$ and variances ${\vec {w}}^{T}\Sigma _{i}{\vec {w}}$ for $i=0,1$. Fisher defined the separation between these two distributions to be the ratio of the variance between the classes to the variance within the classes:

$S={\frac {\sigma _{\text{between}}^{2}}{\sigma _{\text{within}}^{2}}}={\frac {({\vec {w}}\cdot {\vec {\mu }}_{1}-{\vec {w}}\cdot {\vec {\mu }}_{0})^{2}}{{\vec {w}}^{T}\Sigma _{1}{\vec {w}}+{\vec {w}}^{T}\Sigma _{0}{\vec {w}}}}={\frac {({\vec {w}}\cdot ({\vec {\mu }}_{1}-{\vec {\mu }}_{0}))^{2}}{{\vec {w}}^{T}(\Sigma _{0}+\Sigma _{1}){\vec {w}}}}$

This measure is, in some sense, a measure of the signal-to-noise ratio for the class labelling. It can be shown that the maximum separation occurs when

${\vec {w}}\propto (\Sigma _{0}+\Sigma _{1})^{-1}({\vec {\mu }}_{1}-{\vec {\mu }}_{0})$

When the assumptions of LDA are satisfied, the above equation is equivalent to LDA.

Be sure to note that the vector ${\vec {w}}$ is the normal to the discriminant hyperplane. As an example, in a two dimensional problem, the line that best divides the two groups is perpendicular to ${\vec {w}}$.

Generally, the data points to be discriminated are projected onto ${\vec {w}}$; then the threshold that best separates the data is chosen from analysis of the one-dimensional distribution. There is no general rule for the threshold. However, if projections of points from both classes exhibit approximately the same distributions, a good choice would be the hyperplane between projections of the two means, ${\vec {w}}\cdot {\vec {\mu }}_{0}$ and ${\vec {w}}\cdot {\vec {\mu }}_{1}$. In this case the parameter c in threshold condition ${\vec {w}}\cdot {\vec {x}}>c$ can be found explicitly:

$c={\vec {w}}\cdot {\frac {1}{2}}({\vec {\mu }}_{0}+{\vec {\mu }}_{1})={\frac {1}{2}}{\vec {\mu }}_{1}^{T}\Sigma ^{-1}{\vec {\mu }}_{1}-{\frac {1}{2}}{\vec {\mu }}_{0}^{T}\Sigma ^{-1}{\vec {\mu }}_{0}$.

Otsu's method is related to Fisher's linear discriminant, and was created to binarize the histogram of pixels in a grayscale image by optimally picking the black/white threshold that minimizes intra-class variance and maximizes inter-class variance within/between grayscales assigned to black and white pixel classes.

Multiclass LDA

In the case where there are more than two classes, the analysis used in the derivation of the Fisher discriminant can be extended to find a subspace which appears to contain all of the class variability. This generalization is due to C. R. Rao. Suppose that each of C classes has a mean $\mu _{i}$ and the same covariance $\Sigma$. Then the scatter between class variability may be defined by the sample covariance of the class means

$\Sigma _{b}={\frac {1}{C}}\sum _{i=1}^{C}(\mu _{i}-\mu )(\mu _{i}-\mu )^{T}$

where $\mu$ is the mean of the class means. The class separation in a direction ${\vec {w}}$ in this case will be given by

$S={\frac {{\vec {w}}^{T}\Sigma _{b}{\vec {w}}}{{\vec {w}}^{T}\Sigma {\vec {w}}}}$

This means that when ${\vec {w}}$ is an eigenvector of $\Sigma ^{-1}\Sigma _{b}$ the separation will be equal to the corresponding eigenvalue.

If $\Sigma ^{-1}\Sigma _{b}$ is diagonalizable, the variability between features will be contained in the subspace spanned by the eigenvectors corresponding to the C − 1 largest eigenvalues (since $\Sigma _{b}$ is of rank C − 1 at most). These eigenvectors are primarily used in feature reduction, as in PCA. The eigenvectors corresponding to the smaller eigenvalues will tend to be very sensitive to the exact choice of training data, and it is often necessary to use regularisation as described in the next section.

If classification is required, instead of dimension reduction, there are a number of alternative techniques available. For instance, the classes may be partitioned, and a standard Fisher discriminant or LDA used to classify each partition. A common example of this is "one against the rest" where the points from one class are put in one group, and everything else in the other, and then LDA applied. This will result in C classifiers, whose results are combined. Another common method is pairwise classification, where a new classifier is created for each pair of classes (giving C(C − 1)/2 classifiers in total), with the individual classifiers combined to produce a final classification.

Incremental LDA

The typical implementation of the LDA technique requires that all the samples are available in advance. However, there are situations where the entire data set is not available and the input data are observed as a stream. In this case, it is desirable for the LDA feature extraction to have the ability to update the computed LDA features by observing the new samples without running the algorithm on the whole data set. For example, in many real-time applications such as mobile robotics or on-line face recognition, it is important to update the extracted LDA features as soon as new observations are available. An LDA feature extraction technique that can update the LDA features by simply observing new samples is an incremental LDA algorithm, and this idea has been extensively studied over the last two decades. Catterjee and Roychowdhury proposed an incremental self-organized LDA algorithm for updating the LDA features. In other work, Demir and Ozmehmet proposed online local learning algorithms for updating LDA features incrementally using error-correcting and the Hebbian learning rules. Later, Aliyari et al. derived fast incremental algorithms to update the LDA features by observing the new samples.

Practical use

In practice, the class means and covariances are not known. They can, however, be estimated from the training set. Either the maximum likelihood estimate or the maximum a posteriori estimate may be used in place of the exact value in the above equations. Although the estimates of the covariance may be considered optimal in some sense, this does not mean that the resulting discriminant obtained by substituting these values is optimal in any sense, even if the assumption of normally distributed classes is correct.

Another complication in applying LDA and Fisher's discriminant to real data occurs when the number of measurements of each sample exceeds the number of samples in each class. In this case, the covariance estimates do not have full rank, and so cannot be inverted. There are a number of ways to deal with this. One is to use a pseudo inverse instead of the usual matrix inverse in the above formulae. However, better numeric stability may be achieved by first projecting the problem onto the subspace spanned by $\Sigma _{b}$. Another strategy to deal with small sample size is to use a shrinkage estimator of the covariance matrix, which can be expressed mathematically as

$\Sigma =(1-\lambda )\Sigma +\lambda I\,$

where $I$ is the identity matrix, and $\lambda$ is the shrinkage intensity or regularisation parameter. This leads to the framework of regularized discriminant analysis or shrinkage discriminant analysis.

Also, in many practical cases linear discriminants are not suitable. LDA and Fisher's discriminant can be extended for use in non-linear classification via the kernel trick. Here, the original observations are effectively mapped into a higher dimensional non-linear space. Linear classification in this non-linear space is then equivalent to non-linear classification in the original space. The most commonly used example of this is the kernel Fisher discriminant.

LDA can be generalized to multiple discriminant analysis, where c becomes a categorical variable with N possible states, instead of only two. Analogously, if the class-conditional densities $p({\vec {x}}\mid c=i)$ are normal with shared covariances, the sufficient statistic for $P(c\mid {\vec {x}})$ are the values of N projections, which are the subspace spanned by the N means, affine projected by the inverse covariance matrix. These projections can be found by solving a generalized eigenvalue problem, where the numerator is the covariance matrix formed by treating the means as the samples, and the denominator is the shared covariance matrix.

Applications

In addition to the examples given below, LDA is applied in positioning and product management.

Bankruptcy prediction

In bankruptcy prediction based on accounting ratios and other financial variables, linear discriminant analysis was the first statistical method applied to systematically explain which firms entered bankruptcy vs. survived. Despite limitations including known nonconformance of accounting ratios to the normal distribution assumptions of LDA, Edward Altman's 1968 model is still a leading model in practical applications.

Face recognition

In computerised face recognition, each face is represented by a large number of pixel values. Linear discriminant analysis is primarily used here to reduce the number of features to a more manageable number before classification. Each of the new dimensions is a linear combination of pixel values, which form a template. The linear combinations obtained using Fisher's linear discriminant are called Fisher faces, while those obtained using the related principal component analysis are called eigenfaces.

Marketing

In marketing, discriminant analysis was once often used to determine the factors which distinguish different types of customers and/or products on the basis of surveys or other forms of collected data. Logistic regression or other methods are now more commonly used. The use of discriminant analysis in marketing can be described by the following steps:

1. Formulate the problem and gather data — Identify the salient attributes consumers use to evaluate products in this category — Use quantitative marketing research techniques (such as surveys) to collect data from a sample of potential customers concerning their ratings of all the product attributes. The data collection stage is usually done by marketing research professionals. Survey questions ask the respondent to rate a product from one to five (or 1 to 7, or 1 to 10) on a range of attributes chosen by the researcher. Anywhere from five to twenty attributes are chosen. They could include things like: ease of use, weight, accuracy, durability, colourfulness, price, or size. The attributes chosen will vary depending on the product being studied. The same question is asked about all the products in the study. The data for multiple products is codified and input into a statistical program such as R, SPSS or SAS. (This step is the same as in Factor analysis).
2. Estimate the Discriminant Function Coefficients and determine the statistical significance and validity — Choose the appropriate discriminant analysis method. The direct method involves estimating the discriminant function so that all the predictors are assessed simultaneously. The stepwise method enters the predictors sequentially. The two-group method should be used when the dependent variable has two categories or states. The multiple discriminant method is used when the dependent variable has three or more categorical states. Use Wilks’s Lambda to test for significance in SPSS or F stat in SAS. The most common method used to test validity is to split the sample into an estimation or analysis sample, and a validation or holdout sample. The estimation sample is used in constructing the discriminant function. The validation sample is used to construct a classification matrix which contains the number of correctly classified and incorrectly classified cases. The percentage of correctly classified cases is called the hit ratio.
3. Plot the results on a two dimensional map, define the dimensions, and interpret the results. The statistical program (or a related module) will map the results. The map will plot each product (usually in two-dimensional space). The distance of products to each other indicate either how different they are. The dimensions must be labelled by the researcher. This requires subjective judgement and is often very challenging. See perceptual mapping.

Biomedical studies

The main application of discriminant analysis in medicine is the assessment of severity state of a patient and prognosis of disease outcome. For example, during retrospective analysis, patients are divided into groups according to severity of disease – mild, moderate and severe form. Then results of clinical and laboratory analyses are studied in order to reveal variables which are statistically different in studied groups. Using these variables, discriminant functions are built which help to objectively classify disease in a future patient into mild, moderate or severe form.

In biology, similar principles are used in order to classify and define groups of different biological objects, for example, to define phage types of Salmonella enteritidis based on Fourier transform infrared spectra, to detect animal source of Escherichia coli studying its virulence factors etc.

Earth science

This method can be used to separate the alteration zones. For example, when different data from various zones are available, discriminate analysis can find the pattern within the data and classify it effectively.

See also

• Data mining
• Decision tree learning
• Factor analysis
• Kernel Fisher discriminant analysis
• Logit (for logistic regression)
• Multidimensional scaling
• Pattern recognition
• Perceptron
• Preference regression
• Quadratic classifier

References

Further reading

• Duda, R. O.; Hart, P. E.; Stork, D. H. (2000). Pattern Classification (2nd ed.). Wiley Interscience. ISBN 0-471-05669-3. MR 1802993. 
• Hilbe, J. M. (2009). Logistic Regression Models. Chapman & Hall/CRC Press. ISBN 978-1-4200-7575-5. 
• Mika, S.; et al. (1999). "Fisher Discriminant Analysis with Kernels". IEEE Conference on Neural Networks for Signal Processing IX: 41–48. doi:10.1109/NNSP.1999.788121. 
• Mark Burdon and Paul Harpur, ‘Re-Conceptualising Privacy and Discrimination in an Age of Talent Analytics’ (2014) 37 University of New South Wales Law Journal, 2, 679–712.1
• Miranda Terry and Paul Harpur, ‘The New Era of Segmenting Society on Ability Lines: Workplace Analytics and Disability Discrimination’ (Society for Disability Studies, Atlanta USA, 10–13 June 2015).

External links

• ALGLIB contains open-source LDA implementation in C# / C++ / Pascal / VBA.
• Psychometrica.de open-source LDA implementation in Java
• LDA tutorial using MS Excel
• Biomedical statistics. Discriminant analysis

UpToDate Contents

全文を閲覧するには購読必要です。 To read the full text you will need to subscribe.

• 1. アルコール性肝炎：自然経過およびマネージメント alcoholic hepatitis natural history and management
• 2. Model for End-stage Liver Disease (MELD)
• 3. 計算ツール：アルコール性肝炎におけるコルチコステロイド治療の肝炎判別関数 calculator hepatitis discriminant function for corticosteroid therapy in alcoholic hepatitis
• 4. アルコール性脂肪性肝疾患およびアルコール性肝硬変の予後および治療 prognosis and management of alcoholic fatty liver disease and alcoholic cirrhosis
• 5. 消化器学および肝臓学の最新情報 whats new in gastroenterology and hepatology

English Journal

• Quality control of fruit juices by using organic acids determined by capillary zone electrophoresis with poly(vinyl alcohol)-coated bubble cell capillaries.

• Navarro-Pascual-Ahuir M1, Lerma-García MJ1, Simó-Alfonso EF1, Herrero-Martínez JM2.
• Food chemistry.Food Chem.2015 Dec 1;188:596-603. doi: 10.1016/j.foodchem.2015.05.057. Epub 2015 May 14.
• An enhanced method for the determination of organic acids in several fruit juices by capillary zone electrophoresis (CZE) with direct UV-Vis detection has been developed in this work. First, a study with simulated real juice samples was done to find the best separation conditions. Next, several comm
• PMID 26041236

• Multi-element, multi-compound isotope profiling as a means to distinguish the geographical and varietal origin of fermented cocoa (Theobroma cacao L.) beans.

• Diomande D1, Antheaume I1, Leroux M1, Lalande J1, Balayssac S2, Remaud GS3, Tea I4.
• Food chemistry.Food Chem.2015 Dec 1;188:576-82. doi: 10.1016/j.foodchem.2015.05.040. Epub 2015 May 12.
• Multi-element stable isotope ratios have been assessed as a means to distinguish between fermented cocoa beans from different geographical and varietal origins. Isotope ratios and percentage composition for C and N were measured in different tissues (cotyledons, shells) and extracts (pure theobromin
• PMID 26041233

• Physicochemical characterization of cactus pads from Opuntia dillenii and Opuntia ficus indica.

• Méndez LP1, Flores FT1, Martín JD1, Rodríguez Rodríguez EM1, Díaz Romero C2.
• Food chemistry.Food Chem.2015 Dec 1;188:393-8. doi: 10.1016/j.foodchem.2015.05.011. Epub 2015 May 5.
• Physicochemical characteristics (weight, length, width, thickness, moisture, Brix degree, total fiber, protein, ash, pH, acidity, ascorbic acid, total phenolic compounds, P, Na, K, Ca, Mg, Fe, Cu, Zn, Mn and Cr) were determined in cactus pads from Opuntia dillenii and Opuntia ficus indica. The physi
• PMID 26041209

Japanese Journal

• 高次元データの統計的方法論(日本統計学会研究業績賞受賞者特別寄稿論文)

• 青嶋 誠,矢田 和善
• 日本統計学会誌. シリーズJ 43(1), 123-150, 2013-09
• NAID 110009649320

• Application of FTIR-ATR Spectroscopy Coupled with Multivariate Analysis for Rapid Estimation of Butter Adulteration

• Fadzlillah Nurrulhidayah Ahmad,Rohman Abdul,Ismail Amin [他]
• Journal of Oleo Science 62(8), 555-562, 2013-08
• NAID 40019711244

• 各種ハチミツの微量元素分析による採蜜年度の判別

• 吉垣 茂,中村 純,内田 治,永島 俊夫,安藤 達彦
• 日本健康医学会雑誌 22(2), 107-113, 2013-07-31
• 近年,ハチミツはミネラルやアミノ酸など豊富な栄養素が注目され,健康食品として需要が高まっている。一方,ハチミツは医療への利用という点でも注目されており,医療機関で創傷治療に利用されている。ハチミツでは品質を維持するための純度評価方法が幾つか確立されており,現在では地域判別や種類判別が可能になってきている。本研究では,ハチミツでは今まで報告されてない,同じ蜜源植物で,同じ地域で採蜜されたハチミツの採 …
• NAID 110009634820

Related Links

• discriminant analysis

判別分析: Discriminant Analysis 兵庫教育大学 成田 滋 Updated January 11, 1998 [心理・教育のための統計の初歩][論文の書き方] [偏相関分析について][多変量解析は強力で難しくない] [ノンパラメトリック法による分散分析 ...

• 判別分析 Discriminant Analysis - Murata Lab.

判別分析Discriminant Analysis 担当塚本 平成17 年6 月15 日 1 判別問題とは 複数のデータが、グループごとに分けられていることを考える。新たなデー タを観測したときに、そのデータがどのグループに属するかを決定するよう

Related Pictures

---

★リンクテーブル★

リンク元	「判別分析」「discriminant analyses」
関連記事	「analysis」「discriminant」

「判別分析」

　　[★]

英

discriminant analysis、discriminant analyses

「discriminant analyses」

　　[★]

判別分析

関

discriminant analysis

「analysis」

　　[★]

• n.

• 解析、分析、解析法、分析法

関

anal、analyse、analyses、analytical、analyze、assay、dissect、-metry、solve

　 　 　 　 　 　

「discriminant」

　　[★]

• adj.

• 弁別の

関

discriminate、discrimination、discriminative