R 卡方檢定

適合度檢定（Goodness of fit test）

投擲一顆骰子一萬次，檢定六個面出現的機率是否相同

x<-c(sample.int(6,size=10000,replace=T)) #產生1~6的隨機整數一萬個
table(x) #計算次數
chisq.test(table(x)) #檢定出現的機率是否相同
chisq.test(table(x),p=c(1/6,1/6,1/6,1/6,1/6,1/6)) #結果跟上面一樣

529株豌豆分別有AA、Aa和aa各138、265、126株，檢定比例是否符合1：2：1

豌豆基因  AA  Aa   aa
    株數     138  256 126

x<-c(138,265,126)
chisq.test(x,p=c(1/4,1/2,1/4))






獨立性檢定（Test of independence）

某次米國總統大選的男女投票資料如下，檢定性別是否與投票請向有關

           Bush   Perot   Clinton
male     315     152       337
female  346     126       571

male<-c(Bush=315,Perot=152,Clinton=337)
female<-c(Bush=346,Perot=126,Clinton=571)
rbind(male,female) #輸入男女資料後，排成矩陣
chisq.test(rbind(male,female))