Proovi nii!
suguF<-factor(sugu,labels=c("mees","naine"))
terve<-as.numeric(I(tervis>0))
a<-barplot(prop.table(table(terve,suguF),2)[1,]*100,ylim=c(0,60),
main="Meeste ja naiste võrdlus", ylab="Tõenäosus mitte haigestuda, 95% usalduspiiridega",col=c("lightblue","orange"))
lim1 <- binom.test(table(terve[suguF=="mees"]))$conf.int*100
lim2 <- binom.test(table(terve[suguF=="naine"]))$conf.int*100
arrows(a, c(lim1[1], lim2[1]),
a, c(lim1[2], lim2[2]),
code=3, angle=90, col="gray40", lwd=2)
Tagasi