在Python中读取文件并绘制CDF

Question

我需要以秒为单位读取带有时间戳的长文件，并使用numpy或scipy绘制CDF图。我确实尝试过numpy，但似乎输出不是应该的。下面的代码：任何建议表示赞赏。

import numpy as np
import matplotlib.pyplot as plt

data = np.loadtxt('Filename.txt')
sorted_data = np.sort(data)
cumulative = np.cumsum(sorted_data)

plt.plot(cumulative)
plt.show()

Answer 1

您有两种选择：

1：您可以先对数据进行装箱。使用以下numpy.histogram函数可以轻松完成此操作：

将numpy导入为np
导入matplotlib.pyplot作为plt

数据= np.loadtxt（'Filename.txt'）

＃在这里选择你想要多少个垃圾箱
num_bins = 20

＃使用直方图功能对数据进行分箱
计数，bin_edges = np.histogram（数据，bins = num_bins，normed = True）

＃现在找到CDF
cdf = np.cumsum（counts）

＃最后绘制CDF
plt.plot（bin_edges [1：]，cdf）

plt.show（）

2：而不是使用numpy.cumsum，只需sorted_data针对小于数组中每个元素的项目数绘制数组（请参阅此答案以获取更多详细信息https://stackoverflow.com/a/11692365/588071）：

将numpy导入为np

导入matplotlib.pyplot作为plt

数据= np.loadtxt（'Filename.txt'）

sorted_data = np.sort（数据）

yvals = np.arange（len（sorted_data））/ float（len（sorted_data）-1）

plt.plot（sorted_data，yvals）

plt.show（）