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.
file<-"http://cc.oulu.fi/~jklemela/econometrics/USMacro_Quarterly.csv" data<-read.table(file,skip=1,sep=",")
FILENAME myurl URL 'http://cc.oulu.fi/~jklemela/econometrics/USMacro_Quarterly.txt'; DATA USmacro; INFILE myurl firstobs=2; INPUT time $ gdpq tbill; RUN;
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")