> neska=c(0,1,1,1,0,2,2,1,2,1,1,2,1,1,1,1,0,0,0,0,0,0,1,0,0)
>mutil=c(1,2,1,2,1,0,0,1,2,3,4,0,3,2,1,1,1,2,3,3,2,3,3,1,1,1,2,3,0,2,3,1,1,1)
> fneska=prop.table(table(neska))*100
> fmutil=prop.table(table(mutil))*100
> fneska
neska
 0  1  2 
40 44 16 
> fneska=c(40,44,16,0,0)
> dena=cbind(fneska,fmutil)
> bp=barplot(dena,beside=T,names.arg=c("Nesken bikotekideak","Mutilen bikotekideak"),col=c("green","green","green","green","green","lightblue","lightblue","lightblue","lightblue","lightblue"),ylim=c(0,50),las=1,ylab="Sexu bakoitzeko inkestatuen portzentajeak",main="Bikotekide kopurua sexuaren arabera")
> text(bp,3,paste("(%",c(fneska,round(fmutil,digits=1)),")"),cex=0.6,pos=3)
> text(bp,-0.5,c(0,1,2,3,4),cex=1,pos=3)