Suppose $x$ is a random vector in $\mathbb{R}^n$ which is distributed according to $D$.
What is the unbiased estimator of covariance matrix of an N-dimensional random variable?
When $y$ is a i.i.d. random variable and we have access to $(y_1,y_2,\cdots,y_n)$, the sample mean is an unbiased estimator of $\hat{\mu}=\frac{\sum_{i=1}^N}{N}$ and $\hat{\sigma}^2=\frac{1}{N-1}\sum_{i=1}^N(y_i-\hat{\mu})^2$ is an unbiased estimator of variance.
By going to higher dimension in addition to variance we have covariance between each element of the random vector. My question is
$$ \hat{C}=? $$ where $\hat{C}$ is an unbiased estimator of $C = \mathbb{E}[(x-\mu)(x-\mu)^T]$.