用R语言写的EP检验的代码，可以运行但是答案不对，有大佬可以来看看是哪出了问题吗

我是按照书上这个步骤套上去的，但是最后T值算出来特别大，跟答案也对不上。 （这个是我们数理统计的书，不是软件的）

data1=read.csv("C:\\User\\Desktop\\1.csv")
#数据预处理
data2=as.matrix(data1)
data2
data3=as.numeric(data2)
data3
#计算均值和二阶中心矩
mean3=mean(data3)
sum1=0
for (i in data3)
{
 a=(i-mean3)^2
 sum1=sum1+a
 print(sum1)
}
n=99
m2=sum1/n
m2

#存A
sum2=0
for (i in data3)
{
  S1=(-1)*(i-mean3)^2/(4*m2)
  A=exp(S1)
  sum2=sum2+A
  print(sum2)
}
A=sum2
A
#求b
k=head(data3,-1)
j=tail(data3,-1)
sum3=0
for (J in j){
for (K in k){
  S2=(-1)*(J-K)^2/(2*m2)
  B=exp(S2)
  sum3=sum3+B
  print(sum3)
  
}
}
sum3
#输出T
T=1+(n/sqrt(3))+((2/n)*sum3)-(sqrt(2)*sum2)