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