matlab - 获取正态分布变量的 CDF 值的 NaN：Matlab

Question

我似乎无法理解为什么下面给出的简单代码不起作用。我得到NaN了P_j价值观。

% Plotting the CDF for normally distributed variable
mu_j = 0.008;
simga_j = -0.002;
x= 0.008 -0.002*randn(100000,1);
X_j = min(x) : (max(x) - min(x))/100000 : max(x);
P_j = normcdf(X_j,0.008,-0.002); 

figure(1)

plot(X_j,P_j);

谢谢

Answer 1

标准偏差应该是正值。这是因为在统计中，它被计算为sigma = sqrt(E[(X-mu)^2]))。期望中的平方项永远是正的。

如果您调整sigma_j为 be0.002而不是-0.002，并在代码中的其他地方更正该值（例如当您计算xand时P_j），您应该会得到很好的结果。

具体问题是因为normcdf()没有正确处理负sigma，并返回NaN。