natfeedback_code_eLife/README.md

Data and code for Spacek et. al., 2022, "Robust effects of corticothalamic feedback and behavioral state on movie responses in mouse dLGN", submitted to eLife.

https://doi.org/10.1101/776237

Use main.py to load data for all LGN opto experiments as a function of stimulus type (movies, gratings), run state, and opto state. Data are loaded from pandas DataFrames, which are saved as Python pickles in the pickles folder.

Run main.py at the IPython command line, or in Jupyter; the -i runs scripts in IPython/Jupyter's namespace instead of an empty one:

run -i main.py

pvmvis (default) or ntsrmvis can be provided as an --exptype= argument to this main.py, to load and then plot PV-Cre (most of the paper) or Ntsr1-Cre (Fig S4) data, respectively.

Results can be plotted by executing the desired figx.py file in the same namespace:

run -i fig1.py
run -i fig2.py
run -i fig3.py
run -i fig4.py
run -i fig5.py
run -i fig6.py

Not all supplementary panels/figures have their own script, some are integrated into their parent figure's script for convenience:

run -i fig1S33S1.py
run -i fig1S4.py
run -i fig1S5.py
run -i fig1S6.py
run -i fig4S1.py
run -i fig5S2.py
run -i ipos.py

Note that depending on your system, you may have to call plt.show() to actually display figure windows after running one or more of the above scripts. Use cf() to quickly close all figures.