爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 12053|回复: 0

[求助] PYTHON的KS正态检验的问题

[复制链接]

新浪微博达人勋

发表于 2017-7-29 22:55:41 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
    小弟碰到一个KS正态检验的问题,拿PYTHON算月平均气温通不过正态检验,但是SPSS是能通过的,奇怪的很。
     在这之前我一直认为月季平均气温都是服从正态的,但是用PYTHON的KSTEST函数计算了,
并没有通过0.05信度的检验,具体如下图,是我什么地方算错了吗?代码中最后输出TXT的两列就是KS检验的序列
obs_region和obs_region1,将其代入SPSS算了竟然都是通过检验的,该信哪个,已经彻底晕了?是我的PYTHON语句
有问题吗,但其实KSTEST语句检验的也就是这两个变量序列的呀。(附代码和输入数据,1月和7月平均气温)

2.png PYTHON检验结果

SPSS输出

SPSS输出
SPSS检验结果


  1. import numpy as np
  2. from matplotlib import pyplot as plt
  3. from scipy.stats import kstest

  4. obs0 = np.loadtxt("7.txt", dtype=np.str)
  5. year= np.int_(obs0[11:67,0])    #读取年1961-2016
  6. obs= np.float64(np.array(obs0[11:67,1::]))   #读取数据
  7. obs[obs>30000]=np.nan
  8. obs_region=np.nanmean(obs, axis = 1)

  9. y1,y2=kstest(obs_region, 'norm')
  10. print y1,y2

  11. obs_region1=np.nanmean(obs[0:30,0:5], axis = 1)
  12. z1,z2=kstest(obs_region1, 'norm')
  13. print z1,z2
  14. np.savetxt("new.txt",obs_region,fmt='%10.3f')
  15. np.savetxt("new1.txt",obs_region1,fmt='%10.3f')
复制代码







1.txt

75.07 KB, 下载次数: 0, 下载积分: 金钱 -5

7.txt

79.35 KB, 下载次数: 0, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表