Hi R-helpers, Can anyone help me giving instructions/package sources to make a simple bar plot with mean and standard error of variables [Orange colour (which has four PCs) and Irid colour (which has two PCs)] that I'd like to plot in x-axis and the values in y-axis? The diet should be as legend (H & L). Dataframe: Diet Orange colour Irid colour PC1 PC2 PC3 PC4 PC1 PC2 H -3.4 -2.4 0.1 0.4 -2.4 -7.0 H -6.1 -0.7 -0.8 -0.5 5.5 -3.2 L 12.2 -1.1 -0.6 0.4 10.6 -4.6 L -3.4 -2.2 0.1 -0.1 -0.6 -3.2 I'll highly appretiate your cordial help in this regard. Cheers, Moshi [[alternative HTML version deleted]]