plotsolutions recipe.jl

@userplot PlotSolutions
@recipe function f(o::PlotSolutions)
  samples = o.args[1]

  species       --> measuredinds
  t             --> 0:1/3:30
  color_palette --> [colorant"steelblue"]
  linewidth     --> 0.1
  label         --> ""
  seriesalpha   --> 0.1
  layout  --> length(d[:species])

  nspecies = length(d[:species])
  nsamples = size(samples, 1)

  subplot :=  repeat(collect(1:nspecies)', outer=[1, nsamples])

  solutions = hcat([gync(samples[i,:] |> vec, d[:t])[:,d[:species]] for i in 1:nsamples]...)

  delete!(d, :species)
  #delete!(d, :t)

  d[:t], solutions
end