Hi NA, From the documentation here (See section boxplotGroup(__,'primaryLabels',c) ) the number of primary lables should match the number of boxplotGroups. i.e. you appear to be plotting three box plots (A,B,C) but trying to give these box plots eight names ('a' 'b' 'c' 'd' 'e' 'f' 'g' 'h'). I beleive that you are trying to show something like the last example in boxplotGroup(__,'secondaryLabels',s) (shown below) in which case swap your primary and secondary labels. Let me know if that helps, Christopher

The ability to plot the mean values using _boxplot_ is not available as of release R2016b. To work around this issue, you can find these values and plot them manually. The example below shows how to plot the mean value of each group: % Generate random data X = rand(10); % Create a new figure and draw a box plot figure; boxplot(X) % Overlay the mean as green diamonds hold on plot(mean(X), 'dg') hold off

First of all, if you use boxchart (released in R2020a) you'll have a bit more flexibility than the old boxplot function. The way I interpreted what you were asking for is, two grouping factors: one should be grouping by which column of the matrix, and the other should be grouping by which matrix. To accomplish this with boxchart it's easier to specify data as vectors rather than matrices, this let's you explicitly specify the xgroupdata argument and the cgroupdata argument. You'll need to choose which one you want to be represented by color and which you'll want for x - I think from your description you're looking for 7 sets of 3 boxes (so color is 'which matrix' and x is 'which column'. But if I got that wrong you can just reverse them. Finally, I added a little multiplier to the x so that there was a little space between each set of 3 boxes. I think that makes it a little easier to see, but it means having to set the XTicks/XTickLabels explitily. You could make the below code much more compact, I spread it out so you could see my logic: % Sample Data: trial1 = rand(5,7); trial2 = rand(10,7); trial3 = rand(15,7); % These grouping matrices label the columns: grp1 = repmat(1:7,size(trial1,1),1); grp2 = repmat(1:7,size(trial2,1),1); grp3 = repmat(1:7,size(trial3,1),1); % These color matrices label the matrix id: clr1 = repmat(1,size(trial1)); clr2 = repmat(2,size(trial2)); clr3 = repmat(3,size(trial3)); % Combine the above matrices into one for x, y, and c: x = [grp1;grp2;grp3]; y = [trial1;trial2;trial3]; c = [clr1;clr2;clr3]; % Convert those matrices to vectors: x = x(:); y = y(:); c = c(:); % Multiply x by 2 so that they're spread out: x = x*2; % Make the boxchart, boxchart(x(:),y(:),'GroupByColor',c(:)) % Set the x ticks and labels, and add a legend xticks(2:2:14); xticklabels(1:7) xlabel('Category') legend(["Trial 1" "Trial 2" "Trial 3"],'Location','NorthOutside')