float skewness(float *p,int count) { float mea,med,sd,s; mea = mean(p,count); med =median(p,count); sd = standarddeviation(p,count); s = (3*(mea-med))/sd; return s; }
float kurtosis(float *p,int count) { float sd,sum=0,mea,k; int i; mea = mean(p,count); for(i=0; i<count; i++) { sum=sum+pow((p[i]-mea),4); } sd =standarddeviation(*p,count); k = sum/(pow(sd,4)*count-1); }
double HrPwWindow::corr(QVector<double> &Xi, QVector<double> &Yi,int n) { double r = covariance(Xi,Yi,n)/(standarddeviation(Xi,n)*standarddeviation(Yi,n)); //double r=slope(Xi,Yi,n)*slope(Xi,Yi,n)*(variance(Xi,n)/variance(Yi,n)); return (r); }