data(standtabCoihue)
df<-standtabCoihue
df<-df[-nrow(df),]
# Diameter distribution plot
barplot(df$nha, legend = TRUE, beside = TRUE, las=1,
names.arg = as.numeric(df$diam.cl),
ylab="Density (trees/ha)",xlab="Diameter class (cm)")
abline(h=0)
Run the code above in your browser using DataLab