Hi Rita,
Here is one example implementation which leverages the boxplot function to achieve the desired layout using skin response data:
skinResponse = randn(100, 1); % Replace with your actual skin response data
% Grouping variables
chooserCategory = [repmat({'SS'}, 50, 1); repmat({'ES'}, 50, 1)];
songType = [repmat({'Arousing'}, 25, 1); repmat({'Relaxing'}, 25, 1); ...
repmat({'Arousing'}, 25, 1); repmat({'Relaxing'}, 25, 1)];
% Combine grouping variables
group = strcat(chooserCategory, '-', songType);
figure;
boxplot(skinResponse, group, 'Labels', {'SS-Arousing', 'SS-Relaxing', 'ES-Arousing', 'ES-Relaxing'});
xlabel('Category');
ylabel('Skin Response');
Hope it helps.
