- 积分
- 693
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-6-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
小弟碰到一个KS正态检验的问题,拿PYTHON算月平均气温通不过正态检验,但是SPSS是能通过的,奇怪的很。
在这之前我一直认为月季平均气温都是服从正态的,但是用PYTHON的KSTEST函数计算了,
并没有通过0.05信度的检验,具体如下图,是我什么地方算错了吗?代码中最后输出TXT的两列就是KS检验的序列
obs_region和obs_region1,将其代入SPSS算了竟然都是通过检验的,该信哪个,已经彻底晕了?是我的PYTHON语句
有问题吗,但其实KSTEST语句检验的也就是这两个变量序列的呀。(附代码和输入数据,1月和7月平均气温)
PYTHON检验结果
SPSS输出
SPSS检验结果
- import numpy as np
- from matplotlib import pyplot as plt
- from scipy.stats import kstest
- obs0 = np.loadtxt("7.txt", dtype=np.str)
- year= np.int_(obs0[11:67,0]) #读取年1961-2016
- obs= np.float64(np.array(obs0[11:67,1::])) #读取数据
- obs[obs>30000]=np.nan
- obs_region=np.nanmean(obs, axis = 1)
- y1,y2=kstest(obs_region, 'norm')
- print y1,y2
- obs_region1=np.nanmean(obs[0:30,0:5], axis = 1)
- z1,z2=kstest(obs_region1, 'norm')
- print z1,z2
- np.savetxt("new.txt",obs_region,fmt='%10.3f')
- np.savetxt("new1.txt",obs_region1,fmt='%10.3f')
复制代码
|
-
-
1.txt
75.07 KB, 下载次数: 0, 下载积分: 金钱 -5
-
-
7.txt
79.35 KB, 下载次数: 0, 下载积分: 金钱 -5
|