作者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