目前分類:R與統計分析 (4)

瀏覽方式: 標題列表 簡短摘要

 

當兩個變項具有線性相關時,可將兩個變項的直線關聯性以一個直線型迴歸函數表示。常見的簡單線性迴歸模型中只有一個解釋變項,其模式可表示為Yi=a+bXi+ei (Yi為依變項、Xi為自變項/解釋變項、a為截距、b為迴歸係數、ei為殘差) ,其中研究者所收集到的資料爲(Yi, Xi) ,其主要的目的即爲以收集到的依變項及解釋變項的資料來估計出未知的參數a, b

 

實例說明:腰圍對空腹血糖值的解釋力為何?

 

[程式]

 

>model=lm(fpg~waist) #建立迴規模式lm(Y~X),並命名為model

>summary(model) #檢視迴歸統計量

>plot(waist,fpg) #製作散佈散佈圖plot(X,Y)

文章標籤

estat 發表在 痞客邦 留言(2) 人氣()

 

獨立樣本t檢定(independent t test)

研究者欲比較兩組資料的平均值是否有統計上的顯著差異,且兩組資料間必須是獨立樣本。例如:比較實驗組與對照組的身高是否相同。

此範例為比較男生與女生的BMI值是否相同,由於男生與女生是獨立的樣本,故將使用獨立樣本t檢定來分析。


[程式]

>sexbmi=split(bmi,sex) #.將資料分割不同性別BMI分佈情形,並加以命名

>var.test(sexbmi$"1",sexbmi$"2") #.檢定兩組變異數是否相同?

>t.test(sexbmi$"1",sexbmi$"2",var.equal=T) #.檢定兩組樣本平均數是否相同(變異數相等)

>t.test(sexbmi$"1",sexbmi$"2",var.equal=F) #.檢定兩組樣本平均數是否相同(變異數不相等)

文章標籤

estat 發表在 痞客邦 留言(0) 人氣()

 
費雪精確檢定(Fisher’s exact test)
 
2x2列聯表(自由度為1)分析計算卡方值時,若細格內的理論期望次數小於5,且樣本數小於20,雖用葉式校正檢定(Yate’s correction)可校正,但校正結果較為保守,也不一定精確,且其只限於2x2列聯表,若遇到r x c列聯表(例如:2x3, 3x5)時,該如何校正較爲恰當?此時,使用Fisher’s exact test來分析是一個較合適的方法。但分析時該特別注意的是,因爲Fisher’s exact test是以計算機率值來檢定是否顯著,所以當有些格子內數字偏大時,會花費較多的計算時間。此外,費雪精確檢定其檢定的虚無假設爲OR=1(Odds Ratio, 勝算比),換句話說,當p< 0.05時所代表的意義爲兩變項之間,存在有相關性。
 

以下的例子來說明如何以R來執行Fisher’s exact test並列出(Odds Ratio, 勝算比)
 
---研究目的:壓力與飲酒是否相關?
文章標籤

estat 發表在 痞客邦 留言(0) 人氣()

 

卡方檢定(Chi-Square test)


主要適用於兩個類別或次序變數的差異分析,常用的以下四種:(1)適合度檢定(Goodness-of-fit test)(2)同質性檢定、(3)獨立性檢定、(4)顯著性檢定,其中較常見的分析為獨立性檢定,其目的在考驗兩個類別變項是否相關?例如:運動與心血管疾病是否相關?



[
程式]


文章標籤

estat 發表在 痞客邦 留言(0) 人氣()