Note that there are some explanatory texts on larger screens.

plurals
  1. POextracting values from an aovp object in Lmperm
    primarykey
    data
    text
    <p>I have the following object <code>M</code>, from which I need to extract the <code>fstatistic</code>. It is a model generated by the function <code>summaryC</code> of a model generated by <code>aovp</code>, both functions from package <code>lmPerm</code>. I have tried hints for extracting values from normal linear models and from the functions in <code>attr</code>, <code>extract</code> and <code>getElement</code>, but without success. </p> <p>Anybody could give me a hint?</p> <pre><code>&gt; str(M) List of 2 $ Error: vegetation: NULL $ Error: Within :List of 11 ..$ NA : NULL ..$ terms :Classes 'terms', 'formula' length 3 Temp ~ depth .. .. ..- attr(*, "variables")= language list(Temp, depth) .. .. ..- attr(*, "factors")= int [1:2, 1] 0 1 .. .. .. ..- attr(*, "dimnames")=List of 2 .. .. .. .. ..$ : chr [1:2] "Temp" "depth" .. .. .. .. ..$ : chr "depth" .. .. ..- attr(*, "term.labels")= chr "depth" .. .. ..- attr(*, "order")= int 1 .. .. ..- attr(*, "intercept")= int 1 .. .. ..- attr(*, "response")= int 1 .. .. ..- attr(*, ".Environment")=&lt;environment: R_GlobalEnv&gt; ..$ residuals : Named num [1:498] -46.9 -43.9 -46.9 -38.9 -41.9 ... .. ..- attr(*, "names")= chr [1:498] "3" "4" "5" "6" ... ..$ coefficients : num [1:4, 1:4] -2.00 -1.00 -1.35e-14 1.00 2.59 ... .. ..- attr(*, "dimnames")=List of 2 .. .. ..$ : chr [1:4] "depth1" "depth2" "depth3" "depth4" .. .. ..$ : chr [1:4] "Estimate" "Std. Error" "t value" "Pr(&gt;|t|)" ..$ aliased : Named logi [1:4] FALSE FALSE FALSE FALSE .. ..- attr(*, "names")= chr [1:4] "depth1" "depth2" "depth3" "depth4" ..$ sigma : num 29 ..$ df : int [1:3] 4 494 4 ..$ r.squared : num 0.00239 ..$ adj.r.squared: num -0.00367 ..$ **fstatistic** : Named num [1:3] 0.395 3 494 .. ..- attr(*, "names")= chr [1:3] "value" "numdf" "dendf" ..$ cov.unscaled : num [1:4, 1:4] 0.008 -0.002 -0.002 -0.002 -0.002 ... .. ..- attr(*, "dimnames")=List of 2 .. .. ..$ : chr [1:4] "depth1" "depth2" "depth3" "depth4" .. .. ..$ : chr [1:4] "depth1" "depth2" "depth3" "depth4" ..- attr(*, "class")= chr "summary.lmp" - attr(*, "class")= chr "listof" </code></pre> <p>there it goes a reproducible example to play with:</p> <pre><code>Temp=1:100 depth&lt;- rep( c("1","2","3","4","5"), 100) vegetation=rep( c("1","2"), 50) df=data.frame(Temp,depth,vegetation) M=summaryC(aovp(Temp~depth+Error(vegetation),df, perm="")) </code></pre>
    singulars
    1. This table or related slice is empty.
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
 

Querying!

 
Guidance

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload