f分布的檢驗 r語言 - CSDN

2021-01-16 CSDN技術社區

統計學中的t檢驗法和F檢驗法的應用條件是樣本都來自正態總體或近似正態總體,只有符合這個條件,才能用它們來檢驗各樣本所屬的總體參數的差異顯著性。


一、非參數檢驗

1、Kolmogorov-Smirnov正態性檢驗(單樣本)

檢驗單一樣本是否來自某一特定分布。比如檢驗一組數據是否為正態分布。它的檢驗方法是以樣本數據的累積頻數分布與特定理論分布比較,若兩者間的差距很小,則推論該樣本取自某特定分布族。即對於假設檢驗問題:

H0:樣本所來自的總體分布服從某特定分布

H1:樣本所來自的總體分布不服從某特定分布


Fn(X)表示一組隨機樣本的累計概率函數,F0(X)表示分布的分布函數。

當原假設為真時,D的值應較小,若過大,則懷疑原假設,從而,拒絕域為      

對於給定的又,

R實現:

1.1 ks.test()

> x <- rnorm(3000)> ks.test(x,"pnorm") One-sample Kolmogorov-Smirnov testdata: xD = 0.0077, p-value = 0.9943alternative hypothesis: two-sided

若不是正態分布函數的檢驗,只需要將"pnorm"改成相應地分布函數

當已知分布函數時,如:

> ks.test(x,"pnorm",mean(x),sd(x)) One-sample Kolmogorov-Smirnov testdata: xD = 0.0087, p-value = 0.9768alternative hypothesis: two-sided

若出現錯誤或警報,可能是由於:

Warning message:

In ks.test(X, "pnorm", 2.488715,3.012176) :

 ties should not be present for the Kolmogorov-Smirnov test

#若已知可能的分布函數

#樣本數據中存在有相同的值,單樣本K-S檢驗要求檢驗分布是連續的,而連續分布出現相同值的概率為0.如果是出現相同的,則連續分布的假設不成立,則該方法無法使用

1.2 nortest包裡的lillie.test()可以實行更精確的Kolmogorov-Smirnov檢驗

> library(nortest)> lillie.test(x) Lilliefors (Kolmogorov-Smirnov) normality testdata: xD = 0.0087, p-value = 0.8408

可以發現此方法的結果和ks.test的第二種方法的結果是一樣的

1.3 fBasics包裡的normaltest()和ksnormtest()

> ksnormTest(x)Title: One-sample Kolmogorov-Smirnov testTest Results: STATISTIC: D: 0.0077 P VALUE: Alternative Two-Sided: 0.9943 Alternative Less: 0.7946 Alternative Greater: 0.7012 Description: Tue Jan 13 21:24:19 2015 by user:



#建議:大樣本、已知總體均數和標準差,選擇非參數檢驗-單樣本KS檢驗號。。但對於樣本大小,有說大於50、500、5000的,各有所不同的理解。


2、Shapiro-Wilk(W檢驗)

統計量為:


當原假設為真時,W的值應接近於1,若值過小,則懷疑原假設,從而拒絕域為R:W<c

在給定的水平下P{w<c}=a

R實現:

2.1 shapiro.test()實現

> shapiro.test(x) Shapiro-Wilk normality testdata: xW = 0.9997, p-value = 0.9494

在R中,shpiro.test()函數限制為3-5000個樣本之間
2.2 fBasics包裡的shapiroTest()

> shapiroTest(x)Title: Shapiro - Wilk Normality TestTest Results: STATISTIC: W: 0.9997 P VALUE: 0.9494



二、圖示法

Q-Q圖:以樣本的分位數作為橫坐標,以按照正態分布計算的相應分位點作為縱坐標,把樣本表現為直角坐標系的散點。如果數據服從正態分布,則樣本點應圍繞第一象限的對角線分布。

w <- c(75.0, 64.0, 47.4, 66.9, 62.2, 62.2, 58.7, 63.5,
       66.6, 64.0, 57.0, 69.0, 56.9, 50.0, 72.0)
qqnorm(w); qqline(w)


可以發現:基本上呈正態分布。



相關焦點

  • f檢驗 matlab專題及常見問題 - CSDN
    Mann-Kendall顯著性檢驗工具:1、MATLAB2、DPS3、示例數據(數據採用的是魏鳳英老師《現代氣候統計診斷預測技術》中的1900-1990年上海市的年平均氣溫數據)Mann-Kendall原理可詳細參考https://wenku.baidu.com
  • 什麼是f檢驗 - CSDN
    方差檢驗是多變量t檢驗的延續,對於超過兩個樣本的對比檢驗就無法直接使用獨立T檢驗了,這個時候就需要使用卡方檢驗。例子:冰淇凌老闆想知道三種口味的冰淇凌的銷售情況是否一樣,他有如下的數據巧克力味草莓味原味233234321233343344等等等F檢驗又叫方差齊性檢驗,目的是判斷兩個樣本的總體方差是否相等,計算雙總體樣本檢驗的前提條件。
  • 正態分布的常用數據 - CSDN
    常用的正態分布檢驗方法有以下幾種:1.基於偏度和峰度的假設檢驗基於偏度-峰度的檢驗是利用了正態分布偏度(3階矩)和峰度(4階矩)都為0的特點。如果樣本數據能滿足偏度和峰度均為0的假設,則可以認為總體服從正態分布。由於該檢驗是基於偏度檢驗和峰度檢驗的,樣本數量需要8個以上。以下normaltest函數就使用該原理進行正態分布檢驗。
  • stata單位根檢驗 - CSDN
    其中,服從獨立同分布的N(0,)分布。然而,由於 的漸近分布不是標準分布,傳統的 t 檢驗在這裡並不適用。此外,在回歸中是否包含常數或時間趨勢項會導致檢驗統計量有不同的漸近分布,這也意味著在做這類檢驗時清楚設定原假設和備擇假設是非常重要的。4.1 ADF 檢驗在單位根檢驗的原假設下,真實的過程要麼是隨機遊走,要麼是帶漂移項的隨機遊走。
  • 多個變量結果的正態分布 - CSDN
    多元正態分布轉自個人微信公眾號【Memo_Cleon】的統計學習筆記:多元正態分布檢驗的R實現方法。程序包MVN【函數mvn】、程序包mvnTest【函數AD.test、CM.test、DH.test、HZ.test、R.test、S2.test等】、程序包mvtnorm【函數pmvnorm】、程序包energy【函數mvnorm.e、mvnorm.test、mvnorm.etest】……我們以mshapiro.test {mvnormtest}和mvn {MVN}為例演示多元正態分布檢驗的
  • R語言從入門到精通:Day12--R語言統計--回歸分析
    回歸作為一個廣義的概念,涵蓋了許多變種,R語言中也為其提供了強大而豐富的函數和選項(但顯然選項越多,對初學者越不友好),早在2005年,R中就有200多種關於回歸分析的函數 (https://cran.r-project.org/doc/contrib/Ricci-refcard-regression.pdf,這個文檔提供了部分回歸分析函數列表,供大家參考)。
  • android啟動頁設計專題及常見問題 - CSDN
    轉載請註明出處:http://blog.csdn.net/wangjihuanghun/article/details/63255144啟動頁幾乎成為了每個app的標配,有些商家在啟動頁中增加了開屏廣告以此帶來更多的收入。
  • 協整檢驗數據分析 - CSDN
    我們按面板數據的長度,可以分為長面板數據和短面板數據,長面板數據的單位根檢驗方法為:LLC檢驗、Breitung檢驗、IPS檢驗、Fisher式檢驗、HadriLM檢驗。短面板數據的單位根檢驗方法為:HT檢驗。
  • 臺灣地區語言分布概況
    在臺灣地區,曾經出現和現在使用的語言包括:南島語系、漢藏語系、阿爾泰語系和日本語系。泰雅語群(Atayalic)泰雅語(Atayal)分布最廣,從烏來到南投。賽考利克語(Squli)分布在烏來、新竹、南投等地澤敖利語(C'uli')新竹汶水語宜蘭萬大語
  • 線性回歸方程的顯著性驗證,總體驗證的F檢驗與個體驗證的t檢驗
    許栩原創專欄《從入門到高手:線性回歸分析詳解》第8章,顯著性驗證,總體驗證的F檢驗與個體驗證的t檢驗。上一章,我講述了回歸方程的精度,在回歸分析中,我們求出回歸方程後,除了確認回歸方程的精度外,我們要需要對回歸方程進行顯著性驗證,以確認回歸方程的有效性。
  • 英語中狗狗竟然不會「汪汪」叫,小動物間也會有語言障礙!?
    woof 英[wʊf] 美[wʊf]  英語中的鴨子:Quack quack!漢語中的鴨子:嘎嘎嘎! quack 英[kwæk] 美[kwæk] 英語中的青蛙:Ribbit ribbit!漢語中的青蛙:呱呱呱!
  • 專欄F|Cora單詞20動物(下):生活中不可或缺的夥伴
    monkey [ˈmʌŋki] 猴 panda [ˈpændə] 熊貓 bear [beə(r)] 熊 lion [ˈlaɪən] 獅子 tiger [ˈtaɪɡə(r)] 老虎 fox [fɒks] 狐狸 zebra [ˈzebrə] 斑馬 deer [dɪə(r)] 鹿 giraffe
  • 高考英語3500詞彙(R):單詞+詞組+句型
    使改變, 使轉移reduce [rɪˈdjuːs; (US) -ˈduːs] vt. 減少, 縮減reduce the price by 10 percent.價格降低了百分之十Hunger reduced him to stealing.飢餓使他開始偷東西。refer [rɪˈfɜː(r)] vi.
  • 機器學習 · R語言包大全(共99個包)
    本文譯自2017年5月30日Torsten Hothorn整理的《CRAN Task View: Machine Learning & Statistical Learning》,英文原文連結:https://cran.r-project.org
  • R語言 馬賽克圖
    馬賽克圖是《R語言實戰》書籍開篇中的彩圖,偶然遇到了學習分享一下~Library(vcd)mosaic