【关于时间序列的ML】项目 1 :使用 Python 进行 Covid-19 病例 预测
创始人
2024-04-27 12:20:09
0

 🔎大家好,我是Sonhhxg_柒,希望你看完之后,能对你有所帮助,不足请指正!共同学习交流🔎

📝个人主页-Sonhhxg_柒的博客_CSDN博客 📃

🎁欢迎各位→点赞👍 + 收藏⭐️ + 留言📝​

📣系列专栏 - 机器学习【ML】 自然语言处理【NLP】  深度学习【DL】

​​

 🖍foreword

✔说明⇢本人讲解主要包括Python、机器学习(ML)、深度学习(DL)、自然语言处理(NLP)等内容。

如果你对这个系列感兴趣的话,可以关注订阅哟👋

文章目录

使用 Python 进行 Covid-19 病例预测的机器学习项目

数据准备

数据可视化

使用 Python 预测未来 30 天的 Covid-19 病例


在本文中,我将向您介绍一个在接下来的 30 天内使用 Python 预测 Covid-19 病例的机器学习项目。这些类型的预测模型有助于准确预测流行病,这对于获取有关传染病可能传播和后果的信息至关重要。

政府和其他立法机构依靠这些机器学习预测模型和想法来提出新政策并评估应用政策的有效性。

使用 Python 进行 Covid-19 病例预测的机器学习项目

在接下来的 30 天内,我将通过导入必要的 Python 库和数据集来开始使用 Python 进行 Covid-19 病例预测的任务:

数据集1:

Kaggle: Your Home for Data Science

数据集2:

Kaggle: Your Home for Data Science

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import plotly.express as pxfrom fbprophet import Prophet
from sklearn.metrics import r2_scoreplt.style.use("ggplot")df0 = pd.read_csv("CONVENIENT_global_confirmed_cases.csv")
df1 = pd.read_csv("CONVENIENT_global_deaths.csv")

数据准备

现在下一步是数据准备,我将通过组合上述数据集来简单地准备新数据,然后我们将可视化数据的地理图以查看我们将要使用的内容:

world = pd.DataFrame({"Country":[],"Cases":[]})
world["Country"] = df0.iloc[:,1:].columns
cases = []
for i in world["Country"]:cases.append(pd.to_numeric(df0[i][1:]).sum())
world["Cases"]=casescountry_list=list(world["Country"].values)
idx = 0
for i in country_list:sayac = 0for j in i:if j==".":i = i[:sayac]country_list[idx]=ielif j=="(":i = i[:sayac-1]country_list[idx]=ielse:sayac += 1idx += 1
world["Country"]=country_list
world = world.groupby("Country")["Cases"].sum().reset_index()
world.head()
continent=pd.read_csv("continents2.csv")
continent["name"]=continent["name"].str.upper()
CountryCases
0Afghanistan45716.0
1Albania35600.0
2Algeria79110.0
3Andorra6534.0
4Angola14920.0

数据可视化

现在在这里我将准备三个可视化。一个将是地理可视化,以可视化 Covid-19 的全球传播。那么下一个可视化将是查看世界上每天发生的 Covid-19 病例。然后最后一个可视化将是查看世界上每天 Covid-19 的死亡病例。

现在让我们通过查看 Covid-19 的全球传播情况来开始数据可视化:

world["Cases Range"]=pd.cut(world["Cases"],[-150000,50000,200000,800000,1500000,15000000],labels=["U50K","50Kto200K","200Kto800K","800Kto1.5M","1.5M+"])
alpha =[]
for i in world["Country"].str.upper().values:if i == "BRUNEI":i="BRUNEI DARUSSALAM"elif  i=="US":i="UNITED STATES" if len(continent[continent["name"]==i]["alpha-3"].values)==0:alpha.append(np.nan)else:alpha.append(continent[continent["name"]==i]["alpha-3"].values[0])
world["Alpha3"]=alphafig = px.choropleth(world.dropna(),locations="Alpha3",color="Cases Range",projection="mercator",color_discrete_sequence=["white","khaki","yellow","orange","red"])
fig.update_geos(fitbounds="locations", visible=False)
fig.update_layout(margin={"r":0,"t":0,"l":0,"b":0})
fig.show()

 现在让我们看看世界各地的日常案例:

count = []
for i in range(1,len(df0)):count.append(sum(pd.to_numeric(df0.iloc[i,1:].values)))df = pd.DataFrame()
df["Date"] = df0["Country/Region"][1:]
df["Cases"] = count
df=df.set_index("Date")count = []
for i in range(1,len(df1)):count.append(sum(pd.to_numeric(df1.iloc[i,1:].values)))df["Deaths"] = countdf.Cases.plot(title="Daily Covid19 Cases in World",marker=".",figsize=(10,5),label="daily cases")
df.Cases.rolling(window=5).mean().plot(figsize=(10,5),label="MA5")
plt.ylabel("Cases")
plt.legend()
plt.show()

 

 现在让我们来看看 Covid-19 的每日死亡病例:

df.Deaths.plot(title="Daily Covid19 Deaths in World",marker=".",figsize=(10,5),label="daily deaths")
df.Deaths.rolling(window=5).mean().plot(figsize=(10,5),label="MA5")
plt.ylabel("Deaths")
plt.legend()
plt.show()

 

使用 Python 预测未来 30 天的 Covid-19 病例

现在,我将使用 Facebook 先知模型在接下来的 30 天内使用 Python 进行 Covid-19 病例预测任务。Facebook 先知模型使用时间序列方法进行预测。 

让我们看看我们如何在接下来的 30 天内使用 Facebook 先知模型通过 Python 进行 Covid-19 病例预测:

class Fbprophet(object):def fit(self,data):self.data  = dataself.model = Prophet(weekly_seasonality=True,daily_seasonality=False,yearly_seasonality=False)self.model.fit(self.data)def forecast(self,periods,freq):self.future = self.model.make_future_dataframe(periods=periods,freq=freq)self.df_forecast = self.model.predict(self.future)def plot(self,xlabel="Years",ylabel="Values"):self.model.plot(self.df_forecast,xlabel=xlabel,ylabel=ylabel,figsize=(9,4))self.model.plot_components(self.df_forecast,figsize=(9,6))def R2(self):return r2_score(self.data.y, self.df_forecast.yhat[:len(df)])df_fb  = pd.DataFrame({"ds":[],"y":[]})
df_fb["ds"] = pd.to_datetime(df.index)
df_fb["y"]  = df.iloc[:,0].valuesmodel = Fbprophet()
model.fit(df_fb)
model.forecast(30,"D")
model.R2()forecast = model.df_forecast[["ds","yhat_lower","yhat_upper","yhat"]].tail(30).reset_index().set_index("ds").drop("index",axis=1)
forecast["yhat"].plot(marker=".",figsize=(10,5))
plt.fill_between(x=forecast.index, y1=forecast["yhat_lower"], y2=forecast["yhat_upper"],color="gray")
plt.legend(["forecast","Bound"],loc="upper left")
plt.title("Forecasting of Next 30 Days Cases")
plt.show()

 

我希望您喜欢这篇关于使用 Python 编程语言预测未来 30 天 Covid-19 病例的文章。请随时在下面的评论部分提出您宝贵的问题。

相关内容

热门资讯

MySQL下载和安装(Wind... 前言:刚换了一台电脑,里面所有东西都需要重新配置,习惯了所...
Android---Banne... 轮播图是一种很常见的UI。Banner框架能够帮助我们快速开发,完成首页轮播图效果的需...
python -- PyQt5... 控件2 本章我们继续介绍PyQt5控件。这次的有 QPixmap , QLineEdi...
Mysql SQL优化跟踪来看... 背景 使用索引字段进行筛选数据时,explain查询语句发现MySQL居然没有使用索...
UG 6.0软件安装教程 UG 6.0软件安装教程 软件简介: UG 6.0是目前网络最好用、使用最为广泛的大型...
HTML静态网页作业——关于我... 家乡旅游景点网页作业制作 网页代码运用了DIV盒子的使用方法,如盒子的嵌套、浮动、ma...
MFC文件操作  MFC提供了一个文件操作的基类CFile,这个类提供了一个没有缓存的二进制格式的磁盘...
NoSQL数据库之Redis2 Redis 事务 事务的基础概念 关于事务最常见的例子就是银行转账,A 账户给 B 账...
Spring Security... 前言 在 Spring Security 中,默认的登陆方式是以表单形式进行提交参数的...
FourIE:一种最新联合事件... Cross-Task Instance Representation Interactions an...