Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Triage October cell processing failures #465

Open
4 tasks
wbwakeman opened this issue Aug 5, 2020 · 0 comments
Open
4 tasks

Triage October cell processing failures #465

wbwakeman opened this issue Aug 5, 2020 · 0 comments

Comments

@wbwakeman
Copy link

wbwakeman commented Aug 5, 2020

As an electrophysiology manager and author on a Patch-Seq paper, I want to understand why potential cells for October 2020 release are not progressing through the processing pipeline, so that I can make the decision to fail them or take corrective action.

4674 potential October 2020 cells have been run through the pipeline.

Four failure modes are listed below for three processing pipelines:
'EPHYS_NWB_STIMULUS_SUMMARY_V3_QUEUE' (4 + 1 failures )
'EPHYS_QC_V3_QUEUE' ( 60 failures )
'EPHYS_FEATURE_EXTRACTION_V3_QUEUE' ( 284 + 52 + 9 + 5 + 22 + 1 = 385 failures )

Acceptance criteria:

  • A very brief analysis on each type of failure
  • Whether modifications to the code are a reasonable strategy to address the failure
  • Whether any action taken by data generators (e.g. patch-clamp team or MIES developer) could address the issue
  • Create new issue to address action that could be taken

For each failure, the following information is listed:
cell_specimen_id ; cell_specimen_name ; ephys_roi_result_id

EPHYS_NWB_STIMULUS_SUMMARY_V3_QUEUE

ValueError: Found multiple stimulus series {[s.name for s in matching_series]} for sweep number {sweep_number}
637823753;"Htr3a-Cre_NO152;Ai14-349020.03.02.01";637823738
637824915;"Htr3a-Cre_NO152;Ai14-349020.03.02.02";637824890
637833144;"Htr3a-Cre_NO152;Ai14-349020.03.02.03";637833126
637839160;"Htr3a-Cre_NO152;Ai14-349020.03.02.04";637838825
Exception: Could not construct VoltageClampSeries object due to: VoltageClampSeries.__init__: missing argument 'gain'
767875138;"Nos1-CreERT2;Ai14-421031.03.01.01";767874754

EPHYS_QC_V3_QUEUE

marshmallow.exceptions.ValidationError: {'sweep_features': {15: {'stimulus_duration': ['Field may not be null.']}}}

643572450;"H17.06.013.12.05.01";643572437
646845860;"Ntsr1-Cre_GN220;Ai14-357985.03.02.02";646845845
649801264;"Oxtr-T2A-Cre;Pvalb-T2A-FlpO;Ai65-360187.05.01.01";649801249
650481841;"Gad2-IRES-Cre;Ai14-362153.04.01.02";650481821
652960405;"Rbp4-Cre_KL100;Ai14-366313.05.02.01";652960390
656993561;"Vip-IRES-Cre;Ai14-366689.04.01.03";656993546
657185162;"Chrna2-Cre_OE25;Ai14-367172.03.01.02";657185147
673133977;"Slc32a1-IRES-Cre;Ai14-379844.03.02.03";673133958
675716304;"Htr3a-Cre_NO152;Ai14-381048.03.01.01";675715872
678666933;"Slc17a8-IRES2-Cre;Ai14-382382.03.01.02";678666873
679496928;"Slc17a7-IRES2-Cre;Ai14-382512.03.02.02";679496685
681804422;"Vip-IRES-Cre;Ai14-383085.05.01.01";681804177
694827589;"Vipr2-IRES2-Cre;Slc32a1-T2A-FlpO;Ai65-388176.02.01.01";694827458
699569823;"Nos1-CreERT2;Sst-IRES-FlpO;Ai65-391509.03.01.02";699569250
710060452;"Vip-IRES-Cre;Ai14-395420.04.02.01";710060314
714419844;"Vip-IRES-Cre;Ai14-398013.05.01.01";714419700
717355259;"Calb2-IRES-Cre;Sst-IRES-FlpO;Ai65-398664.05.02.01";717354808
737933862;"H18.06.363.11.15.02.01";737933847
740675302;"H18.03.315.12.03.02";740675226
760514840;"H18.06.365.11.09.01.05";760514825
766195495;"Htr3a-Cre_NO152;Pvalb-T2A-FlpO;Ai65-419991.06.01.01";766195138
797127309;"Htr3a-Cre_NO152;Ai14-431562.03.01.01";797127049
797969939;"H18.03.323.11.09.01.02";797969402
803400389;"Sst-IRES-Cre;Ai14-434645.05.01.01";803400329
807349990;"Vip-IRES-Cre;Ai14-435637.05.01.02";807349737
807429099;"Vip-IRES-Cre;Ai14-435637.05.01.04";807428786
814308166;"Crh-IRES-Cre_ZJH;Sst-IRES-FlpO;Ai65-439424.03.01.01";814307869
819191383;"h19.03.002.12.01.01.01.09";819191225
829302487;"Chrna6-IRES2-FlpO-neo;Ai65F-443493.04.02.01";829302463
832594616;"H19.03.305.11.10.01.06";832594509
832689550;"H19.03.305.11.08.01.09";832689531
835578804;"H19.03.306.11.03.01.02";835578575
835693681;"H19.03.306.11.15.02.02";835693658
839918826;"H19.06.351.11.16.01.03";839918580
840015014;"H19.06.351.11.10.01.05";840014977
844128819;"Vip-IRES-Cre;Ai14-453251.04.01.01";844128556
847678139;"Sst-IRES-Cre;Ai14-454977.05.01.01";847677977
848629140;"H19.06.352.99.02.03.01.02";848629111
848656903;"H19.06.352.99.02.03.01.03";848656880
848680984;"H19.06.352.99.02.03.01.05";848680960
852652659;"H19.03.307.99.01.05.01.06";852652585
862363963;"Crh-IRES-Cre_ZJH;Sst-IRES-FlpO;Ai65-459713.04.01.02";862363869
868067910;"H19.03.311.14.08.01.01";868067681
870853500;"Nkx2-1-CreERT2;Ai14-465338.03.02.02";870852958
875433862;"H19.06.355.99.01.03.01.03";875433500
883688278;"Slc17a7-IRES2-Cre;Ai14-467901.03.01.01";883688139
893455915;"H19.03.314.1A.01.01";893455729
893605869;"H19.03.314.1A.02.06";893605727
910434386;"H19.03.316.2A.01.06";910434305
914356072;"H19.03.316.2A.08.04";914355964
918932896;"Vip-IRES-Cre;Ai14-479200.04.02.01";918932569
923113598;"H19.03.319.11.04.01.02";923113568
926786832;"Gad2-IRES-Cre;Ai14-481679.04.01.01";926786473
945743049;"Rbp4-Cre_KL100;Ai14-486151.04.02.01";945742596
949987663;"Slc32a1-IRES-Cre;Ai14-486249.04.01.01";949987241
957296740;"Rbp4-Cre_KL100;Ai14-488825.05.01.01";957296477
982703343;"Ctgf-T2A-dgCre;Ai14-496562.03.01.01";982703055
986231454;"Penk-IRES2-Cre-neo;Slc17a6-IRES2-FlpO-497491-Ai65.05.01.01";986231379
991968323;"Rbp4-Cre_KL100;Ai14-499543.05.01.02";991968113
991976495;"Rbp4-Cre_KL100;Ai14-499543.05.01.03";991976427
1002212559;"H20.26.401.11.52.01.02";1002212371
1003951822;"Sst-IRES-Cre;Ai14-508744.06.06.01";1003951732
1005039684;"Penk-IRES2-Cre-neo;Slc17a6-IRES2-FlpO-508719-Ai65.10.09.03";1005039668
1005134890;"H20.03.303.11.03.01.09";1005134874
1005433317;"Sst-IRES-Cre;Ai14-508750.10.06.04";1005433206
1005471491;"Penk-IRES2-Cre-neo;Slc17a6-IRES2-FlpO-508935-Ai65.10.09.02";1005471409
1005680527;"H20.03.303.11.51.01.06";1005680446
1006380866;"Ctgf-T2A-dgCre;Ai14-509691.05.01.02";1006380781

EPHYS_FEATURE_EXTRACTION_V3_QUEUE

if r[0] <= 0: # This happens when stimulus starts less than 0.5 s after test pulse
TypeError: 'NoneType' object is not subscriptable

637824915;"Htr3a-Cre_NO152;Ai14-349020.03.02.02";637824890
637833144;"Htr3a-Cre_NO152;Ai14-349020.03.02.03";637833126
638246905;"Rorb-IRES2-Cre;Ai14-349007.06.01.02";638246857
639395516;"Gad2-IRES-Cre;Ai14-350671.04.02.01";639395433
639719157;"Chrna2-Cre_OE25;Ai14-351067.03.02.02";639719100
648547708;"Pvalb-IRES-Cre;Ai14-359332.03.01.02";648547556
651098469;"Ndnf-IRES2-dgCre;Ai14-361703.04.01.02";651098318
660049716;"Pdyn-T2A-CreERT2;Ai14-373524.05.02.01";660049701
668458490;"Vip-IRES-Cre;Ai14-377300.03.02.03";668458471
710899937;"Sst-IRES-Cre;Ai14-395724.03.02.01";710899857
711464976;"Pvalb-T2A-CreERT2;Ai14-396257.06.02.01";711464898
728523077;"Sst-IRES-Cre;Ai14-405113.05.01.03";728522964
754436948;"Sst-IRES-Cre;Ai167-415614.05.02.02";754436507
754447289;"Sst-IRES-Cre;Ai167-415614.02.01.01";754447182
754487624;"Sst-IRES-Cre;Ai167-415614.06.02.02";754487490
754549487;"Sst-IRES-Cre;Ai167-415614.06.01.02";754549472
754753909;"Sst-IRES-Cre;Ai167-415615.03.01.01";754753737
754759234;"Sst-IRES-Cre;Ai167-415615.03.02.01";754759153
754781874;"Sst-IRES-Cre;Ai167-415615.03.01.02";754781722
754804871;"Sst-IRES-Cre;Ai167-415615.04.01.02";754804720
754807838;"Sst-IRES-Cre;Ai167-415615.03.02.02";754807701
754873245;"Sst-IRES-Cre;Ai167-415615.02.01.01";754873216
754893893;"Sst-IRES-Cre;Ai167-415615.06.02.01";754893877
755571201;"Vip-IRES-Cre;Npy-IRES2-FlpO;Ai65-413072.04.01.02";755571001
755931935;"Sst-IRES-Cre;Ai167-415616.04.02.01";755931299
756412962;"Sst-IRES-Cre;Ai167-415617.03.01.03";756412893
756430849;"Sst-IRES-Cre;Ai167-415617.05.02.01";756430562
756456229;"Sst-IRES-Cre;Ai167-415617.03.02.01";756455876
756499370;"Sst-IRES-Cre;Ai167-415617.03.02.02";756498869
756500119;"Sst-IRES-Cre;Ai167-415617.06.01.02";756499777
756500874;"Sst-IRES-Cre;Ai167-415617.06.02.01";756500301
756549723;"Sst-IRES-Cre;Ai167-415617.03.02.03";756549045
757080816;"Sst-IRES-Cre;Ai167-415618.03.01.01";757080441
757120969;"Sst-IRES-Cre;Ai167-415618.03.01.02";757120688
757139742;"Sst-IRES-Cre;Ai167-415618.04.02.01";757139686
758002124;"Sst-IRES-Cre;Ai167-416171.03.01.02";758001945
759416712;"Sst-IRES-Cre;Ai167-416176.04.01.01";759416508
799929013;"Sst-IRES-Cre;Ai167-429815.04.01.03";799928896
799944803;"Sst-IRES-Cre;Ai167-429815.06.02.01";799944684
800320286;"Sst-IRES-Cre;Ai167-432052.04.02.01";800320262
807349990;"Vip-IRES-Cre;Ai14-435637.05.01.02";807349737
807429099;"Vip-IRES-Cre;Ai14-435637.05.01.04";807428786
809070689;"Pvalb-IRES-Cre;Ai14-436295.05.01.04";809070461
814308166;"Crh-IRES-Cre_ZJH;Sst-IRES-FlpO;Ai65-439424.03.01.01";814307869
817053274;"Sst-IRES-Cre;Ai167-439541.04.01.01";817053079
817060299;"Sst-IRES-Cre;Ai167-439541.05.02.02";817060065
818550823;"Sst-IRES-Cre;Ai167-441787.04.02.01";818550615
844128819;"Vip-IRES-Cre;Ai14-453251.04.01.01";844128556
855981127;"Sncg-IRES2-FlpO-neo;Ai65F-451117.04.01.02";855980964
870853500;"Nkx2-1-CreERT2;Ai14-465338.03.02.02";870852958
873726936;"Nkx2-1-CreERT2;Ai14-465343.05.01.01";873726656
899755423;"Vip-IRES-Cre;Ai14-472422.04.01.02";899755278
902212242;"Gad2-IRES-Cre;Ai14-473603.03.01.01";902211977
918932896;"Vip-IRES-Cre;Ai14-479200.04.02.01";918932569
921879850;"Sncg-IRES2-FlpO-neo;Ai65F-480219.03.02.01";921879446
922877472;"Sncg-IRES2-FlpO-neo;Ai65F-480221.03.02.03";922876968
924739827;"Sncg-IRES2-FlpO-neo;Ai65F-481159.03.01.01";924739798
935084855;"Sncg-IRES2-FlpO-neo;Ai65F-484137.04.02.02";935084659
935760710;"Oxtr-T2A-Cre;Ai14-484067.04.01.03";935760305
941981857;"Sncg-IRES2-FlpO-neo;Ai65F-485422.04.02.01";941981583
978133978;"Ctgf-T2A-dgCre;Ai14-495722.04.02.02";978133728
992272581;"Vipr2-IRES2-Cre;Slc32a1-IRES2-FlpO;Ai65-499514.03.01.01";992272566
992430607;"Sncg-IRES2-FlpO-neo;Ai65F-499190.03.01.01";992430573
1004779144;"Sst-IRES-Cre;Ai14-508746.06.06.01";1004779043

52 with this error: e.g. err.id 1004779043

File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 193, in _run_module_as_main
"__main__", mod_spec)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 85, in _run_code
exec(code, run_globals)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 109, in <module>
if __name__ == "__main__": main()
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 105, in main
module.args["cell_features"])
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 54, in run_feature_extraction
cell_state, feature_states) = dsft.extract_data_set_features(data_set)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/data_set_features.py", line 378, in extract_data_set_features
sweep_features[s['sweep_number']]['peak_deflect'] = s['peak_deflect']
KeyError: 4

9 with this error e.g. id = 1005657758

Traceback (most recent call last):
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 193, in _run_module_as_main
"__main__", mod_spec)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 85, in _run_code
exec(code, run_globals)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 109, in <module>
if __name__ == "__main__": main()
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 105, in main
module.args["cell_features"])
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 41, in run_feature_extraction
raise er.FeatureError("There are no QC-passed sweeps available to analyze")
ipfx.error.FeatureError: There are no QC-passed sweeps available to analyze

5 with this error e.g. ids 1005108832
1005032767
1003299623
1001443326
1001424563

Traceback (most recent call last):
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 193, in _run_module_as_main
"__main__", mod_spec)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 85, in _run_code
exec(code, run_globals)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 109, in <module>
if __name__ == "__main__": main()
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 105, in main
module.args["cell_features"])
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 82, in run_feature_extraction
plotqc.display_features(qc_fig_dir, data_set, feature_data)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 769, in display_features
make_sweep_page(data_set, qc_fig_dir)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 715, in make_sweep_page
sweep_files = plot_sweep_figures(data_set, image_dir, sizes)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 207, in plot_sweep_figures
v, i, t, r, dt = load_sweep(data_set, sweep_number)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 36, in load_sweep
dt = sweep.t[1] - sweep.t[0]
IndexError: index 1 is out of bounds for axis 0 with size 0

22 with this error
1000612453,999319545,920015861,926949081,912201092,870849339,874356674,880976423,874570247,873810658,875020807,880542017,87552,386,875486563,872282934,872248672,875139943,874444852,852974739,839931440,941908038,751935222

Traceback (most recent call last):
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 193, in _run_module_as_main
"__main__", mod_spec)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/runpy.py", line 85, in _run_code
exec(code, run_globals)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 109, in <module>
if __name__ == "__main__": main()
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 105, in main
module.args["cell_features"])
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 82, in run_feature_extraction
plotqc.display_features(qc_fig_dir, data_set, feature_data)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 769, in display_features
make_sweep_page(data_set, qc_fig_dir)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 715, in make_sweep_page
sweep_files = plot_sweep_figures(data_set, image_dir, sizes)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 174, in plot_sweep_figures
if r_init[0] <= 0:
TypeError: 'NoneType' object is not subscriptable

ValueError: Axis limits cannot be NaN or Inf
ephys_roi_result.id = 660049701

Failure [limsweb7.corp.alleninstitute.org]: Exit status is 1 [log file: /allen/programs/celltypes/production/mousecelltypes/prod58/Ephys_Roi_Result_660049701/202008031836_EPHYS_FEATURE_EXTRACTION_V3_QUEUE_660049701.log]: (last 25 lines)...
) File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 105, in main
module.args["cell_features"])
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/bin/run_feature_extraction.py", line 82, in run_feature_extraction
plotqc.display_features(qc_fig_dir, data_set, feature_data)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 770, in display_features
make_cell_page(data_set, feature_data, qc_fig_dir)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 728, in make_cell_page
cell_files = plot_cell_figures(data_set, feature_data, image_dir, sizes)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 620, in plot_cell_figures
plot_ramp_figures(data_set, cell_features, lims_features, sweep_features, image_dir, sizes, cell_image_files)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 455, in plot_ramp_figures
figs = plot_single_ap_values(data_set, ramps_sweeps, lims_features, sweep_features, cell_features, "ramp")
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/ipfx/plot_qc_figures.py", line 147, in plot_single_ap_values
plt.xlim(spikes[0]["threshold_t"] - 0.002, spikes[0]["fast_trough_t"] + 0.01)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/matplotlib/pyplot.py", line 1396, in xlim
ret = ax.set_xlim(*args, **kwargs)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/matplotlib/axes/_base.py", line 3185, in set_xlim
right = self._validate_converted_limits(right, self.convert_xunits)
File "/allen/aibs/technology/conda/production/fx/lib/python3.6/site-packages/matplotlib/axes/_base.py", line 3101, in _validate_converted_limits
raise ValueError("Axis limits cannot be NaN or Inf")
ValueError: Axis limits cannot be NaN or Inf 
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant