os.path.join(home_dir, '200303/2020-03-03_19-57-03/'), 32, os.path.join(home_dir, '200303/2020-03-03trial_names_ventral.txt')) rec5 = br.Binary_recording( os.path.join(home_dir, '200309/2020-03-09_16-20-42/'), 32, os.path.join(home_dir, '200309/2020-03-09trial_name_joined.txt')) rec6 = br.Binary_recording( os.path.join(home_dir, '200311/2020-03-11_16-51-10/'), 32, os.path.join(home_dir, '200311/2020-03-11trial_name_binary_joined.txt')) rec7 = br.Binary_recording( os.path.join(home_dir, '200318/2020-03-18_15-24-43/'), 32, os.path.join(home_dir, '200318/2020-03-18trial_name.txt')) rec8 = br.Binary_recording( os.path.join(home_dir, '200319/2020-03-19_16-08-45/'), 32, os.path.join(home_dir, '200319/2020-03-19_16-08-45_trial_names.txt')) recs = jr.JoinedRecording( recordings=[rec1, rec2, rec3, rec4, rec5, rec6, rec7, rec8]) rec1_sb = br.Binary_recording( os.path.join(home_dir, '200228/2020-02-28_19-56-29/'), 32, os.path.join(home_dir, '200228/2020-02-28trial_names_ventral.txt'), sniff_basis=True) rec2_sb = br.Binary_recording( os.path.join(home_dir, '200228/2020-02-28_16-37-36/'), 32, os.path.join(home_dir, '200228/2020-02-28trial_names_dorsal.txt'), sniff_basis=True) rec3_sb = br.Binary_recording( os.path.join(home_dir, '200303/2020-03-03_16-44-23/'), 32, os.path.join(home_dir, '200303/2020-03-03trial_names_dorsal.txt'),
scaled_data = scaled_data - comp1_project pcs.append(pca.components_) scaled_data = scaler.fit_transform(scaled_data.T).T return pcs, scaled_data available_cpu_count = len(psutil.Process().cpu_affinity()) os.environ["MKL_NUM_THREADS"] = str(available_cpu_count) home_dir = '/home/camp/warnert/working/Recordings/binary_pulses' rec1 = br.Binary_recording(os.path.join(home_dir, '200228/2020-02-28_19-56-29/'), 32, os.path.join(home_dir, '200228/2020-02-28trial_names_ventral.txt'), sniff_basis=True) rec2 = br.Binary_recording(os.path.join(home_dir, '200228/2020-02-28_16-37-36/'), 32, os.path.join(home_dir, '200228/2020-02-28trial_names_dorsal.txt'), sniff_basis=True) rec3 = br.Binary_recording(os.path.join(home_dir, '200303/2020-03-03_16-44-23/'), 32, os.path.join(home_dir, '200303/2020-03-03trial_names_dorsal.txt'), sniff_basis=True) rec4 = br.Binary_recording(os.path.join(home_dir, '200303/2020-03-03_19-57-03/'), 32, os.path.join(home_dir, '200303/2020-03-03trial_names_ventral.txt'), sniff_basis=True) rec5 = br.Binary_recording(os.path.join(home_dir, '200309/2020-03-09_16-20-42/'), 32, os.path.join(home_dir, '200309/2020-03-09trial_name_joined.txt'), sniff_basis=True) rec6 = br.Binary_recording(os.path.join(home_dir, '200311/2020-03-11_16-51-10/'), 32, os.path.join(home_dir, '200311/2020-03-11trial_name_binary_joined.txt'), sniff_basis=True) rec7 = br.Binary_recording(os.path.join(home_dir, '200318/2020-03-18_15-24-43/'), 32, os.path.join(home_dir, '200318/2020-03-18trial_name.txt'), sniff_basis=True) rec8 = br.Binary_recording(os.path.join(home_dir, '200319/2020-03-19_16-08-45/'), 32, os.path.join(home_dir, '200319/2020-03-19_16-08-45_trial_names.txt'), sniff_basis=True) recs = jr.JoinedRecording(recordings=[rec1, rec2, rec3, rec4, rec5, rec6, rec7, rec8]) trial_names = ['%d_1' % i for i in range(32)] resps = recs.get_multi_trial_response(trial_names, pre_trial_window=2, post_trial_window=2, w_bootstrap=True, bootstrap_limit=1000, bin_size=0.05) resps = np.array(resps) scaler = StandardScaler() reshaped_resps = np.rollaxis(resps, axis=2).reshape(145, 32*1000*80) pcs, scaled_data = SCA(reshaped_resps) np.save(os.path.join(home_dir, 'SCA/200909_pcs.npy'), pcs) np.save(os.path.join(home_dir, 'SCA/200909_data.npy'), scaled_data)