小编Oro*_*ngo的帖子

我在一张图中绘制了0-100岁的死亡率.它的颜色非常丰富,但对于读者而言,没有解释谁是0岁,40岁或100岁.我见过彩虹传说是在线图片,但没有得到如何自己写出来的答案.我用来绘制的线是

plot(m, series = "total", datatype="rate", plot.type="time", xlab="Years", main="Total death rate")

我可以改变什么来获得彩虹传奇？是否要求在图形区域中写出网格线？

数据的片段是

structure(c(0.107446, 0.02901, 0.015211, 0.010825, 0.008799, 
0.006939, 0.110958, 0.028299, 0.014543, 0.010939, 0.00893, 0.007335, 
0.091162, 0.026333, 0.013556, 0.009593, 0.007193, 0.005984, 0.097073, 
0.026088, 0.012508, 0.008976, 0.007057, 0.005734, 0.089595, 0.022413, 
0.011188, 0.008234, 0.006767, 0.005522, 0.093491, 0.02269, 0.011209, 
0.009104, 0.00724, 0.006003, 0.08592, 0.019631, 0.009459, 0.006521, 
0.005931, 0.004634, 0.080883, 0.019495, 0.009131, 0.006662, 0.004753, 
0.004184, 0.09061, 0.021713, 0.008999, 0.006332, 0.004717, 0.004221, 
0.076158, 0.016583, 0.008261, 0.005645, 0.004305, 0.003866, 0.078229, 
0.019572, 0.008195, 0.005844, 0.004747, 0.003715, 0.07538, …

我模拟了Irwin Hall分布的1,000,000个重复,每个重复是100个均匀随机变量的总和.我创建了这个程序,我认为它是正确的,但是MATLAB已经整夜运行而没有完成,所以有些(显然)错了,但我看不到哪里.这是我的计划,我希望有人可能会看到这个问题.

n=100; % Nb of samples
N=1000000; %Nb of replicates

for jj=1:N
    for ii=1:n
        x(ii)=rand();
    end
    s(jj)=sum(x);
    mu=mean(s);
    sigma2=var(s);
    sigma=std(s);
    S(jj)=1-normcdf((70-mu)/sigma);
    if mod(jj, 100000)==0
        jj
    end
end

histogram(S)