# popisné statistiky 1 ?read.csv2 help(read.csv2) read.csv2("deti23.csv") Deti23 = read.csv2("deti23.csv") ls() # dim(Deti23) # summary(Deti23) # Deti23$vekMatky # attach(Deti23) c(mean(vekOtce),median(vekOtce),min(vekOtce),max(vekOtce)) quantile(vekOtce,(0:4)/4) # sort(vekOtce) length(vekOtce) sort(vekOtce)[12] # order(vekOtce) # vekOtce[order(vekOtce)] Deti23[order(vekOtce),] # rank(vekOtce) # source("psGrafR.R") psGrafR("../boxVekOtce.ps",width=6.35,height=4.48/1.5) boxplot(vekOtce,horizontal=TRUE) dev.off() # quantile(vekOtce,(0:4)/4) # c(var(vekOtce),sd(vekOtce),sd(vekOtce)^2) quantile(vekOtce,3/4)-quantile(vekOtce,1/4) # # závislost # plot(vekOtce~vekMatky,data=Deti23) psGrafR("../primka1.ps") plot(vekOtce~vekMatky) abline(0,1,lty=3) dev.off() # plot(vekOtce~vekMatky) abline(0,1,lty=3) identify(vekMatky,vekOtce) # cor(vekOtce,vekMatky) # psGrafR("../HCDgender.ps") plot(HCD~Gender) dev.off() # xtabs(~Gender+HCD) t = table(Gender,HCD) t prop.table(t) prop.table(t,1) prop.table(t,2) # psGrafR("../boxy1.ps") plot(vaha~Gender,col="yellow") dev.off() # tapply(vaha,Gender,length) tapply(vaha,Gender,mean) tapply(vaha,Gender,sd) tapply(vaha,Gender,summary)