>datuak=rnorm(40,mean=172,sd=10)
>datupolitak=round(datuak,digits=0)
>datupolitak
[1] 162 177 181 165 161 172 153 177 156 164 178 183 161 159 171 179 169 178 174
[20] 173 171 173 176 175 176 163 173 167 168 175 177 172 165 176 161 170 169 170
[39] 179 182
>library(agricolae)
>tarteak=seq(150,185,by=5)
>png("superhistogramapolita_01.png")
>par(mfrow = c(1,2),width=1100,heigth=700) 2 grafiko 2 zutabetan eta errenkada batean
>plot.graph.freq(datupolitak,col="lightblue",breaks=tarteak,
xlab="Altuerak",ylab="Ehunekoak",axes=FALSE,frequency=2) #'axes=False' bidez bi ardatzak kendu, gero moldatzeko
>axis(1,at=tarteak)
>axis(2,las=1,at=seq(0,0.30,by=0.05)
,labels = paste(seq(0,30,by=5), "%", sep = "")) #Ardatzeko maiztasun erlatiboak jarri portzentaje moduan
>abline(h=(seq(0,0.30,by=0.05)), col="orange", lty=3) #marra lagungarri horizontalak
>title(main="Altueren histograma") #izenburua aparte jar daiteke
>#orain, maiztasun poligonoa
>histogramapolita=plot.graph.freq(datupolitak,border=FALSE,
breaks=tarteak,xlab="Altuerak",ylab="Ehunekoak",
axes=FALSE,frequency=2) 
>title(main="Altueren maiztasun-poligonoa")
>axis(1,at=tarteak)
>polygon.freq(histogramapolita,lwd=5,col="orange",frequency=2)
>axis(2,las=1,at=seq(0,0.30,by=0.05),labels = paste(seq(0,30,by=5), "%", sep = ""))
>dev.off()