实现经验 CDF C++

implement empirical cdf c++

本文关键字:C++ CDF 经验 实现      更新时间:2023-10-16

我想根据收到的 N 个样本固定的间隔在C++中定义经验 cdf,我应该保存 2 个向量(T 和 Y)!我编写了这段代码,但它不起作用:

double cum_dist_funct::real_cdf_function(vec X,double x,int N,int Ndft,vec &T, vec &Z)
{
    Z.set_size(X.length());
    Z.zeros();
    vec V;
    V.set_size(N);
    V.zeros();
    double f0=0;
    T.set_size(X.length());
    T.zeros();
    T=X;
    for (int i=0; i<Ndft;i++)
    {
        V=X.get(i*N,(i+1)*N-1);
        for(int j=0;j<N;j++)
        {
            if(x<V[0]) f0=0;
            else if (V[j]<=x && x<V[j+1]) f0=j/N;
            else if (x>=V[N-1])      f0=1;
            Z[i*N+j]=f0;
        }
    }
    return f0;
}

提前感谢任何帮助!

  1. 目前还不清楚为什么需要NNdftX.length()。一个值还不够吗?

  2. 正如有人提到的那样,在整数除法j/N你得到零。请把它变成double(j)/double(N)