2012-04-11 16 views
0

私はPython respを使い始めました。 Rpy2なので、私は初心者です。 これまでRpy2 でロードした独自のRパッケージを取得できました。入力変数の結果が生成されました(Rのlm()関数に匹敵する )。 私の結果は以下の通りです:rpy2での配列の索引付け

>>> res 
<ListVector - Python:0x1561828/R:0xe6a5c8> 
[Matrix, Array] 
    coef.f: <class 'rpy2.robjects.vectors.Matrix'> 
    [-10.554713, 1.640681, 0.957045, ..., 0.000000, 0.000000, 0.783108] 
    pred.f: <class 'rpy2.robjects.vectors.Array'> 
    [17.189366, 7.116058, 41.526062, ..., 459.952946, 312.088634, 677.908581] 

私もインデックス1使用して配列のサブパート(pred.f)を取得することができています:

>>>res1 = res[1] 
>>> res1 
<Array - Python:0x1561850/R:0x4b299d0> 
[17.189366, 7.116058, 41.526062, ..., 459.952946, 312.088634, 677.908581] 
>>> print res1 
, , L=300, AR=1.5, SO=1, T=30 

    sigma_stat sigma_mob 
fit 17.189366 152.21461 
lwr 7.116058 76.98245 
upr 41.526062 301.01379 

, , L=300, AR=1.5, SO=2, T=30 

    sigma_stat sigma_mob 
fit 27.62310 284.5423 
lwr 14.54639 173.3316 
upr 52.45825 467.1449 

, , L=300, AR=1.5, SO=3, T=30 

    sigma_stat sigma_mob 
fit 39.75894 459.9529 
lwr 24.07098 312.0886 
upr 65.67369 677.9086 

を今、私はインテリジェントな方法を使用したいですこの配列内の特定の番号にアクセスします。 例私はsigma_stat、フィット、L = 300、AR = 1.5、SO = 3、T = 30 ...の番号を希望します。これは私の配列が6次元の を意味します。

これを行う簡単な方法はありますか? res1( 'sigma_stat'、 'fit'、 'L' = 300、 'AR' = 1.5、 'SO' = 3、 'T' = 30)

ありがとうございますあなたの助けに...

答えて

0

この作品のようなものをお探しですか?

res1.rx('sigma_stat', 'fit', 'L=300', 'AR=1.5', 'SO=3', 'T=30') 

(ディメンション名が正しい順序でなければならないことに注意してください - 私は時間がテストすることを、引数リストの最後に行くべき「sigma_stat」と「フィット」をかもしれない)

関連する問題