def trace_nuts(csv, extras='', skip=0, n_col=4):
    """Trace plots of NUTS state, along with other model parameters.
    """
    from pylab import subplot, plot, gca, title, grid, xticks
    if isinstance(csv, dict):
        csv = [csv]
    if isinstance(extras, str):
        extras = extras.split()
    n_nuts_params = 7
    n_subplots = len(extras) + n_nuts_params
    n_row = n_subplots // n_col + 1
    for csvi in csv:
        i = 1
        for key in csvi.keys():
            if key.endswith('__') or key in extras:
                subplot(n_row, n_col, i)
                plot(csvi[key][skip:], alpha=0.5)
                if key in ('stepsize__',):
                    gca().set_yscale('log')
                title(key)
                grid(1)
                if i < n_subplots - n_col:
                    xticks(xticks()[0], [])
                i += 1
