Como usar este documento

O código abaxio nas caixas cinzas pode ser copiado diretamente ao “Rstudio”, na janela para Scripts.

I. Introdução: Ler e visualizar dados

Neste módulo, vamos aprender:
1. Como criar um script em R e rodá-lo.
2. Como instalar pacotes.
3. Como ler dados de vazão.
4. Como visualizar dados de vazão.

1. Criar pastas e baixar dados, scripts

No Windows: 1.1. Criar uma pasta para os dados: ex.: “C:/minha_pasta/Rcurso2021/dados/”
1.2. Criar uma pasta para os códigos em R: C:/minha_pasta/Rcurso2021/Rscripts/"
1.3. Baixar os dados de vazão aqui e colocar na pasta para dados.
* Se nao quer abrir, “right-click” o link e “abrir numa nova janela”. 1.4. Abre o arquivo de vazão com o bloco de notas (Notepad). R não pode ler diretamente.
1.5. Baixar o script “readVazao_corr.R” aqui, copia numa Rscript, salva como “readVazao_corr.R”, e coloca na pasta para scripts.

2. Estabelecer pastas em R

2.1. Abrir o Rstudio. Familiarizar-se com as quatro janelas.
2.2. Criar um novo script: Arquivo –> Novo –> Rscript
2.3. Fazer commentário no script:

#  Meu primeiro script!

2.4. Salvar o script na pasta para scripts (“C:/minha_pasta/Rcurso2021/Rscripts/”). Pode chamar o script “exemplo_vazao.R”.
2.5. Estabelecer a pasta onde R vai encontrar os dados, e executar a linha:

pasta.dados = "C:/minha_pasta/Rcurso2021/dados/"

Ver a variável “pasta.dados” (escrever o nome da varável e apertar ‘enter’… pode fazer isso ao longo de todo exercício, quando quiser verificar uma variável)

pasta.dados
## [1] "C:/minha_pasta/Rcurso2021/dados/"
class(pasta.dados)
## [1] "character"

Ver os arquivos que estão na pasta.dados:

dir(pasta.dados)
## [1] "Q.Jaru.csv"            "vazoes_C_15560000.csv"

Tem “vazoes_C_15560000.csv”? Se não, verifique se salvou o arquivo .csv no lugar certo.

2.6. Estabelecer a pasta onde o R vai encontrar os scripts:

pasta.R = "C:/minha_pasta/Rcurso2021/Rscripts/"
pasta.R
## [1] "C:/minha_pasta/Rcurso2021/Rscripts/"
class(pasta.R)
## [1] "character"

Ver os arquivos que estão nesta pasta:

dir(pasta.R)
## [1] "exemplo_vazao.R"         "exemplo_vazao_ao_vivo.R"
## [3] "readVazao.R"             "readVazao_b.R"          
## [5] "readVazao_corr.R"

Tem “readVazao.R” e “exemplo_vazao.R”? Se não, verifique se salvou os arquivos no lugar certo.

Note que este código faz a mesma coisa:

dir("C:/minha_pasta/Rcurso2021/Rscripts/")
## [1] "exemplo_vazao.R"         "exemplo_vazao_ao_vivo.R"
## [3] "readVazao.R"             "readVazao_b.R"          
## [5] "readVazao_corr.R"

3. Instalar pacotes e carregar

3.1. Instalar. Só precisa fazar isso uma vez no seu computador.

install.packages("Hmisc")  # Para manipular objetos de clase caracter
install.packages("stringr") # Tambem para manipular objectos de clase caracter

3.2. Carregar pacotes. Tem que fazer isso todas as vezes que iniciar uma nova sessão do R:

library(Hmisc)
## Warning: package 'Hmisc' was built under R version 4.1.1
## Loading required package: lattice
## Loading required package: survival
## Loading required package: Formula
## Warning: package 'Formula' was built under R version 4.1.1
## Loading required package: ggplot2
## Warning: package 'ggplot2' was built under R version 4.1.1
## 
## Attaching package: 'Hmisc'
## The following objects are masked from 'package:base':
## 
##     format.pval, units
library(stringr)
## Warning: package 'stringr' was built under R version 4.1.1

3.3. Executar a função “readVazao.R”:

setwd(pasta.R)
source("readVazao_corr.R")
class(readVazao)
## [1] "function"

4. Ler e formatar o arquivo de vazão:

setwd(pasta.dados)
x = readVazao(pasta.dados)
class(x)
## [1] "data.frame"
head(x)
##         Date Vazao.cms Status
## 1 1977-12-01        NA     NA
## 2 1977-12-02        NA     NA
## 3 1977-12-03        NA     NA
## 4 1977-12-04        NA     NA
## 5 1977-12-05        NA     NA
## 6 1977-12-06        NA     NA

Mudar os nomes das colunas:

names(x) = c("Date","Q.cms","Status")

4.1. Ver o classe dos dados, e formatar o data usando “as.Date”:

class(x$Date)
## [1] "Date"
x$Date = as.Date(x$Date)

Confirmar que funcionou:

class(x$Date)
## [1] "Date"
head(x)
##         Date Q.cms Status
## 1 1977-12-01    NA     NA
## 2 1977-12-02    NA     NA
## 3 1977-12-03    NA     NA
## 4 1977-12-04    NA     NA
## 5 1977-12-05    NA     NA
## 6 1977-12-06    NA     NA

4.2. Plotar os dados:

plot(x$Date,x$Q.cms,type="l",ylab="Vazao cms",xlab="")

5. Salvar o arquivo como csv.

setwd(pasta.dados)
write.csv(x,"Q.Jaru.csv",row.names=FALSE)

Abrir o arquivo em Notepad para verificar que funcionou.