| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2762 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
monitor2885至尊木蟲(chóng) (職業(yè)作家)
隊(duì)長(zhǎng)
|
[求助]
數(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 百度上說(shuō),峰度系數(shù)與其標(biāo)準(zhǔn)誤的比值用來(lái)檢驗(yàn)正態(tài)性。如果該比值絕對(duì)值大于2,將拒絕正態(tài)性。但是經(jīng)過(guò)kstest測(cè)試,數(shù)組s符合正態(tài)分布呀。請(qǐng)問(wèn),怎么回事? |

木蟲(chóng) (小有名氣)
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ù)本來(lái)就不適合做正態(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. |

木蟲(chóng) (小有名氣)
Matlab

至尊木蟲(chóng) (職業(yè)作家)
隊(duì)長(zhǎng)

木蟲(chóng) (小有名氣)
Matlab
|
呃 兄弟不是吧 偶的英語(yǔ)啊 慘不忍睹啊 help kstest 就可以看見(jiàn)具體的了 更進(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的用法就明白了 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600,材料與化工321分求調(diào)劑 +11 | 大饞小子 2026-03-28 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 材料專碩322分 +9 | 哈哈哈吼吼吼哈 2026-04-02 | 9/450 |
|
|
[考研] 321求調(diào)劑 一志愿 浙江工業(yè)大學(xué)生物醫(yī)藥 +5 | 嘿嘿HC 2026-04-01 | 6/300 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +12 | dongfang59 2026-03-27 | 12/600 |
|
|
[考研] 292分,材料與化工,申請(qǐng)調(diào)劑 +19 | 程晴之 2026-04-01 | 23/1150 |
|
|
[考研] 計(jì)算機(jī)265跨調(diào)環(huán)境 +5 | Yumeng_6 2026-03-27 | 5/250 |
|
|
[考研] 一志愿北交大材料工程總分358 +8 | cs0106 2026-04-01 | 9/450 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +13 | Koxui 2026-03-28 | 13/650 |
|
|
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7 | czhcz 2026-03-31 | 9/450 |
|
|
[考研] 材料科學(xué)與工程339求調(diào)劑 +11 | hyz0119 2026-03-31 | 12/600 |
|
|
[考研] 309求調(diào)劑 +19 | 誰(shuí)不是少年 2026-03-29 | 19/950 |
|
|
[考研] 322求調(diào)劑 +8 | 三水sss 2026-04-01 | 8/400 |
|
|
[考研] 309分085801求調(diào)劑 +7 | 學(xué)員Gtwj7W 2026-03-31 | 7/350 |
|
|
[考研] 一志愿a區(qū)211,085601-307分求調(diào)劑 +10 | 黨嘉豪 2026-03-31 | 23/1150 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過(guò)六級(jí)。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +8 | 哇呼哼呼哼 2026-03-27 | 10/500 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|