double x[14] =
{99,100,90,94,135,108,111,119,104,127,109,117,105,125};
Univar X(14,x); X.print(0);
double TP, TN, T, Z; Univar RP, RN;
Wilcoxon_signed_rank( X, 107, RP, RN, TP, TN, T, Z);
RP.print(1); RN.print(1);
cout.precision(4);
cout << TP << " " << TN << " " << T << " " << Z << " " << endl;
double P_value = 2.0*S_nor_CDF(Z);
cout << "P_value = " << P_value << endl;
cout << "Parametric Approach" << endl;
double t = X.t_stat(107); cout << "t_stat = " << t << endl;
P_value =t_2_tailed_prob(t, 13);
cout << "P_value for t test = " << P_value << endl <<endl;
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.136.161.82