Data and Code for "Learning dynamically regulates stimulus discrimination of ventral striatal D1 receptor expressing neurons"

This code is used to perform analyses for the manuscript entitled “Learning dynamically regulates stimulus discrimination of ventral stratal D1 receptor expressing neurons.” CueDiscrimination.m script calculates neural activity and behavioral responses. Set parameters before running (specified below for each figure). Fig. 1D,E Behavioral parameters group = 'All'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; Fig. 3B-E Neural activity and behavioral parameters group = 'D1'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'selectivity'; cue = ''; Fig. 3G-J Neural activity and behavioral parameters group = 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'selectivity'; cue = ''; Fig. 4C,E Neural activity parameters group = 'D1'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS+'; Fig. 4D,F Neural activity parameters group = 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS+'; Fig. 4I Neural activity parameters group = 'D1'; % or 'D2' trialBlocks = [1:25]; bins = 12:35; output = 'modulation'; cue = 'CS+'; Fig. 5C,E Neural activity parameters group = 'D1'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS-'; Fig. 5D,F Neural activity parameters group = 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS-'; Fig. 5I Neural activity parameters group = 'D1'; % or 'D2' trialBlocks = [1:25]; bins = 12:35; output = 'modulation'; cue = 'CS-'; Fig. 6F-H Behavioral parameters group = 'ChR2 (laser)'; % or 'ChR2 (no laser)', 'Control (laser)', or 'Control (no laser)' trialBlocks = [1:100]; bins = 12:35; Fig. S1A-C Behavioral parameters group = 'D1'; % or 'D2' trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; Fig. S1D-E Behavioral parameters group = 'All'; trialBlocks = [1:5; 6:10; 11:15; 16:20; 21:25]; bins = 12:35; Fig. S2A-D ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity and behavioral parameters group = 'D1'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'selectivity'; cue = ''; Fig. S2E-H ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity and behavioral parameters group = 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'selectivity'; cue = ''; Fig. S2I ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity parameters group = 'D1'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS+'; Fig. S2J ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity parameters group = 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS+'; Fig. S2L ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity parameters group = 'D1'; % or 'D2' trialBlocks = [1:25]; bins = 12:35; output = 'modulation'; cue = 'CS+'; Fig. S2M ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity parameters group = 'D1'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS-'; Fig. S2N ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity parameters group = 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:35; output = 'modulation'; cue = 'CS-'; Fig. S2P ***In AnalyzeNeuralActivity, load data file Data_3msTag.mat Neural activity parameters group = 'D1'; % or 'D2' trialBlocks = [1:25]; bins = 12:35; output = 'modulation'; cue = 'CS-'; Fig. S3 Neural activity parameters group = 'D1'; % or 'D2'; trialBlocks = [1:100]; bins = 12:35; output = 'modulation'; cue = 'CS+'; % or 'CS-' Fig. S4 Neural activity parameters group = 'D1'; % or 'D2'; trialBlocks = [1:5; 6:10; 11:15; 16:20; 21:25]; bins = 12:35; output = 'selectivity'; cue = 'CS+'; % or 'CS-' Fig. S6A-D Neural activity parameters group = 'D1'; % or 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 36:40; output = 'modulation'; cue = 'CS+'; Fig. S6F Neural activity parameters group = 'D1'; % or 'D2'; trialBlocks = [1:25]; bins = 36:40; output = 'modulation'; cue = 'CS+'; Fig. S7C,D,G,H Neural activity and behavioral parameters group = 'D1'; % or 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:20; output = 'selectivity'; cue = ''; Fig. S7E,F,I,J Neural activity parameters group = 'D1'; % or 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 12:20; output = 'modulation'; cue = 'CS+'; % or 'CS-' Fig. S7K,L,O,P Neural activity and behavioral parameters group = 'D1'; % or 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 22:35; output = 'selectivity'; cue = ''; Fig. S7M,N,Q,R Neural activity parameters group = 'D1'; % or 'D2'; trialBlocks = [1:25; 26:50; 51:75; 76:100]; bins = 22:35; output = 'modulation'; cue = 'CS+'; % or 'CS-'