-
Notifications
You must be signed in to change notification settings - Fork 1
/
OSLD_Stat_Summary.m
83 lines (68 loc) · 1.89 KB
/
OSLD_Stat_Summary.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
74
75
76
77
78
79
80
81
82
83
clear all;
close all;
clc;
load('OSLD_Accum.mat');
%DIMENSIONAL Organization:
%Patients, Trials, Files, Positions
%File Order:
%IMRT MASK
%IMRT NO MASK
%VMATSS MASK
%VMATSS NO MASK
%VMAT MASK
%VMAT NO MASK
%Every other file is a mask file
%Average across trials for each position, only mask:
mT = squeeze(mean(squeeze(dose(1,:,1:2:end,:)))).*35/100;
sT = squeeze(std(squeeze(dose(1,:,1:2:end,:)))).*35/100;
difIMRT = mT(3,:) - mT(1,:);
errIMRT = sqrt(sT(3,:).^2 + sT(1,:).^2);
difVMATSS = mT(3,:) - mT(2,:);
errVMATSS = sqrt(sT(3,:).^2 + sT(2,:).^2);
figure(1);
plot(difIMRT,'*'); hold all;
plot(errIMRT,'*');
legend('VMAT - IMRT','Error');
xlabel('Position');
ylabel('Measured Skin Dose (Gy)');
set(gca,'XTick',1:4,'XTickLabel',['R Inf';'L Inf';'R Sup';'L Sup']);
xlim([0.5 4.5]);
title('MASK');
ylim([-5 10]);
figure(2);
plot(difVMATSS,'*'); hold all;
plot(errVMATSS,'*');
legend('VMAT - VMATSS','Error');
xlabel('Position');
ylabel('Measured Skin Dose (Gy)');
set(gca,'XTick',1:4,'XTickLabel',['R Inf';'L Inf';'R Sup';'L Sup']);
xlim([0.5 4.5]);
title('MASK');
ylim([0 25]);
%Average across trials for each position, only no mask:
mT = squeeze(mean(squeeze(dose(1,:,2:2:end,:)))).*35/100;
sT = squeeze(std(squeeze(dose(1,:,2:2:end,:)))).*35/100;
difIMRT = mT(3,:) - mT(1,:);
errIMRT = sqrt(sT(3,:).^2 + sT(1,:).^2);
difVMATSS = mT(3,:) - mT(2,:);
errVMATSS = sqrt(sT(3,:).^2 + sT(2,:).^2);
figure(3);
plot(difIMRT,'*'); hold all;
plot(errIMRT,'*');
legend('VMAT - IMRT','Error');
xlabel('Position');
ylabel('Measured Skin Dose (Gy)');
set(gca,'XTick',1:4,'XTickLabel',['R Inf';'L Inf';'R Sup';'L Sup']);
xlim([0.5 4.5]);
title('NO MASK');
ylim([-5 10]);
figure(4);
plot(difVMATSS,'*'); hold all;
plot(errVMATSS,'*');
legend('VMAT - VMATSS','Error');
xlabel('Position');
ylabel('Measured Skin Dose (Gy)');
set(gca,'XTick',1:4,'XTickLabel',['R Inf';'L Inf';'R Sup';'L Sup']);
xlim([0.5 4.5]);
title('NO MASK');
ylim([0 25]);