Scipy randint 与 Numpy randint 的区别

主要区别是scipy.stats.randint可以定义lower/upper tail的概率,并且可以定义随机变量所服从的分布(可以定义概率质量函数PMF - Probability mass function)。
如果只需要生成在特定区间的随机整数,使用numpy.random.randint更加简单,随机数服从离散均匀分布。

scipy.stats.randint:
https://docs.scipy.org/doc/scipy-0.15.1/reference/generated/scipy.stats.randint.html
numpy.random.randint:
https://docs.scipy.org/doc/numpy/reference/generated/numpy.random.randint.html

Reference:
https://stackoverflow.com/questions/49036853/scipy-randint-vs-numpy-randint

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,我会尽力回答您的问题。 首先,使用 numpy 方法产生分布于[37, 55)之间的随机整数 50 个,并将它们转换成顺序统计量,绘制出对应的散点图和直方图的代码如下: ```python import numpy as np import matplotlib.pyplot as plt # 产生 50 个分布于[37, 55)之间的随机整数 np.random.seed(0) data = np.random.randint(37, 55, 50) # 将随机整数转换成顺序统计量 order_stat = np.argsort(data) # 绘制散点图和直方图 plt.figure(figsize=(12, 4)) plt.subplot(1, 2, 1) plt.scatter(np.arange(50), order_stat) plt.xlabel('Index') plt.ylabel('Order Statistics') plt.subplot(1, 2, 2) plt.hist(data, bins=10) plt.xlabel('Value') plt.ylabel('Frequency') plt.show() ``` 运行上述代码得到的散点图和直方图如下: ![numpy_random.png](https://img-blog.csdnimg.cn/20211011112732647.png) 由于是随机产生的数据,每次运行程序得到的图形可能会有所不同。 接下来,使用种子方法产生相同的随机数,观察数据以及直方图的效果。种子方法可以通过 `np.random.seed()` 函数来实现。下面是修改好的代码: ```python import numpy as np import matplotlib.pyplot as plt # 产生相同的随机数 np.random.seed(0) data1 = np.random.randint(37, 55, 50) # 再次产生相同的随机数 np.random.seed(0) data2 = np.random.randint(37, 55, 50) # 将随机整数转换成顺序统计量 order_stat = np.argsort(data1) # 绘制散点图和直方图 plt.figure(figsize=(12, 4)) plt.subplot(1, 2, 1) plt.scatter(np.arange(50), order_stat) plt.xlabel('Index') plt.ylabel('Order Statistics') plt.subplot(1, 2, 2) plt.hist(data1, bins=10, alpha=0.5, label='data1') plt.hist(data2, bins=10, alpha=0.5, label='data2') plt.xlabel('Value') plt.ylabel('Frequency') plt.legend() plt.show() ``` 运行上述代码得到的散点图和直方图如下: ![numpy_seed.png](https://img-blog.csdnimg.cn/20211011113243992.png) 可以看到,两次产生的随机数完全一样,因此得到的顺序统计量和直方图也完全一样。 最后,使用 scipy.stats 方法产生以上随机数,并绘制出对应的散点图和直方图。scipy.stats 提供了许多随机数生成函数,其中包括了产生分布于[37, 55)之间的随机整数的函数 randint()。下面是修改好的代码: ```python import numpy as np import scipy.stats as stats import matplotlib.pyplot as plt # 产生 50 个分布于[37, 55)之间的随机整数 np.random.seed(0) data = stats.randint.rvs(37, 55, size=50) # 将随机整数转换成顺序统计量 order_stat = np.argsort(data) # 绘制散点图和直方图 plt.figure(figsize=(12, 4)) plt.subplot(1, 2, 1) plt.scatter(np.arange(50), order_stat) plt.xlabel('Index') plt.ylabel('Order Statistics') plt.subplot(1, 2, 2) plt.hist(data, bins=10) plt.xlabel('Value') plt.ylabel('Frequency') plt.show() ``` 运行上述代码得到的散点图和直方图如下: ![scipy_random.png](https://img-blog.csdnimg.cn/20211011113851118.png) 如果要让 scipy.stats 每次产生相同的随机整数,可以使用种子方法,即在调用函数之前使用 `np.random.seed()` 函数。例如,可以将上述代码中的第 6 行改为: ```python np.random.seed(0) data = stats.randint.rvs(37, 55, size=50, random_state=0) ``` 其中的 `random_state` 参数就是用来指定随机数产生器的种子的。这样,每次运行程序得到的随机数就会完全一样。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值