In [1]:
options(OutDec = ",")
#==============================================================================
# Descricao da serie de vendas e precos de um certo tipo de doce
#==============================================================================
# Precos
# Vendas
# EUA, mensal,
# Janeiro de 1976 a dezembro de 1981
#
#==============================================================================
# Lendo os dados do arquivo doce.txt
#==============================================================================
dados <- read.table("../dados/doce.txt",header=TRUE,sep=" ")
print(head(dados))
print(tail(dados))
Preco Vendas
1 0,5366 9,535
2 0,5366 9,245
3 1,1220 9,355
4 1,3210 9,805
5 1,1220 8,821
6 1,1220 8,321
Preco Vendas
67 -0,9569 10,140
68 -1,6690 9,806
69 -1,6690 10,270
70 -0,9569 8,804
71 -1,1370 8,620
72 -0,9569 9,459
In [2]:
#==============================================================================
# Definindo como objeto series temporais do R
#==============================================================================
vendas <- ts(dados[,2],frequency=12,start=c(1976,1))
precos <- ts(dados[,1],frequency=12,start=c(1976,1))
In [3]:
#==============================================================================
# Grafico da serie de vendas
#==============================================================================
par(mfrow=c(1,1),lwd=2.0,cex.lab=1.5,cex.axis=1.5,lab=c(6,6,5),
mar=c(5,5,2,2.5),xpd=T,cex.main=2.0,bty="n")
ts.plot(vendas,xlab="Mês",ylab="Vendas",xlim=c(1976,1982),ylim=c(6,12))
points(vendas,pch=15)
In [4]:
#==============================================================================
# Grafico da serie de precos
#==============================================================================
par(mfrow=c(1,1),lwd=2.0,cex.lab=1.5,cex.axis=1.5,lab=c(6,8,5),
mar=c(5,5,2,2.5),xpd=T,cex.main=2.0,bty="n")
ts.plot(precos,xlab="Mês",ylab="Preços",xlim=c(1976,1982),ylim=c(-3,3))
points(precos,pch=15)
