我似乎无法理解为什么我下面给出的简单代码不起作用。我获得了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);
由于
答案 0 :(得分:2)
标准偏差应为正值。这是因为在统计中,它计算为sigma = sqrt(E[(X-mu)^2]))
。期望中的平方词只会是积极的。
如果您将sigma_j
调整为0.002
而不是-0.002
,并更正代码中其他位置的值(例如,当您计算x
和{{1}时}),你应该得到很好的结果。
由于P_j
无法正确处理否定sigma,并返回normcdf()
,因此导致了特定问题。