function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
    % create tiledlayout (R2023a and newer)
    tl = tiledlayout(panel,"vertical");

    % plot 1: speed
    ax1 = nexttile(tl);
    hold(ax1, "on")
    grid(ax1, "on")
    title(ax1, "speed [km/h]")
    for i = 1:length(selected_laps)
        plot(ax1,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).speed_kph)
    end
    % plot 2: accelerator pedal position [%]
    ax2 = nexttile(tl);
    hold(ax2, "on")
    grid(ax2, "on")
    title(ax2, "APP [%]")
    for i = 1:length(selected_laps)
        plot(ax2,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).app_percent)
    end
    linkaxes([ax1, ax2],"x")
    % plot 3: accelerator pedal position histogram
    ax3 = nexttile(tl);
    hold(ax3, "on")
    grid(ax3, "on")
    title(ax3, "APP distribution")
    ylabel(ax3, "\sigma [%]")
    xlabel(ax3, "APP [%]")
    for i = 1:length(selected_laps)
        histogram(ax3,pltData(selected_laps(i)).app_percent,"Normalization","percentage")
    end
    % plot 4: brake pressure histogram
    ax4 = nexttile(tl);
    hold(ax4, "on")
    grid(ax4, "on")
    title(ax4, "BrakeP distribution")
    ylabel(ax4, "\sigma [%]")
    xlabel(ax4, "Brake Pressure Front [bar]")
    for i = 1:length(selected_laps)
        histogram(ax4,pltData(selected_laps(i)).brakePFront_bar,"Normalization","percentage")
    end

    % return null (not relevant for plots!)
    outputArg = [];
end