MeteoInfoLab脚本示例：线性拟合

陈孝君

好贴，共享。

黑色青春

谢谢王老师的脚本

Linda.W

王老师：
      您好。我在画pm2.5与rh的线性拟合的时候，出现这个错误
Traceback (most recent call last):
  File "<iostream>", line 6, in <module>
ValueError: too many values to unpack
数据是只用了两列求相关的数据
RH        PM2.5
81.8        171.3
69.3        215.3
63.3        239.5
83        253.4
80        154
87.3        188.1
87.8        204.3
91.5        193.7
89.8        119.5
96.3        79.9
92        110.7
92.3        100.1
90.5        123.5
94.5        124.8
76.5        55.5
79.5        92.3
不知道是什么原因，希望老师可以指点迷津，感谢老师。

AX32

ldepn 发表于 2015-7-17 11:28
运行的时候出现一个错误
错误出现在：tdata = readtable(fn, delimiter=',', format='%{yyyyMMddhh}D%f%f% ...

no viable alternative at input 'label'，老师命令栏出现这种是什么原因呢。
fn=os.path.join('E:\MeteoInfo\sample\ASCIIL\skincancer.txt')
df=DataFrame.read_table(fn,format='%s%f%2i%f')
lat=df['lat'].values
mort=df['Mort'].values

slope,intercept,r,p,std_eer=stats.linregress(lat,mort)
scatter(lat,mort,label='original data',edge=False)
plot(lat,intercept+slope*lat,'r'label='fitted line')
text(29,100,r'$\hat{y}='+'%.2f'% slope+'x+'+\
    '%.1f'% intercept+'$',fontsize=16)
text(29,88,r'$R^2='+'%.4f'%(r**2)+'$',fontsize=16)
legend()
xlim(27,50)
ylim(75,250)
xticks(arange(30,51,5))
yticks(arange(100,226,25))
title('Skin Cancer Mortality versus State Latitude')
xlabel('Latitude(at center of State)'
ylabel('Mortality（Deaths per 10 million)')

MeteoInfo

AX32 发表于 2022-3-23 17:13
no viable alternative at input 'label'，老师命令栏出现这种是什么原因呢。
fn=os.path.join('E:\Met ...

错误信息里会提示是哪一行代码出错，你仔细检查出错的代码

AX32

谢谢老师解答

AX32

老师您好，按上述的代码运行完，没有报错，但是不出图是什么原因呢

MeteoInfo

AX32 发表于 2022-4-28 18:02
老师您好，按上述的代码运行完，没有报错，但是不出图是什么原因呢

应该是数据文件不存在。

更新的例子参考这里：http://www.meteothink.org/docs/m ... ic.stats.linregress