信度检验的问题，求助

sosbbxc

最近在折腾信度检验的问题，我选取了28次事件，截取事件开始前60天和后60天，一共121天的距平，再将这28次事件的121天距平进行合成，然后选取合成的121天距平中，过程前31天，也就是-40到-10进行平均，作为前期，再选取过程后31天，就是10到40天，作为后期，接着用后期的平均减去前期的平均，得出差异，所以我就想做一下差异的显著性检验。

1. program xgxsimplicit noneinteger ix,iy,it,nx,ny,nt,n1,n2parameter(nx=144,ny=37,nt=121,n1=31,n2=31)real hgt(nx,ny,nt)     real ::aves(nx,ny)=0,avew(nx,ny)=0,fcs(nx,ny)=0,fcw(nx,ny)=0,res(nx,ny)=0open(11,file='D:\monsoon\jianyan\1000_u.dat',form='binary')do it=1,ntdo iy=1,nydo ix=1,nxread(11) hgt(ix,iy,it)end doend doend doclose(11)do it=21,51do iy=1,nydo ix=1,nxaves(ix,iy)=aves(ix,iy)+hgt(ix,iy,it)end doend doend dodo it=71,101do iy=1,nydo ix=1,nxavew(ix,iy)=avew(ix,iy)+hgt(ix,iy,it)end doend doend dodo iy=1,nydo ix=1,nxaves(ix,iy)=aves(ix,iy)/n1avew(ix,iy)=avew(ix,iy)/n2end do end dodo it=21,51do iy=1,nydo ix=1,nxfcs(ix,iy)=fcs(ix,iy)+(hgt(ix,iy,it)-aves(ix,iy))**2end doend doend dodo it=71,101do iy=1,nydo ix=1,nxfcw(ix,iy)=fcw(ix,iy)+(hgt(ix,iy,it)-avew(ix,iy))**2end doend doend dodo iy=1,nydo ix=1,nxfcs(ix,iy)=fcs(ix,iy)/n1fcw(ix,iy)=fcw(ix,iy)/n2end doend dodo iy=1,nydo ix=1,nxres(ix,iy)=(avew(ix,iy)-aves(ix,iy))/sqrt((((n1-1)*fcs(ix,iy)+(n2-1)*fcw(ix,iy))/(n1+n2-2))*(1.0/real(n1)+1.0/real(n2)))end doend doopen(22,file='1k_z.dat',form='binary')do iy=1,nydo ix=1,nxwrite(22)res(ix,iy)end doend doclose(22)end

program xgxsimplicit noneinteger ix,iy,it,nx,ny,nt,n1,n2parameter(nx=144,ny=37,nt=121,n1=31,n2=31)real hgt(nx,ny,nt)     real ::aves(nx,ny)=0,avew(nx,ny)=0,fcs(nx,ny)=0,fcw(nx,ny)=0,res(nx,ny)=0open(11,file='D:\monsoon\jianyan\1000_u.dat',form='binary')do it=1,ntdo iy=1,nydo ix=1,nxread(11) hgt(ix,iy,it)end doend doend doclose(11)do it=21,51do iy=1,nydo ix=1,nxaves(ix,iy)=aves(ix,iy)+hgt(ix,iy,it)end doend doend dodo it=71,101do iy=1,nydo ix=1,nxavew(ix,iy)=avew(ix,iy)+hgt(ix,iy,it)end doend doend dodo iy=1,nydo ix=1,nxaves(ix,iy)=aves(ix,iy)/n1avew(ix,iy)=avew(ix,iy)/n2end do end dodo it=21,51do iy=1,nydo ix=1,nxfcs(ix,iy)=fcs(ix,iy)+(hgt(ix,iy,it)-aves(ix,iy))**2end doend doend dodo it=71,101do iy=1,nydo ix=1,nxfcw(ix,iy)=fcw(ix,iy)+(hgt(ix,iy,it)-avew(ix,iy))**2end doend doend dodo iy=1,nydo ix=1,nxfcs(ix,iy)=fcs(ix,iy)/n1fcw(ix,iy)=fcw(ix,iy)/n2end doend dodo iy=1,nydo ix=1,nxres(ix,iy)=(avew(ix,iy)-aves(ix,iy))/sqrt((((n1-1)*fcs(ix,iy)+(n2-1)*fcw(ix,iy))/(n1+n2-2))*(1.0/real(n1)+1.0/real(n2)))end doend doopen(22,file='1k_z.dat',form='binary')do iy=1,nydo ix=1,nxwrite(22)res(ix,iy)end doend doclose(22)end

程序自己编的。有点粗犷。莫见怪
程序是根据下面这个公式写 的
avew就是过程后期，aves就是过程前期。算出来的结果怎么看都不对。t检验值能到十几吗。。按照书上查的60个样本量95%置信水平t检验量不过2.0啊。。不知道是不是方法用的不对，请高人指点。。一般文献上不会解释这种细致的小问题。所以很费解。。

sosbbxc

program xgxs
implicit none
integer ix,iy,it,nx,ny,nt,n1,n2
parameter(nx=144,ny=37,nt=121,n1=31,n2=31)
real hgt(nx,ny,nt)     
real ::aves(nx,ny)=0,avew(nx,ny)=0,fcs(nx,ny)=0,fcw(nx,ny)=0,res(nx,ny)=0
open(11,file='D:\monsoon\jianyan\1000_u.dat',form='binary')
do it=1,nt
do iy=1,ny
do ix=1,nx
read(11) hgt(ix,iy,it)
end do
end do
end do
close(11)

do it=21,51
do iy=1,ny
do ix=1,nx
aves(ix,iy)=aves(ix,iy)+hgt(ix,iy,it)
end do
end do
end do

do it=71,101
do iy=1,ny
do ix=1,nx
avew(ix,iy)=avew(ix,iy)+hgt(ix,iy,it)
end do
end do
end do

do iy=1,ny
do ix=1,nx
aves(ix,iy)=aves(ix,iy)/n1
avew(ix,iy)=avew(ix,iy)/n2
end do
end do

do it=21,51
do iy=1,ny
do ix=1,nx
fcs(ix,iy)=fcs(ix,iy)+(hgt(ix,iy,it)-aves(ix,iy))**2
end do
end do
end do

do it=71,101
do iy=1,ny
do ix=1,nx
fcw(ix,iy)=fcw(ix,iy)+(hgt(ix,iy,it)-avew(ix,iy))**2
end do
end do
end do

do iy=1,ny
do ix=1,nx
fcs(ix,iy)=fcs(ix,iy)/n1
fcw(ix,iy)=fcw(ix,iy)/n2
end do
end do

do iy=1,ny
do ix=1,nx
res(ix,iy)=(avew(ix,iy)-aves(ix,iy))/sqrt((((n1-1)*fcs(ix,iy)+(n2-1)*fcw(ix,iy))/(n1+n2-2))*(1.0/real(n1)+1.0/real(n2)))
end do
end do

open(22,file='1k_z.dat',form='binary')
do iy=1,ny
do ix=1,nx
write(22)res(ix,iy)
end do
end do
close(22)

end
代码怎么成这样了。。

river

求和你都不赋初值吗？

sosbbxc

river 发表于 2014-11-5 07:28
求和你都不赋初值吗？

前面不是赋了一堆初值为0的数组了吗？

river

sosbbxc 发表于 2014-11-5 11:14
前面不是赋了一堆初值为0的数组了吗？

你自己把前几步累加计算一下，看看是怎么个循环的。。。初值不是直接赋值在最前面就行了的，这个是要根据循环来的好伐