= read.csv2(
dados "IMB201009.txt",
header=FALSE,
sep="")
IBGE1
Farei aqui reproduções das análises do Guto usando o notebook do R.
Pois é, o R também tem notebook que podem facilmente criar arquivos em PDF, HTML e DOC.
Adeus editor de texto. Se o seu texto é um relatório que envolve análise de dados, é muito melhor fazer por aqui.
Aqui vai a reprodução da análise IBGE 1.
Primeiramente abrimos os dados baixados do site do IBGE, conforme explicado pelo Guto. Pode ser qualquer arquivo. O meu foi IMB201009.txt. Vamos guardar o arquivo na variável dados
, que é do tipo dataframe
.
Adapte o caminho completo do arquivo para o seu computador. A função read.csv2
lê dados separados por vírgulas.
Para ver os dados use View(dados)
ou print(dados)
.
Os nomes das colunas (V1 a V4) foram dados automaticamente.
Agora vamos executar os códigos para gerar um gráfico:
=length(dados$V4) #Número de linhas de V4
nplot(1:n, #Plota x e y
$V4) dados
Vamos plotar do novo mudando algumas características do gráfico.
plot(1:n, #Plota novamente em linha
$V4,
dadostype='line',
xlab='',ylab='')
Agora vamos ajeitar o eixo x (tempo), que ainda não é do tipo datetime. Primeiro concatenamos as duas colunas de Data e Hora e depois transformamos para o tipo datetime.
=paste(dados$V1,dados$V2)
tempo=as.POSIXct(tempo,format = "%d/%m/%Y %H:%M:%S") tempo
Cria o plot completo contendo a série V3 em azul e V4 em vermelho.
plot(tempo, #Eixo x
$V3, #Eixo y
dadostype='line', #Tipo de linha
col='blue', #Cor da linha
main='IMB201009.txt', #Título do gráfico
ylab='Nível (m)', #Título do eixo y
xlab='Tempo') #Título do eixo x
lines(tempo, #Adiciona nova série (V4) em vermelho
$V4,
dadostype='line',
col='red')
legend('topright', #Adiciona legenda
legend=c('Var1','Var2'),
lty=c(1,1),
col=c('blue','red'),
cex=0.6)