一维的时间指数和降水回归后做检验时怎么求t值

感觉自己萌萌哒

首先谢谢大家的帮助
我是用一个事件的爆发时间指数标准化后和32年间的5月全球月平均降水做回归，ncl有说用student_t做检验，但是看到很多说用公式存在一些问题，建议直接求t值查表，我自己也试了一下student_t结果不太对
然后尝试直接算t，可是时间指数是一个一维数组并没有经纬度和降水进行匹配，直接套了图中的公式进去算结果也不是正常范围结果达到了几百，想问一下这种情况怎么算t值啊，谢谢

求t值我是这么写的，SFW是32年标准化的时间指数，precip是格点144*72的32年5月月平均降水，

画图画的t>1.96和t<-1.96

n1=32*144*72
n2=32
X1 = dim_avg_n_Wrap(precip, 0)
S1 = dim_variance_n_Wrap(precip,0)
X2 = dim_avg_Wrap(SFW)
S2 = dim_variance_Wrap(SFW)
t = (X1-X2)/sqrt((1.0/n1+1.0/n2)*(S1*(n1-1.0)+S2*(n2-1.0))/(n1+n2-2.0))
t!0 ="lat"
t&lat = lat
t!1 ="lon"
t&lon = lon

Sybil11111

请问你解决了吗，可否指教一下~