Databases
StackOverflow2013
Postsdbo
POUsing lapply on a list of models
 

Note that there are some explanatory texts on larger screens.

plurals
  1. POUsing lapply on a list of models
    primarykey
    data
    text
    <p>I have generated a list of models, and would like to create a summary table.</p> <p>As and example, here are two models:</p> <pre><code>x &lt;- seq(1:10) y &lt;- sin(x)^2 model1 &lt;- lm(y ~ x) model2 &lt;- lm(y ~ x + I(x^2) + I(x^3)) </code></pre> <p>and two formulas, the first generating the equation from components of formula</p> <pre><code>get.model.equation &lt;- function(x) { x &lt;- as.character((x$call)$formula) x &lt;- paste(x[2],x[1],x[3]) } </code></pre> <p>and the second generating the name of model as a string</p> <pre><code>get.model.name &lt;- function(x) { x &lt;- deparse(substitute(x)) } </code></pre> <p>With these, I create a summary table</p> <pre><code>model.list &lt;- list(model1, model2) AIC.data &lt;- lapply(X = model.list, FUN = AIC) AIC.data &lt;- as.numeric(AIC.data) model.models &lt;- lapply(X = model.list, FUN = get.model) model.summary &lt;- cbind(model.models, AIC.data) model.summary &lt;- as.data.frame(model.summary) names(model.summary) &lt;- c("Model", "AIC") model.summary$AIC &lt;- unlist(model.summary$AIC) rm(AIC.data) model.summary[order(model.summary$AIC),] </code></pre> <p>Which all works fine. I'd like to add the model name to the table using get.model.name</p> <pre><code>x &lt;- get.model.name(model1) </code></pre> <p>Which gives me "model1" as I want.</p> <p>So now I apply the function to the list of models</p> <pre><code>model.names &lt;- lapply(X = model.list, FUN = get.model.name) </code></pre> <p>but now instead of <em>model1</em> I get <em>X[[1L]]</em></p> <p>How do I get <em>model1</em> rather than <em>X[[1L]]</em>?</p> <p>I'm after a table that looks like this:</p> <pre><code> Model Formula AIC model1 y ~ x 11.89136 model2 y ~ x + I(x^2) + I(x^3) 15.03888 </code></pre>
    singulars
    1. PO
      singulars
      1. PTAnswer
    1. This table or related slice is empty.
    1. PTQuestion
    1. USArun
    1. USIsaiah
    plurals
    1. PL
      singulars
      1. LTLinked
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. PO
      singulars
      1. PTAnswer
    2. PO
      singulars
      1. PTAnswer
    3. PO
      singulars
      1. PTAnswer
    1. VO
      singulars
      1. POUsing lapply on a list of models
      1. USBen
      1. VTFavorite
    1. COyou've created an unnamed list in `model.list` and you're passing each element of this list to `get.model.name`. So, `X[[1]]` is indeed passed the first time, and it fetches what you've asked for. @baptiste's overcomes this issue by creating simply a named list and avoiding the complications. If his solution is not what you're looking for, you'll have to explain, perhaps, why you're doing it this way.
      singulars
      1. POUsing lapply on a list of models
      1. USArun
    2. COThanks Arun and baptiste. I've added what I'm trying to produce.
      singulars
      1. POUsing lapply on a list of models
      1. USIsaiah
    3. COThanks again. I've just replaced `model.list <- list(model1, model2)` with `model.list <- list(model1=model1, model2=model2)` after reading your comments and it works well now.
      singulars
      1. POUsing lapply on a list of models
      1. USIsaiah
 

Querying!

 
Guidance
A row detail

Detail views are divided into sections. All the information in the data section comes from columns in the selected row. The other sections display data from other, related rows.

Related data can be related in a to-one or a to-many fashion. Captions of data related in a to-many fashion link to a list view showing a filtered view of the table.

Try moving around until you find a non-empty to-many entry and click on the label to get to one. You can move back to the root by clicking on the database name in the header.

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