estimate the parameters of t-distribution and fitting with MLE Here is Fitting the t-Distribution by Maximum Likelihood t-method in book Statistics and Data Analysis for Financial Engineering with R examples page 113 and 168. 
But I cannot understand
For the univariate case (first picture):


*

*How does the fitting work? Is it the method called expectation–maximization (EM) algorithm? 

*And I understood as first we obtain $\mu,\sigma$ by fitting, then use MLE to obtain the parameter $\nu,$ is it correct? 
For the multivariate case (second picture):


*Does that mean we first compute $\nu,$ then use MLE to compute $\mu$ and $\Lambda?$ 

*If so, then how could we compute $\nu$ with unknown $\mu$ and $\Lambda?$ And the logic totally reverses compared with the univariate case. 

*I think the parameter estimation should be consistent between univariate and multivariate cases and I must have a big misunderstanding.


 A: This is just a partial answer in that it addresses the univariate questions.
The method in that textbook tells you to maximize the log of the likelihood so it is using maximum likelihood (using odd notation in my opinion as it seems to assume unjustifiably that $\mu$ and $\sigma$ are known and that it displays those two parameters that have nothing to do with a "standardized" t-distribution - maybe that's a typo in the book).
One can also use the method of moments where one equates the sample moments to the associated expected moments.  For this distribution one could use $E(X)=\mu$, $Var(X)=\nu \sigma^2/(\nu-2)$ (for $\nu>2$), and $E[(X-\mu)^4]=\frac{3 \nu ^2 \sigma ^4}{(\nu -4) (\nu -2)}$ (for $\nu>2$).  If $\bar{x}$ is the sample mean, $s^2$ is the sample variance, and $m_4$ is the sample fourth central moment, then $\hat{\mu}=\bar{x}$, $\hat{\sigma}^2=\frac{m_4 s^2}{2 m_4-3 s^4}$, and $\hat{\nu}=\frac{2 \left(2 m_4-3 s^4\right)}{m_4-3 s^4}$.  (The distribution is symmetric about $\mu$ so the third central moment doesn't help us estimate the parameters.  In other words, the expected third central moment is zero and therefore not dependent on any of the parameters.)
These values can be used as starting values for an iterative maximum likelihood procedure (although I think that $\bar{x}$ will be the maximum likelihood estimator for $\mu$).
