Comment stocker des "tableaux" de modèles statistiques?

Est-il un R structure de données dans laquelle je peux stocker un certain nombre de lm ou lmer ou gam objets? J a enfermé les tableaux, et on peut mettre beaucoup de choses dans les cellules d'une telle boîte de tableau. Je pense que c'est ce que je cherche dans R.

J'ai essayé les listes et les trames de données, en vain; j'ai pensé listes.

> testlist <- list()
> testlist[1] <- subject1.2008.gam
Warning message:
In testlist[1] <- subject1.2008.gam :
  number of items to replace is not a multiple of replacement length
>

Sinon, est-il un moyen de créer et d'utiliser un nom de variable sur la GAUCHE de <-?

Enfin, peut-être que vous avez une meilleure idiome pour moi à prendre en compte. Je suis en train de créer une collection de GAM modèles sur un ensemble de sujets et des années, par exemple. Plus tard, je veux être capable de tracer ou de prédire à partir de ces modèles, donc je pense que j'ai besoin de garder le modèle complet autour. Parce que je veux être en mesure d'utiliser ce code avec différents ensembles de données plus tard, j'aimerais ne pas coder en dur le nom de la gam objets, ni leur nombre.

Moment, j'ai commencé par mettre le gam() appel dans une boucle, je pense que l'un des apply() fonctions pourrait fonctionner mieux, mais j'ai encore besoin d'un endroit pour stocker le résultat.

source d'informationauteur Bill

  1. 18

    Vous avez besoin de la [[ opérateur pour les listes, essayez

    testlist[[1]] <- subject1.2008.gam

    Les autres d'habitude astuce est que vous pouvez pré-allouer si vous savez comment beaucoup d'éléments que vous pourriez avoir, je le fais souvent

    testlist <- vector(mode="list", length=N)

    pour un N.

  2. 3

    Utilisation [[ pour accéder à la liste des éléments:

    library(mgcv)
set.seed(0) ## simulate some data... 
dat <- gamSim(1,n=400,dist="normal",scale=2)

mods <- vector(mode = "list", length = 3)
for(i in seq_along(mods)) {
    mods[[i]] <- gam(y ~ s(x0) + s(x1) + s(x2) + s(x3), data = dat)
}

    :

    > str(mods, max = 1)
List of 3
 $ :List of 43
  ..- attr(*, "class")= chr [1:3] "gam" "glm" "lm"
 $ :List of 43
  ..- attr(*, "class")= chr [1:3] "gam" "glm" "lm"
 $ :List of 43
  ..- attr(*, "class")= chr [1:3] "gam" "glm" "lm"
  3. 2

    Les autres réponses montrent comment utiliser un index et [[ ]] mais vous pouvez aussi faire quelque chose comme 

    x1  <- 1:10  ; y1  <-  30*x1 + rnorm(10)
x2  <- rnorm(20)  ; y2  <- 30*x2 + 100 + rnorm(20)
lm1 <- lm(y1 ~ x1); lm2 <- lm(y2 ~ x2) 

testlist <- list( A = lm1, Z = lm2 ) 
testlist$Z
testlist$Z$model$y2