| 10 | 1/1 | 返回列表 |
| 查看: 2749 | 回復(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)性。如果該比值絕對值大于2,將拒絕正態(tài)性。但是經(jīng)過kstest測試,數(shù)組s符合正態(tài)分布呀。請問,怎么回事? |

至尊木蟲 (職業(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ā)表 | |
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +7 | 紙魚ly 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 081700 調(diào)劑 267分 +7 | 迷人的哈哈 2026-03-23 | 7/350 |
|
|
[考研] 材料調(diào)劑 +3 | 匹克i 2026-03-23 | 3/150 |
|
|
[考研] 265求調(diào)劑 +10 | 梁梁校校 2026-03-17 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|