Repository navigation
Expand file tree
/
Copy pathplot_baseline_intervention.m
More file actions
57 lines (52 loc) · 1.95 KB
/
Copy pathplot_baseline_intervention.m
File metadata and controls
57 lines (52 loc) · 1.95 KB
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
subplot(2,2,1)
load('ParamVectorsDokanTofa.mat');
plot((1:length(mfPrevArray(:,1)))/12,100*(mfPrevArray(:,:)),'Color',[0.6 0.6 0.6]);
hold on;
aa = ConstructBinomialErrorBars(MfData);
p1= errorbar(aa(:,1),aa(:,2),aa(:,3),aa(:,4),'or','MarkerFaceColor','r','MarkerEdgeColor','r');
hold on;
xlabel('Age (years)');
ylabel('Mf prevalence (%)');
title('Low prevalence site (DokanTofa)');
hold on;
subplot(2,2,2)
load('ParamVectorsPeneng.mat');
plot((1:length(mfPrevArray(:,1)))/12,100*(mfPrevArray(:,:)),'Color',[0.6 0.6 0.6]);
hold on;
aa = ConstructBinomialErrorBars(MfData1);
p1= errorbar(aa(:,1),aa(:,2),aa(:,3),aa(:,4),'-or','MarkerFaceColor','r','MarkerEdgeColor','r');
hold on;
aa = ConstructBinomialErrorBars(MfData2);
p1= errorbar(aa(:,1),aa(:,2),aa(:,3),aa(:,4),'-ob','MarkerFaceColor','b','MarkerEdgeColor','b');
hold on;
aa = ConstructBinomialErrorBars(MfData3);
p1= errorbar(aa(:,1),aa(:,2),aa(:,3),aa(:,4),'-ok','MarkerFaceColor','k','MarkerEdgeColor','k');
hold on;
xlabel('Age (years)');
ylabel('Mf prevalence (%)');
title('High prevalence site (Peneng)');
hold on;
subplot(2,2,3)
load('Intv_WK_test_new_DokanTofa.mat')
plot((1:length(mfPrevIntv1(:,1)))/12,(mfPrevIntv1(:,:)),'Color',[0.6 0.6 0.6]);
hold on;
MfData1(:,1)=MfData1(:,1)/12;
aa = ConstructBinomialErrorBars(MfData1);
p1= errorbar(aa(:,1),aa(:,2),aa(:,3),aa(:,4),'or','MarkerFaceColor','r','MarkerEdgeColor','r');
hold on;
set(gca,'XLim',[0 7]);
xlabel('Time (years)');
ylabel('Mf prevalence (%)');
hold on;
subplot(2,2,4)
load('Intv_WK_test_new_Peneng.mat')
plot((1:length(mfPrevIntv1(:,1)))/12,(mfPrevIntv1(:,:)),'Color',[0.6 0.6 0.6]);
hold on;
MfData0(:,1)=MfData0(:,1)/12;
aa = ConstructBinomialErrorBars(MfData0);
p1= errorbar(aa(:,1),aa(:,2),aa(:,3),aa(:,4),'or','MarkerFaceColor','r','MarkerEdgeColor','r');
hold on;
set(gca,'XLim',[0 7]);
xlabel('Time (years)');
ylabel('Mf prevalence (%)');
hold on;