StdVisualModel/s1_get_ROImean.m at master · PumpkinPop/StdVisualModel · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
clear all; close all; clc

% This script is create to collect the target data from raw dataset

%%
% Add path of data to the script
addpath(genpath(fullfile(pwd,'fMRIdata')));
save_address = fullfile(pwd, 'ROImean' );

%% Load Dataset 1 (Catherine 06-19)

load roiBetas

% v1 mean ( voxel _ mean _ dataset)
v_mean_69( 1 , : ) = roiBetas.roiBetamn{11};

% v2 mean
v_mean_69( 2 , : ) = roiBetas.roiBetamn{12};

% v3 mean
v_mean_69( 3 , : ) = roiBetas.roiBetamn{13};

% save
save([save_address , '\v_mean_69'] , 'v_mean_69')

% delete the variable to save up memory
clear ( 'roiBetas')

%% Load Dataset 2 (Catherine 10-05)

load roiBetas2

% v1 mean
v_mean_05( 1 , : ) = (roiBetas.roiBetamn{1} + roiBetas.roiBetamn{5})./2;

% v2 mean
v_mean_05( 2 , : ) = (roiBetas.roiBetamn{2} + roiBetas.roiBetamn{6})./2;

% v3 mean
v_mean_05( 3 , : ) = (roiBetas.roiBetamn{3} + roiBetas.roiBetamn{7})./2;

% save
save([save_address , '\v_mean_05'] , 'v_mean_05')

% delete the variable to save up memory
clear ( 'roiBetas')

%% Load Dataset 3 (Kendrick dataset 3)

% Choose the right stimuli we need
labelVec_K = [[70:1:99],[107:1:115]];

% load data
load dataset03

% v1 mean
v_mean_K1( 1 , : ) = mean(betamn(roi==2 , labelVec_K));

% v2 mean
v_mean_K1( 2 , : ) = mean(betamn(roi==3 , labelVec_K));

% v3 mean
v_mean_K1( 3 , : ) = mean(betamn(roi==4 , labelVec_K));

% save
save([save_address , '\v_mean_K1'] , 'v_mean_K1')

%% Load Dataset 4 (Kendrick dataset 4)

% Choose the right stimuli we need
labelVec_K = [[70:1:99],[107:1:115]];

% load data
load dataset04

% v1 mean
v_mean_K2( 1 , : ) = mean(betamn(roi==2 , labelVec_K));

% v2 mean
v_mean_K2( 2 , : ) = mean(betamn(roi==3 , labelVec_K));

% v3 mean
v_mean_K2( 3 , : ) = mean(betamn(roi==4 , labelVec_K));

% save
save([save_address , '\v_mean_K2'] , 'v_mean_K2')