Rの出力からのすべての情報を含む係数を抽出する

Question

係数を、グラスの出力からのSE、t値、およびp値と一緒にRに取得する必要があります。

library(nlme) 
fm1 <- gls(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), Ovary, 
     correlation = corAR1(form = ~ 1 | Mare)) 

summary(fm1) 
Generalized least squares fit by REML 
    Model: follicles ~ sin(2 * pi * Time) + cos(2 * pi * Time) 
    Data: Ovary 
    AIC  BIC logLik 
    1571.455 1590.056 -780.7273 

Correlation Structure: AR(1) 
Formula: ~1 | Mare 
Parameter estimate(s): 
    Phi 
0.7532079 

Coefficients: 
        Value Std.Error t-value p-value 
(Intercept)  12.216398 0.6646437 18.380373 0.0000 
sin(2 * pi * Time) -2.774712 0.6450478 -4.301561 0.0000 
cos(2 * pi * Time) -0.899605 0.6975383 -1.289685 0.1981 

Correlation: 
       (Intr) s(*p*T 
sin(2 * pi * Time) 0.000  
cos(2 * pi * Time) -0.294 0.000 

Standardized residuals: 
     Min   Q1   Med   Q3   Max 
-2.41180365 -0.75405234 -0.02923628 0.63156880 3.16247697 

Residual standard error: 4.616172 
Degrees of freedom: 308 total; 305 residual 

誰かが私にこれを理解させるのを助けてくれたら大変感謝します。前もって感謝します。

Answer 1

はこれを試してみてください。

> cs <- as.data.frame(summary(fm1)$tTable) 
> cs 
         Value Std.Error t-value  p-value 
(Intercept)  12.2163982 0.6646437 18.380373 2.618737e-51 
sin(2 * pi * Time) -2.7747122 0.6450478 -4.301561 2.286284e-05 
cos(2 * pi * Time) -0.8996047 0.6975383 -1.289685 1.981371e-01 
> cs$t 
[1] 18.380373 -4.301561 -1.289685 
> cs$p 
[1] 2.618737e-51 2.286284e-05 1.981371e-01 

Answer 2

は、あなたがテーブルのいくつかの種類のために上記の値をしたいと仮定すると、ソリューションは比較的簡単です。より一般的に

sumfm1 <- summary(fm1) 

sumfm1$tTable 

         Value Std.Error t-value  p-value 
(Intercept)  12.2163982 0.6646437 18.380373 2.618737e-51 
sin(2 * pi * Time) -2.7747122 0.6450478 -4.301561 2.286284e-05 
cos(2 * pi * Time) -0.8996047 0.6975383 -1.289685 1.981371e-01 

あなたが任意のRオブジェクトのstr()関数を呼び出した場合、あなたは通常、（試行錯誤のビットで）あなたが必要とする結果を抽出する方法を見つけ出すことができます。

EDIT：スウェーブファイルに移動する必要がある場合は、上記のオブジェクトのxtableを呼び出すとすべてが有効になります。