作者fr407041 (so sad so blue)
看板NCTU-STAT96G
標題[心得] imsl chi-squar檢定
時間Wed Aug 1 22:16:41 2007
痾 我不知道我這樣跑對不對
但我讓它complier過了
所以想跟大家討論一下
int N_CATEGORIES = 10;//我猜這應該是把分配分幾塊的意思
IMSLS_F_CHI_SQUARED_TEST imsls_f_chi_squared_test;
IMSLS_F_NORMAL_CDF imsls_f_normal_cdf;
imsls_f_normal_cdf =
(IMSLS_F_NORMAL_CDF)GetProcAddress(hStatLib,"imsls_f_normal_cdf");
imsls_f_chi_squared_test =
(IMSLS_F_CHI_SQUARED_TEST)GetProcAddress(hStatLib,"imsls_f_chi_squared_test");
float p_value ;
p_value = imsls_f_chi_squared_test ((float *)imsls_f_normal_cdf ,
N_OBSERVATIONS , N_CATEGORIES ,y , 0);
printf("pvalue = %f\n", p_value);
我出來的pvalue是0
--
一直以為幸福在遠方，在可以追逐的未來。
我的雙眼保持著眺望，我的雙耳仔細聆聽，
唯恐疏忽錯過。後來才發現，那些握過的手，
唱過的歌，流過的淚，愛過的人ꐊ
............所謂的曾經，就是幸福。☆席慕蓉☆
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.249.227
※ 編輯: fr407041 來自: 140.113.249.227 (08/01 22:20)
1F：推 uso1626:一樣0... 08/01 22:24
2F：推 fr407041:這跟Shapiro-Wilk TEST結果一樣 所以我認為有可能對 08/01 22:24
3F：推 pael:我也是 0 數字好像沒錯 因為我將資料拿去其它統計軟體檢驗 08/01 22:50
4F：推 pael:資料本身好像真的不服從 normal @@ 08/01 22:52
5F：推 fr407041:庭偉兄已經測出這可以跑 08/01 23:28
6F：推 pael:我剛用 IMSL 本身產生 Normal 分配並測試卡方檢定的功能 08/01 23:30
7F：推 pael:有不是 "0" 的 p-value 跑出來 所以上面應該真的是 "0" 08/01 23:31
9F：推 pael:當然你可也以任意的改 #define 裡的 seed 跟 number 試試 ^^" 08/01 23:34
10F：推 fr407041:感謝版主辛苦測試 08/01 23:36
11F：推 pael:別這麼說啦 只是小小的測試而已 也不見得就沒問題 @@ 08/02 00:56
12F：推 pd0320:感謝分享!! 終於知道錯在哪了 XD 08/02 02:16