Investment | R Documentation |

## US Investment Data

### Description

US data for fitting an investment equation.

### Usage

data(Investment)

### Format

An annual time series from 1963 to 1982 with 7 variables.

- GNP
nominal gross national product (in billion USD),

- Investment
nominal gross private domestic investment (in billion USD),

- Price
price index, implicit price deflator for GNP,

- Interest
interest rate, average yearly discount rate
charged by the New York Federal Reserve Bank,

- RealGNP
real GNP (= GNP/Price),

- RealInv
real investment (= Investment/Price),

- RealInt
approximation to the real interest rate
(= Interest - 100 * diff(Price)/Price).

### Source

Table 15.1 in Greene (1993)

### References

Greene W.H. (1993), *Econometric Analysis*, 2nd edition.
Macmillan Publishing Company, New York.

Executive Office of the President (1984), *Economic Report of the
President*. US Government Printing Office, Washington, DC.

### Examples

## Willam H. Greene, Econometric Analysis, 2nd Ed.
## Chapter 15
## load data set, p. 411, Table 15.1
data(Investment)
## fit linear model, p. 412, Table 15.2
fm <- lm(RealInv ~ RealGNP + RealInt, data = Investment)
summary(fm)
## visualize residuals, p. 412, Figure 15.1
plot(ts(residuals(fm), start = 1964),
type = "b", pch = 19, ylim = c(-35, 35), ylab = "Residuals")
sigma <- sqrt(sum(residuals(fm)^2)/fm$df.residual) ## maybe used df = 26 instead of 16 ??
abline(h = c(-2, 0, 2) * sigma, lty = 2)
if(require(lmtest)) {
## Newey-West covariances, Example 15.3
coeftest(fm, vcov = NeweyWest(fm, lag = 4))
## Note, that the following is equivalent:
coeftest(fm, vcov = kernHAC(fm, kernel = "Bartlett", bw = 5, prewhite = FALSE, adjust = FALSE))
## Durbin-Watson test, p. 424, Example 15.4
dwtest(fm)
## Breusch-Godfrey test, p. 427, Example 15.6
bgtest(fm, order = 4)
}
## visualize fitted series
plot(Investment[, "RealInv"], type = "b", pch = 19, ylab = "Real investment")
lines(ts(fitted(fm), start = 1964), col = 4)
## 3-d visualization of fitted model
if(require(scatterplot3d)) {
s3d <- scatterplot3d(Investment[,c(5,7,6)],
type = "b", angle = 65, scale.y = 1, pch = 16)
s3d$plane3d(fm, lty.box = "solid", col = 4)
}