analysis.plot(
        model=mdl,
        time=time,
        yDatas={'RNA': totalRna},
        fileName=os.path.join(output_directory, 'Total-RNA.png')