%Guide for performing a linear-nonlinear model (LN-model)
%analysis from voltage recordings of bipolar cells (BC) under natural
%stimuli
%Check out the Documentation_Data.pdf, which contains relevant information
%Check BC_mainWhiteNoise.m -> from here you get the filter and output
%nonlinearity for performing the prediction to the movie.
clear;
close all;

%% Preprocessing of the natural movie contrast signal
loadData='Z:\naturalMovies\'; %path of the BC files for natural movies (has to be adapted to your own path)
filename='16d200016Comp.mat'; %BC file to analyse: for example here we use file 16d200016Comp from cell 5, recorded in retina 2 on 161219
Hz=25; % given in the excel file LogInfo_NaturalMovies.xlsx for the corresponding file
movienbr=4; %given in the excel file LogInfo_NaturalMovies.xlsx for the corresponding file
pixelSize=9; %given in the excel file LogInfo_SpatioTemporalWhiteNoise.xlsx from the white noise stimulus of the corresponding cell
[movieSignal]=BC_getNaturalMovies(pixelSize);%example code to preprocess the natural movies

%% Prediction of the response to the movie from the LN-model with spatiotemporal white noise
[naturalMovie_output]=BC_predNaturalMovies(loadData,filename,movieSignal,STA_output,NL_output,Hz);%example code to get the prediction to natural movies from the spatiotemporal white noise LN-model