Tietokoneharjoitus 1

Data sisaltaa 2 muuttujaa:

RealGDP: Quarterly values of Real GDP for the United States in Billions of Chained (2000) Dollars Seasonally Adjusted, Annual Rate.

TBillRate: Quarterly values of the rate on 3-month Treasury Bills. Quaterterly averages of daily rates in percentage points at an annual rate.

Datan lukeminen R:aan

file<-"http://cc.oulu.fi/~jklemela/econometrics/USMacro_Quarterly.csv"
data<-read.table(file,skip=1,sep=",")

Datan lukeminen SAS:iin

FILENAME myurl URL 'http://cc.oulu.fi/~jklemela/econometrics/USMacro_Quarterly.txt';

DATA USmacro;
   INFILE myurl firstobs=2;
   INPUT time $ gdpq tbill;
RUN;

Tehtävä 5

Laske muuttujan RealGDP keskiarvo ja keskihajonta ja piirrä aikasarjasta kuvio.

gdp<-data[,2]
plot(gdp)

mean(gdp)
sd(gdp)

#[1] 5041.707
#[1] 2657.933

# muita tulostuksia

start<-c(1947,1)
end<-c(2004,4)
gdp.ts<-ts(gdp,start=start,end=end,frequency=4)
plot(gdp.ts,cex.axis=2)

time<-matrix(0,length(gdp),1)
delta<-0.25
for (i in 1:length(gdp)) time[i]<-1947+i*delta
plot(time,gdp,type="l",cex.axis=1.5)

t<-data[,1]
plot(t,gdp)

Laske muuttujan TbillRate keskiarvo ja keskihajonta ja piirrä aikasarjasta kuvio.

tbill<-data[,3]
plot(tbill,type="l")

mean(tbill)
sd(tbill)

#[1] 4.749526
#[1] 2.949203

SAS:

PROC means data = USmacro;
run;

PROC sgplot data = usmacro;
series x=time y= gdpq;
series x=time y= tbill / y2axis;
xaxis display = none;
run;

Ohjeita muuttujien muuntamiseen:

gdp<-data[,2]
tbill<-data[,3]
n<-length(tbill)

loggdp<-log(gdp)
plot(loggdp,type="l")

y<-log(gdp[2:n])-log(gdp[1:(n-1)])
plot(y,type="l")

x<-tbill[2:n]-tbill[1:(n-1)]
plot(x,type="l")