0

R でリッジ回帰パッケージを使用していますが、係数リストの T 統計、p 値を抽出するのに問題があります。

hereのサンプルコードを次に示します。

data(GenCont)
mod <- linearRidge(Phenotypes ~ ., data = as.data.frame(GenCont))
summary(mod)

これは戻ります

Call:
linearRidge(formula = Phenotypes ~ ., data = as.data.frame(GenCont))

Coefficients:
             Estimate Scaled estimate Std. Error (scaled) t value (scaled) Pr(>|t|)    
(Intercept)  1.533386              NA                  NA               NA       NA    
SNP1         0.277296        4.045409            0.266120           15.201  < 2e-16 ***
SNP2        -0.110458       -1.256154            0.216332            5.807 6.38e-09 ***

しかし、プログラムで「係数」のデータにアクセスしたいのです。私がlmを使用していた場合、次のようなことを理解しています

coef(summary(mod))

動作するはずですが、これは私が求めていることをしていないようです。私は単にNULLを取得します。

このデータにアクセスする方法について何かアイデアはありますか?

4

1 に答える 1

3

リッジ パラメーターが自動的に選択される状況で係数テーブルを取得するには、summary()結果を何らかのオブジェクトとして保存してから、リスト要素にアクセスする必要がありますsummaries$summary1$coefficients。summary()function で見ることができるオブジェクトの全体構造str(sumar)。

sumar<-summary(mod)
sumar$summaries$summary1$coefficients
                Estimate Scaled estimate Std. Error (scaled) t value (scaled)     Pr(>|t|)
(Intercept)  1.533385893              NA                  NA               NA           NA
SNP1         0.277296215     4.045408706           0.2661197      15.20146202 0.000000e+00
SNP2        -0.110457822    -1.256153622           0.2163319       5.80660491 6.375233e-09
SNP3        -0.110457822    -1.256153622           0.2163319       5.80660491 6.375233e-09
SNP4         0.005229639     0.011635212           0.3716925       0.03130332 9.750276e-01
SNP5         0.531172545     6.323006229           0.3153685      20.04958196 0.000000e+00
SNP6        -0.119163534    -1.373227248           0.2230470       6.15667175 7.428960e-10
SNP7         0.113843942     0.113730041           0.3721807       0.30557749 7.599264e-01
SNP8        -0.099148877    -1.028580596           0.3558067       2.89084074 3.842128e-03
SNP9        -0.008320553    -0.008312229           0.3723863       0.02232152 9.821915e-01
SNP10        0.058562323     0.101128163           0.3715670       0.27216670 7.854938e-01
SNP11       -0.096526424    -1.495698673           0.3292496       4.54275034 5.552500e-06
SNP12       -0.334279101    -0.333944654           0.3722483       0.89710186 3.696646e-01

t 値のみを取得するには、このテーブルの 4 番目の列を選択します。

> sumar$summaries$summary1$coefficients[,4]
(Intercept)        SNP1        SNP2        SNP3        SNP4        SNP5        SNP6        SNP7        SNP8        SNP9 
         NA 15.20146202  5.80660491  5.80660491  0.03130332 20.04958196  6.15667175  0.30557749  2.89084074  0.02232152 
      SNP10       SNP11       SNP12 
 0.27216670  4.54275034  0.89710186

summary() 関数によって生成された他の要素にアクセスするには、適切なリスト要素を選択するだけです。サマリー オブジェクトの構造は、選択する要素を示します。

str(sumar$summaries$summary1)
List of 4
 $ coefficients: num [1:13, 1:5] 1.53339 0.2773 -0.11046 -0.11046 0.00523 ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:13] "(Intercept)" "SNP1" "SNP2" "SNP3" ...
  .. ..$ : chr [1:5] "Estimate" "Scaled estimate" "Std. Error (scaled)" "t value (scaled)" ...
 $ df          : Named num [1:3] 3.12 1.21 5.04
  ..- attr(*, "names")= chr [1:3] "model" "variance" "residual"
 $ nPCs        : int 1
 $ lambda      : num 2.21

たとえば、取得するにはdf

sumar$summaries$summary1$df
   model variance residual 
3.120934 1.205412 5.036457

結果は、新しいオブジェクトを作成せずに選択することもできます。summary(mod)オブジェクト名の代わりに使用するだけsumarです。

summary(mod)$summaries$summary1$df
   model variance residual 
3.120934 1.205412 5.036457
于 2013-08-19T10:14:05.403 に答える