Problema proposto:

Este é um estudo sobre pessoas que cometeram suicídio. Dentre estas pessoas, existe o interesse de saber o que leva indivíduos a cometerem suicídio com arma de fogo.

As variáveis do os estudo eram:

aux.idade = idades aux.id_legal = 1, tinha id. para porte de arma , o c.c. trab.armado = 1, trabalha com porte de arma , 0 c.c. sexo = 1=masculino, 2=feminino aux.raca = 1=branco, 2=preto, 3=amarelo, 4=pardo, 5=indigena, 6=NA aux.estcivil = 1= solteiro, 2=uniao estavel, 3= casado, 4= viuvo, 5=divorciado, 6=NA aux.esc = 1 = 1 a 3 anos, 2= 4 a 7 anos, 3= 8 a 11 anos, 5= 12+, 6=NA

Objetivos:

-O que leva um individuo a cometer suicidio com arma de fogo?

-usar estado (arrumar microregioes) para explicar (se houver problema pode usar grandes regioes)

-proporcao de PAF por UF

-Faca o grafico, mapas

-pode dividir por faixa etaria

Foi utilizado:

RStudio , pacotes ggplot2, mapas, modelos generalizados, regressão logística

Neste relatório serão apresentados os resultados do estudo sobre o perfil dos indivíduos que cometem suicídio por arma de fogo nos estados brasileiros. Para esse estudo foram coletados uma amostra de 1500 indivíduos que cometeram ou não suicídio PAF em 20 estados brasileiros.

Iremos apresentar uma visão geral, através de estatística descritiva, de como estão distribuidos os suicídios nos estados brasileiros por: Genero, raça, escolariedade, estado civil, trabalha armado e idade.

Identificaremos também, através de um modelo de regressão linear, quais fatores levam um indivíduo a cometer suicídio por arma de fogo.

Analise exploratoria dos dados

Perfil da Amostra:

##  suicidio_paf   aux.idade     aux.id_legal trab.armado3        SEXO     
##  Nao:1293     Min.   : 1.00   Nao: 211     Nao:1288     Masculino:1189  
##  Sim: 207     1st Qu.:16.00   Sim:1289     Sim: 212     Feminino : 311  
##               Median :27.00                                             
##               Mean   :29.71                                             
##               3rd Qu.:41.00                                             
##               Max.   :81.00                                             
##      aux.raca              aux.estciv               aux.esc    regiao  
##  branco  :835   Menos que 3 anos:1276   Menos que 3 anos:637   CO:362  
##  preto   : 70   Mais que 3 anos : 224   Mais que 3 anos :863   NE:177  
##  amarelo :  6                                                  NO: 76  
##  pardo   :490                                                  SE:487  
##  indigena: 15                                                  SU:398  
##  Outros  : 84

Quantidade e proporção de suicídios PAF:

##           Nao     Sim
## [1,] 1293.000 207.000
## [2,]    0.862   0.138

Quantidade e Idade dos Indivíduos de acordo com o Fato de se Cometer ou Não Cometer Suicídio

  • Percebemos que 207.000 indivíduos, 13,8% da nossa amostra, cometeram suicídio PAF.

  • De acordo com que a idade do indivíduo aumenta, a quantidade de suicídios diminui.

Genero e Raça de acordo com o Fato de se Cometer ou Não Suicídio

  • A amostra é formada por 20% de indivíduos do sexo feminino e 80% so sexo masculino. Dentre os indivíduos do sexo masculino 15,5% cometeram suicídio, e dentre o feminino 7% cometeram suicídio.

  • Os percentuais de suicídios entre as raças ‘branca’, ‘preta’, ‘amarela’, ‘parda’, ‘indigêna’ e outras foram de 15,6%, 11,4%, 16,7%, 10,2%, 0%, 0,2%, respectivamente.

Estado Civil e Região de acordo com o Fato de se Cometer ou Não Suicídio

  • 85% da nossa amostra é classificada como ‘Menos que 3 anos’ no Estado Civil, e a taxa de suicídio nesse quesito é de 13,9%. Já para os 15% que possuem relacionamentos com mais de 3 anos a taxa de suicídio é de 12,9%.
  • Os percentuais de suicídios para as regiões CO, NE, NO, SE, SU foram de 12,4%, 10,1%, 11,8%, 12,3% e 18,8%, respectivamente.

Regioes

Interessante notar que existe grande variedade quando se leva em consideração as regiões, vejamos:

## dados$regiao: CO
## [1] 407
## -------------------------------------------------------- 
## dados$regiao: NE
## [1] 195
## -------------------------------------------------------- 
## dados$regiao: NO
## [1] 85
## -------------------------------------------------------- 
## dados$regiao: SE
## [1] 547
## -------------------------------------------------------- 
## dados$regiao: SU
## [1] 473

Baixas correlações, correlação alta apenas para aux.id.legal


Modelo para descrever PAF:

Os parâmetros do modelo são estimados pelo método de máxima verossimilhança (MV)

## 
## Call:  glm(formula = suicidio_paf ~ aux.idade + trab.armado3 + SEXO, 
##     family = binomial("logit"), data = dados)
## 
## Coefficients:
##  (Intercept)     aux.idade  trab.armado3          SEXO  
##     -2.46093      -0.01222       0.53156       1.07127  
## 
## Degrees of Freedom: 1499 Total (i.e. Null);  1496 Residual
## Null Deviance:       1204 
## Residual Deviance: 1175  AIC: 1183

Intervalos de 95% de confiança

Os intervalos de 95% de confiança para os parâmetros do modelo, com base na estatística de Wald:

##                    2.5 %       97.5 %
## (Intercept)  -2.97463618 -1.947228902
## aux.idade    -0.02133009 -0.003104303
## trab.armado3  0.10841820  0.954696180
## SEXO          0.59056864  1.551972265

Portando, o modelo que poderá ser utilizado para descrever os dados será o seguinte:

\[ ln \Big( \dfrac{(p_i) }{(1-(p_i))} \Big) = \beta_0 + \beta_1 X_{i1} +\beta_2 X_{21} +\beta_3 X_{31} \]

  • As variáveis:
    • \(X_{i1}\): Idade do indivíduo i
    • \(X_{i2}\): 1=se trabalha armado , 0 se não trabalha armado
    • \(X_{i3}\): 1=se sexo masculino , 0 sexo feminino
  • Os parâmetros :
    • \(p_i\): a probabilidade de cometer suicídio por arma de fogo
    • \(\beta_0\): intercepto do modelo (sem interpretação prática)
    • \(\beta_1\): o quanto varia a chance de suicídio por arma de fogo (em escala logarítimica) a aumentar a idade em 1 ano
    • \(\beta_2\): o efeito do nível " trabalha armado" em comparação a “não trabalhar armado”
    • \(\beta_3\): o efeito do nível " sexo masculino" em comparação ao “sexo feminino”

Medidas de associação (razões de chance)

Pode-se demonstrar matematicamente que a razão de chance é o exponencial da estimativa pontual


##  (Intercept)    aux.idade trab.armado3         SEXO 
##   0.08535532   0.98785713   1.70157993   2.91908571
## [1] 5.693878

Os intervalos de confiança para as razões de chance (odds ratio - OR), fixando o nível de confiança de 95%:

##                   2.5 %    97.5 %
## (Intercept)  0.05106601 0.1426689
## aux.idade    0.97889579 0.9969005
## trab.armado3 1.11451374 2.5978812
## SEXO         1.80501453 4.7207716

Razões de chance e os seus respectivos intervalos de confiança (com 3 casas decimais).

##                OR1 2.5 % 97.5 %
## (Intercept)  0.085 0.051  0.143
## aux.idade    0.988 0.979  0.997
## trab.armado3 1.702 1.115  2.598
## SEXO         2.919 1.805  4.721

Interpretação das razões de chance (odds ratio):

Tanto a idade quanto, a idade legal para porte de armas, quanto o sexo estão significativamente relacionadas com a chance de suicídio com arma de fogo entre as pessoas que cometeram suicídio (OBS: Note que o p-valor foi menor que o nível de significância de 5% e o IC para OR não inclui a unidade).

A chance do indivídio cometer suicídio com arma de fogo diminui em 98,8% ao aumentar em 1 ano de idade.

Indivíduos com idade legal para porte de armas tem chance de cometer suicídio com arma de fogo 70.2% maior do que os indivíduos que não tem idade legal para porte de armas.

Indivíduos do sexo masculino tem chance de cometer suicídio com arma de fogo quase 2 vezes mais do que as do sexo feminino.


Mapas:

Para a confecção dos mapas, foram consideradas apenas as variaveis dicotomicas: PAF, identidade com porte de armas, trabalho armado, sexo, estado civil e escolaridade.

Portanto, foi feita a contagem dessas variáveis levando em conta apenas quem cometeu PAF, avaliadas por regioes e por estados

Número de indivíduos que cometeram PAF:

Número de indivíduos que cometeram PAF e trabalhavam armado:


Número de indivíduos que cometeram PAF e são do sexo masculino:

Número de indivíduos que cometeram PAF e possuem mais de 3 anos de estudo: