-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathwrapper_pointingError.m
73 lines (52 loc) · 2.01 KB
/
wrapper_pointingError.m
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
%% ==== 'pointingError' ANALYSIS ==== %%
%% Takes X error and Y error to compute pointing error
absList = [0,1]; %absolute = 1, no absolute = 0
for a = 1:length(absList)
absSwitch = absList(a);
if ~absSwitch
thisVarStr = 'pointingError';
axisVal_allTargets_y_ttest = [-10 50];
axisVal_sideOfSpace_y_ttest = [-10 30];
axisVal_allTargets_y_BDST = [-10 40];
axisVal_sideOfSpace_y_BDST = [-10 40];
%% Grab X/Y Error
[dataX,tmpH] = getData('XError',rawData,sNames);
[dataY,tmpH] = getData('YError',rawData,sNames);
elseif absSwitch
thisVarStr = 'pointingError_ABSOLUTE';
axisVal_allTargets_y_ttest = [-5 50];
axisVal_sideOfSpace_y_ttest = [-5 25];
axisVal_allTargets_y_BDST = [-40 40];
axisVal_sideOfSpace_y_BDST = [-25 25];
%% Grab X/Y Error
[dataX,tmpH] = getData_AbsoluteVersion('XError',rawData,sNames);
[dataY,tmpH] = getData_AbsoluteVersion('YError',rawData,sNames);
end
%% Compute distance error (from Rob):
% sqroot( squ(x error) + squ(y error) )
data = [];
for c = 1:4
tmpX = dataX.targetMean(:,:,c);
tmpY = dataY.targetMean(:,:,c);
tmpX = tmpX.^2;
tmpY = tmpY.^2;
tmpD = (tmpX + tmpY);
data.targetMean(:,:,c) = sqrt(tmpD);
end
%% Run main deficit tests
outDir = fullfile('results',thisVarStr);
try
rmdir(outDir,'s');
catch
%NOOP
end
mkdir(outDir)
diary(fullfile(outDir,[thisVarStr,'_inputForCrawford-SingleBayes_ES.txt']));
stats = doCrawfordTtest_allTargets(thisVarStr,data,conditionNames,outDir,axisVal_allTargets_y_ttest);
stats = doCrawfordTtest_sideOfSpace(thisVarStr,data,conditionNames,outDir,axisVal_sideOfSpace_y_ttest);
diary OFF
diary(fullfile(outDir,[thisVarStr,'_inputForCrawford-DissocsBayes_ES.txt']));
plotBDST_allTargets(thisVarStr,data,conditionNames,outDir,axisVal_allTargets_y_BDST);
plotBDST_sideOfSpace(thisVarStr,data,conditionNames,outDir,axisVal_sideOfSpace_y_BDST);
diary OFF
end