#exampleArgs="--fi_expr ../exprData/impute-t10_libNorm10K_cells.foreskin_genesGeq1Pct.all.kcyte.pkl \
#--fi_stageIDs ../clusterCells/kasp.ka10.k30_impute.t10_gene.geq5UMIgeq100Cell.all.csv \
#--cells ./cellNames.passFilter.tmp \
#--superstages progenitor:stage1,stage2,stage3;\
#progenitor.stage4:stage1,stage2,stage3,stage4;\
#differentiated:stage5,stage6,stage7;\
#differentiated.stage4:stage4,stage5,stage6,stage7;\
#all.noStage8:stage1,stage2,stage3,stage4,stage5,stage6,stage7 \
#--corrMethod logTpm \
#--fo corr_tmpA"
#args = parser.parse_args(exampleArgs.split())
if args.superstages is not None:
superstage_dict = IOutils.parseStrToDict(args.superstages,
valueType="str",
pairSep=";")
try:
superstage_dict = OrderedDict([(k, [int(x) for x in v.split(",")])
for k, v in superstage_dict.items()])
except:
superstage_dict = OrderedDict([
(k, [int(re.search(r'(\d+)$', x).group(1)) for x in v.split(",")])
for k, v in superstage_dict.items()
])
############################################################################################################
## Load data
cellData = IOutils.loadCellData(
OrderedDict([("expr", args.fi_expr), ("pcComps", args.fi_stageIDs)]))
cellData = pd.concat([
