# rm(list=ls()) library(RcmdrPlugin.TeachingDemos) # (nn = c(11,9,13,11,8,5,10,6,13,8,8,9)) sum(nn) chisq.test(nn,p=c(31,28,31,30,31,30,31,31,30,31,30,31)/365) # # reprezentativnost výběru # (nn = c(56,72,54,18)) (n = sum(nn)) chisq.test(nn,p=c(35,35,20,10)/100) (ocek = n*c(0.35,0.35,0.20,0.10)) nn-ocek (nn-ocek)^2 (nn-ocek)^2/ocek sum((nn-ocek)^2/ocek) qchisq(0.95,3) # kritická hodnota # co když četnosti poloviční: (nn2 = nn/2) (n2 = sum(nn2)) (ocek2 = n2*c(0.35,0.35,0.20,0.10)) nn2-ocek2 (nn2-ocek2)^2 (nn2-ocek2)^2/ocek2 sum((nn2-ocek2)^2/ocek2) chisq.test(nn2,p=c(35,35,20,10)/100) # už NENÍ VÝZNAMNÉ!!! # # barva květů a tvar pylových zrnek (nn = c(296,27,19,85)) (n = sum(nn)) (ocek = n*c(9,3,3,1)/16) (nn-ocek)^2/ocek sum((nn-ocek)^2/ocek) chisq.test(nn,p=c(9,3,3,1)/16) # (tab = matrix(nn,nrow=2,ncol=2,byrow=TRUE)) rownames(tab) = c("oválný tvar","kulatý tvar") colnames(tab) = c("purpurová barva","červená barva") addmargins(tab) (nn1 = addmargins(tab,1)[3,]) # číslo rozměru, který ZŮSTANE chisq.test(nn1,p=c(3,1)/4) (nn2 = addmargins(tab,2)[,3]) # číslo rozměru, který ZŮSTANE chisq.test(nn2,p=c(3,1)/4) # # příklad antigen # (nn = c(18,17,6)) (thHat <- (2*nn[1]+nn[2])/(2*sum(nn))) (piHat = c(thHat^2,2*thHat*(1-thHat),(1-thHat)^2)) chisq.test(c(18,17,6),p=piHat) # POZOR, p-hodnota špatně!! 1-pchisq(chisq.test(c(18,17,6),p=piHat)$statistic,1) # # Ichs (t = matrix(c(14,11,14,78,55,28,24,189,55,44,24,175,73,42,17,106),4,4)) chisq.test(t) # Baden (t = matrix(c(1768,946,115,807,1387,438,189,746,288,47,53,16),3,4)) chisq.test(t) # krevní skupiny (t=matrix(c(121,118,120,95,79,121,33,30),2,4)) chisq.test(t) # McNemarův test t = matrix(c(4,7,1,3,21,15,3,11,35),3,3) mcnemar.test(t) # hraboš t = matrix(c(4,11,27,473),2,2) chisq.test(t) chisq.test(t,correct=FALSE) fisher.test(t) mcnemar.test(matrix(c(4,11,27,473),2,2),correct=FALSE)