Hi,

I am using Stata16 and would like to produce a multicoloured coefplot from a matrix that contains means, p25 and p75 from a tabstat command:
Code:
tabstat dom1-dom9, by(cluster_num) stat(mean p25 p75) save
matrix c1 = r(Stat1)
My coefplot code is:
Code:
    coefplot matrix(c1[1]), ///
        ci((2 3)) ///
        msymbol(x) msize(large) ///
        ylabel(, labsize(small)) xlabel(, labsize(small)) xlab(0(20)100) xsca(range(0(20)100)) ///
        graphregion(color(white) fcolor(white) lcolor(white) ilcolor(white) ifcolor(white)) ///
        title("Cluster 1", ///
        size(small)) name("scores_cluster_1", replace)
I would like to plot each interquartile range in a different colour. Unfortunately, I am using a martix for the coefplot, so I wonder if there is a workaround for this?

Thanks,
Jane.