fn=websave('MATLAB_exampledata.xlsx','https://www.mathworks.com/matlabcentral/answers/uploaded_files/1136975/MATLAB_exampledata.xlsx');
iscat=vartype('cellstr');
tData=convertvars(tData,@iscellstr,'categorical');
tData.Mean=mean(tData{:,vartype('double')},2)
tData =
Group Range percentage percentage_1 percentage_2 Mean
_______ _____ __________ ____________ ____________ ______
Control <30 98.276 100 93.103 97.126
Control 30-39 1.7241 0 6.8966 2.8736
Control 40< 0 0 0 0
Drug 1 <30 89.394 87.692 90 89.029
Drug 1 30-39 10.606 12.308 10 10.971
Drug 1 40< 0 0 0 0
Drug 2 <30 61.765 75.61 100 79.125
Drug 2 30-39 36.765 24.39 0 20.385
Drug 2 40< 1.4706 0 0 0.4902
Drug 3 <30 33.962 1.6129 88.043 41.206
Drug 3 30-39 52.83 43.548 11.957 36.112
Drug 3 40< 13.208 54.839 0 22.682
R=categories(tData.Range);
if i==1, posn=hT(i).Position; L=posn(1);B=posn(2);H=posn(4);end
hSC=scatter(tData(ix,:), 'Group', {'percentage', 'percentage_1', 'percentage_2'},'filled');
hL=plot(double(tData.Group(ix)).'+0.5*[-1;1],double(tData.Mean(ix)).'+[-0;0], 'r-');
hT(i).XAxisLocation='origin';
if i>1, hT(i).YAxis.Visible='off';end
hLg=legend([hSC(:);hL(1)],'Pct 0','Pct 1','Pct 2','Mean', ...
'Location','eastoutside','Orientation','vertical');
posn=hT(i).Position; W=posn(1)+posn(3)-L;
hA=annotation('rectangle',[L B W H]);
Would be my starting point with a traditional scatter can't run the cigarScatter on online platform; I'll see if I can find some time to look at it some later...that's most of what you wanted on the axes; it isn't so easy to draw the horizontal line with categorical variables; yline works but goes across the whole axis range so that doesn't help and there's only the one point on a categorical axes at which can plot. One way is to introduce "fake" categories into the variable to introduce more granularity in the axis. Another is to lay another axis on top the the existing that is numeric and plot into it; alternatively, just do as the cigar plot guy did; do it all numerically and the just label the ticks; the tick labels and the ticks are independent properties.
Mayhaps that will give some ideas on approaches...
ADDENDUM: "Stealing" from @Chunru using the table variable names and the new feature (introduced R2022b) of being able to use other than a categorical variable in a categoricalRuler axis, just a few enhancements. Also to add the overall figure box is not too bad with the annotation function; just have to compute the outer boundaries of the axes from first to last to find the overall coordintes.
To add the dividing lines between is another application of annotation with a line; one needs to add the computation of the x-location between the first two and the last two axes to come up with the position at which to draw it; I've left that as "exercise for student" -- it would follow the idea illustrated above in the computation of the overall width -- one has to find the RH position of the first axis by adding the Left and Width position values and then add halve the distance between that location and the Left postion of the second.
With the advent of being able to use numeric data on categorical axis, all it would take would be to add a little jitter to the x position vector to be able to create the "cigarscater" figure without having to fool with it. I looked at it a little last night but it had so much extra baggage coming along with it to handle all the "what ifs" that it does that I didn't feel like trying to delve into just how it actually worked to try to do anything more with it.
