I don't understand what you are looking for. The line barplot(t(ld)) is placing the four treatments side by side, and is also stacking the three values of typ. Please indicate what you would like the figure to look like. Rich