2

私は次のコードを持っています：

test <- ca.jo(x, type='trace', ecdet='const', K=2)

私が要約（テスト）を書いているとき、次のことが起こります：

最初の列に正規化された固有ベクトル: (これらは共和分関係です)

              gld.l2     gdx.l2
    gld.l2  1.000000  1.0000000
    gdx.l2 -1.488325 -0.1993057

これらの正規化された固有ベクトルをどのように呼び出すことができますか?

私が書いているとき

slot(test, "Vorg")

次のデータのみを取得します

              gld.l2       gdx.l2
    gld.l2 -0.01346063 -0.012380092
    gdx.l2  0.02003378  0.002467422

しかし、私は正規化されたものを呼びたいです。

4

1 に答える 1

3 
data(denmark)
sjd <- denmark[, c("LRM", "LRY", "IBO", "IDE")]
sjd.vecm <- ca.jo(sjd, ecdet = "const", type="eigen", K=2, spec="longrun",
                  season=4)
sm <- summary(sjd.vecm)
sm@V
            LRM.l2     LRY.l2     IBO.l2     IDE.l2   constant
LRM.l2    1.000000  1.0000000  1.0000000   1.000000  1.0000000
LRY.l2   -1.032949 -1.3681031 -3.2266580  -1.883625 -0.6336946
IBO.l2    5.206919  0.2429825  0.5382847  24.399487  1.6965828
IDE.l2   -4.215879  6.8411103 -5.6473903 -14.298037 -1.8951589
constant -6.059932 -4.2708474  7.8963696  -2.263224 -8.0330127

詳細を確認したい場合str(sm)があります。

于 2012-07-13T22:15:12.473 に答える