Wikipedia says "the relationship between the standard error and the standard deviation is such that, for a given sample size, the standard error equals the standard deviation divided by the square root of the sample size." At Standard errors for maximum likelihood estimates in FindDistributionParameters, there is this code:
standarderrors[data_, dist_, paramlist_, mleRule_] :=
Block[{len, infmat, cov},
len = Length[data];
(* compute negative of expected Fisher information *)
infmat = -D[LogLikelihood[dist, data], {paramlist, 2}]/len /. mleRule;
(* invert to get asymptotic covariance for Sqrt[n](theta-theta0) *)
cov = Inverse[infmat];
(* standard errors are the Sqrt of diagonal elements divided by sample size *)
Sqrt[Diagonal[cov]/len]
]
which suggests it gives the Standard Error. But in the same link, the above code is later referred to as giving the Standard Deviation. Which is it?
len
in that code can be completely removed and you'll get the same results.) $\endgroup$ – JimB Jan 29 '18 at 18:20