我想以对数-对数比例创建图,但是文件中的某些列包含-nun值,这会给我一条错误消息:
“ ValueError:数据没有正值,因此无法进行日志缩放。”
我的问题是如何解决这个问题?
请注意,我尝试不使用参数"nonposy="clip"
来解决它,但是它不起作用。
文件包含的示例
tau DDM DDMerr n_pairs
3.00 -nan -nan 0
3.75 -nan -nan 0
4.69 -nan -nan 0
5.86 4.997519e-09 7.082515e-09 1
7.32 3.943581e-07 9.197634e-09 59
9.16 1.466947e-06 3.455545e-07 1
代码:
Data = numpy.loadtxt(sf_file)
tau = Data[:,0]
DDM = Data[:,1]
DDM_err = Data[:,2]
n_pairs = Data[:,3]
#Kolmogorov
kolm= (tau)**(5./3.)
fig=plt.figure()
ax = plt.axes()
ax.set_xscale('log', nonposy="clip")
ax.set_yscale('log', nonposy="clip")
plt.errorbar(tau, DDM, DDM_err, label= "SF", marker='o', ms=7, color='magenta')
plt.plot(tau, kolm, color='g',linewidth=2,label='Kolmogorov')
plt.show()