Homec4science
Diffusion ShoulderDB_ImportSegmentMeasureAnalyze (aterrier-plotUpdate) Browse

R8218/statSCase.m
c003a89f3c06aterrier-plotUpdate

Last Change
Raw File

statSCase.m
View Options

% This script load the entire SCaseDB and performs statisctical analyses
%% Initialisation of paths
% Should be replaced by a function
dataDir = '../../../data'; % location of data
logDir = 'log'; % log file in xls folder
xlsDir = [dataDir '/Excel/xlsFromMatlab'];
matlabDir = [dataDir '/matlab'];
tic;
%% Load ScaseDB
matlabFile = 'SCaseDB.mat';
matlabFile = [matlabDir '/' matlabFile];
if exist(matlabFile, 'file')
load(matlabFile); % Load SCase dat file
else
error('Not valid matlab file');
end
%% Select SCase with valid anatomical measurements
SCase = ShoulderCase; % Create new object ShoulderCase
iSCase = 0; % Index for SCase
for iSCaseDB = 1:length(SCaseDB) % Loop on all SCase
if ~isempty(SCaseDB(iSCaseDB).shoulder.side)
iSCase = iSCase + 1;
SCase(iSCase) = SCaseDB(iSCaseDB);
end
end
%% Select normal and pathological SCase
SCaseN = ShoulderCase; % Create new object ShoulderCase for normal
SCaseP = ShoulderCase; % Create new object ShoulderCase for pahtological
iSCaseN = 0; % Index for SCaseN
iSCaseP = 0; % Index for SCaseP
for iSCase = 1:length(SCase) % Loop on all valide SCase
if strcmp(SCase(iSCase).id(1), 'N') % SCase is normal
iSCaseN = iSCaseN + 1;
SCaseN(iSCaseN) = SCase(iSCase);
else
iSCaseP = iSCaseP + 1;
SCaseP(iSCaseP) = SCase(iSCase);
end
end
%% Compare normal to pathological
% Check way to compare histogram with different number of patients
for iSCaseN = 1:length(SCaseN)
varN(iSCaseN) = SCaseN(iSCaseN).shoulder.scapula.acromion.IE;
end
for iSCaseP = 1:length(SCaseP)
varP(iSCaseP) = SCaseP(iSCaseP).shoulder.scapula.acromion.IE;
end
hold off
histogram(varN);
hold on;
histogram(varP);
%%
toc
c4science · Help