比较最大似然估计与克拉美-冯·米塞斯估计
采用最大似然估计法和优化克拉美-冯·米塞斯检验统计量法对正态分布进行参数估计,绘图比较这两种方法.
data = BlockRandom[SeedRandom[2];
RandomVariate[NormalDistribution[], 100]];cvmFit[μ_ ? NumericQ, σ_ ? NumericQ, dist_] := CramerVonMisesTest[data, dist[μ, σ], "TestStatistic"];
cvmPar = Quiet[FindMinimum[cvmFit[μ, σ, NormalDistribution], {{μ}, {σ}}]];mlePar = FindDistributionParameters[data, NormalDistribution[μ, σ]];Show[Plot3D[cvmFit[μ, σ, NormalDistribution], {μ, -3, 3}, {σ, 0, 6}, ColorFunction -> ColorData["SouthwestColors"], PlotRange -> {{-2, 2}, {0, 3}, {-10, 35}}, MeshFunctions -> {#3&}, MeshStyle -> Gray, Mesh -> 30, BoxRatios -> 1], Graphics3D[{Darker@Blue, Thickness[.01], Line[{{μ, σ, 40}, {μ, σ, -20}} /. cvmPar[[2]]]}], Graphics3D[{Darker@Red, Thickness[.01], Line[{{μ, σ, 40}, {μ, σ, -20}} /. mlePar]}]]