澳洲抛砖引玉, 简单计算NSW的病例/病死率/疫苗相关性

在澳大利亚医疗保健





proxies = {
  'http': 'xxxxx',
  'https': 'xxxxx'
}
import requests
import pandas as pd
from scipy.stats import linregress
import matplotlib.pyplot as plt
from sklearn.preprocessing import MinMaxScaler
import numpy as np

urls = ['https://covidlive.com.au/report/daily-cases/nsw',
        'https://covidlive.com.au/report/daily-deaths/nsw',
       'https://covidlive.com.au/report/daily-vaccinations/nsw']

T = 13
T2 = 42
d = pd.DataFrame()
for url in urls:
    html = requests.get(url,proxies=proxies).content
    tab_list = pd.read_html(html)
    data = tab_list[1]
    data.set_index("DATE", inplace=True)
    d = pd.concat([d, data], axis=1)

d.index = pd.to_datetime(d.index)
data = d[['DEATHS','DOSES','CASES','NEW']]
data = data.dropna()
data = data.assign(DAILYDEATH = data["DEATHS"] - data["DEATHS"].shift(-1),
                  CASE_T = data["NEW"].shift(-1*T),
                  CASELOG = data["NEW"]/data["NEW"].shift(-1),
                  DEATHLOG = data["DEATHS"]/data["DEATHS"].shift(-1),
                  DOSE_T = data["DOSES"].shift(-1*T2))
data = data.dropna()
data.fillna(0, inplace=True)
data["CASELOG"] = data["CASELOG"].apply(np.log)
data["DEATHLOG"] = data["DEATHLOG"].apply(np.log)
data['CFR'] = data.DAILYDEATH/data.CASE_T
data = data.sort_index()
data = data.loc["2021-07-01":]
#[['NEW','DAILYDEATH','CASE_T','INCREASELOG','DOSE_T','CFR']]
standardizer = MinMaxScaler().fit(data)
data = pd.DataFrame(data, columns=['DEATHS', 'DOSES', 'CASES', 'NEW', 'DAILYDEATH', 'CASE_T',
       'CASELOG', 'DEATHLOG', 'DOSE_T', 'CFR'])


data[["CFR", "DOSE_T"]].corr()

        CFR        DOSE_T
CFR        1.000000        -0.108629
DOSE_T        -0.108629        1.000000
#完全疫苗率上升,死亡率下降

data[["CASELOG", "DOSE_T"]].corr()
CASELOG        DOSE_T
CASELOG        1.000000        -0.197275
DOSE_T        -0.197275        1.000000
#完全疫苗率上升,病例扩散速度下降

data[["CASELOG", "CFR"]].corr()
CASELOG        CFR
CASELOG        1.000000        0.091137
CFR        0.091137        1.000000

#病例扩散速度增加,死亡增加

评论
公共卫生/疫苗/病毒和统计学都是有相当高的专业门槛的。粗糙计算一下欢迎指正。

评论

太粗糙了

评论

欢迎指正。哪些方面应该提高精度、如何提高精度等等。谢谢。

评论

你的态度还是正确的。要提准确及精确度需要长时间的认真学习。就像“疫苗”一样,研发不到一年,那就不会知道一年后发生什么。需要足够长的时间。

评论

敢问一下你是不是有在中国国企当领导的经验?说话高屋建瓴,洋洋洒洒,永远正确-任何人都需要提高、学习啊。
我们做技术的,一般说这种人就是:虽然心里没数,不妨碍到处放屁。
当然如果你并没有到处放屁的习惯,请不要自己认领这个头衔。有建设性意见就提吧。

评论
安慰一下楼主,新闻版有不少ID的贴是不能回的。个人觉得a picture is worth a thousand words, 对于我们这些不是数学专业的人来说,我们要看图。

评论

赞同一图胜过千言万语。比如下面这张 E1B12C86-B51A-4080-98F2-C721E8D957BF.jpeg

评论


很多前两排的 ,尤其是develop too fast和own research。

评论

真不文明

评论

Who stands for gain的也不少哈哈。近期犹多。

评论

有个问题请务必回答一下:
你对疫苗的看法,都是基于同样的自信发表的吗?这种自信,是天生的呢,还是习得的呢?

评论

真可怜
澳洲中文论坛热点
悉尼部份城铁将封闭一年,华人区受影响!只能乘巴士(组图)
据《逐日电讯报》报导,从明年年中开始,因为从Bankstown和Sydenham的城铁将因Metro South West革新名目而
联邦政客们具有多少房产?
据本月早些时分报导,绿党副首领、参议员Mehreen Faruqi已获准在Port Macquarie联系其房产并建造三栋投资联

医疗,保健,保险,澳洲华人,澳大利亚论坛

医疗保健

Covid 又开始了。。

澳大利亚娃学校老师得了,好几个学生也中招,我也中了,朋友也有几个。 这快夏天了怎么又开始了。 评论 常態化呀,沒有多少人關心了 评论 现在没人在意了,就让他过去吧 评论 留个鼻涕而 ...

医疗保健

牙齿矫正后嘴歪了

澳大利亚女儿牙齿矫正,嘴巴歪了 图一是矫正前学校拍的证件照 图二是昨天我用手机拍的,明显可以看到人中歪了 图三也是昨天拍的,笑的时候更歪了 几个月前跟orthodontist说过这个问题,当时 ...

医疗保健

NSW看儿童的呼吸科专家

澳大利亚我的GP给我推荐的呼吸科专科医生不看小孩。 他说,只能推荐儿科医生。 我还是想问问,有没有看儿童的呼吸科专科医生。 我们在麦考利,Epping Eastwood 一带。远点也行啊。 谢谢 评论 呼 ...

医疗保健

好贵! 牙科挂号要收266刀

澳大利亚第一次去医院看牙病, 啥都不懂, 请各位大神指点: 前几天牙疼, 可能是智齿发炎或上火了.我看到公司附近有个公立牙科医院叫The Royal Dental Hospital of Melbourne, 就拿着medicare卡去看急诊, 谁知 ...

医疗保健

问个买保险的问题

澳大利亚父母移民过来了,想买保险。可是怎么才能免LHC呢?网上买有loading。哪位同学有经验指导一下。谢谢 评论 拿medicare的第一年内必须买。超过了就没办法了。 如果一年内,和Medicare要封信 ...

医疗保健

40岁后,左侧鼻孔总是堵着的,怎么办

澳大利亚40岁后,左侧鼻孔总是堵着的,不是全堵住,但是呼吸很费力,好像鼻腔变细了似得。 如果用手堵住右边的鼻孔,左边的需要大力呼吸,而且呼吸声音很大。 于是总是用右侧鼻孔呼吸, ...

医疗保健

花粉症治疗新方法

澳大利亚最近越发严重了墨尔本的花粉真的是 眼睛巨痒有时候连嗓子和耳朵都是痒的telfast zyrtec claratine 各种吃还是不行 听说中医可以治疗 请问大家有人试过么 效果怎么样 是针灸还是拔罐?? ...

医疗保健

感冒需要去医院吗?

澳大利亚上周女儿花粉症发作,我买了空气清新机,结果第二天起女儿发烧38度,才发现是感冒了(上周日),还有干咳,吃 Codral Day Night + Dry Cough 。烧退了,但还是咳嗽,周三开始有痰,痰的 ...

医疗保健

分分急问,13岁的孩子,头疼用什么药?

澳大利亚如题,昨晚说头疼的睡不着。有什么办法缓解吗?或者吃药。她不要去看医生 评论 Panadol 评论 panadol,吃了就好就没问题,我女儿头痛已经拍了MRI的程度,后来被GP查出来咬合牙齿有问题 ...

医疗保健

请问澳洲治疗心脏疾病 水平怎么样

澳大利亚确诊为心衰 正在考虑回澳洲看还是在国内看 请熟悉这方面的朋友给点建议 谢谢 评论 以前不知道在哪里看到个news 间接提到澳洲心脏治疗水平 世界领先。。。 评论 楼主多大岁数了? ...

医疗保健

我父亲在国内确诊为前列腺癌

澳大利亚我父亲八十几岁了,在国内确诊为前列腺癌,他有澳洲PR,现在要回澳洲治疗,不知道这边是否承认中国的活检结果,是不是还要让再做活检。另外看专科医生,再等医院治疗是不是要排 ...

医疗保健

请教,老人夜里常抽筋怎么办?

澳大利亚家里近70岁老人,身体尚好. 近几年来夜里常抽筋,白天常未发生过.已看了多次GP,没有多少改善.谢谢! 评论 不知道别人抽筋原因 我知道自己的原因:怀孕的时候抽筋过,缺钙。现在 ...

医疗保健

BUPA extra到底有买的必要吗

澳大利亚目前一家三口买了BUPA私保,Silver Plus Essential Hospital with Your Choice Extras 60,一个月大约320。 有同学知道,大约这个EXTRA占了多少钱吗? 另外,这个一年核算下来,是不是如果只是洗洗牙, ...

医疗保健

关于小腿和脚发凉

澳大利亚我最近几年时不时的感觉左小腿和脚凉,右边没有,即便是夏天也有感觉。有时睡觉会感觉左膝盖后侧和床接触的部分冷。去看了2个GP,都只是看看,说没有肿,外面也看不出问题,说 ...