def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list, args.select)
if not args.static:
logging.debug('Remove zero velocity samples')
data = ntp.remove_non_positive_velocity_samples(data)
if args.ue == 'e398':
# Rename MAC downlink throughput in Application downlink throughput if need be
ntp.process_data(data, ntp.process_lte_rename_mac_to_app)
# Get basic data
ntp.process_data(data, ntp.process_lte_app_throughput)
ntp.process_data(data, ntp.process_lte_pdcp_throughput)
column_list = [
'PDCP downlink throughput', 'Application throughput downlink'
]
if args.select is None:
df = tl.concat_pandas_data([df[column_list] for df in data])
else:
df = data
print(df['PDCP downlink throughput'].describe())
print(df['Application throughput downlink'].describe())
# Normalize
f_norm = lambda x: x / 1e6
plt.ion()
plt.figure()
plt.subplot2grid((2, 2), (0, 0), colspan=2)
x = np.arange(0, 120, 1)
dpl.plot_ecdf_pair(
df['Application throughput downlink'].dropna().apply(f_norm),
df['PDCP downlink throughput'].dropna().apply(f_norm), x,
'Application th.', 'PDCP th.', 'Mbit/s')
plt.legend(loc=0)
plt.subplot2grid((2, 2), (1, 0), colspan=1)
dpl.plot_density(
df['Application throughput downlink'].dropna().apply(f_norm), x,
'App. th.', 'Mbit/s')
plt.subplot2grid((2, 2), (1, 1), colspan=1)
dpl.plot_density(df['PDCP downlink throughput'].dropna().apply(f_norm), x,
'PDCP th.', 'Mbit/s')
if args.print:
plt.savefig(args.print, dpi=300, bbox_inches='tight')
input('Press any key')
def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list,args.select)
if args.ue == 'e398':
# Rename MAC downlink throughput in Application downlink throughput if need be
ntp.process_data(data,ntp.process_lte_rename_mac_to_app)
# Get basic data
ntp.process_data(data,ntp.process_lte_app_bw_prb_util)
ntp.process_data(data,ntp.process_lte_app_bw_prb_util_bw20)
ntp.process_data(data,ntp.process_lte_app_bw_prb_util_bw10)
ntp.process_data(data,ntp.process_lte_app_bw_prb_util_bw15)
# Spectral efficiency
ntp.process_data(data,ntp.process_se_bw_norm)
ntp.process_data(data,ntp.process_lte_se_rb)
column_list = ['DL bandwidth','DL bandwidth full',
'Application throughput downlink',
'PRB utilization DL','PRB utilization DL full',
'SE norm','SE RB norm']
if args.select is None:
df = dl.concat_pandas_data([df[column_list] for df in data ])
else:
df = data
print(df['SE RB norm'].describe())
print(df['SE norm'].describe())
print(df['Application throughput downlink'].describe())
# Normalize
f_norm = lambda x: x/1e6
plt.ion()
plt.figure()
plt.subplot2grid((2,1), (0,0))
dpl.plot_ts(df['Application throughput downlink'].dropna(),
'MAC DL throughput' if args.ue == 'e398' else 'App. throughput',
'Mbit/s',marker_size=10,ylim=None)
plt.subplot2grid((2,1), (1,0))
x = np.arange(0,155,1)
dpl.plot_ecdf(df['Application throughput downlink'].dropna().apply(f_norm),x,
'App. throughput','Mbit/s')
plt.figure()
dpl.plot_ecdf_pair(df['SE norm'].dropna(),
df['SE RB norm'].dropna(),np.linspace(0,8.5,86),
'SE norm','SE RB norm',
'bit/s/Hz')
input('Press any key')
def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list, args.select)
if not args.static:
logging.debug("Remove zero velocity samples")
data = ntp.remove_non_positive_velocity_samples(data)
if args.ue == "e398":
# Rename MAC downlink throughput in Application downlink throughput if need be
ntp.process_data(data, ntp.process_lte_rename_mac_to_app)
# Get basic data
ntp.process_data(data, ntp.process_lte_app_throughput)
ntp.process_data(data, ntp.process_lte_pdcp_throughput)
column_list = ["PDCP downlink throughput", "Application throughput downlink"]
if args.select is None:
df = tl.concat_pandas_data([df[column_list] for df in data])
else:
df = data
print(df["PDCP downlink throughput"].describe())
print(df["Application throughput downlink"].describe())
# Normalize
f_norm = lambda x: x / 1e6
plt.ion()
plt.figure()
plt.subplot2grid((2, 2), (0, 0), colspan=2)
x = np.arange(0, 120, 1)
dpl.plot_ecdf_pair(
df["Application throughput downlink"].dropna().apply(f_norm),
df["PDCP downlink throughput"].dropna().apply(f_norm),
x,
"Application th.",
"PDCP th.",
"Mbit/s",
)
plt.legend(loc=0)
plt.subplot2grid((2, 2), (1, 0), colspan=1)
dpl.plot_density(df["Application throughput downlink"].dropna().apply(f_norm), x, "App. th.", "Mbit/s")
plt.subplot2grid((2, 2), (1, 1), colspan=1)
dpl.plot_density(df["PDCP downlink throughput"].dropna().apply(f_norm), x, "PDCP th.", "Mbit/s")
if args.print:
plt.savefig(args.print, dpi=300, bbox_inches="tight")
input("Press any key")
def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list, args.select)
if not args.static:
logging.debug('Remove zero velocity samples')
data = ntp.remove_non_positive_velocity_samples(data)
if args.ue == 'e398':
# Rename MAC downlink throughput in Application downlink throughput if need be
ntp.process_data(data, ntp.process_lte_rename_mac_to_app)
# Get basic data
ntp.process_data(data, ntp.process_lte_app_bw_prb_util)
ntp.process_data(data, ntp.process_lte_app_bw_prb_util_bw20)
ntp.process_data(data, ntp.process_lte_app_bw_prb_util_bw10)
ntp.process_data(data, ntp.process_lte_app_bw_prb_util_bw15)
# Spectral efficiency
ntp.process_data(data, ntp.process_se_bw_norm)
column_list = [
'DL bandwidth', 'DL bandwidth full', 'PRB utilization DL',
'PRB utilization DL full', 'PRB utilization DL 10',
'PRB utilization DL 15', 'PRB utilization DL 20', 'SE', 'SE norm',
'SE 10 norm', 'SE 15 norm', 'SE 20 norm'
]
if args.select is None:
df = tl.concat_pandas_data([df[column_list] for df in data])
else:
df = data
print(df['SE'].describe())
print(df['SE norm'].describe())
# Remove outliers because of bandwidth normalization issues
df['SE norm'][df['SE norm'] > 7.5] = np.nan
print('90% percentile for SE norm: {:.1f} bit/s/Hz'.format(
np.percentile(df['SE norm'].dropna(), 90)))
print('95% percentile for SE norm: {:.1f} bit/s/Hz'.format(
np.percentile(df['SE norm'].dropna(), 95)))
print('99% percentile for SE norm: {:.1f} bit/s/Hz'.format(
np.percentile(df['SE norm'].dropna(), 99)))
print(df['SE 10 norm'].describe())
print(df['SE 15 norm'].describe())
print(df['SE 20 norm'].describe())
plt.ion()
plt.figure()
plt.subplot2grid((2, 2), (0, 0), colspan=2)
x = np.arange(0, 8.1, 0.1)
dpl.plot_ecdf_pair(df['SE norm'].dropna(), df['SE'].dropna(), x,
'Spectral efficiency (PRB norm.)',
'Spectral efficiency', 'bit/s/Hz')
plt.legend(loc=0)
plt.subplot2grid((2, 2), (1, 0), colspan=1)
dpl.plot_density(df['SE'].replace([np.inf, -np.inf], np.nan).dropna(), x,
'Spectral efficiency', 'bit/s/Hz')
plt.subplot2grid((2, 2), (1, 1), colspan=1)
dpl.plot_density(
df['SE norm'].where(df['SE norm'] < 8).replace([np.inf, -np.inf],
np.nan).dropna(), x,
'Spectral efficiency (PRB norm.)', 'bit/s/Hz')
if args.print:
plt.savefig(args.print[0], dpi=300, bbox_inches='tight')
plt.figure()
if len(df['SE 15 norm'].dropna()) > 0:
dpl.plot_ecdf_triplet(df['SE 10 norm'].dropna(),
df['SE norm'].dropna(),
df['SE 15 norm'].dropna(), x,
'Spectral efficiency 10 MHz (PRB norm.)\n',
'Spectral efficiency (PRB norm.)\n',
'Spectral efficiency 15 MHz (PRB norm.)',
'bit/s/Hz')
else:
dpl.plot_ecdf_pair(df['SE 10 norm'].dropna(), df['SE norm'].dropna(),
x, 'Spectral efficiency 10 MHz (PRB norm.)\n',
'Spectral efficiency (PRB norm.)\n', 'bit/s/Hz')
plt.legend(loc=0)
if args.print:
plt.savefig(args.print[1], dpi=300, bbox_inches='tight')
plt.figure()
if len(df['SE 15 norm'].dropna()) > 0:
dpl.plot_ecdf_quad(df['SE 10 norm'].dropna(), df['SE norm'].dropna(),
df['SE 15 norm'].dropna(), df['SE'].dropna(), x,
'SE 10 MHz (PRB norm.)', 'SE (PRB norm.)',
'SE 15 MHz (PRB norm.)', 'SE', 'bit/s/Hz')
plt.xlabel('[bit/s/Hz]')
else:
dpl.plot_ecdf_triplet(df['SE 10 norm'].dropna(),
df['SE norm'].dropna(), df['SE'].dropna(), x,
'Spectral efficiency 10 MHz (PRB norm.)\n',
'Spectral efficiency (PRB norm.)\n',
'Spectral efficiency', 'bit/s/Hz')
plt.legend(loc='best')
plt.tight_layout()
if args.print:
plt.savefig(args.print[2], dpi=300, bbox_inches='tight')
input('Press any key')
def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list,args.select)
if not args.static:
logging.debug('Remove zero velocity samples')
data = ntp.remove_non_positive_velocity_samples(data)
# Get basic data
ntp.process_data(data,ntp.process_lte_bw)
# Now check the RS-SNR difference between 10 and 15 MHz
ntp.process_data(data,ntp.process_lte_rs_snr)
ntp.process_data(data,ntp.process_lte_rs_snr_bw)
column_list = ['DL bandwidth','DL bandwidth full',
'RS SNR/Antenna port - 1','RS SNR/Antenna port - 1 10','RS SNR/Antenna port - 1 15',
'RS SNR/Antenna port - 2','RS SNR/Antenna port - 2 10','RS SNR/Antenna port - 2 15']
if args.select is None:
df = tl.concat_pandas_data([df[column_list] for df in data ])
else:
df = data
print(df['RS SNR/Antenna port - 1'].describe())
print(df['RS SNR/Antenna port - 2'].describe())
print('95% percentile for AP 1 and 2: {:.1f}, {:.1f} dB'.format(np.percentile(df['RS SNR/Antenna port - 1'].dropna(),95),np.percentile(df['RS SNR/Antenna port - 2'].dropna(),95)))
print('99% percentile for AP 1 and 2: {:.1f}, {:.1f} dB'.format(np.percentile(df['RS SNR/Antenna port - 1'].dropna(),99),np.percentile(df['RS SNR/Antenna port - 2'].dropna(),99)))
plt.ion()
plt.figure()
plt.subplot2grid((2,1), (0,0))
x = np.linspace(-15,37,53)
if len(df['RS SNR/Antenna port - 1 15'].dropna()) > 0:
dpl.plot_ecdf_triplet(df['RS SNR/Antenna port - 1 10'].dropna(),
df['RS SNR/Antenna port - 1'].dropna(),
df['RS SNR/Antenna port - 1 15'].dropna(),x,
'RS-SNR AP1 10 MHz','RS-SNR AP1','RS-SNR AP1 15 MHz','dB')
else:
dpl.plot_ecdf_pair(df['RS SNR/Antenna port - 1 10'].dropna(),
df['RS SNR/Antenna port - 1'].dropna(),x,
'RS-SNR AP1 10 MHz','RS-SNR AP1','dB')
plt.xlim([-12,33])
if args.blind:
plt.xticks([])
plt.subplot2grid((2,1), (1,0))
if len(df['RS SNR/Antenna port - 2 15'].dropna()) > 0:
dpl.plot_ecdf_triplet(df['RS SNR/Antenna port - 2 10'].dropna(),
df['RS SNR/Antenna port - 2'].dropna(),
df['RS SNR/Antenna port - 2 15'].dropna(),x,
'RS-SNR AP2 10 MHz','RS-SNR AP2','RS-SNR AP2 15 MHz','dB')
else:
dpl.plot_ecdf_pair(df['RS SNR/Antenna port - 2 10'].dropna(),
df['RS SNR/Antenna port - 2'].dropna(),x,
'RS-SNR AP2 10 MHz','RS-SNR AP2','dB')
plt.xlim([-12,33])
if args.blind:
plt.xticks([])
if args.print:
plt.savefig(args.print,dpi=300,bbox_inches='tight')
input('Press any key')
def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list,args.select)
if not args.static:
logging.debug('Remove zero velocity samples')
data = ntp.remove_non_positive_velocity_samples(data)
if args.ue == 'e398':
# Rename MAC downlink throughput in Application downlink throughput if need be
ntp.process_data(data,ntp.process_lte_rename_mac_to_app)
# Get basic data
ntp.process_data(data,ntp.process_lte_app_bw_prb_util)
ntp.process_data(data,ntp.process_lte_app_bw_prb_util_bw20)
ntp.process_data(data,ntp.process_lte_app_bw_prb_util_bw10)
ntp.process_data(data,ntp.process_lte_app_bw_prb_util_bw15)
# Spectral efficiency
ntp.process_data(data,ntp.process_se_bw_norm)
column_list = ['DL bandwidth','DL bandwidth full',
'PRB utilization DL','PRB utilization DL full',
'PRB utilization DL 10','PRB utilization DL 15','PRB utilization DL 20',
'SE','SE norm','SE 10 norm','SE 15 norm','SE 20 norm']
if args.select is None:
df = tl.concat_pandas_data([df[column_list] for df in data ])
else:
df = data
print(df['SE'].describe())
print(df['SE norm'].describe())
# Remove outliers because of bandwidth normalization issues
df['SE norm'][df['SE norm'] > 7.5] = np.nan
print('90% percentile for SE norm: {:.1f} bit/s/Hz'.format(np.percentile(df['SE norm'].dropna(),90)))
print('95% percentile for SE norm: {:.1f} bit/s/Hz'.format(np.percentile(df['SE norm'].dropna(),95)))
print('99% percentile for SE norm: {:.1f} bit/s/Hz'.format(np.percentile(df['SE norm'].dropna(),99)))
print(df['SE 10 norm'].describe())
print(df['SE 15 norm'].describe())
print(df['SE 20 norm'].describe())
plt.ion()
plt.figure()
plt.subplot2grid((2,2), (0,0),colspan=2)
x = np.arange(0,8.1,0.1)
dpl.plot_ecdf_pair(df['SE norm'].dropna(),
df['SE'].dropna(),x,
'Spectral efficiency (PRB norm.)',
'Spectral efficiency',
'bit/s/Hz')
plt.legend(loc=0)
plt.subplot2grid((2,2), (1,0),colspan=1)
dpl.plot_density(df['SE'].replace([np.inf, -np.inf], np.nan).dropna(),x,
'Spectral efficiency','bit/s/Hz')
plt.subplot2grid((2,2), (1,1),colspan=1)
dpl.plot_density(df['SE norm'].where(df['SE norm'] < 8).replace([np.inf, -np.inf], np.nan).dropna(),x,
'Spectral efficiency (PRB norm.)','bit/s/Hz')
if args.print:
plt.savefig(args.print[0],dpi=300,bbox_inches='tight')
plt.figure()
if len(df['SE 15 norm'].dropna()) > 0:
dpl.plot_ecdf_triplet(df['SE 10 norm'].dropna(),
df['SE norm'].dropna(),
df['SE 15 norm'].dropna(),x,
'Spectral efficiency 10 MHz (PRB norm.)\n',
'Spectral efficiency (PRB norm.)\n',
'Spectral efficiency 15 MHz (PRB norm.)',
'bit/s/Hz')
else:
dpl.plot_ecdf_pair(df['SE 10 norm'].dropna(),
df['SE norm'].dropna(),x,
'Spectral efficiency 10 MHz (PRB norm.)\n',
'Spectral efficiency (PRB norm.)\n',
'bit/s/Hz')
plt.legend(loc=0)
if args.print:
plt.savefig(args.print[1],dpi=300,bbox_inches='tight')
plt.figure()
if len(df['SE 15 norm'].dropna()) > 0:
dpl.plot_ecdf_quad(df['SE 10 norm'].dropna(),
df['SE norm'].dropna(),
df['SE 15 norm'].dropna(),
df['SE'].dropna(),x,
'SE 10 MHz (PRB norm.)',
'SE (PRB norm.)',
'SE 15 MHz (PRB norm.)','SE',
'bit/s/Hz')
plt.xlabel('[bit/s/Hz]')
else:
dpl.plot_ecdf_triplet(df['SE 10 norm'].dropna(),
df['SE norm'].dropna(),df['SE'].dropna(),x,
'Spectral efficiency 10 MHz (PRB norm.)\n',
'Spectral efficiency (PRB norm.)\n','Spectral efficiency',
'bit/s/Hz')
plt.legend (loc='best')
plt.tight_layout ()
if args.print:
plt.savefig(args.print[2],dpi=300,bbox_inches='tight')
input('Press any key')
def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list, args.select)
if not args.static:
logging.debug("Remove zero velocity samples")
data = ntp.remove_non_positive_velocity_samples(data)
if args.ue == "e398":
# Rename MAC downlink throughput in Application downlink throughput if need be
ntp.process_data(data, ntp.process_lte_rename_mac_to_app)
# Get basic data
ntp.process_data(data, ntp.process_lte_app_bw_prb_util)
ntp.process_data(data, ntp.process_lte_app_bw_prb_util_bw20)
ntp.process_data(data, ntp.process_lte_app_bw_prb_util_bw10)
ntp.process_data(data, ntp.process_lte_app_bw_prb_util_bw15)
# Spectral efficiency
ntp.process_data(data, ntp.process_se_bw_norm)
# Alternate version
ntp.process_data(data, ntp.process_lte_se_rb)
column_list = [
"DL bandwidth",
"DL bandwidth full",
"PRB utilization DL",
"PRB utilization DL full",
"PRB utilization DL 10",
"PRB utilization DL 15",
"PRB utilization DL 20",
"SE",
"SE norm",
"SE 10 norm",
"SE 15 norm",
"SE 20 norm",
"SE RB",
"SE RB norm",
]
if args.select is None:
df = tl.concat_pandas_data([df[column_list] for df in data])
else:
df = data
print(df["SE"].describe())
print(df["SE norm"].describe())
print(df["SE 10 norm"].describe())
print(df["SE 15 norm"].describe())
print(df["SE 20 norm"].describe())
print(df["SE RB"].describe())
print(df["SE RB norm"].describe())
plt.ion()
plt.figure()
plt.subplot2grid((2, 2), (0, 0), colspan=2)
x = np.arange(0, 8.1, 0.1)
dpl.plot_ecdf_pair(
df["SE norm"].dropna(),
df["SE"].dropna(),
x,
"Spectral efficiency (PRB norm.)",
"Spectral efficiency",
"bit/s/Hz",
)
plt.legend(loc=0)
plt.subplot2grid((2, 2), (1, 0), colspan=1)
dpl.plot_density(df["SE"].replace([np.inf, -np.inf], np.nan).dropna(), x, "Spectral efficiency", "bit/s/Hz")
plt.subplot2grid((2, 2), (1, 1), colspan=1)
dpl.plot_density(
df["SE norm"].where(df["SE norm"] < 8).replace([np.inf, -np.inf], np.nan).dropna(),
x,
"Spectral efficiency (PRB norm.)",
"bit/s/Hz",
)
if args.print:
plt.savefig(args.print[0], dpi=300, bbox_inches="tight")
plt.figure()
if len(df["SE 15 norm"].dropna()) > 0:
dpl.plot_ecdf_triplet(
df["SE 10 norm"].dropna(),
df["SE norm"].dropna(),
df["SE 15 norm"].dropna(),
x,
"Spectral efficiency 10 MHz (PRB util. norm.)\n",
"Spectral efficiency (PRB util. norm.)\n",
"Spectral efficiency 15 MHz (PRB util. norm.)",
"bit/s/Hz",
)
else:
dpl.plot_ecdf_pair(
df["SE 10 norm"].dropna(),
df["SE norm"].dropna(),
x,
"Spectral efficiency 10 MHz (PRB util. norm.)\n",
"Spectral efficiency (PRB util. norm.)\n",
"bit/s/Hz",
)
plt.legend(loc=0)
if args.print:
plt.savefig(args.print[1], dpi=300, bbox_inches="tight")
plt.figure()
plt.subplot2grid((2, 2), (0, 0), colspan=2)
dpl.plot_ecdf_pair(
df["SE norm"].dropna(),
df["SE RB norm"].dropna(),
x,
"Spectral efficiency (PRB util. norm.)\n",
"Spectral efficiency (RB norm.)",
"bit/s/Hz",
)
plt.legend(loc=0)
plt.subplot2grid((2, 2), (1, 0), colspan=1)
dpl.plot_density(df["SE"].replace([np.inf, -np.inf], np.nan).dropna(), x, "Spectral efficiency", "bit/s/Hz")
plt.subplot2grid((2, 2), (1, 1), colspan=1)
dpl.plot_density(df["SE RB norm"].dropna(), x, "Spectral efficiency (RB norm.)", "bit/s/Hz")
input("Press any key")
def main(args):
data_file_list = tl.get_data_file_list(args.library)
if args.list:
tl.print_list(data_file_list)
sys.exit(0)
data = tl.load_data_file(data_file_list, args.select)
if not args.static:
logging.debug('Remove zero velocity samples')
data = ntp.remove_non_positive_velocity_samples(data)
# Get basic data
ntp.process_data(data, ntp.process_lte_bw)
# Now check the RS-SNR difference between 10 and 15 MHz
ntp.process_data(data, ntp.process_lte_rs_snr)
ntp.process_data(data, ntp.process_lte_rs_snr_bw)
column_list = [
'DL bandwidth', 'DL bandwidth full', 'RS SNR/Antenna port - 1',
'RS SNR/Antenna port - 1 10', 'RS SNR/Antenna port - 1 15',
'RS SNR/Antenna port - 2', 'RS SNR/Antenna port - 2 10',
'RS SNR/Antenna port - 2 15'
]
if args.select is None:
df = tl.concat_pandas_data([df[column_list] for df in data])
else:
df = data
print(df['RS SNR/Antenna port - 1'].describe())
print(df['RS SNR/Antenna port - 2'].describe())
print('95% percentile for AP 1 and 2: {:.1f}, {:.1f} dB'.format(
np.percentile(df['RS SNR/Antenna port - 1'].dropna(), 95),
np.percentile(df['RS SNR/Antenna port - 2'].dropna(), 95)))
print('99% percentile for AP 1 and 2: {:.1f}, {:.1f} dB'.format(
np.percentile(df['RS SNR/Antenna port - 1'].dropna(), 99),
np.percentile(df['RS SNR/Antenna port - 2'].dropna(), 99)))
plt.ion()
plt.figure()
plt.subplot2grid((2, 1), (0, 0))
x = np.linspace(-15, 37, 53)
if len(df['RS SNR/Antenna port - 1 15'].dropna()) > 0:
dpl.plot_ecdf_triplet(df['RS SNR/Antenna port - 1 10'].dropna(),
df['RS SNR/Antenna port - 1'].dropna(),
df['RS SNR/Antenna port - 1 15'].dropna(), x,
'RS-SNR AP1 10 MHz', 'RS-SNR AP1',
'RS-SNR AP1 15 MHz', 'dB')
else:
dpl.plot_ecdf_pair(df['RS SNR/Antenna port - 1 10'].dropna(),
df['RS SNR/Antenna port - 1'].dropna(), x,
'RS-SNR AP1 10 MHz', 'RS-SNR AP1', 'dB')
plt.xlim([-12, 33])
if args.blind:
plt.xticks([])
plt.subplot2grid((2, 1), (1, 0))
if len(df['RS SNR/Antenna port - 2 15'].dropna()) > 0:
dpl.plot_ecdf_triplet(df['RS SNR/Antenna port - 2 10'].dropna(),
df['RS SNR/Antenna port - 2'].dropna(),
df['RS SNR/Antenna port - 2 15'].dropna(), x,
'RS-SNR AP2 10 MHz', 'RS-SNR AP2',
'RS-SNR AP2 15 MHz', 'dB')
else:
dpl.plot_ecdf_pair(df['RS SNR/Antenna port - 2 10'].dropna(),
df['RS SNR/Antenna port - 2'].dropna(), x,
'RS-SNR AP2 10 MHz', 'RS-SNR AP2', 'dB')
plt.xlim([-12, 33])
if args.blind:
plt.xticks([])
if args.print:
plt.savefig(args.print, dpi=300, bbox_inches='tight')
input('Press any key')
