# 8. přednáška # library(RcmdrPlugin.TeachingDemos) # # dvouvýběrový t-test a intervaly spolehlivosti data(Kojeni) # omezit se na maturantky attach(Kojeni[Kojeni$Vzdelani=="maturita",]) # nyní přímo uvidíme jen informace o maturantkách search() tapply(hmotnost,Hoch,shapiro.test) var.test(hmotnost~Hoch) t.test(hmotnost~Hoch,var.equal=TRUE) (tDivka<-t.test(hmotnost[Hoch=="dívka"])) (tHoch<-t.test(hmotnost[Hoch=="hoch"])) source("errorBars.R") errorBars(hmotnost,Hoch,type="CI") # # wilcox.test(hosi,divky) wilcox.test(hosi,divky,corr=FALSE) (srovnej = outer(hosi,divky,"-")) (W = sum(srovnej>0)+sum(srovnej==0)/2) ks.test(hosi,divky) # plot(c(min(hosi),sort(hosi)),(0:15)/15,type="s",xlim=range(Vysky$vyska), ylim=0:1,xlab="výšky",ylab="",col=4) lines(c(min(divky),sort(divky)),(0:12)/12,type="s",col=2) abline(h=0:1) arrows(140.5,10/15,140.5,3/12,length=0.05,code=3) text(139.75,0.5,"D") # # analýza rozptylu # data(Med) plot(cu~Misto,data=Med,ylab="Cu",col="yellow") plot(lnCu~Misto,data=Med,ylab="log(Cu)",col="yellow") anova(lm(lnCu~Misto,data=Med)) summary(aov(lnCu~Misto,data=Med)) bartlett.test(lnCu~Misto,data=Med) library(car) levene.test(lnCu~Misto,data=Med) shapiro.test(resid(lm(lnCu~Misto,data=Med))) shapiro.test(resid(aov(lnCu~Misto,data=Med))) TukeyHSD(aov(lnCu~Misto,data=Med)) plot(TukeyHSD(aov(lnCu~Misto,data=Med))) with(Med,round(tapply(lnCu,Misto,mean),3)) # plot(vek.m~Vzdelani,data=Kojeni,xlab="vzdě lání",ylab="věk",col="yellow") kruskal.test(vek.m~Vzdelani,data=Kojeni)# s přihlédnutím ke shodám