- 积分
- 48
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-8
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我在做一个臭氧对应紫外辐射高低值年(已知2001为高值年,2008为低值年)的求平均,然后进行T检验的工作。现有2001和2008年臭氧nc数据,其中o3
信息为
编了求平均和T检验的程序如下: %读取高值年2001年的数据 o3_01=ncread('o3_2001_150-1hPa.nc','o3');
level_01=ncread('o3_2001_150-1hPa.nc','level');
long_01=ncread('o3_2001_150-1hPa.nc','longitude');
lat_01=ncread('o3_2001_150-1hPa.nc','latitude');
time_01=ncread('o3_2001_150-1hPa.nc','time');
%读取低值年2008年的数据
o3_08=ncread('o3_2008_150-1hPa.nc','o3');
level_08=ncread('o3_2008_150-1hPa.nc','level');
long_08=ncread('o3_2008_150-1hPa.nc','longitude');
lat_08=ncread('o3_2008_150-1hPa.nc','latitude');
time_08=ncread('o3_2008_150-1hPa.nc','time');
x=nanmean(o3_01,4);
y=nanmean(o3_08,4);
%作图
for k=1:13
figure;
contour(lat_01,long_01,o3_01(:,:,k));
contour(lat_08,long_08,o3_08(:,:,k));
end
for i=1:241
for j=1:480
x1=squeeze(x(j,i,:));
y1=squeeze(y(j,i,:));
s=ttest2(x1,y1);
ss(j,i)=s;
end
end
%作T检验图
figure;
contour(ss)
然后只出了13张臭氧图,最后的T检验图上什么都没有
并且显示
而且打开变量ss后发现全是0,T检验该怎么弄?
|
|