c1.jl

using Gadfly
using Distributions

function plot_theme()
  Theme(panel_fill = colorant"#3B444B", default_color = colorant"#DE3163",
    grid_color = colorant"#848482")
end

x = rand(Dirichlet([10; 1; 1]), 3000)

# View the samples in yz plane
myplot = plot(x = x1[2, :], y = x[3, :], Guide.xlabel("y"), Guide.ylabel("z"), plot_theme)

draw(SVG("myplot.svg", 6inch, 4inch), myplot)