| 10 | 1/1 | 返回列表 |
| 查看: 2748 | 回復(fù): 9 | ||
monitor2885至尊木蟲 (職業(yè)作家)
隊(duì)長
|
[求助]
數(shù)據(jù)是否符合正態(tài)分布
|
|
s=[1.43;1.48;1.15;1.57;1.96;3.46;3.49;1.35;2.17;2.4;2.43;2.68;2.38;1.43;1.42;1;1.24;1.59;1.67;1.32;2.43]; kurtosis(s)/std(s) ans = 4.0153 百度上說,峰度系數(shù)與其標(biāo)準(zhǔn)誤的比值用來檢驗(yàn)正態(tài)性。如果該比值絕對(duì)值大于2,將拒絕正態(tài)性。但是經(jīng)過kstest測試,數(shù)組s符合正態(tài)分布呀。請(qǐng)問,怎么回事? |

至尊木蟲 (職業(yè)作家)
隊(duì)長

木蟲 (小有名氣)
Matlab
|
>> s=[1.43;1.48;1.15;1.57;1.96;3.46;3.49;1.35;2.17;2.4;2.43;2.68;2.38;1.43;1.42;1;1.24;1.59;1.67;1.32;2.43]; >> p_judge(s,0.05) H1 = 0 s1 = 0.3186 該數(shù)據(jù)源服從正態(tài)分布。 H2 = 0 s2 = 0.4848 該數(shù)據(jù)源服從γ分布。 H3 = 1 s3 = 4.4963e-004 該數(shù)據(jù)源不服從泊松分布。 H4 = 1 s4 = 0.0011 該數(shù)據(jù)源不服從指數(shù)分布。 phat = 1.4361 pci = 1.1841 1.8253 H5 = 0 s5 = 0.1979 該數(shù)據(jù)源服從rayleigh分布。 >> 進(jìn)一步的看看 >> [h,p,istat,cv]=lillietest(s,0.05) h = 1 p = 0.0253 istat = 0.2013 cv = 0.1877 很明顯拒絕原假設(shè);即 數(shù)據(jù)在0.05的置信水平下不是正態(tài)分布;承認(rèn)原假設(shè)的概率為2.53%左右 >> subplot(211) >> normplot(s) >> subplot(212) >> hist(s) 由圖形可以明顯看出,不服從分布檢驗(yàn)的 所以我去查看了kstest 的幫助文檔 終于知道為什么了:這個(gè)函數(shù)本來就不適合做正態(tài)分布的檢驗(yàn).原文是: The Kolmogorov-Smirnov test requires that CDF be predetermined. It is not accurate if CDF is estimated from the data. To test x against a normal distribution without specifying the parameters, use lillietest instead. |

木蟲 (小有名氣)
Matlab

至尊木蟲 (職業(yè)作家)
隊(duì)長

木蟲 (小有名氣)
Matlab
|
呃 兄弟不是吧 偶的英語啊 慘不忍睹啊 help kstest 就可以看見具體的了 更進(jìn)一步的解釋是 The function KSTEST by default, performs a two-tailed test. In case of such tests, if the significance level is alpha (0.05 by default), the null hypothesis is rejected if the P-value is less than alpha/2. To account for the tail in the computations, you have to specify the 4th input to the KSTEST function ('tail'). 默認(rèn)的是5%水平的雙邊檢驗(yàn) 在缺省條件下, 如果P值小于α/ 2,零假設(shè)被拒絕...CDF啊,那個(gè)是一個(gè)分布概率啊 你看看kstest的用法就明白了 |

至尊木蟲 (職業(yè)作家)
隊(duì)長

木蟲 (小有名氣)
Matlab

至尊木蟲 (職業(yè)作家)
隊(duì)長

至尊木蟲 (職業(yè)作家)
隊(duì)長

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專業(yè)求調(diào)劑 +9 | hanamiko 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 335分 | 材料與化工專碩 | GPA 4.07 | 有科研經(jīng)歷 +3 | cccchenso 2026-03-23 | 3/150 |
|
|
[考研] 081700 調(diào)劑 267分 +5 | 迷人的哈哈 2026-03-23 | 5/250 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|