1.hist Դ?源码?
2.声线年龄层怎么划分？

hist Դ??

       项目内容

       案例选择商品类目：沙发；数量：共页个商品；筛选条件：天猫、销量从高到低、源码价格元以上。源码

       以下是源码分析，源码点击文末链接

       项目目的源码

       1. 对商品标题进行文本分析，词云可视化。源码idea导入tomcat源码

       2. 不同关键词word对应的源码sales统计分析。

       3. 商品的源码价格分布情况分析。

       4. 商品的源码销量分布情况分析。

       5. 不同价格区间的源码商品的平均销量分布。

       6. 商品价格对销量的源码影响分析。

       7. 商品价格对销售额的源码影响分析。

       8. 不同省份或城市的源码商品数量分布。

       9. 不同省份的源码商品平均销量分布。

       注：本项目仅以以上几项分析为例。源码

       项目步骤

       1. 数据采集：Python爬取淘宝网商品数据。

       2. 数据清洗和处理。

       3. 文本分析：jieba分词、wordcloud可视化。

       4. 数据柱形图可视化barh。

       5. 数据直方图可视化hist。

       6. 数据散点图可视化scatter。

       7. 数据回归分析可视化regplot。

       工具&模块：

       工具：本案例代码编辑工具Anaconda的Spyder。

       模块：requests、retrying、missingno、jieba、matplotlib、wordcloud、imread、seaborn等。

       原代码和相关文档后台回复“淘宝”下载。

       一、爬取数据

       因淘宝网是反爬虫的，虽然使用多线程、修改headers参数，但仍然不能保证每次%爬取，所以，我增加了循环爬取，直至所有页爬取成功停止。

       说明：淘宝商品页为JSON格式，人脉社区app源码这里使用正则表达式进行解析。

       代码如下：

       二、数据清洗、处理：

       (此步骤也可以在Excel中完成，再读入数据)

       代码如下：

       说明：根据需求，本案例中只取了item_loc、raw_title、view_price、view_sales这4列数据，主要对标题、区域、价格、销量进行分析。

       代码如下:

       三、数据挖掘与分析：

       1. 对raw_title列标题进行文本分析：

       使用结巴分词器，安装模块pip install jieba。

       对title_s（list of list格式）中的每个list的元素（str）进行过滤，剔除不需要的词语，即把停用词表stopwords中有的词语都剔除掉：

       为了准确性，这里对过滤后的数据title_clean中的每个list的元素进行去重，即每个标题被分割后的词语唯一。

       观察word_count表中的词语，发现jieba默认的词典无法满足需求。

       有的词语（如可拆洗、不可拆洗等）却被cut，这里根据需求对词典加入新词（也可以直接在词典dict.txt里面增删，然后载入修改过的dict.txt）。

       词云可视化：

       安装模块wordcloud。

       方法1：pip install wordcloud。

       方法2：下载Packages安装：pip install 软件包名称。

       软件包下载地址：lfd.uci.edu/~gohlke/pyt...

       注意：要把下载的软件包放在Python安装路径下。

       代码如下：

       分析

       1. 组合、整装商品占比很高；

       2. 从沙发材质看：布艺沙发占比很高，比皮艺沙发多；

       3. 从沙发风格看：简约风格最多，北欧风次之，其他风格排名依次是美式、中式、日式、法式等；

       4. 从户型看：小户型占比最高、大小户型次之，大户型最少。草莓熊代码源码

       2. 不同关键词word对应的sales之和的统计分析：

       （说明：例如词语‘简约’，则统计商品标题中含有‘简约’一词的商品的销量之和，即求出具有‘简约’风格的商品销量之和）

       代码如下：

       对表df_word_sum中的word和w_s_sum两列数据进行可视化。

       （本例中取销量排名前的词语进行绘图）

       由图表可知：

       1. 组合商品销量最高；

       2. 从品类看：布艺沙发销量很高，远超过皮艺沙发；

       3. 从户型看：小户型沙发销量最高，大小户型次之，大户型销量最少；

       4. 从风格看：简约风销量最高，北欧风次之，其他依次是中式、美式、日式等；

       5. 可拆洗、转角类沙发销量可观，也是颇受消费者青睐的。

       3. 商品的价格分布情况分析：

       分析发现，有一些值太大，为了使可视化效果更加直观，这里我们选择价格小于的商品。

       代码如下：

       由图表可知：

       1. 商品数量随着价格总体呈现下降阶梯形势，价格越高，在售的商品越少；

       2. 低价位商品居多，价格在-之间的商品最多，-之间的次之，价格1万以上的商品较少；

       3. 价格1万元以上的商品，在售商品数量差异不大。

       4. 商品的销量分布情况分析：

       同样，为了使可视化效果更加直观，这里我们选择销量大于的商品。

       代码如下：

       由图表及数据可知：

       1. 销量以上的商品仅占3.4%，其中销量-之间的商品最多，-之间的次之；

       2. 销量-之间，商品的数量随着销量呈现下降趋势，且趋势陡峭，低销量商品居多；

       3. 销量以上的商品很少。

       5. 不同价格区间的商品的平均销量分布：

       代码如下：

       由图表可知：

       1. 价格在-之间的商品平均销量最高，-之间的次之，元以上的最低；

       2. 总体呈现先增后减的趋势，但最高峰处于相对低价位阶段；

       3. 说明广大消费者对购买沙发的需求更多处于低价位阶段，在元以上价位越高平均销量基本是越少。

       6. 商品价格对销量的影响分析：

       同上，为了使可视化效果更加直观，这里我们选择价格小于的蚌埠专业商城源码商品。

       代码如下：

       由图表可知：

       1. 总体趋势：随着商品价格增多其销量减少，商品价格对其销量影响很大；

       2. 价格-之间的少数商品销量冲的很高，价格-之间的商品多数销量偏低，少数相对较高，但价格以上的商品销量均很低，没有销量突出的商品。

       7. 商品价格对销售额的影响分析：

       代码如下：

       由图表可知：

       1. 总体趋势：由线性回归拟合线可以看出，商品销售额随着价格增长呈现上升趋势；

       2. 多数商品的价格偏低，销售额也偏低；

       3. 价格在0-的商品只有少数销售额较高，价格2万-6万的商品只有3个销售额较高，价格6-万的商品有1个销售额很高，而且是最大值。

       8. 不同省份的商品数量分布：

       代码如下：

       由图表可知：

       1. 广东的最多，上海次之，江苏第三，尤其是广东的数量远超过江苏、浙江、上海等地，说明在沙发这个子类目，广东的店铺占主导地位；

       2. 江浙沪等地的数量差异不大，基本相当。

       9. 不同省份的商品平均销量分布：

       代码如下：

       热力型地图

       源码：Python爬取淘宝商品数据挖掘分析实战

声线年龄层怎么划分？

       从EDA、音频预处理到特征工程和数据建模的完整源代码演示

       大多数人都熟悉如何在图像、文本或表格数据上运行数据科学项目。 但处理音频数据的样例非常的少见。 在本文中，将介绍如何在机器学习的帮助下准备、探索和分析音频数据。 简而言之：与其他的形式（例如文本或图像）类似我们需要将音频数据转换为机器可识别的格式。

       音频数据的有趣之处在于您可以将其视为多种不同的模式：

       · 可以提取高级特征并分析表格数据等数据。

       · 可以计算频率图并分析图像数据等数据。

       · 可以使用时间敏感模型并分析时间序列数据等数据。

       · 可以使用语音到文本模型并像文本数据一样分析数据。

       在本文中，我们将介绍前三种方法。 首先看看音频数据的实际样子。

       音频数据的格式

       虽然有多个 Python 库可以处理音频数据，但我们推荐使用 librosa。 让我们加载一个 MP3 文件并绘制它的内容。

       # Import librosa

       import librosa

       # Loads mp3 file with a specific sampling rate, here kHz

       y, sr = librosa.load("c4_sample-1.mp3", sr=_)

       # Plot the signal stored in 'y'

       from matplotlib import pyplot as plt

       import librosa.display

       plt.figure(figsize=(, 3))

       plt.title("Audio signal as waveform")

       librosa.display.waveplot(y, sr=sr);

       这里看到的是句子的波形表示。

       1、象棋源码小程序波形 - 信号的时域表示

       之前称它为时间序列数据，但现在我们称它为波形？ 当只看这个音频文件的一小部分时，这一点变得更加清晰。 下图显示了与上面相同的内容，但这次只有 .5 毫秒。

       我们看到的是一个时间信号，它以不同的频率和幅度在值 0 附近振荡。该信号表示气压随时间的变化，或扬声器膜（或耳膜）的物理位移 . 这就是为什么这种对音频数据的描述也称为波形的原因。

       频率是该信号振荡的速度。 低频例如 Hz 可能是低音吉他的声音，而鸟儿的歌声可能是 Hz 的更高频率。 我们人类语言通常介于两者之间。

       要知道这个信号在单位时间内从连续信号中提取并组成离散信号的采样个数，我们使用赫兹（Hz）来表示每秒的采样个数。 ' 或 k Hz表示美标采集了次。 我们在上图中可以看到的 1' 个时间点代表了 .5 毫秒（/ = 0.）的音频信号。

       2、傅里叶变换——信号的频域表示

       虽然之前的可视化可以告诉我们什么时候发生了（即 2 秒左右似乎有很多波形信号），但它不能真正告诉我们它发生的频率。 因为波形向我们显示了有关时间的信息，所以该信号也被称为信号的时域表示。

       可以使用快速傅立叶变换，反转这个问题并获得关于存在哪些频率的信息，同时丢弃掉关于时间的信息。 在这种情况下，信号表示被称为信号的频域表示。

       让我们看看之前的句子在频域中的表现。

       import scipy

       import numpy as np

       # Applies fast fourier transformation to the signal and takes absolute values

       y_freq = np.abs(scipy.fftpack.fft(y))

       # Establishes all possible frequency

       # (dependent on the sampling rate and the length of the signal)

       f = np.linspace(0, sr, len(y_freq))

       # Plot audio signal as frequency information.

       plt.figure(figsize=(, 3))

       plt.semilogx(f[: len(f) // 2], y_freq[: len(f) // 2])

       plt.xlabel("Frequency (Hz)")

       plt.show();

       可以在此处看到大部分信号在 ~ 到 ~ Hz 之间（即 ² 到 ³ 之间）。 另外，似乎还有一些从 1' 到 ' Hz 的内容。

       3、频谱图

       我们并不总是需要决定时域或频域。 使用频谱图同时表示这两个领域中的信息，同时将它们的大部差别保持在最低限度。 有多种方法可以创建频谱图，但在本文中将介绍常见的三种。

       3a 短时傅里叶变换 (STFT)

       这时是之前的快速傅立叶变换的小型改编版本，即短时傅立叶变换 (STFT)， 这种方式是以滑动窗口的方式计算多个小时间窗口（因此称为“短时傅立叶”）的 FFT。

       import librosa.display

       # Compute short-time Fourier Transform

       x_stft = np.abs(librosa.stft(y))

       # Apply logarithmic dB-scale to spectrogram and set maximum to 0 dB

       x_stft = librosa.amplitude_to_db(x_stft, ref=np.max)

       # Plot STFT spectrogram

       plt.figure(figsize=(, 4))

       librosa.display.specshow(x_stft, sr=sr, x_axis="time", y_axis="log")

       plt.colorbar(format="%+2.0f dB")

       plt.show();

       与所有频谱图一样，颜色代表在给定时间点给定频率的量（响度/音量）。 +0dB 是最响亮的，-dB 接近静音。 在水平 x 轴上我们可以看到时间，而在垂直 y 轴上我们可以看到不同的频率。

       3b 梅尔谱图

       作为 STFT 的替代方案，还可以计算基于 mel 标度的梅尔频谱图。 这个尺度解释了我们人类感知声音音高的方式。 计算 mel 标度，以便人类将由 mel 标度中的 delta 隔开的两对频率感知为具有相同的感知差异。

       梅尔谱图的计算与 STFT 非常相似，主要区别在于 y 轴使用不同的刻度。

       # Compute the mel spectrogram

       x_mel = librosa.feature.melspectrogram(y=y, sr=sr)

       # Apply logarithmic dB-scale to spectrogram and set maximum to 0 dB

       x_mel = librosa.power_to_db(x_mel, ref=np.max)

       # Plot mel spectrogram

       plt.figure(figsize=(, 4))

       librosa.display.specshow(x_mel, sr=sr, x_axis="time", y_axis="mel")

       plt.colorbar(format="%+2.0f dB")

       plt.show();

       与 STFT 的区别可能不太明显，但如果仔细观察，就会发现在 STFT 图中，从 0 到 Hz 的频率在 y 轴上占用的空间比在 mel 图中要大得多 .

       3c 梅尔频率倒谱系数 (MFCC)

       梅尔频率倒谱系数 (MFCC) 是上面梅尔频谱图的替代表示。 MFCC 相对于 梅尔谱图的优势在于特征数量相当少（即独特的水平线标度），通常约为 。

       由于梅尔频谱图更接近我们人类感知音高的方式，并且 MFCC 只有少数几个分量特征，所以大多数机器学习从业者更喜欢 使用MFCC 以“图像方式”表示音频数据。 但是对于某些问题，STFT、mel 或波形表示可能会更好。

       让我们继续计算 MFCC 并绘制它们。

       # Extract 'n_mfcc' numbers of MFCCs components (here )

       x_mfccs = librosa.feature.mfcc(y, sr=sr, n_mfcc=)

       # Plot MFCCs

       plt.figure(figsize=(, 4))

       librosa.display.specshow(x_mfccs, sr=sr, x_axis="time")

       plt.colorbar()

       plt.show();

       数据清洗

       现在我们更好地理解了音频数据的样子，让我们可视化更多示例。

       在这四个示例中，我们可以收集到有关此音频数据集的更多问题：

       · 大多数录音在录音的开头和结尾都有一段较长的静默期（示例 1 和示例 2）。 这是我们在“修剪”时应该注意的事情。

       · 在某些情况下，由于按下和释放录制按钮，这些静音期会被“点击”中断（参见示例 2）。

       · 一些录音没有这样的静音阶段，即一条直线（示例 3 和 4）。

       · 在收听这些录音时，有大量背景噪音。

       为了更好地理解这在频域中是如何表示的，让我们看一下相应的 STFT 频谱图。

       当听录音时，可以观察到样本 3 具有覆盖多个频率的不同背景噪声，而样本 4 中的背景噪声相当恒定。 这也是我们在上图中看到的。 样本 3 在整个过程中都非常嘈杂，而样本 4 仅在几个频率上（即粗水平线）有噪声。 我们不会详细讨论如何消除这种噪音，因为这超出了本文的范围。

       但是让我们研究一下如何消除此类噪音并修剪音频样本的“捷径”。 虽然使用自定义过滤函数的更手动的方法可能是从音频数据中去除噪声的最佳方法，但在我们的例子中，将推荐使用实用的 python 包 noisereduce。

       import noisereduce as nr

       from scipy.io import wavfile

       # Loop through all four samples

       for i in range(4):

       # Load audio file

       fname = "c4_sample-%d.mp3" % (i + 1)

       y, sr = librosa.load(fname, sr=_)

       # Remove noise from audio sample

       reduced_noise = nr.reduce_noise(y=y, sr=sr, stationary=False)

       # Save output in a wav file as mp3 cannot be saved to directly

       wavfile.write(fname.replace(".mp3", ".wav"), sr, reduced_noise)

       聆听创建的 wav 文件，可以听到噪音几乎完全消失了。 虽然我们还引入了更多的代码，但总的来说我们的去噪方法利大于弊。

       对于修剪步骤，可以使用 librosa 的 .effects.trim() 函数。每个数据集可能需要一个不同的 top_db 参数来进行修剪，所以最好进行测试，看看哪个参数值好用。 在这个的例子中，它是 top_db=。

       # Loop through all four samples

       for i in range(4):

       # Load audio file

       fname = "c4_sample-%d.wav" % (i + 1)

       y, sr = librosa.load(fname, sr=_)

       # Trim signal

       y_trim, _ = librosa.effects.trim(y, top_db=)

       # Overwrite previous wav file

       wavfile.write(fname.replace(".mp3", ".wav"), sr, y_trim)

       现在让我们再看一下清理后的数据。

       看样子好多了

       特征提取

       数据是干净的，应该继续研究可以提取的特定于音频的特征了。

       1、开始检测

       通过观察一个信号的波形，librosa可以很好地识别一个新口语单词的开始。

       # Extract onset timestamps of words

       onsets = librosa.onset.onset_detect(

       y=y, sr=sr, units="time", hop_length=, backtrack=False)

       # Plot onsets together with waveform plot

       plt.figure(figsize=(8, 3))

       librosa.display.waveplot(y, sr=sr, alpha=0.2, x_axis="time")

       for o in onsets:

       plt.vlines(o, -0.5, 0.5, colors="r")

       plt.show()

       # Return number of onsets

       number_of_words = len(onsets)

       print(f"{ number_of_words} onsets were detected in this audio signal.")

       >>> 7 onsets were detected in this audio signal

       2、录音的长度

       与此密切相关的是录音的长度。录音越长，能说的单词就越多。所以计算一下录音的长度和单词被说出的速度。

       duration = len(y) / sr

       words_per_second = number_of_words / duration

       print(f"""The audio signal is { duration:.2f} seconds long,

       with an average of { words_per_second:.2f} words per seconds.""")

       >>> The audio signal is 1. seconds long,

       >>> with an average of 4. words per seconds.

       3、节奏

       语言是一种非常悦耳的信号，每个人都有自己独特的说话方式和语速。因此，可以提取的另一个特征是说话的节奏，即在音频信号中可以检测到的节拍数。

       # Computes the tempo of a audio recording

       tempo = librosa.beat.tempo(y, sr, start_bpm=)[0]

       print(f"The audio signal has a speed of { tempo:.2f} bpm.")

       >>> The audio signal has a speed of . bpm.

       4、基频

       基频是周期声音出现时的最低频率。在音乐中也被称为音高。在之前看到的谱图图中，基频(也称为f0)是图像中最低的亮水平条带。而在这个基本音之上的带状图案的重复称为谐波。

       为了更好地说明确切意思，下面提取基频，并在谱图中画出它们。

       # Extract fundamental frequency using a probabilistic approach

       f0, _, _ = librosa.pyin(y, sr=sr, fmin=, fmax=, frame_length=)

       # Establish timepoint of f0 signal

       timepoints = np.linspace(0, duration, num=len(f0), endpoint=False)

       # Plot fundamental frequency in spectrogram plot

       plt.figure(figsize=(8, 3))

       x_stft = np.abs(librosa.stft(y))

       x_stft = librosa.amplitude_to_db(x_stft, ref=np.max)

       librosa.display.specshow(x_stft, sr=sr, x_axis="time", y_axis="log")

       plt.plot(timepoints, f0, color="cyan", linewidth=4)

       plt.show();

       在 Hz 附近看到的绿线是基本频率。 但是如何将其用于特征工程呢？ 可以做的是计算这个 f0 的具体特征。

       # Computes mean, median, 5%- and %-percentile value of fundamental frequency

       f0_values = [

       np.nanmean(f0),

       np.nanmedian(f0),

       np.nanstd(f0),

       np.nanpercentile(f0, 5),

       np.nanpercentile(f0, ),

       ]

       print("""This audio signal has a mean of { :.2f}, a median of { :.2f}, a

       std of { :.2f}, a 5-percentile at { :.2f} and a -percentile at { :.2f}.""".format(*f0_values))

       >>> This audio signal has a mean of ., a median of ., a

       >>> std of 4., a 5-percentile at . and a -percentile at ..

       除以上说的技术意外，还有更多可以探索的音频特征提取技术，这里就不详细说明了。

       音频数据集的探索性数据分析 (EDA)

       现在我们知道了音频数据是什么样子以及如何处理它，让我们对它进行适当的 EDA。 首先下载一个数据集Kaggle 的 Common Voice 。 这个 GB 的大数据集只是来自 Mozilla 的 + GB 大数据集的一个小的快照。 对于本文这里的示例，将只使用这个数据集的大约 9' 个音频文件的子样本。

       看看这个数据集和一些已经提取的特征。

       1、特征分布调查

       目标类别年龄和性别的类别分布。

       目标类别分布是不平衡的

       下一步，让我们仔细看看提取的特征的值分布。

       除了 words_per_second，这些特征分布中的大多数都是右偏的，因此可以从对数转换中获益。

       import numpy as np

       # Applies log1p on features that are not age, gender, filename or words_per_second

       df = df.apply(

       lambda x: np.log1p(x)

       if x.name not in ["age", "gender", "filename", "words_per_second"]

       else x)

       # Let's look at the distribution once more

       df.drop(columns=["age", "gender", "filename"]).hist(

       bins=, figsize=(, ))

       plt.show();

       好多了，但有趣的是 f0 特征似乎都具有双峰分布。 让我们绘制与以前相同的内容，但这次按性别分开。

       正如怀疑的那样，这里似乎存在性别效应！ 但也可以看到，一些 f0 分数（这里特别是男性）比应有的低和高得多。 由于特征提取不良，这些可能是异常值。 仔细看看下图的所有数据点。

       # Plot sample points for each feature individually

       df.plot(lw=0, marker=".", subplots=True, layout=(-1, 3),

       figsize=(, 7.5), markersize=2)

       plt.tight_layout()

       plt.show();

       鉴于特征的数量很少，而且有相当漂亮的带有明显尾部的分布，可以遍历它们中的每一个，并逐个特征地确定异常值截止阈值。

       2、特征的相关性

       下一步，看看所有特征之间的相关性。 但在这样做之前需要对非数字目标特征进行编码。 可以使用 scikit-learn 的 OrdinalEncoder 来执行此操作，但这可能会破坏年龄特征中的正确顺序。 因此在这里手动进行映射。

       # Map age to appropriate numerical value

       df.loc[:, "age"] = df["age"].map({

       "teens": 0,

       "twenties": 1,

       "thirties": 2,

       "fourties": 3,

       "fifties": 4,

       "sixties": 5})

       # Map gender to corresponding numerical value

       df.loc[:, "gender"] = df["gender"].map({ "male": 0, "female": 1})

       现在可以使用 pandas 的 .corr() 函数和 seaborn 的 heatmap() 来更深入地了解特征相关性。

       import seaborn as sns

       plt.figure(figsize=(8, 8))

       df_corr = df.corr() *

       sns.heatmap(df_corr, square=True, annot=True, fmt=".0f",

       mask=np.eye(len(df_corr)), center=0)

       plt.show();

       非常有趣！提取的 f0 特征似乎与性别目标有相当强的关系，而年龄似乎与任何其他的特征都没有太大的相关性。

       3、频谱图特征

       目前还没有查看实际录音。 正如之前看到的，有很多选择（即波形或 STFT、mel 或 mfccs 频谱图）。

       音频样本的长度都不同，这意味着频谱图也会有不同的长度。 因此为了标准化所有录音，首先要将它们剪切到正好 3 秒的长度：太短的样本会被填充，而太长的样本会被剪掉。

       一旦计算了所有这些频谱图，我们就可以继续对它们执行一些 EDA！ 而且因为看到“性别”似乎与录音有特殊的关系，所以分别可视化两种性别的平均梅尔谱图，以及它们的差异。

       男性说话者的平均声音低于女性。 这可以通过差异图中的较低频率（在红色水平区域中看到）的更多强度来看出。

       模型选择

       现在已经可以进行建模了。我们有多种选择。关于模型，我们可以……

       · 训练我们经典（即浅层）机器学习模型，例如 LogisticRegression 或 SVC。

       · 训练深度学习模型，即深度神经网络。

       · 使用 TensorflowHub 的预训练神经网络进行特征提取，然后在这些高级特征上训练浅层或深层模型

       而我们训练的数据是

       · CSV 文件中的数据，将其与频谱图中的“mel 强度”特征相结合，并将数据视为表格数据集

       · 单独的梅尔谱图并将它们视为图像数据集

       · 使用TensorflowHub现有模型提取的高级特征，将它们与其他表格数据结合起来，并将其视为表格数据集

       当然，有许多不同的方法和其他方法可以为建模部分创建数据集。因为我们没有使用全量的数据，所以在本文我们使用最简单的机器学习模型。

       经典（即浅层）机器学习模型

       这里使用EDA获取数据，与一个简单的 LogisticRegression 模型结合起来，看看我们能在多大程度上预测说话者的年龄。除此以外还使用 GridSearchCV 来探索不同的超参数组合，以及执行交叉验证。

       from sklearn.linear_model import LogisticRegression

       from sklearn.preprocessing import RobustScaler, PowerTransformer, QuantileTransformer

       from sklearn.decomposition import PCA

       from sklearn.pipeline import Pipeline

       from sklearn.model_selection import GridSearchCV

       # Create pipeline

       pipe = Pipeline(

       [

       ("scaler", RobustScaler()),

       ("pca", PCA()),

       ("logreg", LogisticRegression(class_weight="balanced")),

       ]

       )

       # Create grid

       grid = {

       "scaler": [RobustScaler(), PowerTransformer(), QuantileTransformer()],

       "pca": [None, PCA(0.)],

       "logreg__C": np.logspace(-3, 2, num=),

       }

       # Create GridSearchCV

       grid_cv = GridSearchCV(pipe, grid, cv=4, return_train_score=True, verbose=1)

       # Train GridSearchCV

       model = grid_cv.fit(x_tr, y_tr)

       # Collect results in a DataFrame

       cv_results = pd.DataFrame(grid_cv.cv_results_)

       # Select the columns we are interested in

       col_of_interest = [

       "param_scaler",

       "param_pca",

       "param_logreg__C",

       "mean_test_score",

       "mean_train_score",

       "std_test_score",

       "std_train_score",

       ]

       cv_results = cv_results[col_of_interest]

       # Show the dataframe sorted according to our performance metric

       cv_results.sort_values("mean_test_score", ascending=False)

       作为上述 DataFrame 输出的补充，还可以将性能得分绘制为探索的超参数的函数。 但是因为使用了有多个缩放器和 PCA ，所以需要为每个单独的超参数组合创建一个单独的图。

       在图中，可以看到总体而言模型的表现同样出色。 当降低 C 的值时，有些会出现更快的“下降”，而另一些则显示训练和测试（这里实际上是验证）分数之间的差距更大，尤其是当我们不使用 PCA 时。

       下面使用 best_estimator_ 模型，看看它在保留的测试集上的表现如何。

       # Compute score of the best model on the withheld test set

       best_clf = model.best_estimator_

       best_clf.score(x_te, y_te)

       >>> 0.

       这已经是一个很好的成绩了。 但是为了更好地理解分类模型的表现如何，可以打印相应的混淆矩阵。

       虽然该模型能够检测到比其他模型更多的 岁样本（左混淆矩阵），但总体而言，它实际上在对 岁和 岁的条目进行分类方面效果更好（例如，准确率分别为 % 和 %）。

       总结

       在这篇文章中，首先看到了音频数据是什么样的，然后可以将其转换成哪些不同的形式，如何对其进行清理和探索，最后如何将其用于训练一些机器学习模型。如果您有任何问题，请随时发表评论。

       最后本文的源代码在这里下载：

       /post/5cfeda9d3dc

       作者：Michael Notter