求助，关于降水与位涡相关系数的T检验

勇敢的悲伤

下面是我求T值的程序，和公式对照过，感觉没有算错。但是出来的值都是16点多，这个肯定不对了吧，但是看了一下午都没有看出来错在哪里，麻烦各位高手给看一下。谢过了·~~

program main
implicit none
integer,parameter::nt=92,nx=144,ny=37,undef=-9.99e+08
integer i,j,t,num
real r1(nt),r2(nt),c(nt),s1,s2(nx,ny),fc1,fc2(nx,ny)
real r(nt),mpv(nx,ny,nt),tjianyan(nx,ny),nr(nx,ny)
open(1,file='f:\huatu\xgr1961-1991.grd',form='binary')
open(2,file='f:\huatu\850mpv-1961-1991.grd',form='binary')
do t=1,nt
       read(1) r(t)
enddo
      
do t=1,nt
do j=1,ny
do i=1,nx
  read(2) mpv(i,j,t)
  enddo
  enddo  
enddo
!计算降水量平均值
s1=0
do i=1,nt
s1=s1+r(i)
enddo
s1=s1/nt
!计算降水方差
fc1=0
do i=1,nt
fc1=fc1+(r(i)-s1)**2
enddo
fc1=fc1/nt
!计算位涡平均值
do j=1,ny
do i=1,nx
s2(i,j)=0
do t=1,nt
if(mpv(i,j,t).ne.undef)then
s2(i,j)=s2(i,j)+mpv(i,j,t)
else
s2(i,j)=undef
endif
enddo
enddo
enddo
do j=1,ny
do i=1,nx
if(s2(i,j).ne.undef)then
s2(i,j)=s2(i,j)/nt
else
s2(i,j)=undef
endif
enddo
enddo
!计算位涡方差
do j=1,ny
do i=1,nx
fc2(i,j)=0
do t=1,nt
if(mpv(i,j,t).ne.undef.and.s2(i,j).ne.undef)then
fc2(i,j)=fc2(i,j)+(mpv(i,j,t)-s2(i,j))**2
else
fc2(i,j)=undef
endif
enddo
enddo
enddo
do j=1,ny
do i=1,nx
if(fc2(i,j).ne.undef)then
fc2(i,j)=fc2(i,j)/nt
else
fc2(i,j)=undef
endif
enddo
enddo
!计算t值
do j=1,ny
do i=1,nx
if(fc2(i,j).ne.undef.and.s2(i,j).ne.undef)then
tjianyan(i,j)=(s1-s2(i,j))/sqrt(((nt-1)*fc1+(nt-1)*fc2(i,j))/(nt+nt-2)*sqrt(1.0/nt+1.0/nt))
else
tjianyan(i,j)=undef
endif
enddo
enddo

open(3,file='f:\huatu\tjianyan.txt')
do j=1,ny
do i=1,nx
write(3,*) abs(tjianyan(i,j))
enddo
enddo

end

river

这么长的程序，没人会从头到尾看一遍的

wuwei2163

迷茫啊，看了半天都不明白

cris

进行两独立样本t检验的前提条件：
1.两样本应该是相互独立。
2.样本来自的两个总体应该服从正态分布
你的降水和位势涡度是否符合上述两个条件呢？问题可能出在这里，你需要将降水数据以及位涡数据进行标准化处理，从而使处理后的数据呈正态分布。

勇敢的悲伤

wuwei2163 发表于 2013-6-24 23:56

                               
登录/注册后可看大图

迷茫啊，看了半天都不明白

其实就是根据T检验的那个公式算的。。。

勇敢的悲伤

cris 发表于 2013-6-25 10:00

                               
登录/注册后可看大图

进行两独立样本t检验的前提条件：
1.两样本应该是相互独立。
2.样本来自的两个总体应该服从正态分布

如果求的是相关系数的显著性，是不是可以直接用相关系数来作为验证？意思就是相关系数超过哪个数的话就通过验证了？？

cris

勇敢的悲伤 发表于 2013-6-26 22:01

                               
登录/注册后可看大图

如果求的是相关系数的显著性，是不是可以直接用相关系数来作为验证？意思就是相关系数超过哪个数的话就通 ...

可以这么理解。你可以查表，置信度α的选取可以根据需要来，当两者的相关系数超过你所选取的特定置信度对应的某个数值时候，就可以说他们通过了检验，二者显著相关。

勇敢的悲伤

cris 发表于 2013-6-27 10:14

                               
登录/注册后可看大图

可以这么理解。你可以查表，置信度α的选取可以根据需要来，当两者的相关系数超过你所选取的特定置信度对 ...

只要是相关系数的绝对值超过某一值就可以了吧？还是只能正值呢？

cris

勇敢的悲伤 发表于 2013-6-27 11:01

                               
登录/注册后可看大图

只要是相关系数的绝对值超过某一值就可以了吧？还是只能正值呢？

如果是正值，说明二者是正相关的；反之，如果是负值就是说明负相关。

勇敢的悲伤

cris 发表于 2013-6-27 11:20

                               
登录/注册后可看大图

如果是正值，说明二者是正相关的；反之，如果是负值就是说明负相关。

嗯嗯。那就应该对了。哈哈~~谢谢啦~~