\documentclass[11pt,twoside]{article}\makeatletter

\IfFileExists{xcolor.sty}%
  {\RequirePackage{xcolor}}%
  {\RequirePackage{color}}
\usepackage{colortbl}
\usepackage{wrapfig}
\usepackage{ifxetex}
\ifxetex
  \usepackage{fontspec}
  \usepackage{xunicode}
  \catcode`⃥=\active \def⃥{\textbackslash}
  \catcode`❴=\active \def❴{\{}
  \catcode`❵=\active \def❵{\}}
  \def\textJapanese{\fontspec{Noto Sans CJK JP}}
  \def\textChinese{\fontspec{Noto Sans CJK SC}}
  \def\textKorean{\fontspec{Noto Sans CJK KR}}
  \setmonofont{DejaVu Sans Mono}
  
\else
  \IfFileExists{utf8x.def}%
   {\usepackage[utf8x]{inputenc}
      \PrerenderUnicode{–}
    }%
   {\usepackage[utf8]{inputenc}}
  \usepackage[english]{babel}
  \usepackage[T1]{fontenc}
  \usepackage{float}
  \usepackage[]{ucs}
  \uc@dclc{8421}{default}{\textbackslash }
  \uc@dclc{10100}{default}{\{}
  \uc@dclc{10101}{default}{\}}
  \uc@dclc{8491}{default}{\AA{}}
  \uc@dclc{8239}{default}{\,}
  \uc@dclc{20154}{default}{ }
  \uc@dclc{10148}{default}{>}
  \def\textschwa{\rotatebox{-90}{e}}
  \def\textJapanese{}
  \def\textChinese{}
  \IfFileExists{tipa.sty}{\usepackage{tipa}}{}
\fi
\def\exampleFont{\ttfamily\small}
\DeclareTextSymbol{\textpi}{OML}{25}
\usepackage{relsize}
\RequirePackage{array}
\def\@testpach{\@chclass
 \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else
  \ifnum \@lastchclass=7 5 \else
   \ifnum \@lastchclass=8 \tw@ \else
    \ifnum \@lastchclass=9 \thr@@
   \else \z@
   \ifnum \@lastchclass = 10 \else
   \edef\@nextchar{\expandafter\string\@nextchar}%
   \@chnum
   \if \@nextchar c\z@ \else
    \if \@nextchar l\@ne \else
     \if \@nextchar r\tw@ \else
   \z@ \@chclass
   \if\@nextchar |\@ne \else
    \if \@nextchar !6 \else
     \if \@nextchar @7 \else
      \if \@nextchar (8 \else
       \if \@nextchar )9 \else
  10
  \@chnum
  \if \@nextchar m\thr@@\else
   \if \@nextchar p4 \else
    \if \@nextchar b5 \else
   \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi
   \fi \fi  \fi  \fi  \fi  \fi  \fi \fi \fi \fi \fi \fi}
\gdef\arraybackslash{\let\\=\@arraycr}
\def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}}
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\def\abbr{}
\def\corr{}
\def\expan{}
\def\gap{}
\def\orig{}
\def\reg{}
\def\ref{}
\def\sic{}
\def\persName{}\def\name{}
\def\placeName{}
\def\orgName{}
\def\textcal#1{{\fontspec{Lucida Calligraphy}#1}}
\def\textgothic#1{{\fontspec{Lucida Blackletter}#1}}
\def\textlarge#1{{\large #1}}
\def\textoverbar#1{\ensuremath{\overline{#1}}}
\def\textquoted#1{‘#1’}
\def\textsmall#1{{\small #1}}
\def\textsubscript#1{\@textsubscript{\selectfont#1}}
\def\textxi{\ensuremath{\xi}}
\def\titlem{\itshape}
\newenvironment{biblfree}{}{\ifvmode\par\fi }
\newenvironment{bibl}{}{}
\newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par }
\newenvironment{citbibl}{}{\ifvmode\par\fi }
\newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi }
\newenvironment{docDate}{}{\ifvmode\par\fi }
\newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi }
\newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par }
\newenvironment{msHead}{\vskip 6pt}{\par}
\newenvironment{msItem}{\vskip 6pt}{\par}
\newenvironment{rubric}{}{}
\newenvironment{titlePart}{}{\par }

\newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}}
\newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}}
\newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}}
\newcolumntype{P}[1]{){\arraybackslash}p{#1}}
\newcolumntype{B}[1]{){\arraybackslash}b{#1}}
\newcolumntype{M}[1]{){\arraybackslash}m{#1}}
\definecolor{label}{gray}{0.75}
\def\unusedattribute#1{\sout{\textcolor{label}{#1}}}
\DeclareRobustCommand*{\xref}{\hyper@normalise\xref@}
\def\xref@#1#2{\hyper@linkurl{#2}{#1}}
\begingroup
\catcode`\_=\active
\gdef_#1{\ensuremath{\sb{\mathrm{#1}}}}
\endgroup
\mathcode`\_=\string"8000
\catcode`\_=12\relax

\usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry}
\usepackage{framed}

\definecolor{shadecolor}{gray}{0.95}
\usepackage{longtable}
\usepackage[normalem]{ulem}
\usepackage{fancyvrb}
\usepackage{fancyhdr}
\usepackage{graphicx}
\usepackage{marginnote}

\renewcommand{\@cite}[1]{#1}


\renewcommand*{\marginfont}{\itshape\footnotesize}

\def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif}

  \pagestyle{fancy}

\usepackage[pdftitle={Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a Power Supplier},
 pdfauthor={}]{hyperref}
\hyperbaseurl{}

	 \paperwidth210mm
	 \paperheight297mm
              
\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{3}
\clubpenalty=8000
\emergencystretch 3em
\hbadness=4000
\hyphenpenalty=400
\pretolerance=750
\tolerance=2000
\vbadness=4000
\widowpenalty=10000

\renewcommand\section{\@startsection {section}{1}{\z@}%
     {-1.75ex \@plus -0.5ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
     {-1.75ex\@plus -0.5ex \@minus- .2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
     {-1.5ex\@plus -0.35ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\large}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
     {-1ex \@plus-0.35ex \@minus -0.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\normalsize}}
\renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}%
     {1.5ex \@plus1ex \@minus .2ex}%
     {-1em}%
     {\reset@font\normalsize\bfseries}}


\def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt}
 \@tempdima 1.5em \begingroup
 \parindent \z@ \rightskip \@pnumwidth 
 \parfillskip -\@pnumwidth 
 \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par
 \endgroup}
\def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
\def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
\def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
\def\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
\@ifundefined{c@section}{\newcounter{section}}{}
\@ifundefined{c@chapter}{\newcounter{chapter}}{}
\newif\if@mainmatter 
\@mainmattertrue
\def\chaptername{Chapter}
\def\frontmatter{%
  \pagenumbering{roman}
  \def\thechapter{\@roman\c@chapter}
  \def\theHchapter{\roman{chapter}}
  \def\thesection{\@roman\c@section}
  \def\theHsection{\roman{section}}
  \def\@chapapp{}%
}
\def\mainmatter{%
  \cleardoublepage
  \def\thechapter{\@arabic\c@chapter}
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \pagenumbering{arabic}
  \setcounter{secnumdepth}{6}
  \def\@chapapp{\chaptername}%
  \def\theHchapter{\arabic{chapter}}
  \def\thesection{\@arabic\c@section}
  \def\theHsection{\arabic{section}}
}
\def\backmatter{%
  \cleardoublepage
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \setcounter{secnumdepth}{2}
  \def\@chapapp{\appendixname}%
  \def\thechapter{\@Alph\c@chapter}
  \def\theHchapter{\Alph{chapter}}
  \appendix
}
\newenvironment{bibitemlist}[1]{%
   \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}%
	}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty \clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}%
  {\def\@noitemerr
    {\@latex@warning{Empty `bibitemlist' environment}}%
    \endlist}

\def\tableofcontents{\section*{\contentsname}\@starttoc{toc}}
\parskip0pt
\parindent1em
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\newenvironment{reflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemsep}{0pt}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\itshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{sansreflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\itemsep}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\upshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{specHead}[2]%
 {\vspace{20pt}\hrule\vspace{10pt}%
  \phantomsection\label{#1}\markright{#2}%

  \pdfbookmark[2]{#2}{#1}%
  \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}%
  }{}
      \def\TheFullDate{1970-01-01 (revised: 01 January 1970)}
\def\TheID{\makeatother }
\def\TheDate{1970-01-01}
\title{Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a Power Supplier}
\author{}\makeatletter 
\makeatletter
\newcommand*{\cleartoleftpage}{%
  \clearpage
    \if@twoside
    \ifodd\c@page
      \hbox{}\newpage
      \if@twocolumn
        \hbox{}\newpage
      \fi
    \fi
  \fi
}
\makeatother
\makeatletter
\thispagestyle{empty}
\markright{\@title}\markboth{\@title}{\@author}
\renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2\p@
\belowdisplayshortskip 4\p@ plus2\p@ minus2\p@
\def\@listi{\leftmargin\leftmargini
               \topsep 2\p@ plus1\p@ minus1\p@
               \parsep 2\p@ plus\p@ minus\p@
               \itemsep 1pt}
}
\makeatother
\fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm}
\fancyhf{} 
\setlength{\headheight}{14pt}
\fancyhead[LE]{\bfseries\leftmark} 
\fancyhead[RO]{\bfseries\rightmark} 
\fancyfoot[RO]{}
\fancyfoot[CO]{\thepage}
\fancyfoot[LO]{\TheID}
\fancyfoot[LE]{}
\fancyfoot[CE]{\thepage}
\fancyfoot[RE]{\TheID}
\hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true}
\fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}}

\date{}
\usepackage{authblk}

\providecommand{\keywords}[1]
{
\footnotesize
  \textbf{\textit{Index terms---}} #1
}

\usepackage{graphicx,xcolor}
\definecolor{GJBlue}{HTML}{273B81}
\definecolor{GJLightBlue}{HTML}{0A9DD9}
\definecolor{GJMediumGrey}{HTML}{6D6E70}
\definecolor{GJLightGrey}{HTML}{929497} 

\renewenvironment{abstract}{%
   \setlength{\parindent}{0pt}\raggedright
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
   \textcolor{GJBlue}{\large\bfseries\abstractname\space}
}{%   
   \vskip8pt
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
}

\usepackage[absolute,overlay]{textpos}

\makeatother 
      \usepackage{lineno}
      \linenumbers
      
\begin{document}

             \author[1]{Marcos Antonio Lopes Freixo  Filho}

             \author[2]{Ana Julia  Righetto}

\renewcommand\Authands{ and }

\date{\small \em Received: 1 January 1970 Accepted: 1 January 1970 Published: 1 January 1970}

\maketitle


\begin{abstract}
        


There are several technologies in the area of engineering, new and old, currently used with the objective of reducing the impacts of meteorological phenomena on electrical energy distribution networks. However, they still constitute one of the biggest challenges for companies around the world, the mitigation of the number of interruptions in the energy supply, in the face of the weather. This research work seeks something similar, which is to offer another analysis, from the point of view of data science, studying the correlation between the meteorological phenomena recorded by the Climatempo Monitoring and Alert System and the interruptions in the energy supply in two main cities of Maranhão State, with unknown causes, recorded by the control system used by the integrated operations center of the local power supplier, in an attempt to find some kind of fragility in the electrical system, not yet identified, that allows its improvement.

\end{abstract}


\keywords{DEC, FEC, clima, transformadores de distribuição.}

\begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) 
\textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark}
\end{textblock*}

\begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) 
\textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}}
\end{textblock*}


\begin{textblock*}{10cm}(1.05cm,3cm)
{{\textit{CrossRef DOI of original article:}} \underline{}}
\end{textblock*}\let\tabcellsep& 	 	 		 \par
expurgado dos resultados das distribuidoras de energia nos cálculos da agência reguladora, e os fenômenos climáticos, ou meteorológicos que por ventura derrubem o sistema, como descargas elétricas, ventanias, entre outros fatores.\par
Segundo \hyperref[b7]{Silva (2013)}, em nosso país, aproximadamente 70\% dos desligamentos na transmissão e 40\% na distribuição, são provocados por raios, sendo o número de transformadores queimados por descargas elétricas atmosféricas em torno de 40\%. Tais números causam um impacto considerável no fornecimento de energia, o que pode ser constatado pela alta correlação entre a incidência de descargas elétricas e os índices de qualidade da maioria das empresas do setor elétrico.\par
No trabalho \hyperref[b13]{de Wronscki et al (2003)}, após a utilização de modelos estatísticos, foram encontradas evidências de maior número de atendimentos e maior tempo médio de atendimento nas chamadas emergenciais de interrupção no fornecimento para a condições meteorológicas de vento e chuva, enquanto houve um menor número para a condição temperatura.\par
A distribuidora de energia do Maranhão, que durante anos permaneceu entre as três melhores empresas do país no quesito qualidade no fornecimento de energia, despencou em 2021 para os últimos lugares, de acordo com o último ranking divulgado pela [ANEEL] (2022). A empresa acredita que isto não se deve aos investimentos realizados na rede, que na verdade, ano após ano são cada vez maiores, mas se deve a um ano atípico de muitos registros de fortes precipitações, enchentes, ventanias, descargas elétricas atmosféricas, entre outros fenômenos climáticos.\par
Este trabalho tem como objetivo verificar e quantificar a correlação entre os fenômenos meteorológicos e a quantidade de ocorrências de interrupção no fornecimento de energia elétrica registradas pelo centro de operações integradas da distribuidora energética do Estado do Maranhão, para os municípios de São Luís e Imperatriz. 
\section[{II. Material e Métodos}]{II. Material e Métodos}\par
O princípio que irá nortear este trabalho de pesquisa, a formar a base de dados final para estudos, ou seja, após a fase de "data wrangling", é que esta base de dados reflita muito bem o fenômeno a ser estudado, que é a interrupção no fornecimento de energia, que seja capaz de explicar o preceito principal de que, apesar das dificuldades enfrentadas particularmente no país, em suma, o fornecimento de energia elétrica passa mais tempo em continuidade do fornecimento, do que interrompido, ou suspenso.\par
Para se ter uma dimensão sobre tal fato, o indicador do [DEC] em 2021, valor apurado na média para todo o Brasil, segundo o Governo Federal (2022), ficou em 11,84 horas. Em outras palavras, de um total de 8.760 horas do ano de 2021, em 99,86\% do tempo a energia esteve em pleno fornecimento, sem interrupção, no Brasil inteiro, como pode ser constatado na Figura 1. 
\section[{Fonte: Governo Federal (2022) Figura 1: Recorte do Resultado do [DEC] em 2021}]{Fonte: Governo Federal (2022) Figura 1: Recorte do Resultado do [DEC] em 2021}\par
Ainda que o resultado para o Maranhão tenha ficado acima deste valor, registrando 29 horas, estando 99,67\% do tempo em pleno fornecimento, conforme informado pelo distribuidora de energia em seu site de relacionamento com os investidores, em março de 2022, este número continua apontando para um fato supra relevante: a base de dados deve dar alicerce para realizar o estudo de um fenômeno raro, ou seja, que além dos dados sobre as interrupções no fornecimento de energia, a base deve possuir cada hora em que o fenômeno não aconteceu, o período em que o fornecimento de energia não foi interrompido.\par
Seguindo estra premissa, este trabalho de pesquisa experimental, vem trabalhar a hipótese da existência de uma correlação entre os fenômenos meteorológicos, e a quantidade de vezes em que o sistema elétrico teve sua distribuição interrompida, porém somente dos dados sinalizados como causa desconhecida, atendida por transformadores de distribuição, em dois municípios do Estado do Maranhão.\par
O centro de operações integradas[COI], setor responsável por manter o sistema de distribuição de energia elétrica em pleno atendimento para a população maranhense, possui em seus controles, uma base de registros de interrupções do fornecimento de energia, e seus equipamentos afetados, chamados de componentes, com informações que caracterizam em detalhe a sua tipificação, informações a princípio muito relevantes para esta pesquisa para verificarmos se há um padrão de componentes mais suscetíveis, ou não, às intemperes climáticas.\par
Com autorização da empresa, foram colhidos dados dos últimos doze meses, compreendidos entre agosto de 2021 e julho de 2022, somente dos equipamentos de transformadores de distribuição, já que toda interrupção de energia, que afete o bem-estar dos consumidores, acaba passando por este tipo de equipamento, segundo o especialista do próprio setor.\par
Outro fato importante é que, com a ajuda deste especialista, foram excluídas desta amostra de dados, todas as interrupções de energia derivadas da intervenção humana, como acidentes de trânsito, desligamentos programados, sobrecargas do equipamento causadas pelo consumo, não realização de poda de árvores abaixo da rede elétrica, para citar alguns, o que deixa a base de informações mais limpa para a realização deste estudo, já que todas as interrupções citadas têm a sua causa conhecida, e registrada no sistema.\par
Coletada a primeira parte dos dados, a segunda parte foi fornecida junto com a área conhecida como geoprocessamento, e lá foram coletadas as informações a respeito do clima, para as mesmas cidades. A companhia de energia possui um contrato de prestação de serviços com a empresa Climatempo, e através do Sistema de Monitoramento e Alertas do Climatempo, com a autorização da distribuidora, foram coletadas também informações do mesmo período da base coletada anteriormente, referentes a pluviometria, descargas elétricas atmosféricas, rajadas de vento, velocidade média do vento ao longo do dia, temperatura média, máxima e mínima.\par
As variáveis possuem precisões diferentes para relatar o período em que aconteceram, variando entre a precisão da hora, minuto e segundo em que ocorreram, como os dados de interrupções de energia, e os dados das descargas elétricas atmosféricas, e até informação que trata de resultado válido pelo dia inteiro, como velocidade do vento, e as informações sobre temperatura. Todas foram colocadas na mesma escala de tempo, demonstrando os fatos ocorridos a cada uma hora do dia, durante 365 dias. Assim, foram criadas quatro variáveis para identificar o tempo: ano, mês, dia e hora.\par
Um fato importante é que tal base neste momento contém apenas informações relativas às interrupções, e como estava o clima naquele exato momento, entretanto, como já mencionado, ainda que uma distribuidora no território nacional tenha os piores indicadores de continuidade no fornecimento de energia, mesmo assim o sistema elétrico desta distribuidora hipotética estaria em pleno funcionamentona esmagadora parte do tempo. Desta forma, foram reproduzidos os dados que representam o período em que a energia estava emplenadistribuição, para cada tipo diferente de transformador de distribuição, trazido na amostra de interrupções, com o mesmo período de tempo da amostra, em dias e horas, e devidamente cruzada com a base climatológica, para se conhecer as exatas condições do tempo neste mesmo exato período. Para diferenciar as duas informações na base de dados, foi criada uma variável chamada de "interrupções", do tipo numérica, atribuindo a quantidade "1" para as observações relativas a uma interrupção no fornecimento de energia, e "0" para as observações que representam o período em que a energia estava em pleno funcionamento.\par
Cabe salientar que foram excluídas desta base de dados de observações em que a energia estava em plena distribuição, as observações que coincidem temporalmente com a base de interrupções, afinal, exatamente neste exato período, a energia fortuitamente estava suspensa em seu fornecimento. Também, no intuito de explorar uma outra possibilidade que pode ser relevante para a pesquisa, foi criada uma variável chamada "raio", do tipo binária, para indicar a presença de raios naquela exata hora relada atribuindo o valor 1, e 0 quando o fato não havia ocorrido, o que é diferente da variável "descargas\textunderscore atm", que apresenta a quantidade de descargas elétricas atmosféricas dentro daquela mesma observação.\par
Finalmente com as duas bases de dados unidas (com e sem interrupções), encerrando a etapa de"data wrangling" da coleta de dados, para que se garanta a aleatoriedade da informação, quanto ao objeto desta pesquisa experimental, e permitindo analisar se este estudo se trata de um caso em que o melhor modelo estatístico deva levar em conta o estudo de um fenômeno raro, considerando o fato de que energia elétrica passa mais tempo em plena distribuição, do que suspensa.\par
Por fim, em razão do tamanho da base de dados para este estudo, se fossem considerados todos os municípios do Estado do Maranhão, ainda que utilizando dados relativos ao período de apenas um ano, foram escolhidos apenas dois dos principais municípios do Estado: São Luís, e Imperatriz. O estudo poderia se limitar a capital do Estado (São Luís), por seu tamanho e importância. Todavia, a cidade de Imperatriz possui características completamente diferentes, no que se diz respeito ao clima, relevo, economia, ou outras não citadas, o que torna ainda mais interessante checar neste estudo, se há semelhanças nos resultados obtidos entre os dois municípios, o que poderia levantar outras hipóteses, para outros trabalhos futuros. 
\section[{III. Modelo de Regressão de Poisson}]{III. Modelo de Regressão de Poisson}\par
A metodologia estatística que foiadotada neste estudo foia aplicação do modelo de regressão de Poisson. Os modelos de regressão Poisson e o modelo de regressão binomial negativo fazem parte do que éconhecido por modelos de regressão para dados de contagem, que têm porobjetivo analisar o comportamento, em função de variáveis preditoras, de umadeterminada variável dependente que se apresenta na forma quantitativa, comvalores discretos e não negativos, exatamente como o objeto de estudo deste trabalho, onde a princípio não é possível fracionar o evento "interrupção do sistema elétrico" (não considerase neste estudo os casos em que a energia esteve em oscilação quanto ao fornecimento), tendo sempre sua quantificação como números discretos e positivos.\par
Considerando o que afirma Fávero e Belfiore (2022, p.149), ao mencionar que no caso de ocorrência de fenômenos raros, com baixa probabilidade de sucesso (p ? 0), sob determinada exposição (unidade temporal, espacial, social, etc.)., como no caso desta pesquisa experimental, em determinado intervalo de tempo, o modelo de distribuição Poisson se torna o prioritário a ser explorado. 
\section[{Expressão do modelo de Poisson:}]{Expressão do modelo de Poisson:}ln?? i ? = ln?? poisson i ? = ? + ? 1 . X 1i + ? 2 . X 2i + . . . + ? k . X ki form. (\textbf{1})\par
Em que ln?? i ? ou ln?? poisson i ? é definido a partir de uma função de ligação canônica que é proposta com base na definição dos logaritmos das funções de verossimilhança oriunda da função densidade/ probabilidade da distribuição Poisson, ?representa os termos do intercepto, ? k são os coeficientes de cada variável explicativa e correspondem aos parâmetros a serem estimados, eX ki são as variáveis explicativas. Fórmula da probabilidade da distribuição Poisson:p(Y i = m) = e -? i .? i m m! form. (2)\par
Em que p(Y i = m) é afunção da probabilidade de um evento, no modelo de regressão de Poisson, determinada por uma observação i (i = 1, 2, ..., n), "n" é o tamanho da amostra, e possui a probabilidade de ocorrência de uma contagem "m", em uma determinada exposição (período, área, região, entre outros exemplos), em que "?" é o número esperado de ocorrências ou a taxa média estimada de incidência do fenômeno em estudo para uma dada exposição A título de comparação de melhor modelo, também foi aplicado neste estudo o modelo de regressão binomial negativo, também conhecido como Poisson Gama, outro modelo de contagem aqui citado. A diferença básica entre os dois modelos,é que o modelo binomial negativo é mais utilizado para os casos de superdispersão, ou inflação de zeros da amostra, onde a variância da variável dependente é maior que sua própria média. Segundo Fávero e Belfiore (2022, p.696), a superdispersão é comumente gerada pela presença de maior heterogeneidade nos dados entre observações da amostra. 
\section[{Expressão do modelo Binomial Negativo:}]{Expressão do modelo Binomial Negativo:}ln?? i ? = ln?? bneg i ? = ? + ? 1 . X 1i + ? 2 . X 2i + . . . + ? k . X ki form. (\textbf{3})\par
Em que ln?? i ? ou ln?? bneg i ? é definido a partir de uma função de ligação canônica que é proposta com base na definição dos logaritmos das funções de verossimilhança oriunda da função densidade/ probabilidade da distribuição Poisson Gama, ou Binomial Negativa, ?representa os termos do intercepto, ? k são os coeficientes de cada variável explicativa e correspondem aos parâmetros a serem estimados, e X ki são as variáveis explicativas.\par
Fórmula da probabilidade da distribuição Binomial Negativa:p(Y i = m) = ? ? .m i ? -1 .e -m i .? (?-1)! form. (\textbf{4})\par
Em que p(Y i = m) é a função de probabilidade no modelo binomial negativo determinada por uma observação i (i = 1, 2, ..., n), "n" é o tamanho da amostra, onde ? é chamado de parâmetro de forma e deve ser maior que zero, e ? é chamado de parâmetro de taxa de decaimento, também maior que zero.\par
A hipótese nula para este estudo será a baixa ou nenhuma correlação entre os fenômenos meteorológicos e as interrupções no fornecimento de energia elétrica em uma distribuidora, abrindo a possibilidade de estudar outros fatores, como a fadiga dos equipamentos, vida útil da rede de distribuição, ou mais fatores, que possam ter preponderância maior na suspensão ocasional da distribuição energética de eletricidade.\par
Ainda, foram feitas análises de correspondências múltiplas, que se trata de uma técnica multivariada que possibilita a investigação de associação com mais de duas variáveis categóricas, conformeFávero e Belfiore (2022, p.463), no intuito de perceber a relação de interdependência entre as variáveis utilizadas neste estudo, trabalhando somente com as observaçõesonde ocorreram as interrupções do fornecimento de energia elétrica, nas duas amostras os dados.\par
As técnicas de análise de correspondência são métodos de representação de linhas e colunas de tabelas cruzadas de dados, como coordenadas em um gráfico, chamado de mapa perceptual, a partirdo qual se podem interpretar as similaridades e diferenças de comportamento entre variáveis, e entre categorias.\par
Essa técnica tem como principal objetivo avaliar a significância dessas similaridades, reduzindo suas dimensões, de modo a determinar coordenadasdas categorias com base na distribuição dos dados em tabelas cruzadaspara,a partir dessas coordenadas, construirmapas perceptuais, que nada mais são que diagramas de dispersão que representam as categorias das variáveisna forma de pontos em relação a eixos de coordenadas ortogonais, fazendo delas,na realidade, mapas de categorias.\par
O método para realizar a análise consiste em transformar as variáveis qualitativas em variáveis binárias, obtendo daí uma matriz binária. Com base nesta matriz binária, é obtido a inércia principal total na [ACM]. Supondo que a matriz binária seja semelhante a uma tabela de contingência da análise de correspondência, é possível obter a inércia principal parcial das dimensões, seus autovalores, autovetores, e posterior coordenadas desta matriz.\par
Outro método alternativo é a combinação em uma única matriz, as tabelas de contingências com os cruzamentos de todos os pares de variáveis. Essa matriz resultante, quadrada e simétrica, é conhecida por matriz de Burt, como explica \hyperref[b0]{Baltar (2005)}.\par
Só participaram da [ACM] as variáveis que apresentaram associação estaticamente significativa, com pelo menos uma variável contida na análise, por meio do teste do X 2 . Caso alguma delas não apresentasse associação com outras, estas ficaram excluídas da análise de correspondência. Ao elaborar o mapa perceptual, são representadas no mapa as coordenadas dasdimensões que apresentam a inércia principal parcial maior do que a média dainércia principal total por dimensão.\par
A aplicação dos modelos estatísticos se deu através da linguagem R, utilizando o software RStudio versão 1.4.1106. A base final foi exportada para o formato de planilha, e tratada em banco de dados SQL Server. 
\section[{IV. Resultados e Discussão a) São Luís}]{IV. Resultados e Discussão a) São Luís}\par
Foram carregados os dados da amostra de São Luís, coletada na distribuidora de energia no R com um total de 272.151 observações e 12 variáveis. Os dados podem ser vistos através da função "glimpse()", conforme apêndice.\par
Para reforçar essa possibilidade conceitual, de que o objeto de estudo deste trabalho de pesquisa, se trata de um fenômeno raro, foi realizado o diagnóstico preliminar para observação de eventual igualdade entre a média e a variância da variável dependente INTERRUPCAO, cujos resultados foram respectivamente 0,0180268 e 0,0177019, valores muito próximos, o que indica fortemente que de fato, o modelo estatístico mais indicado para explicar a contagem de interrupções de energia elétrica da amostra de dados de São Luís, seja um modelo de regressão Poisson.\par
Utilizando a função "glm()", fez-se a estimação deste modelo, colocando a variável INTERRUPCAO como a variável dependente, e as outras onze variáveis, já citadas, como variáveis explicativas.\par
O primeiro resultado apurado com base no modelo de regressão de Poisson, confirma a hipótese desta pesquisa experimental, de que existe de fato uma correlação entre os fenômenos meteorológicos e as vezes em que o fornecimento de energia elétrica foi interrompido aos consumidores de São Luís, pois pelo menos um beta de uma variável com informações da meteorologia no momento em que acontece a interrupção, teve o seu p-valormenor que 0.05, ou seja, podemos admitir o modelo estatísticos perfeitamente correlacionado, para fins preditivos, como vê-se na Tabela 1, através da função "summary()". Para fins de comparação, e se aproximar da certeza de que o modelo estatístico não se trata de uma distribuição Binomial Negativa, corroborando com o diagnóstico preliminar que compara a média e a variância da base coletada de dados de São Luís, fezse necessário carregar o modelo binomial negativo desta mesma base, utilizando a função "glm.nb()".\par
Segundo Bussab e Morettin (2005), o princípio da verossimilhança explica que se deve escolher aquele valor do parâmetro desconhecido que maximiza a probabilidade de obter a amostra particular observada, ou seja, o valor que torna aquela amostra a "mais provável". Esse princípio foi enunciado por Ronald Aylmer Fisher pela primeira vez em 1912 e, em 1922, deu-lhe forma mais completa, introduzindo a expressão Log-Likelihood(verossimilhança).\par
Em outras palavras, o método de máxima verossimilhança consiste em estimar o parâmetro ? por uma amostra aleatória X1, ..., Xn independente e igualmente distribuída de X. Na Figura 2 é feita a comparação entre os dois modelos estatísticos, comparando o Log-Likelihoodde ambos os modelos estatísticos aplicados na amostra de São Luís. 
\section[{Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a}]{Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a}\par
Power Supplier Apenas para exemplificar a importância desta pesquisa experimental, tomando o modelo estatístico fundamentadamente adotado, através da função "predict", é possível estimar em aproximadamente 36,6\%, a probabilidade de ocorrer uma interrupção do fornecimento de energia elétrica de um transformador de distribuição com três fases conectadas (ABC), de carga de 112,5Kva, sob uma chuva de 15 milímetros, com uma temperatura mínima de 22ºC, na presença de raios, rajadas de vento de 30Km/h,conforme pode ser visto no apêndice deste trabalho. 
\section[{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version}]{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version}\par
Ainda, cabe neste trabalho explorar a relação que cada uma destas variáveis que explicam as interrupções no fornecimento de energia, na amostra utilizada no modelo de Poisson, se relaciona-se entre si. Para isto, foi gerada um segundo "dataframe", a partir da base denominada de final.\par
Em um novo "wrangling" dos dados, para que se possa realizar uma análise de correspondência, foram criadas novas variáveis para minimizar a quantidade de rótulos a serem exibidos num mapa perceptual que será feito com a análise de correspondência múltipla  {\ref [ACM]}.\par
A variável PLUVIOMETRIA, foi substituída pela variável PLUV\textunderscore TIPO, agora como uma variável categórica, contendo uma classificação da chuva baseada na sua intensidade medida na quantidade de milímetros registrados naquela hora. A classificação da intensidade da chuva obedeceu aos critérios adotados pela Organização Meteorológica Mundial -OMM (2018), e para nossa amostra, se dividiu entre fraca, moderada e forte.\par
A variável RAJADA foi substituída pela variávelRAJADA\textunderscore TIPO, e suas categorias foram classificadas segundo a escala Beaufort, também disciplinada pela OMM, escala baseada na velocidade do vento registrada conforme Tabela 2.A variável VENTO foi excluída por possíveis relações de multicolinearidade na amostra com outras variáveis, como já mostradoneste trabalho.\par
A variável RAIO foi mantida, mas para fins de exibição dos rótulos de dados no mapa perceptual da [ACM], os valores "0" e "1", foram substituídas respectivamente pelas categorias "RAIO\textunderscore NAO" e "RAIO\textunderscore SIM". A variável DESCARGAS\textunderscore ATM foi excluída, sem substituta, por não ter sido encontrado quaisquer classificações advindas de instituições oficiais ou norma científica, nas pesquisas deste trabalho, que definissem a intensidade de descargas elétricas atmosféricas. 
\section[{Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a}]{Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a}\par
Power Supplier 
\section[{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version I Year 2023}]{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version I Year 2023}\par
Tabela De forma semelhante, as três variáveis de temperatura foram substituídas por três novas variáveis categóricas, TEMP\textunderscore MIN\textunderscore FAIXA, TEMP\textunderscore MAX\textunderscore FAIXA e TEMP\textunderscore MED\textunderscore FAIXA, onde de forma arbitrária, por também não ter encontrado literatura semelhante às outras variáveis anteriormente citadas, que agrupam as temperaturas por faixa.\par
Por último, para fins demonstrativos, destinados a orientar estudos posteriores a este trabalho, foram selecionados apenas três transformadores de distribuição na amostra, em razão de serem os mais utilizados pela Companhia, da igual forma, a fase escolhida a foi a "ABC".\par
Através do comando "mca()", obtemos a análise de correspondência múltipla da amostra. Foram obtidas 15 dimensões, e a inércia total explicada se deu conforme exposto no Figura 3. 
\section[{Fonte: Resultados originais da pesquisa Figura 3: Inércia Total Explicada da [ACM] da Amostra}]{Fonte: Resultados originais da pesquisa Figura 3: Inércia Total Explicada da [ACM] da Amostra}\par
Por fim, para fins de exemplo de estudos que possam ser explorados pelas áreas de engenharia, com outros equipamentos, em outras regiões do Estado, ou mesmo em subestações, linhas de transmissão, tem-se o mapa perceptual da Figura 4. 
\section[{Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a}]{Correlations between Meteorological Phenomena and Interruptions in the Supply of Electric Energy in a}\par
Power Supplier 
\section[{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version I Year 2023}]{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version I Year 2023}\par
Fonte: Resultados originais da pesquisa Figura 4: Mapa Perceptual da [ACM] da Amostra Para Exemplo de Estudo Em uma breve análise do mapa, é possível perceber relações importantes entre algumas variáveis, que configuram as situações mais comuns para interrupções no fornecimento de energia elétrica. Vê-se no mapa que os transformadores de distribuição do exemplo da amostra, estão mais vulneráveis a um dia sem descargas elétricas atmosféricas, com rajada de vento calma (velocidade do vento próxima a zero), com uma temperatura média entre 28ºC e 30ºC, e uma chuva fraca (menos de 2,5 mm/h). Mas também interessante notar o quanto uma rajada de vento classificado como Brisa Forte, na Escala de Beaufort, com velocidade do vento acima dos 38Km/h, pode ser mais preponderante a derrubar o sistema elétrico, do que uma chuva moderada.\par
Entretanto, antes de quaisquer conclusões antecipadas, cabe verificar se as mesmas distâncias se mantêm quando vistas em gráficos de três dimensões.\par
que pode se tornar um trabalho de rotina a se adotar em uma distribuidora de energia elétrica. 
\section[{b) Imperatriz}]{b) Imperatriz}\par
Após a carga dos dados da cidade de Imperatriz, segunda principal cidade do Estado do Maranhão em importância econômica, para efeitos de comparação com os resultados encontrados para a capital do Estado, no início deste trabalho, serão seguidos exatos mesmos procedimentos tomados no estudo anterior. A base de dados correspondente a um ano de informações, sobre todos os transformadores de distribuição instalados, contém 280.442 observações.\par
Do mesmo modo como aconteceu com os dados da primeira amostra, para o estudo com os dados do município de Imperatriz, a média e a variância da variável dependente INTERRUPCAO encontradas foram respectivamente 0.0056233, e 0.0055917, números muito próximos, que mostram que o diagnóstico preliminar para observação de eventual igualdade,também aponta para um estudo de um fenômeno raro, mais uma vez, corroborando a tese da Distribuição de Poisson como melhor modelo estatístico para esta nova amostra. Também, para fins de comparação com o estudo anterior, o mesmo exemplo de "predict" feito para a mostra de São Luís, foi feito para a amostra de Imperatriz, com as mesmas variáveis, e o resultado foi bem diferente entre as duas cidades. Enquanto a probabilidade de uma interrupção no fornecimento de energia elétrica na primeira amostra foi de 36,6\%, nas mesmas condições do tempo, e mesmas características de transformadores, o resultado foi uma probabilidade de apenas 5,16\% do sistema elétrico ser interrompido pelos mesmos fatores (sete vezes menor). 
\section[{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version}]{Global Journal of Research in Engineering ( F ) XXIII Issue IV Version}\par
Ao realizar a análise de correspondências múltiplas, mais diferenças são observáveis entre os dois resultados. Em razão das disparidades climáticas entre as duas cidades, mais categorias foram adicionadas para os cálculos, gerando 19dimensões para a análise, três dimensões a mais que na amostra de São Luís, conforme o Figura 5, com a inércia total explicada da análise de correspondência múltipla. De acordo com Lima e Ferreira (2013), separadas por aproximadamente 600km, uma da outra, as duas cidades de fato possuem climas bastante divergentes. São Luís, que fica ao norte do Estado, é uma cidade litorânea, mais sujeita a ventos fortes, e temperaturas mais equilibradas ao longo do ano, do que a cidade de Imperatriz, ao sul, com um tipo de clima misto, ora seco como o cerrado e ao sertão, ora chuvoso como o clima amazônico, por fazer fronteira entre estes dois biomas, além de ventos mais fracos, e temperaturas mais extremas do que na capital. 
\section[{Fonte: Resultados originais da pesquisa}]{Fonte: Resultados originais da pesquisa}\par
Ao analisar o mapa perceptual da amostra de Imperatriz na Figura 6, é interessante perceber que a chuva violenta (com mais de 50mm/h) tem pouquíssima relação com as interrupções no fornecimento de energia de Imperatriz. Isto pode mostrar o quanto o sistema elétrico na região foi bem preparado pela distribuidora para suportar as fortes intempéries do período chuvoso.\par
Já a chuva moderada (entre 2,5 a 10 mm/h) se mostra muito mais relevante para a interrupção do fornecimento de energia do que em São Luís. A rajada de brisa forte, que em São Luís se mostra importante, em Imperatriz é a rajada de brisa fraca (ventos repentinos de 12 a 19km/h) que se mostra mais preponderante.\par
Tal mapa pode auxiliar as áreas de engenharias a encontrar soluções que mitiguem a quantidade de suspensões do fornecimento de energia elétrica, segundo as condições apresentadas. V. Conclusão É possível afirmar, através deste trabalho de pesquisa realizado com dados reais de uma distribuidora de energia elétrica do Estado do Maranhão, com duas de suas principais cidades, que há de fato, com um nível de significância de 5\%, fortes correlações associadas aos fenômenos meteorológicos e as interrupções no fornecimento de energia, com causas desconhecidas, para ambas as amostras, com ressalvas que indicam que o estudo precisa ser de fato específico por região, e que, para fins preditivos, o modelo de distribuição de Poisson se encaixa perfeitamente. 
\section[{Fonte: Resultados originais da pesquisa}]{Fonte: Resultados originais da pesquisa}\par
Ainda, o estudo possibilita às distribuidoras, que possam trabalhar de modo preventivo, identificando as fragilidades do sistema elétrico quanto ao tema, e explorar ações que possam reduzir a frequência com que acontecem as interrupções, analisando não só transformadores de distribuição, mas equipamento por equipamento, ou até mesmo estruturas mais robustas, como subestações, através modelo de Poisson, e de mapas perceptuais de análise de correspondência múltipla, de modo a melhor o indicador do  {\ref [FEC]}.\par
Por fim, também é possível criar um sistema de alertas que envie para o Centro de Operações Integradas, que tipo de equipamento, onde, e a que probabilidade, este possa ter sua distribuição de energia interrompida, de acordo com a previsão do tempo, combinado ao Sistema de Monitoramento de Alertas do ClimaTempo, com previsões em tempo real, podendo antecipadamente ao fenômeno,deslocar equipes de plantão a lugares estratégicos, para que estas possam tomar ações que reestabeleçam o fornecimento de energiacom mais celeridade. Desta forma, melhorando também o indicador do  {\ref [DEC]}.\par
Por fim, este estudo sugere que outras variáveis sejam acrescentadas para um estudo semelhante, como a vida útil da rede, transformador de distribuição, ou o que queira ser estudado. Também seria interessante acrescentar o fabricante, e lote adquirido, para identificar se uma linha de materiais comprados pela distribuidora não possui mais fragilidade à fenômenos meteorológicos. 
\section[{Agradecimento}]{Agradecimento}\par
A Deus por ter me dado toda a força deste mundo para chegar até aqui, com tantas pedras pelo caminho, incertezas que quase me levaram a desistir, dificuldades financeiras, mas principalmente a difícil arte de revezar o tempo com a família, amigos e trabalho, tarefa das mais complicadas que já passei nos últimos tempos. Só eu e Ele sabíamos o que estávamos construindo. Agradecer em especial aos amigos Evandro Lima Meireles, da Gerência Corporativa de Geoprocessamento, e Eliezer Silva de Oliveira, da Gerência Corporativa de Operações, que com seus esforços na fase de mineração de dados, contribuíram enormemente com este trabalho, além também das pessoas que de alguma forma compreenderam este projeto de vida, e minimamente colaboraram para a conclusão deste trabalho.\begin{figure}[htbp]
\noindent\textbf{62}\includegraphics[]{image-2.png}
\caption{\label{fig_0}I Year 2023 6 ©Figura 2 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{5}\includegraphics[]{image-3.png}
\caption{\label{fig_1}Fonte: Resultados originais da pesquisaFigura 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6}\includegraphics[]{image-4.png}
\caption{\label{fig_2}Figura 6 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-5.png}
\caption{\label{figure5}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-6.png}
\caption{\label{figure6}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-7.png}
\caption{\label{figure7}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-8.png}
\caption{\label{figure8}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-9.png}
\caption{\label{figure9}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-10.png}
\caption{\label{figure10}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-11.png}
\caption{\label{figure11}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1} \par 
\begin{longtable}{P{0.23631664853101197\textwidth}P{0.07121871599564744\textwidth}P{0.07676822633297062\textwidth}P{0.008324265505984765\textwidth}P{0.1294885745375408\textwidth}P{0.0018498367791077257\textwidth}P{0.10405331882480957\textwidth}P{0.12023939064200218\textwidth}P{0.10174102285092491\textwidth}}
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA15}\tabcellsep \multicolumn{2}{l}{-3,0510}\tabcellsep 0,1446\tabcellsep -21,0960\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA150}\tabcellsep \multicolumn{2}{l}{-0,4834}\tabcellsep 0,0484\tabcellsep -9,9800\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA1500}\tabcellsep \multicolumn{2}{l}{-7,1696}\tabcellsep 1,0024\tabcellsep -7,1520\tabcellsep 0,0000\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA1750}\tabcellsep \multicolumn{2}{l}{-6,0710}\tabcellsep 0,5833\tabcellsep -10,4080\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA2000}\tabcellsep \multicolumn{2}{l}{-6,4765}\tabcellsep 0,7120\tabcellsep -9,0960\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA225}\tabcellsep \multicolumn{2}{l}{-2,9794}\tabcellsep 0,1396\tabcellsep -21,3490\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA25}\tabcellsep \multicolumn{2}{l}{-7,1520}\tabcellsep 1,0015\tabcellsep -7,1410\tabcellsep 0,0000\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA250}\tabcellsep \multicolumn{2}{l}{-7,1690}\tabcellsep 1,0024\tabcellsep -7,1520\tabcellsep 0,0000\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA2600}\tabcellsep \multicolumn{2}{l}{-7,1696}\tabcellsep 1,0024\tabcellsep -7,1520\tabcellsep 0,0000\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA30}\tabcellsep \multicolumn{2}{l}{-2,9908}\tabcellsep 0,1431\tabcellsep -20,9010\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA300}\tabcellsep \multicolumn{2}{l}{-3,9313}\tabcellsep 0,2119\tabcellsep -18,5510\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA3000}\tabcellsep \multicolumn{2}{l}{-7,1696}\tabcellsep 1,0024\tabcellsep -7,1520\tabcellsep 0,0000\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA350}\tabcellsep \multicolumn{2}{l}{-6,4765}\tabcellsep 0,7120\tabcellsep -9,0960\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA3500}\tabcellsep \multicolumn{2}{l}{-7,1690}\tabcellsep 1,0024\tabcellsep -7,1520\tabcellsep 0,0000\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA45}\tabcellsep \multicolumn{2}{l}{-1,2695}\tabcellsep 0,0760\tabcellsep -16,7130\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \tabcellsep KVA5\tabcellsep \multicolumn{2}{l}{-5,9503}\tabcellsep 0,4556\tabcellsep -13,0590\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA500}\tabcellsep \multicolumn{2}{l}{-4,6847}\tabcellsep 0,3004\tabcellsep -15,5930\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA612,5}\tabcellsep \multicolumn{2}{l}{-6,4765}\tabcellsep 0,7120\tabcellsep -9,0960\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA725}\tabcellsep \multicolumn{2}{l}{-6,4765}\tabcellsep 0,7120\tabcellsep -9,0960\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA75}\tabcellsep \multicolumn{2}{l}{-0,1176}\tabcellsep 0,0463\tabcellsep -2,5430\tabcellsep 0,0110\\
\multicolumn{4}{l}{MEDIA\textunderscore CLIENTES}\tabcellsep \multicolumn{2}{l}{-0,0006}\tabcellsep 0,0004\tabcellsep -1,3650\tabcellsep 0,1722\\
\tabcellsep \multicolumn{3}{l}{PLUVIOMETRIA}\tabcellsep 0,0154\tabcellsep \tabcellsep 0,0025\tabcellsep 6,0680\tabcellsep 0,0000\\
\tabcellsep \multicolumn{3}{l}{TEMP\textunderscore MIN}\tabcellsep \multicolumn{2}{l}{-0,0720}\tabcellsep 0,0149\tabcellsep -4,8420\tabcellsep 0,0000\\
\tabcellsep \multicolumn{3}{l}{TEMP\textunderscore MAX}\tabcellsep \multicolumn{2}{l}{-0,0200}\tabcellsep 0,0132\tabcellsep -1,5130\tabcellsep 0,1303\\
\tabcellsep \multicolumn{3}{l}{TEMP\textunderscore MEDIA}\tabcellsep NA\tabcellsep \tabcellsep NA\tabcellsep NA\tabcellsep NA\\
\multicolumn{4}{l}{DESCARGAS\textunderscore ATM}\tabcellsep 0,0015\tabcellsep \tabcellsep 0,0007\tabcellsep 2,0870\tabcellsep 0,0369\\
\tabcellsep \tabcellsep \tabcellsep RAIO\tabcellsep 0,3236\tabcellsep \tabcellsep 0,0633\tabcellsep 5,1090\tabcellsep 0,0000\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{RAJADA}\tabcellsep 0,0042\tabcellsep \tabcellsep 0,0010\tabcellsep 4,0830\tabcellsep 0,0000\\
\tabcellsep \multicolumn{3}{l}{VENTOS\textunderscore KMH}\tabcellsep 0,0021\tabcellsep \tabcellsep 0,0015\tabcellsep 1,4400\tabcellsep 0,1499\\
\multicolumn{6}{l}{Fonte: Resultados originais da pesquisa}\\
\multicolumn{7}{l}{No intuito de evitar efeitos prováveis de}\\
\multicolumn{7}{l}{multicolinearidade entre as variáveis meteorológicas}\\
\multicolumn{7}{l}{que naturalmente tendem a ter resultados próximos, foi}\\
\multicolumn{7}{l}{feito uma nova estimação do modelo, onde foram}\\
\multicolumn{7}{l}{consideradas apenas as variáveis que foram}\\
\multicolumn{7}{l}{estatisticamente significantes para o objeto deste}\\
\multicolumn{7}{l}{estudo. Desta forma, foram retiradas as variáveis}\\
TEMP\textunderscore MAX,\tabcellsep \multicolumn{3}{l}{TEMP\textunderscore MEDIA,}\tabcellsep \multicolumn{3}{l}{DESCARGAS\textunderscore ATM,}\\
VENTOS\textunderscore KMH,\tabcellsep \tabcellsep e\tabcellsep atribuindo\tabcellsep o\tabcellsep nome\tabcellsep de\\
\multicolumn{2}{l}{modelo\textunderscore poisson2.}\tabcellsep \multicolumn{2}{l}{Variáveis}\tabcellsep \multicolumn{2}{l}{Estimativa}\tabcellsep Erro padrão\tabcellsep Valor de Z\tabcellsep Pr(>|z|)\\
\tabcellsep \multicolumn{3}{l}{(Intercepto)}\tabcellsep 0,7762\tabcellsep \tabcellsep 0,6353\tabcellsep 1,2220\tabcellsep 0,2218\\
\tabcellsep \multicolumn{3}{l}{FASESABC}\tabcellsep \multicolumn{2}{l}{-0,3258}\tabcellsep 0,4756\tabcellsep -0,6850\tabcellsep 0,4934\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{FASESAC}\tabcellsep \multicolumn{2}{l}{-2,0067}\tabcellsep 0,5547\tabcellsep -3,6170\tabcellsep 0,0003\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{FASESB}\tabcellsep \multicolumn{2}{l}{-4,9644}\tabcellsep 0,5544\tabcellsep -8,9540\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA10}\tabcellsep \multicolumn{2}{l}{-5,8039}\tabcellsep 0,4141\tabcellsep -14,0140\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA1000}\tabcellsep \multicolumn{2}{l}{-5,3779}\tabcellsep 0,4167\tabcellsep -12,9070\tabcellsep <2e-16\\
\tabcellsep \tabcellsep \multicolumn{2}{l}{KVA112,5}\tabcellsep 0,0339\tabcellsep \tabcellsep 0,0417\tabcellsep 0,8140\tabcellsep 0,4158\end{longtable} \par
 
\caption{\label{tab_0}Tabela 1 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{} \par 
\begin{longtable}{P{0.0429140832595217\textwidth}P{0.15132860938883969\textwidth}P{0.15057573073516387\textwidth}P{0.10766164747564215\textwidth}P{0.39751992914083256\textwidth}}
Grau\tabcellsep Designação\tabcellsep m/s\tabcellsep km/h\tabcellsep Efeitos em terra\\
0\tabcellsep Calmo\tabcellsep <0,3\tabcellsep <1\tabcellsep Fumaça sobe na vertical\\
1\tabcellsep Aragem\tabcellsep 0,3 a 1,5\tabcellsep 1 a 5\tabcellsep Fumaça indica direção do vento\\
2\tabcellsep Brisa leve\tabcellsep 1,6 a 3,3\tabcellsep 6 a 11\tabcellsep As folhas das árvores movem; os moinhos começam a trabalhar\\
3\tabcellsep Brisa fraca\tabcellsep 3,4 a 5,4\tabcellsep \multicolumn{2}{l}{12 a 19 As folhas agitam-se e as bandeiras desfraldam ao vento}\\
4\tabcellsep Brisa moderada\tabcellsep 5,5 a 7,9\tabcellsep 20 a 28\tabcellsep Poeira e pequenos papéis levantados; movem-se os galhos das árvores\\
5\tabcellsep Brisa forte\tabcellsep 8 a 10,7\tabcellsep \multicolumn{2}{l}{29 a 38 Movimentação de grandes galhos e árvores pequenas}\\
6\tabcellsep Vento fresco\tabcellsep \multicolumn{2}{l}{10,8 a 13,8 39 a 49}\tabcellsep Movem-se os ramos das árvores; dificuldade em manter um guarda-chuva aberto; assobio em fios de postes\\
7\tabcellsep Vento forte\tabcellsep \multicolumn{2}{l}{13,9 a 17,1 50 a 61}\tabcellsep Movem-se as árvores grandes; dificuldade em andar contra o vento\\
8\tabcellsep Ventania\tabcellsep \multicolumn{2}{l}{17,2 a 20,7 62 a 74}\tabcellsep Quebram-se galhos de árvores; dificuldade em andar contra o vento; barcos permanecem nos portos\\
9\tabcellsep Ventania forte\tabcellsep \multicolumn{2}{l}{20,8 a 24,4 75 a 88}\tabcellsep Danos em árvores e pequenas construções; impossível andar contra o vento\\
10\tabcellsep Tempestade\tabcellsep \multicolumn{3}{l}{24,5 a 28,4 89 a 102 Árvores arrancadas; danos estruturais em construções}\\
11\tabcellsep \multicolumn{4}{l}{Tempestade violenta 28,5 a 32,6 103 a 117 Estragos generalizados em construções}\\
\multicolumn{3}{l}{Fonte: Resultados originais da pesquisa}\tabcellsep \tabcellsep \end{longtable} \par
 
\caption{\label{tab_1}}\end{figure}
 			\label{foot_0}\footnote{\label{foot_0} © 2023 Global Journals} 			\label{foot_1}\footnote{\label{foot_1} Então, como pode ser visto na Tabela 3, também foi aplicado o modelo de Poisson, e já é possível verificar em termos quantitativos, que há menos variáveis de fenômenos meteorológicos se correspondendo com a variável dependente, interrupção no fornecimento de energia, do que na cidade de São Luís, sendo elas as variáveis PLUVIOMETRIA e RAJADA. A primeira se referindo ao volume de chuvas, e a segunda com a força repentina do vento, ambas com p-valor abaixo de 0.05, o que significa são estatisticamente significantes para explicar as interrupções no sistema elétrico, há um nível de significância de 5\%.Outra diferença em relação aos dados de São Luís, é que enquanto no estudo da capital, em que dois tipos de fases da amostra (AC, e B) se mostraram estatisticamente significante para fins preditivos, a um intervalo de confiança de 95\%, nenhuma fase da amostra de Imperatriz teve resultado semelhante. Entretanto, vale pontuar que a fase ABC possui uma significância de 10\%, enquanto na amostra de São Luís, a mesma não teve significância importante para o modelo estatístico de Poisson.Já em relação a potência dos transformadores (variável KVA), resultado muito semelhante se obteve ao estudo da amostra de dados de São Luís. Em Imperatriz, praticamente todos os transformadores de distribuição mostraram p-valor abaixo de 0.05, quanto a sua potência, ficando de fora apenas os transformadores com potência de 45kva, apresentando um p-valor igual a 0.16979.} 		 		\backmatter  			  				\begin{bibitemlist}{1}
\bibitem[Bussab et al. ()]{b1}\label{b1} 	 		\textit{},  		 			W O Bussab 		,  		 			P A Morettin 		,  		 			Estatística Básica 		.  		2005. São Paulo: Saraiva. 5.  	 
\bibitem[Lima and Ferreira ()]{b6}\label{b6} 	 		\textit{},  		 			R M B F Lima 		,  		 			A J A Ferreira 		.  	 	 		\textit{Estudos de Geografia do Maranhão. Maranhão. EDUFMA}  		2013.  	 
\bibitem[Silva ()]{b7}\label{b7} 	 		\textit{A relação entre raios e chuvas nos arredores de Manaus: análise de casos extremos. Dissertação (Mestrado) -INPA/UEA},  		 			J C Silva 		.  		2013. Manaus, AM.  	 
\bibitem[Baltar ()]{b0}\label{b0} 	 		\textit{Análise fatorial múltipla para tabelas de contingência},  		 			V T Baltar 		.  		2005.  		 			Tese de Doutorado. Universidade de São Paulo 		 	 
\bibitem[Wronscki ()]{b13}\label{b13} 	 		‘Apêndice: Operações utilizadas no R Studio para a realização deste estudo’.  		 			V R Wronscki 		.  	 	 		\textit{Seminário Brasileiro sobre Qualidade da Energia Elétrica},  		 (Aracaju, SE)  		2003.  	 
\bibitem[Canalenergia ()]{b2}\label{b2} 	 		\textit{Cemig é multada por indicadores de qualidade},  		 			Canalenergia 		.  		 \url{<https://www.canalenergia.com.br/noticias/53066819/cemig-e-multada-por-indicadores-de}  		2022.  	 
\bibitem[Weathers]{b8}\label{b8} 	 		\textit{Clima e condições meteorológicas médias em Imperatriz no ano todo},  		 			Park Weathers 		.  		 \url{<https://pt.weatherspark.com/y/302}  		 	 	 (Disponível em) 
\bibitem[Weathers]{b10}\label{b10} 	 		\textit{Clima e condições meteorológicas médias em São Luís no ano todo},  		 			Park Weathers 		.  		 \url{<https://pt.weatherspark.com/y/305}  		 	 	 (Disponível em) 
\bibitem[Clima-caracter%C3%ADstico-em-Imperatriz-Bra sil-durante-o-ano> Acesso em 15 de out ()]{b9}\label{b9} 	 		\textit{Clima-caracter\%C3\%ADstico-em-Imperatriz-Bra sil-durante-o-ano> Acesso em 15 de out},  		2022.  	 
\bibitem[Clima-caracter%C3%ADstico-em-S%C3%A3o-Lu %C3%ADs-Brasil-durante-o-ano> Acesso em 15 de out ()]{b11}\label{b11} 	 		\textit{Clima-caracter\%C3\%ADstico-em-S\%C3\%A3o-Lu \%C3\%ADs-Brasil-durante-o-ano> Acesso em 15 de out},  		2022.  	 
\bibitem[Grupo Equatorial Energia ()]{b4}\label{b4} 	 		 			Grupo Equatorial Energia 		.  		 \url{<https://ri.equatorialenergia.com.br/pt-br/divulgacao-e-resultados/central-de-resultados}  		\textit{Central de Resultados. Disponível em},  		2022.  	 
\bibitem[Guide to Instruments and Methods of Observation: Volume I -Measurement of Meteorological Variables Genebra, Suíça: Organização Meteorológica Mundial]{b12}\label{b12} 	 		‘Guide to Instruments and Methods of Observation: Volume I -Measurement of Meteorological Variables’.  		 \url{<https://library.wmo.int/doc\textunderscore num.php?explnum\textunderscore id=10616}  	 	 		\textit{Genebra, Suíça: Organização Meteorológica Mundial}  		World Meteorological Organization.  (8) .  	 	 (7 ed. 2018. 618p. Disponível em. Acesso em 04 de out. de 2022) 
\bibitem[Imprensa Nacional [IM]. [ANEEL] divulga desempenho e ranking das distribuidoras sobre fornecimento de energia em 2021 ()]{b5}\label{b5} 	 		\textit{Imprensa Nacional [IM]. [ANEEL] divulga desempenho e ranking das distribuidoras sobre fornecimento de energia em 2021},  		 \url{<https://www.gov.br/aneel/pt-br/assuntos/noticias/2022/aneel-divulga-desempenho-e-ranking-das-distribuidoras-sobre-fornecimento-de-energia-em-2021}  		2022.  	 
\bibitem[Fávero and Belfiore ()]{b3}\label{b3} 	 		\textit{Manual de Análise de Dados},  		 			L P Fávero 		,  		 			P Belfiore 		.  		2022. Rio de Janeiro: LTC.  	 
\end{bibitemlist}
 			 		 	 
\end{document}
