顯示具有 統計 標籤的文章。 顯示所有文章
顯示具有 統計 標籤的文章。 顯示所有文章

2016-02-22

用 R 算眾數

今天心血來潮,想重新學一下 R 來跑統計(大概重新學過 3 次了 ...)。用金門大學 陳鍾誠 老師的教材當範例:
習題 1 : 請模擬從班上隨機抽學生一萬次,看看誰與你最有緣 (抽到次數最多)
解答:假設班上有 59 人,那麼編為 1 到 59 號,於是我們可以用下列程式,進行 1 萬次抽樣, 並繪出統計圖。
 
> x=sample(1:59, 10000, TRUE)> hist(x, breaks=0.5:60)執行結果
想說來跑一下眾數,沒想到基本 package 沒有這個功能,冏。只好開始 google。最後大概有兩種公式,簡單整理一下。
第一個
as.numeric(names(table(x)))[table(x) == max(table(x))]

連結:來源一來源二
第二個
which.max(table(x))
連結:來源一來源二來源三
第三個
mode_dat=NULL
for(i in 1:1000){
dat=rnorm(1000)
ds=density(dat)
mode_dat=c(mode_dat,ds$x[which(ds$y==max(ds$y))])
}
連結:來源(看不懂
打完收工
借放:ANOVA


2016-01-24

R 資源整理



又到了說想學 R 的時候(大概說了三次以上了吧)
想先將 R 的相關資源記錄下來,之後再慢慢整理出自己學習 R 的歷程

目前

  • 有四本,先記錄一下

過去整理


部落格




2015-11-21

統計概念整理筆記

動機

我一直知道自己的統計能力滿薄弱的(汗顏)。但也因此,我遇到有疑惑的觀念,就會去問人或是去google。前陣子因為想嘗試WinGanova,卻發現出來的報表我不太會解讀(沒看到事後比較的部份)。想說自己對於事後比較使用的時間好像還不太明白,這幾天就又開始google...。總之,我想要替自己的統計觀念做個記錄,透過電腦玩物上介紹的筆記3層空間。一方面我可以動態的整理,也可以將資料放在網路上。或許會對將來的某些人有幫助...吧

↓以下是筆記


目前 or 未來想釐清的事情

  • 熟悉Excel、R、SPSS、WinGanova的操作
    • Excel:語法、圖表、VB
    • R:客製化出R的anova報表、圖形
    • SPSS:大概多跑幾次就記住了。可能多了解其他統計方法吧
    • WinGanova:大概多跑幾次就記住了

  • 整理變異數分析的流程、報表與圖形(單/雙、獨立/相依 的流程不太一樣...)
    1. 樣本特性描述
      • 數據:平均數、標準差
      • 作圖:通常不需要
    2. 預設的假設檢定
      • 概念
        • 書上看到的(行為科學統計學,p506)
        • 報表整理出來的
          • 變異數同質檢定(Levene統計量、顯著性)
          • 球面檢定(Mauchly'sW 係數、卡方分配值、P值)
      • 數據:痾,顯著性P值?
      • 作圖:不需要
    3. 整體考驗
      • 概念
        • 依據因子多寡,分成單因子、多因子實驗設計
        • 依據因子各水準的相關性,分成受試者內因子(相依樣本)、受試者間因子(獨立樣本)
        • 綜合以上,可分成完全受試者間、完全受試者內、以及混合設計三種實驗設計方式。
      • 數據:Df ,Sum Sq,Mean Sq, F value,Pr(>F)
      • 作圖:做表(變異數摘要表)、做交互作用圖(bar or line)
    4. 事後考驗
      • 單純主要效果(如果有的話,和變異數摘要表類似)
        • 數據:和整體考驗類似
        • 作圖:做表(變異數摘要表)、作圖(看狀況)
      • 事後檢定(如果有的話,應該會兩兩比較)
        • 數據:兩兩比較的顯著性(?)
        • 作圖:做表(通常和單純主要效果放在一起)


概念整理


描述統計


推論統計(研究生2.0:哪種統計方法好體育學院PPT

  • T檢定(R演習室
    • 單一樣本T檢定
    • 雙樣本 相依樣本T檢定
    • 雙樣本 獨立樣本T檢定
  • 變異數分析(網友的課堂筆記
    • 型1型4平方和問題(量化研究與統計分析 p8-14)
      • 型1平方和:階層化拆解,後進的因子不受之前的因子影響,多用在共變數分析、多項式回歸模式
      • 型1平方和:適用於只有主要效果的模型
      • 型1平方和:最嚴格,適用各組人數不等時
      • 型1平方和:適用有空白細格
    • 因子多寡
      • 單因子
      • 多因子
    • 樣本特性
      • 獨立
      • 相關
        • 單因子相關設計,變異數 SSt = SSb.subject + SSw
          SSw = SSa + SSr
        • SSt = 總變異、SSb = 受試者間(組間)的變異、SSw = 受試者內的變異、SSa = 受試者內因子造成的變異、SSr = 受試者內個人的誤差
    • 多重比較(研究生2.0中興生物系統工程 陳加忠
      • 因為F檢定是看整個模型的顯著程度,沒辦法看到像是A因子裡頭,A1 A2 A3 彼此之間是否有差異,因此必須進行事後檢定(事前檢定比較少用到)。
      • 除了one way需要,我覺得two way檢驗單純主要效果後,為了了解固定B因子情況下,A因子各個水準差異如何,應該也會用到事後檢定...(p9-7 提到的決策數,提到確實會有多重比較...)


相關資料封存

統計影片

統計書本

其他

  • 當初高統想記錄高統的網誌 ... XD(高統之記錄
  • http://www.minitab.com.tw/product/minitab/quick-start_weather.php


心得
  • 151121 不曉得這次熱情可以撐多久 ... orz

2012-12-21

集中趨勢:平均數中位數眾數

定義

平均數 (mean) 所有數字相加,除以個數 集中趨勢的指標
中位數 (median) 將數字從小到大排序,取中間數 集中趨勢的指標
眾數 (mode) 集合中出現頻率最多的數字 集中趨勢的指標

 

例子

一數字集合 U 有八個數字

U:{ 23, 29, 20, 32, 23, 21, 33, 25 }

平均數 (mean) 相加除個數
(23 + 29 + 20 + ... +25 ) / 8 = 25.75
中位數 (median) 排序後,中間兩值相除
U':{20 ... 23, 25, ... 33}
(23+25) /2=24
眾數 (mode) 出現頻率最高的數為 23