title: Default Linear Discriminant Analysis solver in R MASS
body: <p>What's the default solver of Linear Discriminant Analysis (LDA) in R MASS? For example when running something along the lines</p>
  
  <pre><code>mylda &lt;- lda(dep~ind1+ind2+ind3)
  </code></pre>
  
  <p><a href="https://stat.ethz.ch/R-manual/R-devel/library/MASS/html/lda.html" rel="nofollow">R MASS docs on LDA</a> are not really clear. There is a mention of parameter <code>method</code>, which can be "momentum" or "mle" (maximum likelihood estimation), but can't see how this relates to scikit-learn, where <a href="http://scikit-learn.org/stable/modules/generated/sklearn.discriminant_analysis.LinearDiscriminantAnalysis.html#sklearn.discriminant_analysis.LinearDiscriminantAnalysis" rel="nofollow">scikit-learn API docs on LDA</a> clearly state that the following can be uses for the solver: Singular value decomposition ("svd"), Least squares ("lsqr"), or Eigenvalue decomposition ("eigen").</p>