def METHOD_NAME(statistics, output_file):
    figsize=(4,10)
    plt.rc('text', usetex=True)
    plt.rcParams['text.latex.preamble']='\\usepackage{relsize} \\usepackage{amsmath}'
    plt.rc('font', family='sanserif', size="16")
    plt.figure(figsize=figsize)
    for i_error in range(3):
        ax = plt.subplot(3,1,i_error+1)
        ax.set_ylabel(ylabels[i_error])
        plt.xlim(2e-3,1.1)
        for model,label,marker in zip(models,labels,markers):
            ax.loglog(statistics[model]["Time (seconds)"],statistics[model][errors[i_error]], label=label)
        if (i_error == 2):
            ax.set_ybound(lower=1e-5,upper=None)
    plt.xlabel("Time")
    ax.legend(bbox_to_anchor=(1.05, 1), loc='upper left', borderaxespad=0.)
    plt.savefig(output_file, bbox_inches='tight',dpi=200)
    return None