“Numpy数据分析与挖掘:高效学习重点技能“

目录

# 开篇

# 补充

zeros & ones

eye

1. numpy数组的创建

1.1 array

1.2 range

1.3 arange

1.4 常见的数据类型

1.5 astype

1.6 random.random() & round

2. numpy数组计算和数组计算

2.1 reshape

2.2 shape

2.3 将一维数组变成多维数组

2.4 指定一维数组成为多维数组

2.5 使用代表24行,仅1列

2.6 将数组变成一维数组

2.7 flatten

2.8 数组和数组计算

2.9 广播原则

2.9.1 数组t5 和数组t6 相加

2.9.2 数组t5 和数组t8 相减

2.10 NumPy 中轴( axis )的概念及其应用

2.11 transpose

2.12 swapaxes

3. numpy读取本地数据

3.1 如何使用 NumPy 读取 CSV 文件

3.2 示例:读取 YouTube 视频数据

3.3 loadtxt

4. numpy中的索引和切片

4.1 行切片

4.1.1 t1[2:]

4.1.2 t1[[2, 8, 10]]

4.1.3 取不连续的多行

4.2 列切片

4.2.1 t1[:, 0]

4.2.2 t1[:, 2:]

4.2.3 t1[:, [0, 2]]

4.2.4 t1[2, 3]

4.2.5 t1[2:5, 1:4]

4.2.6 t1[[0, 2], [0, 1]]

4.2.7 numpy高级索引方式

4.2.8 t1[[0, 2, 2], [0, 1, 3]]

5. numpy中更多的索引方式

5.1 t1[:, 2:4] = 0

5.2 t1 < 15000

5.3 t1[t1 < 15000] = 10000

5.4 numpy中的三元运算符

5.4.1 np.where(t1 < 15000, 0, 10)

5.5 numpy中的clip裁剪

5.5.1 clip

6. numpy数据的拼接

6.1 vstack & hstack

6.2 数组的行列交换

6.3 将两组数据拼接在一起进行科学计算

7. numoy中的随机方法

7.1 argmax

7.2 获取最小值的位置

7.3 np.zeros((3, 4))

7.4 np.ones((3, 4))

7.5 eye

7.6 案例

7.6.1 获取最大值的位置

7.6.2 获取最小值的位置

7.7 numpy生成随机数

7.8 生成二维均匀分布的随机数组

7.9 生成二维标准正态分布随机数组

7.10 生成二维整数随机数组

7.11 生成二维均匀分布的随机数组

7.12 生成二维正态分布随机数组

7.13 设置随机数种子

8. numpy中的nan和inf

8.1 np.nan 和np.inf

8.2 count_nonzero & isnan

9. numpy常用统计方法

9.1 sum

9.2 nan与任何值计算都是nan

9.3 numpy中的nan的注意点

9.4 axis 参数

9.5 sum 求和

9.6 mean 均值

9.7 median 中值

9.8 max 最大值

9.9 min 最小值

9.10 ptp 极值范围

9.11 std 标准差

10. numpy中填充nan和youtube数据的练习

10.1 需求:统计行均值时,出现nan,如何处理nan的值?

10.1.1 模拟nan数据

10.1.2 实现

10.1.3 分析

10.2 需求:绘制出youybe的数据,制出评论数量的直方图

10.2.1 数据获取

10.2.2 直方图


# 开篇

参考:

  1. numpy:用来处理数值型的数据;

  2. pandas:用来处理数值,以及字符串,字典、列表等数据的;

# 补充

zeros & ones

np.zeros((t2.shape[0], 1)) & ones_data = np.ones((t2.shape[0], 1)) :构造全为0的数据,构造全为1的数据:

  • zeros_data = np.zeros((t2.shape[0], 1)) # 创建一个3行1列的数组,全为0
  • ones_data = np.ones((t2.shape[0], 1)) # 创建一个3行1列的数组,全为1
t2 = np.arange(12, 24).reshape(3, 4)
print(t2)
print("*" * 100)
zeros_data = np.zeros((t2.shape[0], 1))  # 创建一个3行1列的数组,全为0
print(zeros_data)
print("*" * 100)
ones_data = np.ones((t2.shape[0], 1))  # 创建一个3行1列的数组,全为1
print(ones_data)

eye

np.eye(3) # 创建一个3行3列的对角矩阵:

t2 = np.eye(3)  # 创建一个3行3列的对角矩阵
print(t2)

1. numpy数组的创建

1.1 array

np.array([1,2,3]) :创建numpy数组

import numpy as np

t1 = np.array([1, 2, 3])
print(t1)
print(type(t1))

1.2 range

range(10) :生成0-9的数组

# 默认生成0到9,十个数字
t2 = np.array(range(10))
print(t2)

1.3 arange

np.arange(4, 10, 2) :取4-10区间,步长为2

t3 = np.arange(4, 10, 2)
print(t3)

np.array(range(1, 4), dtype=float) :指定返回的类型为浮点数

t4 = np.array(range(1, 4), dtype=float)
print(t4)
print(t4.dtype)

1.4 常见的数据类型

类型类型代码说明
int8, uint8i1, u1有符号和无符号的8位(1个字节)整数
int16, uint16i2, u2有符号和无符号的16位(2个字节)整数
int32, uint32i4, u4有符号和无符号的32位(4个字节)整数
int64, uint64i8, u8有符号和无符号的64位(8个字节)整数
float16f2半精度浮点数
float32f4或f标准的单精度浮点数。与C的float兼容
float64f8或d标准的双精度浮点数。与C的double和Python的float对象兼容
float128f16或g扩展精度浮点数
complex64, complex128c8, c16分别用两个32位、64位或128位浮点数表示的复数
complex256c32
bool?存储True和False值的布尔类型

1.5 astype

np.astype("int8") :调整数据类型

# numpy中的bool类型
t5 = np.array([1, 1, 0, 1, 0, 1], dtype=bool)
print(t5)
print(t5.dtype)

# 调整数据类型
t6 = t5.astype("int8")
print(t6)
print(t6.dtype)

1.6 random.random() & round

random.random() & np.round(t7, 2) :随机生成小数 & 保留后两位小数

# 随机生成10个小数
t7 = np.array([random.random() for i in range(10)])
print(t7)
print(t7.dtype)

# 执行后两位小数
t8 = np.round(t7, 2)
print("*" * 100)
print(t8)

2. numpy数组计算和数组计算

2.1 reshape

np.arange(100, 124).reshape(4, 6) :默认生成100, 124区间数字,形成一个4行6列的数组

t6 = np.arange(100, 124).reshape(4, 6)
print(t6)

2.2 shape

np.shape :返回数组的元组

t1 = np.arange(12)
print(t1)
print(t1.shape)  # (12,) 用元组返回数组的行数
print("*" * 100)

t2 = np.array([[1, 2, 3], [4, 5, 6]])
print(t2)
print(t2.shape)  # (2, 3) 用元组返回数组的行数2,返回3数组的列数
print("*" * 100)

t3 = np.array([[[1, 2], [3, 4]], [[5, 6], [7, 8]]])
print(t3)
print(t3.shape)  # (2, 2, 2) 第一个2代表数组有两个块,每块有2行,每行有2列

2.3 将一维数组变成多维数组

np.reshape(3, 4) :将一维数组变成多维数组

t4 = np.arange(12)
print(t4)
print("*" * 100)
print(t4.reshape((3, 4)))  # 将t4一维数组变成3行4列的数组
print("*" * 100)
print(t4.reshape((3, 5)))  # 3行5列的数组,报错

2.4 指定一维数组成为多维数组

np.arange(24).reshape((2, 3, 4)) :指定一维数组成为多维数组

# 指定一维数组成为多维数组
t5 = np.arange(24).reshape((2, 3, 4))  # 2个块,每块3行,每行4列
print(t5)

注:有返回值,但是不会对自身数据进行修改;

可以使用:变量 = 变量 方式赋值;

# 指定一维数组成为多维数组
t5 = np.arange(24).reshape((2, 3, 4))  # 2个块,每块3行,每行4列
print(t5)
print("*" * 100)
print(t5.reshape((4, 6)))
print("*" * 100)
print(t5)  # 改变数组形状,不会改变数组内容

2.5 使用代表24行,仅1列

np.reshape(24, 1) :使用代表24行,仅有1列

注意:np.reshape(x, xx) 前面参数代表行,后面参数代表列;

t5 = np.arange(24)
print(t5.reshape(24, 1))  # 24行1列

2.6 将数组变成一维数组

np.reshape((t5.shape[0] * t5.shape[1],)) :将数组变成一维数组

t5.shape[0] 代表行;

t5.shape[1] 代表列;

相乘就是一个首元素等于24元组;

t5 = np.arange(24).reshape((4, 6))
print(t5)
print("*" * 100)
t6 = t5.reshape((t5.shape[0] * t5.shape[1],))  # 将数组变成一维数组
print(t6)

2.7 flatten

np.flatten() :快速的帮我们把数组变成一维数组;

t5 = np.arange(24).reshape((4, 6))
print(t5)
print("*" * 100)
print(t5.flatten())  # 将数组变成一维数组

2.8 数组和数组计算

np+2 :2加到数组中每个元素上

就是当numpy用于与树脂进行计算的时候,这个计算会被用于每一个数字;

t5 = np.arange(24).reshape((4, 6))
print(t5)
print("*" * 100)
print(t5 + 2)  # 2加到数组中每个元素上
print("*" * 100)
print(t5 / 2)  # 2除到数组中每个元素上
print("*" * 100)
print(t5 / 0)  # 0除到数组中每个元素上

2.9 广播原则

问题:

  • 数组形状为 (3, 3, 3) 能否与形状为 (3, 2) 的数组进行计算?不可以
  • 数组形状为 (3, 3, 2) 能否与形状为 (3, 2) 的数组进行计算?可以

关键点:

  • 维度对齐:可以将维度指的是 shape 所对应的数字个数。
  • 例子:每列的数据减去每列的平均值的结果。

理解:

  • 广播机制:NumPy 的广播机制允许形状不同的数组进行某些操作,但要求在维度和形状上能够相互兼容。一般情况下,计算中维度要么完全匹配,要么其中一个数组的维度为 1,以便进行广播。

举例:

  • 形状为 (3, 3, 3) 的数组与形状为 (3, 2) 的数组进行操作时,必须要满足特定的维度对齐规则。
  • 在这个例子中,可以考虑每列的数据减去每列的平均值的操作,即通过对齐和广播,使得形状不同的数组能够进行兼容操作。

总结:

  • 检查数组的形状是否兼容,利用广播机制处理不同形状的数组。
  • 实际操作时,需要考虑每个维度是否可以通过扩展或对齐使得操作合法。
2.9.1 数组t5 和数组t6 相加
t5 = np.arange(24).reshape((4, 6))
t6 = np.arange(100, 124).reshape((4, 6))
print(t5)
print("*" * 100)
print(t6)
print("*" * 100)
print(t5 + t6)  # 两个数组相加

2.9.2 数组t5 和数组t8 相减

在数组运算中,不同的列,但是相同的行,是可以进行运算的;

t5 = np.arange(24).reshape((4, 6))
print(t5)
print("*" * 100)
t8 = np.arange(4).reshape(4, 1)
print(t8)
print("*" * 100)
print(t5 - t8)  # 两个数组相减, 会在每一列上都减去相同维度的数组

t5 = np.arange(24).reshape((4, 6))
print(t5)
print("*" * 100)
t9 = np.arange(10)
print(t9)
print("*" * 100)
print(t5 - t9)  # 在数组元素对应不上的时候,是无法进行运算

在元素对应不上的时候,是无法进行运算;

2.10 NumPy 中轴( axis )的概念及其应用

轴(axis)

在 NumPy 中,轴可以理解为方向,使用 0, 1, 2, ... 数字表示:

  • 对于一维数组,只有一个 0 轴。
  • 对于二维数组(如 shape=(2,2)),有 0 轴和 1 轴。
  • 对于三维数组(如 shape=(2,2,3)),有 0, 1, 2 轴。

计算应用

有了轴的概念后,计算会更加方便。例如:

  • 计算一个二维数组的平均值,必须指定是计算哪个方向上的数字的平均值。

实例问题

在前面的知识中,轴在哪里?

import numpy as np
np.arange(0, 10).reshape((2, 5))
  • reshape 中,2 表示 0 轴长度(包含数据的条数),5 表示 1 轴长度(包含数据的条数)。因此,数组变成 2x5,共 10 个数据。

关键点总结

  • 轴的编号:在 NumPy 中,从 0 开始编号,0 轴表示第一个维度,以此类推。
  • 应用场景:在进行计算(如求和、平均等)时,需要明确指定计算的轴。
  • 多维数组:多维数组中的轴编号根据数组的维度递增。

2.11 transpose

np.transpose(t2) :数组转置

行列转置,将行转置为列;

t2 = np.arange(12).reshape(3, 4)
print(t2)
print("*" * 100)
t3 = np.transpose(t2)  # 转置
print(t3)

2.12 swapaxes

np.swapaxes(x, 0, 1) :0, 1 表示交换第一维和第二维

x = np.array([[1, 2, 3]])
print(x)
print("*" * 100)
y = np.swapaxes(x, 0, 1)  # 0, 1 表示交换第一维和第二维
print(y)

np.swapaxes(x, 0, 2) :0, 2 表示交换第一维和第三维

x = np.array([[[0, 1], [2, 3]], [[4, 5], [6, 7]]])
print(x)
print("*" * 100)
y = np.swapaxes(x, 0, 2)  # 0, 2 表示交换第一维和第三维
print(y)

3. numpy读取本地数据

预备知识:

热门YouTube视频统计和评论:Trending YouTube Video Statistics and Comments (kaggle.com)

  • 数据来源:

    Kaggle 是一个国际知名的数据科学竞赛平台

    数据来自 Kaggle:https://www.kaggle.com/datasnaek/youtube/data

    • 处理后的数据:

3.1 如何使用 NumPy 读取 CSV 文件

NumPy 提供了 np.loadtxt 函数用于读取 CSV 文件。它是一个多功能工具,可以读取文本文件中的数据,并将其转换为 NumPy 数组。

使用 np.loadtxt 的基本语法:

np.loadtxt(
    frame,
    dtype=float,
    delimiter=',',
    skiprows=0,
    usecols=None,
    unpack=False
)
参数解释示例
fname文件名或文件路径。可以是字符串、生成器、.gz 或 .bz2 压缩文件'data.csv'
dtype数据类型,决定读取的数据类型。默认是 np.floatdtype=int
delimiter分隔符,用于分隔 CSV 文件中的列。默认是空白字符delimiter=','
skiprows跳过文件开头的几行。通常用来跳过表头。skiprows=1
usecols指定要读取的列。可以是索引或元组类型。usecols=(0, 1, 2, 3)
unpack如果为 True,数据将会分解到不同的数组变量中。为 False 时,数据会读入一个数组unpack=True

3.2 示例:读取 YouTube 视频数据

假设你有一个包含 YouTube 视频点击、喜欢、不喜欢、评论数量的 CSV 文件 youtube_data.csv,其数据如下:

views,likes,dislikes,comment_count
1234567,12345,234,5678
2345678,23456,345,6789
...

可以按照以下步骤读取并处理这个数据:

import numpy as np

# 读取CSV数据
data = np.loadtxt(
    'youtube_data.csv',  # 文件名
    dtype=float,         # 数据类型
    delimiter=',',       # 分隔符
    skiprows=1,          # 跳过第一行(表头)
    usecols=(0, 1, 2, 3) # 读取指定的列
)

# 分别获取不同列的数据
views, likes, dislikes, comment_count = data.T  # 转置使数据分别赋值给变量

# 打印数据
print("Views:", views)
print("Likes:", likes)
print("Dislikes:", dislikes)
print("Comment Count:", comment_count)

重要提示: 

  1. 数据来源:要确保你的 CSV 文件路径正确,并且格式与代码中的 delimiter 一致。
  2. 数据类型:根据你的数据情况选择合适的 dtype
  3. 列选择:如果你只需要特定列,可以使用 usecols 进行选择。

3.3 loadtxt

np.loadtxt() :文件读取

# coding: utf-8
import numpy as np

us_file_path = "../youyube_video_data/GBvideos-numpy.csv"

t1 = np.loadtxt(us_file_path, delimiter=',', dtype="int", encoding='utf-8')
print(t1)

4. numpy中的索引和切片

4.1 行切片

4.1.1 t1[2:]

t1[2:] :从第二行开始到最后一行进行切片

print(t1[2:]) # 2:表示从第3行开始

4.1.2 t1[[2, 8, 10]]

t1[[2, 8, 10]] :指定行切片

print(t1[[2, 8, 10]]) # 2,8,10表示第3,9,11行

4.1.3 取不连续的多行

t1[1, :] & t1[2:, :] & t1[[2, 10, 3], :] :取不连续的多行

  • t1[x, :] :x代表行,在逗号后面, 代表着列,x, : 代表返回所有的列;
    1. t1[1, :] :指定行取值,并展示所有的列;
    2. t1[2:, :] :从第三行开始取值,并展示所有的列;
    3. t1[[2, 10, 3], :] :指定行取值,并展示所有的列;
print(t1)
print("*" * 100)
print(t1[1, :])  # 1表示第2行
print("*" * 100)
print(t1[2:, :])  # 2:表示从第3行开始
print("*" * 100)
print(t1[[2, 10, 3], :])  # 2,10,3表示第3,11,4行

4.2 列切片

4.2.1 t1[:, 0]

t1[:, 0] :指定第0列切片

print(t1[:, 0]) # 0表示第1列

4.2.2 t1[:, 2:]

t1[:, 2:] :取连续的多列

print(t1[:, 2:]) # 2:表示从第3列开始

4.2.3 t1[:, [0, 2]]

t1[:, [0, 2]] :指定第0,2列切片

print(t1[:, [0, 2]]) # 0,2表示第1,3列

4.2.4 t1[2, 3]

t2 = t1[2, 3] :指定行、指定列切片

print(t1)
print("*" * 100)
t2 = t1[2, 3]  # 2表示第3行,3表示第4列
print(t2)
print(type(t2))

4.2.5 t1[2:5, 1:4]

t2 = t1[2:5, 1:4] :指定多行、指定多列切片

取第3行到第5行,第2列到第4列的结果;

print(t1)
print("*" * 100)
t2 = t1[2:5, 1:4]  # 2:5表示从第3行开始,1:4表示从第2列开始
print(t2)

4.2.6 t1[[0, 2], [0, 1]]

t2 = t1[[0, 2], [0, 1]] :高级索引取多个不同相邻的点

print(t1)
print("*" * 100)
t2 = t1[[0, 2], [0, 1]]
print(t2)

4.2.7 numpy高级索引方式

        在 t2 = t1[[0, 2], [0, 1]] 中,[0, 1] 表示索引的位置,并不是代表具体的列,而是用于与 [0, 2] 组合来定位具体的元素。具体来说,这里 numpy 是采用了高级索引方式来处理这个表达式。

numpy** 的高级索引会将索引对 (i, j) 视为独立的一对对待,而不是作为行列索引的组合。**

让我们详细解释这段代码的含义:

t2 = t1[[0, 2], [0, 1]]
  • t1[[0, 2], [0, 1]] 表示从 t1 中获取索引 (0, 0)(2, 1) 的元素。
  • t1[[0, 2]] 选取了第1行和第3行:即 [7426393, 78240, 13548, 705][142819, 13119, 151, 1141]
  • t1[[0, 2], [0, 1]] 会将上述行中对应索引的元素取出,即 (0, 0)(2, 1) 的元素。

简而言之,这段代码是取 t1(0, 0) 位置的元素和 (2, 1) 位置的元素。因此,结果是 [7426393, 13119]

4.2.8 t1[[0, 2, 2], [0, 1, 3]]

t2 = t1[[0, 2, 2], [0, 1, 3]] :高级索引取多个不同相邻的点

print(t1)
print("*" * 100)
t2 = t1[[0, 2, 2], [0, 1, 3]]
print(t2)

5. numpy中更多的索引方式

5.1 t1[:, 2:4] = 0

t1[:, 2:4] = 0 :将所有行的3到4列都赋值为0

print(t1)
print("*" * 100)
t1[:, 2:4] = 0  # 2:4表示所有行从第3列开始都是0
print(t1)

5.2 t1 < 15000

t1 < 15000 :返回数组中小于15000的布尔值

print(t1)
print("*" * 100)
print(t1 < 15000)  # 返回小于15000的布尔值

5.3 t1[t1 < 15000] = 10000

t1[t1 < 15000] = 10000 :小于15000的赋值为10000

print(t1)
print("*" * 100)
t1[t1 < 15000] = 10000  # 小于15000的赋值为10000
print(t1)

5.4 numpy中的三元运算符

5.4.1 np.where(t1 < 15000, 0, 10)

t2 = np.where(t1 < 15000, 0, 10) :小于15000的赋值为0,大于等于15000的赋值为10

(numpy的三元运算符)

print(t1)
print("*" * 100)
t2 = np.where(t1 < 15000, 0, 10)  # 小于15000的赋值为0,大于等于15000的赋值为10
print(t2)

5.5 numpy中的clip裁剪

5.5.1 clip

t2 = t1.clip(15000, 50000) :小于15000的赋值为15000,大于等于50000的赋值为50000

print(t1)
print("*" * 100)
t2 = t1.clip(15000, 50000)  # 小于15000的赋值为15000,大于等于50000的赋值为50000
print(t2)

6. numpy数据的拼接

  • 数组的拼装:

    • 注意:

      竖直拼接的时候:每一列代表的意义都必须相同;

      如果每一列的意义不同,这时候应该交换某一组数的列,让其和另外一类相同;

    np.vstack((t1,t2)) :竖直拼接;将第二个数组拼接在第一个数组的最后行后面

    np.hstack((t1,t2)) :水平拼接;将第二个数组的列拼接到一个数组每一个最后列后面

6.1 vstack & hstack

np.vstack((t1,t2)) & np.hstack((t1,t2))

t2 = np.arange(12).reshape(2, 6)
print(t2)
print("*" * 100)
t3 = np.arange(12, 24).reshape(2, 6)
print(t3)
print("*" * 100)
t4 = np.vstack((t2, t3))  # 纵向堆叠
print(t4)
print("*" * 100)
t5 = np.hstack((t2, t3))  # 横向堆叠
print(t5)

6.2 数组的行列交换

  • 数组的行列交换:

    t[[1, 2], :] = t[[2, 1], :] # 行1和2交换

    t[:, [0, 2]] = t[:, [2, 0]] # 列0和2交换

  • 示例:
    t = np.arange(12, 24).reshape(3, 4)
    print(t)
    print("*" * 100)
    
    t[[1, 2], :] = t[[2, 1], :]  # 行1和2交换
    print(t)
    print("*" * 100)
    
    t[:, [0, 2]] = t[:, [2, 0]]  # 列0和2交换
    print(t)
    

6.3 将两组数据拼接在一起进行科学计算

t2 = np.arange(12, 24).reshape(3, 4)
print(t2)
print("*" * 100)
zeros_data = np.zeros((t2.shape[0], 1))  # 创建一个3行1列的数组,全为0
print(zeros_data)
print("*" * 100)
ones_data = np.ones((t2.shape[0], 1))  # 创建一个3行1列的数组,全为1
print(ones_data)
print("*" * 100)

us = np.hstack((t2, zeros_data))  # 分别添加0和1的列
uk = np.hstack((t2, ones_data))  # 分别添加0和1的列
print(us)
print("*" * 100)
print(uk)
print("*" * 100)

final_data = np.vstack((us, uk))
print(final_data)

7. numoy中的随机方法

7.1 argmax

np.argmax(t, axis=0) :获取最大值的位置

  • np.argmax 用于返回数组中最大值的索引。
  • axis=0 表示沿数组的第0轴(即列)进行操作。
  • 获取数组 t 沿指定轴的最大值位置:
np.argmax(t, axis=0)

7.2 获取最小值的位置

np.argmin(t, axis=1) :获取最小值的位置

  • np.argmin 用于返回数组中最小值的索引。
  • axis=1 表示沿数组的第1轴(即行)进行操作。
  • 获取数组 t 沿指定轴的最小值位置:
np.argmin(t, axis=1)

7.3 np.zeros((3, 4))

np.zeros((3, 4)) :创建一个全0的数组

  • np.zeros 创建一个指定形状的全0数组。
  • 这里创建的数组形状为 (3, 4),即3行4列的数组。
  • 创建一个形状为 (3, 4) 的全0数组:
np.zeros((3, 4))

7.4 np.ones((3, 4))

np.ones((3, 4)) :创建一个全1的数组

  • np.ones 创建一个指定形状的全1数组。
  • 这里创建的数组形状为 (3, 4),即3行4列的数组。
  • 创建一个形状为 (3, 4) 的全1数组:
np.ones((3, 4))

7.5 eye

np.eye(3):创建一个对角线为1的正方形数组(方阵)

  • np.eye 创建一个单位矩阵,即对角线为1,其他元素为0的矩阵。
  • 这里创建的单位矩阵是一个3x3的方阵。
  • 创建一个形状为 (3, 3) 的单位矩阵:
np.eye(3)

7.6 案例

7.6.1 获取最大值的位置

np.argmax(t, axis=0) :获取最大值的位置

t = np.eye(4)
print(t)
print("*" * 100)
print(np.argmax(t, axis=0))

7.6.2 获取最小值的位置

np.argmax(t, axis=1) :获取最小值的位置

t = np.eye(4)
print(t)
print("*" * 100)
print(np.argmax(t, axis=0))  # 求最大下标的索引值

t[t == 1] = -1  # 将1替换为-1
print("*" * 100)
print(t)
print("*" * 100)
print(np.argmin(t, axis=1))  # 求最小下标的索引值

7.7 numpy生成随机数

  • 分布的补充:

    1.均匀分布:在相同的大小范围内的出现概率是等可能的;

    2.正态分布:呈钟型,两头低,中间高,左右对称;

7.8 生成二维均匀分布的随机数组

np.random.rand(2, 3) :生成二维均匀分布的随机数组

函数: np.random.rand(d0, d1, ..., dn)

参数描述: d0, d1, ..., dn 维度的均匀分布的随机数组;

解释:创建 d0dn 维度的浮点数随机数组,范围从 0 到 1;

import numpy as np

array_rand = np.random.rand(2, 3)
print(array_rand)
# 输出示例:[[0.17678197 0.83879827 0.59548066]
#           [0.00811262 0.05331716 0.16584186]]

7.9 生成二维标准正态分布随机数组

np.random.randn(2, 3) :生成二维标准正态分布随机数组

函数:np.random.randn(d0, d1, ..., dn)

参数描述: d0, d1, ..., dn 维度的标准正态分布随机数组

解释:创建 d0dn 维度的标准正态分布随机数组,均值 0,标准差 1;

array_randn = np.random.randn(2, 3)
print(array_randn)
# 输出示例:[[-0.56922926 -0.46477167  1.53166647]
#           [ 0.51248294  1.61693618 -1.09743841]]

7.10 生成二维整数随机数组

np.random.randint(0, 10, (2, 3)) :生成二维整数随机数组

函数:np.random.randint(low, high, size)

参数描述: low: 起始值、 high: 结束值、 size: 形状

解释:从 lowhigh 之间生成随机整数,形状由 size 指定;

array_randint = np.random.randint(0, 10, (2, 3))
print(array_randint)
# 输出示例:[[7 4 2]
#           [0 8 4]]

7.11 生成二维均匀分布的随机数组

np.random.uniform(0, 1, (2, 3)) :生成二维均匀分布的随机数组

函数:np.random.uniform(low, high, size)

参数描述: low: 起始值、 high: 结束值、 size: 形状

解释:产生具有均匀分布的随机数数组,范围从 lowhigh

array_uniform = np.random.uniform(0, 1, (2, 3))
print(array_uniform)
# 输出示例:[[0.38019132 0.8574209  0.34507822]
#           [0.28912108 0.95612227 0.33245865]]

7.12 生成二维正态分布随机数组

np.random.normal(0, 1, (2, 3)) :生成二维正态分布随机数组

函数:np.random.normal(loc, scale, size)

参数描述:loc: 均值、 scale: 标准差、 size: 形状

解释:从正态分布中抽取随机样本,均值为 loc,标准差为 scale

array_normal = np.random.normal(0, 1, (2, 3))
print(array_normal)
# 输出示例:[[ 0.32859216 -0.89933692  0.53515463]
#           [ 0.89322116 -0.13982256 -1.05578878]]

7.13 设置随机数种子

np.random.seed(42) :设置随机数种子

函数:np.random.seed(s)

参数描述:s: 随机数种子

解释:设置随机数种子,可以确保生成的随机数相同;

np.random.seed(42)
array_seed = np.random.rand(2, 3)
print(array_seed)
# 输出示例:[[0.37454012 0.95071431 0.73199394]
#           [0.59865848 0.15601864 0.15599452]]

8. numpy中的nan和inf

8.1 np.nan 和np.inf

  • NaN 用于表示一个未定义或不可表示的值:
    • 从文件读取数据时,如果有缺失值,将其表示为 NaN
    • 进行不合理的计算时,例如 0 除以 0,结果会是 NaN
import numpy as np

nan_value = np.nan
print(nan_value)  # 输出: nan
print(type(nan_value))  # 输出: <class 'float'>
  • Inf 表示正无穷大,-Inf 表示负无穷大。它们通常出现在以下情况中:

    • 进行无穷大的计算时,例如除数为 0。
    • 超过数值范围时,例如超出浮点数表示范围的结果。

    生成 Inf-Inf

import numpy as np

inf_value = np.inf
neg_inf_value = -np.inf
print(inf_value)  # 输出: inf
print(type(inf_value))  # 输出: <class 'float'>

print(neg_inf_value)  # 输出: -inf
print(type(neg_inf_value))  # 输出: <class 'float'>
  • 使用 NaNInf 的案例:
  1. 读取文件中的缺失值
    import numpy as np
    
    data = np.genfromtxt('data.csv', delimiter=',')
    # 如果文件中某些数据缺失,它们会被读取为 NaN
    
  2. 处理无效计算
    result = np.log(-1)
    print(result)  # 输出: nan
    
  3. 处理无穷大
    result = 1.0 / 0.0
    print(result)  # 输出: inf
    

8.2 count_nonzero & isnan

np.nan 和np.inf 的注意点 & np.count_nonzero(t != t) & np.isnan(t)

两个nan是不相等的:np.nan == np.nan # False

  1. np.nan != np.nan:np.nan != np,nan # True
  2. nan和任何值计算都为nan;

可以利用以上的特性,判断数组中nan的个数:

t = np.array([1., 2., np.nan, np.nan])
print(t)
print("*" * 100)
t2 = np.count_nonzero(t != t)  # 判断t中为nan的个数
print(t2)

由于2,那么如何判断一个数字是否为nan呢:

  • 通过np.isnan(a) 来判断,返回bool类型,比如希望把nan替换为0;
t = np.array([1., 2., np.nan, np.nan])
print(t)
print("*" * 100)
t[np.isnan(t)] = 0  # 将nan替换为0
print(t)

np.isnan(t)

9. numpy常用统计方法

  • 标准差:当值越小,数值越聚集,当值越大,数值越离散;

9.1 sum

np.sum(t3, axis=0) :统计行或列的结果

t3 = np.arange(12).reshape(3, 4)
print(t3)
print("*" * 100)
t4 = np.sum(t3, axis=0)  # 统计每行全部列的结果
print(t4)
print("*" * 100)
t5 = np.sum(t3, axis=1)  # 统计每行全部行的结果
print(t5)

9.2 nan与任何值计算都是nan

t3 = np.arange(12, dtype=float).reshape(3, 4)
print(t3)
print("*" * 100)

t3[2, 3] = np.nan  # 给浮点数组赋值 NaN
print(t3)
print("*" * 100)
print(np.sum(t3, axis=0))
print("*" * 100)
print(np.nansum(t3, axis=0))  # 使用 np.nansum 忽略 NaN 进行求和

9.3 numpy中的nan的注意点

nan处理会存在的情况:

在计算统计中,如果把nan的值替换为0进行统计,这样并不合适,替换之前的平均值如果大于0,替换之后均值肯定会变小;

所以一般的方式是把缺失的数值替换为均值(中值)或者是直接删除缺失值的一行:

9.4 axis 参数

axis=None: 默认值,对整个数组进行操作。

axis=0: 沿垂直方向(列)操作。

axis=1: 沿水平方向(行)操作。

9.5 sum 求和

求和: np.sum(t3, axis=None)

  • 计算数组元素的总和。
  • axis=None 表示对整个数组进行求和;指定轴(0 或 1)时对特定轴进行求和。

9.6 mean 均值

均值: np.mean(t3, axis=None)

  • 计算数组元素的均值(平均值)。
  • axis=None 表示对整个数组计算均值;指定轴(0 或 1)时对特定轴计算均值。
t3 = np.arange(12, dtype=float).reshape(3, 4)
print(t3)
print("*" * 100)
print(np.mean(t3, axis=None))  # 统计全部元素的均值

9.7 median 中值

中值: np.median(t, axis=None)

  • 计算数组元素的中位数。
  • axis=None 表示对整个数组计算中位数;指定轴(0 或 1)时对特定轴计算中位数。
t3 = np.arange(12, dtype=float).reshape(3, 4)
print(t3)
print("*" * 100)
print(np.median(t3, axis=None))  # 统计全部元素的中值

9.8 max 最大值

最大值: np.max(t3, axis=None)

  • 返回数组中的最大值。
  • axis=None 表示返回整个数组中的最大值;指定轴(0 或 1)时对特定轴返回最大值。
t3 = np.arange(12, dtype=float).reshape(3, 4)
print(t3)
print("*" * 100)
print(np.max(t3, axis=None))  # 统计全部元素的最大值

9.9 min 最小值

最小值: np.min(t3, axis=None)

  • 返回数组中的最小值。
  • axis=None 表示返回整个数组中的最小值;指定轴(0 或 1)时对特定轴返回最小值。
t3 = np.arange(12, dtype=float).reshape(3, 4)
print(t3)
print("*" * 100)
print(np.min(t3, axis=None))  # 统计全部元素的最小值

9.10 ptp 极值范围

极值范围: np.ptp(t3, axis=None)

极值 是指数组中最大值与最小值之间的差值,即 max - min。极值的计算函数是 numpy.ptp()(ptp 意为 "peak to peak")。

公式:

函数: t.ptp(axis=None)

  • 计算数组中最大值和最小值之差。
  • axis=None 表示对整个数组计算极值范围;指定轴(0 或 1)时对特定轴计算极值范围。
t3 = np.arange(12, dtype=float).reshape(3, 4)
print(t3)
print("*" * 100)
print(np.ptp(t3, axis=None))  # 统计全部元素的极值范围

9.11 std 标准差

标准差: np.std(t3, axis=None)

  • 标准差用于衡量数据集中的值偏离均值的程度。标准差越大,表示数据点离均值越远;标准差越小,表示数据点更集中于均值附近。
  • 计算数组元素的标准差。
  • axis=None 表示对整个数组计算标准差;指定轴(0 或 1)时对特定轴计算标准差。
t3 = np.arange(12, dtype=float).reshape(3, 4)
print(t3)
print("*" * 100)
print(np.std(t3, axis=None))  # 统计全部元素的标准差

10. numpy中填充nan和youtube数据的练习

10.1 需求:统计行均值时,出现nan,如何处理nan的值?

将统计除nan列的均值,再将均值赋值给当前列的nan值上;

10.1.1 模拟nan数据
t2 = np.arange(12).reshape((3, 4)).astype("float")
print(t2)
print("*" * 100)
t2[1, 2:] = np.nan  # 给浮点数组赋值 NaN
print(t2)

10.1.2 实现
def fill_ndarray(t2):
    for i in range(t2.shape[1]):  # 遍历每一列
        temp_col = t2[:, i]  # 当前列
        nan_num = np.count_nonzero(temp_col != temp_col)  # 当前列中NaN的个数
        if nan_num != 0:
            temp_not_nan_col = temp_col[temp_col == temp_col]  # 当前列中非NaN的列
            temp_col[np.isnan(temp_col)] = temp_not_nan_col.mean()  # 将NaN替换为当前列中非NaN的平均值
    return t2


if __name__ == '__main__':
    t2 = np.arange(12).reshape((3, 4)).astype("float")
    t2[1, 2:] = np.nan  # 给浮点数组赋值 NaN
    print(t2)
    print("*" * 100)
    t3 = fill_ndarray(t2)
    print(t3)

10.1.3 分析

统计存在有nan的列,将nan的列中其它不等于nan的值进行mean求出均值,

在把均值赋值给当前列中所有等于nan的值,

再return返回给t2;

10.2 需求:绘制出youybe的数据,制出评论数量的直方图

10.2.1 数据获取
import numpy as np

us_file_path = "../youyube_video_data/GBvideos-numpy.csv"
t1 = np.loadtxt(us_file_path, delimiter=',', dtype="int", encoding='utf-8')
print(t1)
print("*" * 100)
# 取评论的数据
t_us_comments = t1[:-1]
print(t_us_comments.max(), t_us_comments.min())
print("*" * 100)

# 选择比500万评论数小的数据
t_us_comments = t_us_comments[t_us_comments < 5000000]
print(t_us_comments.max(), t_us_comments.min())

10.2.2 直方图

import numpy as np
from pyecharts.charts import Bar
from pyecharts import options as opts

us_file_path = "../youyube_video_data/GBvideos-numpy.csv"
t1 = np.loadtxt(us_file_path, delimiter=',', dtype="int", encoding='utf-8')
t_us_comments = t1[:-1]  # 获取最后一列评论数

# 过滤掉评论数小于等于0的数据
t_us_comments = t_us_comments[t_us_comments > 0]

# 过滤掉评论数超过5万的数据
t_us_comments = t_us_comments[t_us_comments < 50000]

# 数据分组:选择合适的分组区间
bins = np.linspace(t_us_comments.min(), t_us_comments.max(), 100)
hist, bin_edges = np.histogram(t_us_comments, bins=bins)  # 计算直方图

bar = Bar()

# 设置 x 轴和 y 轴的值
# 构造 x 轴的数据,按照固定间隔
x_labels = [f"{i * 1000}" for i in range(len(bin_edges))]
bar.add_xaxis(x_labels[:-1])  # 只保留到倒数第二个标签,避免过多重复
bar.add_yaxis("评论数", hist.tolist(), label_opts=opts.LabelOpts(is_show=False))  # 设置 label_opts,is_show=False 表示不显示标签

# 设置图表的标题和坐标轴
bar.set_global_opts(
    title_opts=opts.TitleOpts(title="视频评论数直方图"),
    xaxis_opts=opts.AxisOpts(name="评论数区间", axislabel_opts={"rotate": 1000}),
    yaxis_opts=opts.AxisOpts(name="频次", type_="log", max_=50000)  # 使用对数刻度
)

# 渲染图表并保存
bar.render("comments_histogram.html")

print("直方图已保存为 'comments_histogram.html'")

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/788627.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

计算给定数字的阶乘

1 问题 计算给定数字的阶乘. 2 方法 使用while循环。使用for循环。使用函数。 通过实验、实践等证明提出的方法是有效的&#xff0c;是能够解决开头提出的问题。 代码清单 1 使用while循环numberint(input(请输入一个数字:))factorial1i1while i<number: factorialfactor…

【问题记录】VsCode中以管理员权限运行Powershell

问题展示 今天在尝试运行nodemon命令的时候出问题&#xff0c;显示没法识别&#xff0c;经过分析发现是管理员权限的问题&#xff0c;由于是在vscode里面进行开发&#xff0c;因此特此进行配置。 方法一 直接在vscode命令行中输入如下命令&#xff1a; Start-Process powers…

【Linux】vim详解

1.什么是vi/vim? 简单来说&#xff0c;vi是老式的文本编辑器&#xff0c;不过功能已经很齐全了&#xff0c;但是还是有可以进步的地方。vim则可以说是程序开发者的一项很好用的工具&#xff0c;就连 vim的官方网站&#xff08; http://www.vim.org&#xff09;自己也说vim是一…

SpringBoot:SpringBoot中如何实现对Http接口进行监控

一、前言 Spring Boot Actuator是Spring Boot提供的一个模块&#xff0c;用于监控和管理Spring Boot应用程序的运行时信息。它提供了一组监控端点&#xff08;endpoints&#xff09;&#xff0c;用于获取应用程序的健康状态、性能指标、配置信息等&#xff0c;并支持通过 HTTP …

使用OpenCV的absdiff函数报错

1.absdiff用法 absdiff函数用于计算两个输入图像之间每个像素的差异&#xff0c;并返回结果图像。 void cv::absdiff ( InputArray src1,InputArray src2,OutputArray dst ) //eg&#xff1a;比较两图像的差异 /*cv::Mat diff;cv::absdiff(depLeft32, imDepth, diff…

CentOS6禁止锁屏

在电源中设置后还是会锁屏, 原因是有屏幕保护程序 电源管理都 “从不” 一些AI的回答 在CentOS 6系统中&#xff0c;如果你想要禁用锁屏功能&#xff0c;可以编辑/etc/kbd/config文件。这个文件通常包含了键盘相关的设置&#xff0c;包括密码策略和屏幕锁定选项。 首先打开终…

Python | Leetcode Python题解之第226题翻转二叉树

题目&#xff1a; 题解&#xff1a; class Solution:def invertTree(self, root: TreeNode) -> TreeNode:if not root:return rootleft self.invertTree(root.left)right self.invertTree(root.right)root.left, root.right right, leftreturn root

Qt 异步实现事件的定时执行 - QTimer和QThread的联合使用

异步实现事件的定时执行 - QTimer和QThread的联合使用 引言一、核心源码二、其信号和槽函数简述三、定时器及其moveToThread简述 引言 在 Qt 中&#xff0c;如果想要定时执行某些事件或函数&#xff0c;通常会使用 QTimer 类。QTimer 允许设置一个时间间隔&#xff0c;当这个时…

MySQL架构你了解多少?

MySQL是一个服务器-客户端应用&#xff0c;MySQL8.0服务器是由连接池、服务管理工具和公共组件、NoSQL接口、SQL接口、解析器、优化器、缓存、存储引擎、文件系统组成。MySQL还为各种编程语言提供了一套用于外部程序访问服务器的连接器。整体架构图如下所示: MySQLConnectors:为…

基于eBPF的procstat软件追踪等待锁和持有锁的时间

在并发编程中&#xff0c;锁的使用是保证线程安全的重要手段。然而&#xff0c;过度使用锁或者锁竞争可能导致性能瓶颈。为了分析程序中锁的使用情况&#xff0c;我们可以借助procstat软件来追踪程序加锁时间和等待锁的时间。procstat是一个基于eBPF&#xff08;extended Berke…

double和float的区别与使用

double和float类型的区别与使用 在Java中&#xff0c;double和float都是基本数据类型&#xff0c;用于表示浮点数&#xff08;即带有小数点的数&#xff09;。 它们在精度和范围上有所不同&#xff1a; double类型提供了更高的精度和更大的范围&#xff0c;而float类型则精度更…

读人工智能全传08人工智能的今天

1. 人工智能的今天 1.1. 未来&#xff0c;或许有些领域会有非常明显的人工智能痕迹&#xff0c;有些领域则不会 1.2. 2018年&#xff0c;来自计算机视觉处理器公司英伟达的研究人员证明了人工智能软件能够创造出虚假的人物照片&#xff0c;并且能够完全…

旷野之间2 - 如何训练医疗保健小型语言模型(AI-SLM)

​​​​ 在本文中,我们将研究如何针对疾病症状训练一个小型医疗保健语言模型。为此,我们将从HuggingFace获取数据集(用于训练我们的模型):https://huggingface.co/datasets/QuyenAnhDE/Diseases_Symptoms QuyenAnhDE/Diseases_Symptoms 数据集来自 Hugging Face。图片来源…

【初阶数据结构】1.算法复杂度

文章目录 1.数据结构前言1.1 数据结构1.2 算法1.3 如何学好数据结构和算法 2.算法效率2.1 复杂度的概念2.2 复杂度的重要性 3.时间复杂度3.1 大O的渐进表示法3.2 时间复杂度计算示例3.2.1 示例13.2.2 示例23.2.3 示例33.2.4 示例43.2.5 示例53.2.6 示例63.2.7 示例7 4.空间复杂…

C语言编程4:复合赋值,递增递减运算符,局部变量与全局变量,本地变量,转义字符

一篇文章带你玩转C语言基础语法4&#xff1a;复合赋值&#xff0c;递增递减运算符&#xff0c;局部变量与全局变量&#xff0c;本地变量&#xff0c;转义字符 一、复合赋值&#x1f33f; 1.1&#x1f4a0;定义 赋值就是给任意一个变量或者常量赋一个值&#xff0c;这个值可以…

《操作系统真象还原》学习笔记:第2章——编写MBR主引导记录

2.1 计算机的启动过程 载入内存&#xff1a; &#xff08;1&#xff09; 程序被加载器&#xff08;软件或硬件&#xff09;加载到内存某个区域 &#xff08;2&#xff09;CPU 的 cs:ip 寄存器被指向这个程序的起始地址 2.2 软件接力第一棒&#xff0c;BIOS 2.2.1 实模式下的…

多模态:Nougat详解

文章目录 前言一、模型结构1. encoder2. decoder3. set 二、数据增强三、数据splitting the pages 四、实验评估repetitions during inference 五、代码1. 环境安装2. Dataset&#xff08;dataset.py&#xff09;3. Model&#xff08;model.py&#xff09; 总结 前言 科学知识…

UML建模工具Draw.io简介

新书速览|《UML 2.5基础、建模与设计实践 Draw.io是一个非常出色的免费、开源、简洁、方便的绘图软件&#xff0c;利用这款软件可以绘制出生动有趣的图形&#xff0c;包括流程图、地图、网络架构图、UML用例图、流程图等。它支持各种快捷键&#xff0c;免费提供了1000多张画图…

【基于深度学习方法的激光雷达点云配准系列之GeoTransformer】——粗配准

【GeoTransformer系列】——粗配准 1.coarse_matching1.1 概要1.2 功能1.3 超参1.4 input1.5 output2 coarse_target2.1 概要2.2 功能2.3 input2.4 output在模型部分有了初步了解之后, 接下来我们对后续的粗配准、精配准等部分依次进行解读. 本篇主要来看粗配准部分, 代码是在G…

机器人三定律及伦理分析

全世界的机器人定律并没有一个统一的标准或体系&#xff0c;但是在科学文献中&#xff0c;最广为人知的是由科幻小说家阿西莫夫提出的“机器人三定律”。本文将以这些定律为基础&#xff0c;分析现有的机器人伦理和实际应用中的问题&#xff0c;给出若干实例&#xff0c;并对相…