Skip to main content

plotting - How to plot paired smooth histogram/distribution plots?



I've been trying to get paired distribution (aka "violin") plots like those shown below for a few hours, but all my attempts have failed.


enter image description here


The key features here are



  1. paired smooth histograms/distribution plots with a common vertical abscissa and divergent horizontal ordinates;

  2. different fill colors for the two subplots;


The closest Mathematica has are DistributionCharts, but these plots are not paired (i.e. they are always symmetrical).




I first tried SmoothHistogram, but it appears that there's no simple way to get a SmoothHistogram with a vertical abscissa.



Next I tried PairedSmoothHistogram, but I can't manage to assign different fill colors to the two sides.


BlockRandom[SeedRandom[0];
Quiet[
PairedSmoothHistogram[
RandomVariate[NormalDistribution[], 100]
, RandomVariate[NormalDistribution[], 100]
, AspectRatio -> Automatic
, Axes -> {False, True}
, Ticks -> None
, Spacings -> 0

, ImageSize -> 30
, Filling -> Axis
]
, OptionValue::nodef
]
]


Mathematica graphics




Then I tried a combination of SmoothKernelDistribution and either ListPlot, ParametricPlot or ContourPlot, but this won't work because neither ParametricPlot nor ContourPlot accepts a Fill option, and I can't figure out how to get ListPlot to fill the spaces between the curves and the vertical axis.


For example,


violin[data1_, data2_, rest___] := Module[
{ d1 = SmoothKernelDistribution[data1]
, d2 = SmoothKernelDistribution[data2]
, x
, xrange
}, xrange = {x, Sequence @@ (#[data1, data2] & /@ {Min, Max})}
; ParametricPlot[ {{-PDF[d1, x], x}, {PDF[d2, x], x}}
, Evaluate@xrange

, rest
, PlotRange -> All
, Axes -> {None, True}
, Ticks -> None
, PlotRangePadding -> {Automatic, None}
]
];

BlockRandom[SeedRandom[0];
violin[RandomVariate[NormalDistribution[], 100],

RandomVariate[NormalDistribution[], 100],
ImageSize -> 30
]
]


Mathematica graphics





I didn't expect this would be so hard. My brain is now fried.



Does anyone know how one does this?



Answer



Here is something using a custom ChartElementFunction


Module[{c = 0},
half[{{xmin_, xmax_}, {ymin_, ymax_}}, data_, metadata_] := (c++;
Map[Reverse[({0, Mean[{xmin, xmax}]} + # {1, (-1)^c})] &,
First@Cases[
First@Cases[InputForm[SmoothHistogram[data, Filling -> Axis]],
gc_GraphicsComplex :> Normal[gc], ∞],
p_Polygon, ∞], {2}])]


(thanks to @halirutan for reminding me about how to do closures in WL).


data = RandomVariate[NormalDistribution[0, 1], {4, 2, 100}];

DistributionChart[data, BarSpacing -> -1, ChartElementFunction -> half]

Mathematica graphics


Comments

Popular posts from this blog

front end - keyboard shortcut to invoke Insert new matrix

I frequently need to type in some matrices, and the menu command Insert > Table/Matrix > New... allows matrices with lines drawn between columns and rows, which is very helpful. I would like to make a keyboard shortcut for it, but cannot find the relevant frontend token command (4209405) for it. Since the FullForm[] and InputForm[] of matrices with lines drawn between rows and columns is the same as those without lines, it's hard to do this via 3rd party system-wide text expanders (e.g. autohotkey or atext on mac). How does one assign a keyboard shortcut for the menu item Insert > Table/Matrix > New... , preferably using only mathematica? Thanks! Answer In the MenuSetup.tr (for linux located in the $InstallationDirectory/SystemFiles/FrontEnd/TextResources/X/ directory), I changed the line MenuItem["&New...", "CreateGridBoxDialog"] to read MenuItem["&New...", "CreateGridBoxDialog", MenuKey["m", Modifiers-...

How to thread a list

I have data in format data = {{a1, a2}, {b1, b2}, {c1, c2}, {d1, d2}} Tableform: I want to thread it to : tdata = {{{a1, b1}, {a2, b2}}, {{a1, c1}, {a2, c2}}, {{a1, d1}, {a2, d2}}} Tableform: And I would like to do better then pseudofunction[n_] := Transpose[{data2[[1]], data2[[n]]}]; SetAttributes[pseudofunction, Listable]; Range[2, 4] // pseudofunction Here is my benchmark data, where data3 is normal sample of real data. data3 = Drop[ExcelWorkBook[[Column1 ;; Column4]], None, 1]; data2 = {a #, b #, c #, d #} & /@ Range[1, 10^5]; data = RandomReal[{0, 1}, {10^6, 4}]; Here is my benchmark code kptnw[list_] := Transpose[{Table[First@#, {Length@# - 1}], Rest@#}, {3, 1, 2}] &@list kptnw2[list_] := Transpose[{ConstantArray[First@#, Length@# - 1], Rest@#}, {3, 1, 2}] &@list OleksandrR[list_] := Flatten[Outer[List, List@First[list], Rest[list], 1], {{2}, {1, 4}}] paradox2[list_] := Partition[Riffle[list[[1]], #], 2] & /@ Drop[list, 1] RM[list_] := FoldList[Transpose[{First@li...

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...