看板 FJU-STAT95B 關於我們 聯絡資訊
================================以下開始複製================================== #include <C:\devcpp\all.h> int main() { long int i, j; srand((unsigned)time(NULL)); cout.setf(ios::fixed); cout.precision(10); // $$$$$$$$ Start your program here $$$$$$$$$$$$$$$$$$$$$ Matrix Y; Y.read_from_file("c:\\devcpp\\anova.dat",2); Y.print(3); Univar sum_of_rank, sample_size, T_stat; Kruskal_Wallis_Test( Y, sum_of_rank, sample_size, T_stat ); Y.print(3); sum_of_rank.print(1); sample_size.print(0); T_stat.print(4); cout.precision(6); double P_value = 1.0 - X2_CDF(T_stat,3); cout << "P_value = " << P_value << endl; Univar KW; for(i=0; i < 10000; i++) { Matrix Y2(21,1); Y2.uni_ran(10,30); Y2.append_column(Y.col(1)); Kruskal_Wallis_Test( Y2, sum_of_rank, sample_size, T_stat ); KW.push_back(T_stat); } KW.freq_dist(); // $$$$$$$ End of your program $$$$$$$$$$$$$$$$$$$$$$$$$ cin.get(); return 0; ===================================複製結束=================================== -- ▄▄▄▄▄◣ ◢▄▄▄▄▄ KEI 沒穿過 acer 水皎嫂 凱馨食品建大 雅芳食品 Mizuno 這套戰服 雲林尚讚 Mizuno南投最美▇▇ 九泰光電 廣吉食品 ///小美 別說你 inon 參加過 ittt 優 美 興家安速 F-1 ▄▄▄▄▄◣ ◢▄▄▄▄▄ KEI 沒穿過 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.136.161.81
abidog:水拉!! 05/25 14:51
※ 編輯: weiching0323 來自: 140.136.161.81 (05/25 15:21)