計量経済学のためのR環境

> x1 <- c( 1.5, 2.3, 3.8,   4.2,  5.6,   6.3)
> x2 <- c( 1.0, 3.9, 8.6, 15.25, 28.6, 32.68)
> y  <- c(-3.2, 2.0, 2.2,   6.3,  5.8,  12.5)
> eq1 <- lm(y ~ x1 + x2)
> summary (eq1)

Call:
lm(formula = y ~ x1 + x2)

Residuals:
     1      2      3      4      5      6 
-1.125  2.053 -1.517  1.421 -2.858  2.027 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)
(Intercept) -5.73749    5.62573  -1.020    0.383
x1           2.42235    2.82058   0.859    0.454
x2           0.02905    0.39772   0.073    0.946

Residual standard error: 2.713 on 3 degrees of freedom
Multiple R-Squared: 0.8416,     Adjusted R-squared: 0.7361 
F-statistic: 7.972 on 2 and 3 DF,  p-value: 0.06302 

> eq2 <- dwtest(eq1)
> summary (eq2)
            Length Class  Mode     
statistic   1      -none- numeric  
method      1      -none- character
alternative 1      -none- character
p.value     1      -none- numeric  
data.name   1      -none- character
> eq3 <- arma (x1, c(1,0))
> summary (eq3)


Call:
arma(x = x1, order = c(1, 0))

Model:
ARMA(1,0)

Residuals:
    Min      1Q  Median      3Q     Max 
-0.5374 -0.2993 -0.1108  0.4570  0.4907 

Coefficient(s):
           Estimate  Std. Error  t value Pr(>|t|)    
ar1          0.9296      0.1159    8.022 1.11e-15 ***
intercept    1.2049      0.4366    2.759  0.00579 ** 
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Fit:
sigma^2 estimated as 0.2101,  Conditional Sum-of-Squares = 0.84,  AIC = 11.67

>

トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2008-12-11 (木) 22:19:25 (5616d)