#Exercício 1.2
#a)

#Buscar a Localização do Documento:
setwd("C:/Users/HP/Desktop/Aula1eAula2") 

#Importar Dados do Ficheiro em questão:
cotinina = read.csv("C:/Users/HP/Desktop/Aula1eAula2/ex1.2.csv", sep = ';', header = T)

# ou alternativamente: cotinina = read.csv("ex1.2.csv", sep = ';', header = T)

#Retira-se a linha correspondente ao total, pois não é necessária:
cotinina = cotinina[-nrow(cotinina),]

cotinina

  #Classes:
    classes = cotinina$classesCotinina
    classes
    
  #Valores dos Fumadores:
    valFum = cotinina$Fumadores
    valFum
    
  #Valores dos Não Fumadores:
    valNFum = cotinina$naoFumadores
    valNFum

  frFum = valFum/1539
  frNFum = valNFum/3445

data.frame('Frequencia absoluta Fumadores' = valFum,
           'Frequencia relativa Fumadores' = frFum,
           'Frequencia absoluta Não Fumadores' = valNFum,
           'Frequencia relativa Não Fumadores' = frNFum)

  dev.new() #Abre uma janela de visualização. 
  
  #Gráfico das Distribuições Pedidas:
  par(mar=c(5, 4, 4, 8), xpd=TRUE)
  matplot(matrix(c(cumsum(frFum),cumsum(frNFum)),length(frFum),2), 
          type="l", ylab="Distribuições das Frequências de Cotinina",xlab="Classes", xaxt="n")
  legend("bottomright", c("Fumadores", "Não Fumadores"),lty=1, col=1:2)
  axis(side=1, at = c(1:length(frFum)), labels = c('0-13', '14-49', '50-99', '100-149', '150-199', '200-249', '250-299', '300-399'))
  

#b)
  
library('ggplot2')
new_cotinina = data.frame(new_classes = c(classes, classes),
                          all = c(valFum, valNFum),
                          distincao = c( rep('Fumador', 8) , rep('Não Fumador',8)))

#Histograma para Comparação das duas Populações:
ggplot( new_cotinina, aes(x=new_classes, y=all, fill=distincao )) +
  geom_bar( position='dodge', stat = 'identity') +
  scale_fill_manual(values = c('royalblue', 'darkorange'))+
  labs(title='Histograma para Comparação das Duas Populações', x= 'Classes', y='Frequência', fill='Distinção')+
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),
        panel.background = element_blank(), axis.line = element_line(colour = "black")) #Estético


#c)

#O número de pessoas fumadoras para um mesmo nível de cotinina é maior do que o número de pessoas 
#não fumadoras, pelo que se pode concluir que o tabaco leva a um aumento de produção desta substância.