Analise de Dados COVID-19 em Sao Paulo 2020-2023

Fonte httpscovidsaudegovbr rmlistls filechoose setwdCUsersmeiDocuments getwd arquivo1HISTPAINELCOVIDBR2020Parte113set2023csv arquivo2HISTPAINELCOVIDBR2020Parte213set2023csv arquivo3HISTPAINELCOVIDBR2021Parte113set2023csv arquivo4HISTPAINELCOVIDBR2021Parte213set2023csv arquivo5HISTPAINELCOVIDBR2022Parte113set2023csv arquivo6HISTPAINELCOVIDBR2022Parte213set2023csv arquivo7HISTPAINELCOVIDBR2023Parte113set2023csv arquivo8HISTPAINELCOVIDBR2023Parte213set2023csv initSystime dados1readcsv2arquivo1 headerTRUE sep dados2readcsv2arquivo2 headerTRUE sep dados3readcsv2arquivo3 headerTRUE sep dados4readcsv2arquivo4 headerTRUE sep dados5readcsv2arquivo5 headerTRUE sep dados6readcsv2arquivo6 headerTRUE sep dados7readcsv2arquivo7 headerTRUE sep dados8readcsv2arquivo8 headerTRUE sep Final Systime init Final empilhar os dados dadosrbinddados 1dados 2dados3dados4dados5dados6dados7dados8 objectsize dados1000000 rmdados 1dados 2dados3dados4dados5dados6dados7dados8 names dados dim dados head dados Data mining unicosestado unique dadosestado unicosestado unicoscidadesaopaulouniquedadosmunicipiowhichdadosestadoSP unicoscidadesaopaulo vírgula pois quero todas as colunas dadoscidsaopaulo dados which dadosestado SP dadosmunicipio São Paulo dim dadoscidsaopaulo head dadoscidsaopaulo tail dadoscidsaopaulo ydadoscidsaopaulocasosAcumulado2lengthdadoscidsaopaulocasosAcumulado ymdadoscidsaopauloobitosAcumulado2lengthdadoscidsaopaulocasosAcumulado Tlengthy tseq1T plotyttypep plot yttype p xlab t ylab main Casos Acumulados Cid São Paulo ablinev180hNULLcolbluelty3 plotymt typep xlab t ylab mainObitos Acumulados Cid Sao Paulo abline v 180h NULL dadoscidsaopaulo data 230 dimnames dados 2 Crio um vetor y com a variavel dependente date asDate dadoscidsaopaulodata 2T1 headdate plotyasDatedateformatdmYxaxtn ylab xlab main Casos Acumulados COVID Cid São Paulo type p ablinev 180h NULL axisDate para preencher o eixo X axisDate1 atseqdate1 datelengthdate byweek formatdmy las2 ystary1whichdate20230228 ystary1whichdate20210701 ystary1whichdate20201001 ystary1whichdate20230905 ystarprevywhichdate20230906lengthy datelengthy tstar1lengthystar tstarprevlengthy star 1lengthy tstar1lengthystar Tstarlengthystar plotystarasDatedate1Tstar format dmY xaxt n ylab xlab main Casos Acumulados COVID Cid São Paulo typep axisDate1 atseqdate1 dateTstar byweek formatdmy las2 escrevendo a função log da verossimilhança para o modelo logistico funcMLlogistico function theta Defino meus coeficiente a theta 1 a alpha b theta 2 b beta g theta3 g gamma sigma varystarTstar1Tstar05 Criase um vetor de residuos uystara1expbgtstar Escrevo a função de log da verossimilhanca lnL com uma distribuição normal lnL 05lengthy star log2pi05lengthystarlogabssigma 05sumu 2 sigma2 na frente da funcao para maximizala no algoritmo de otimizacao Agora calculo o valor da funcao return lnL Defino os chutes iniciais maxystar10 acmaxystar bclogacminystar1 gc003669434 chutescacbcgc Funcao nlm Non Linear Minimization lnLmaxlog nlm funcMLlogistico chutes steptol 1e6 gradtol 1e6 printlevel 1 iterlim 100000 ndigit 12 fscale 1 tail ystar lnLmaxlogestimate ahlog lnLmaxlog estimate 1 bhlog lnLmaxlog estimate 2 ghloglnLmaxlogestimate3 ahlog bhlog ghlog yhlogahlog1expbhlogghlogtstar yhlogprevahlog1expbhlogghlogtstarprev ystarprev plotystarprev typep linesyhlogprev plot ystartstartypepxlabDia ylabmainCasos COVID Cid São Paulo lines yhlogcol red Medidas de bondade do ajuste para o modelo logístico residlogprevystarprevyhlogprev histresidlogprev residlogystaryhlog histresidlog MSElogsumresidlog2Tstar MSElog RMSElogsqrtMSElog RMSElog MAElogsumabsresidlogTstar MAElog MAPElogsumabsresidlogystarTstar MAPElog Medidas para chute de pha de 0 sobre maxy funcMLgompfunctiontheta Defino meus coeficiente a theta 1 a alpha b theta2 b beta g theta3 g gamma sigma varystarlengthystar1lengthystar05 Criase um vetor de residuos uystaraexpexpbgtstar Escrevo a função de log da verossimilhanca lnL com uma distribuição normal lnL 05lengthy star log2pi05lengthystarlogabssigma 05sumu 2 sigma2 na frente da funcao para maximizala no algoritmo de otimizacao Agora calculo o valor da funcao return lnL Defino os chutes iniciais maxystar10 acmaxystar bcloglogaclogminystar gc2503299e02 chutescacbcgc Funcao nlm Non Linear Minimization lnLmaxlog nlm funcMLgomp chutes steptol 1e6 gradtol 1e6 printlevel 1 iterlim 100000 ndigit 12 fscale 1 tail ystar lnLmaxlogestimate a hgomp lnLmaxlogestimate 1 b hgomp lnLmaxlogestimate 2 ghgomplnLmaxlogestimate3 ahgomp bhgomp ghgomp yhgompahgomp1expbhgompghgomptstar plotystartstartypepxlabDia ylabmainCasos COVID Cid São Paulo linesyhgompcolblue L9232258 0 residgompystaryhgomp histresidgomp MSEgompsumresidgomp2Tstar MSEgomp RMSEgompsqrtMSEgomp RMSEgomp MAEgompsumabsresidgompTstar MAEgomp MAPEgompsumabsresidgompystarTstar MAPEgomp