course-idma2016/matlab/idma_box_mean.m


tmp0=prof_T; tmp0.prof(prof_Tweight.prof==0)=NaN;
kk=2;

nr=prof_T.nr;
tmp_nb=NaN*zeros(240,nr);
tmp_mean=NaN*zeros(240,nr);
tmp_median=NaN*zeros(240,nr);
tmp_std=NaN*zeros(240,nr);
for yy=1992:2011;
for mm=1:12;
  dd=[datenum(yy,mm,0) datenum(yy,mm+1,0)];
  tt=(yy-1992)*12+mm;
  %
  tmp1=MITprof_subset(tmp0,'date',dd);
  tmp1=tmp1.prof;
  %
  tmp_nb(tt,:)=sum(~isnan(tmp1),1);
  tmp_mean(tt,:)=nanmean(tmp1,1);
  tmp_median(tt,:)=nanmedian(tmp1,1);
  tmp_std(tt,:)=nanstd(tmp1,1);
  %
end;
end;

figure;
tim=1992+[0.5:239.5]/12;
subplot(3,1,1); plot(tim,tmp_nb(:,kk));
legend('sample size','Location','NorthWest');
subplot(3,1,2); plot(tim,tmp_mean(:,kk)); ylabel('degree C');
hold on; plot(tim,tmp_median(:,kk),'r'); 
legend('sample mean','sample median','Location','NorthWest');
subplot(3,1,3); plot(tim,tmp_std(:,kk)); ylabel('degree C'); 
legend('sample std','Location','NorthWest');

1	gforget	1.1
2			tmp0=prof_T; tmp0.prof(prof_Tweight.prof==0)=NaN;
3			kk=2;
4
5			nr=prof_T.nr;
6			tmp_nb=NaN*zeros(240,nr);
7			tmp_mean=NaN*zeros(240,nr);
8			tmp_median=NaN*zeros(240,nr);
9			tmp_std=NaN*zeros(240,nr);
10			for yy=1992:2011;
11			for mm=1:12;
12			dd=[datenum(yy,mm,0) datenum(yy,mm+1,0)];
13			tt=(yy-1992)*12+mm;
14			%
15			tmp1=MITprof_subset(tmp0,'date',dd);
16			tmp1=tmp1.prof;
17			%
18			tmp_nb(tt,:)=sum(~isnan(tmp1),1);
19			tmp_mean(tt,:)=nanmean(tmp1,1);
20			tmp_median(tt,:)=nanmedian(tmp1,1);
21			tmp_std(tt,:)=nanstd(tmp1,1);
22			%
23			end;
24			end;
25
26			figure;
27			tim=1992+[0.5:239.5]/12;
28			subplot(3,1,1); plot(tim,tmp_nb(:,kk));
29			legend('sample size','Location','NorthWest');
30			subplot(3,1,2); plot(tim,tmp_mean(:,kk)); ylabel('degree C');
31			hold on; plot(tim,tmp_median(:,kk),'r');
32			legend('sample mean','sample median','Location','NorthWest');
33			subplot(3,1,3); plot(tim,tmp_std(:,kk)); ylabel('degree C');
34			legend('sample std','Location','NorthWest');
35