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

https://doi.org/10.7554/eLife.70469

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.