R: может ли `assign` использоваться для моделей?

Question

Может ли функция assign в R использоваться для назначения имен объектов для моделей (например, моделей gls)? Я предполагаю, что не б/с я получаю предупреждение:

> Warning messages: 
In assign(paste0(deparse(substitute(mod)), "_", i, j), update(mod, : 
    only the first element is used as variable name 

В результате объекты не создаются

Есть ли способ сделать это?

Вот мой fucntion код, если это поможет:

#Choose best corARMA structure for model of choice: 
    corARMA.chooser <- function(mod,min = 0,max = 3) { 
    #This function creates 1. object for each combo of ARMA {0:3} 2. AIC table comparing all of these models 
    mod <- get('mod') 
    aic.arma <- AIC(mod) 
    ps <- 0 
    qs <- 0 
    for(i in min:max) { 
     js <- if(i == 0) c(1:max) else c(min:max) 
     for(j in js) { 
     arma <- corARMA(p = i, q = j) 
     assign(paste0(deparse(substitute(mod)),'_',i,j), update(mod, .~., correlation = arma), envir = .GlobalEnv) 
     aic.arma <- c(aic.arma, AIC(get(paste0(deparse(substitute(mod)),'_',i,j)))) 
     ps <- c(ps, i) 
     qs <- c(qs, i) 
     } 
     aic.arma.out <- data.frame(ps, qs, aic.arma) 
     aic.arma.out 
    } 
    } 

Update:

Я попытался использовать подход списка, но я получаю сообщение об ошибке:

Error in names(mod.list) <- c(names(mod.list), paste0(deparse(substitute(mod)), : 
'names' attribute [1275] must be the same length as the vector [1] 

Answer 1

EDIT: что на самом деле разрывает ваше имя переменной, это строка mod <- get('mod'), где вы перезаписываете свой именованный экземпляр mod, почему вы это сделать? Если изменить функцию этого он ведет себя, как я бы ожидать, что это:

corARMA.chooser <- function(modIn,min = 0,max = 3) { 
#This function creates 1. object for each combo of ARMA {0:3} 2. AIC table comparing all of these models mod <- get('modIn') aic.arma <- AIC(modIn) ps <- 0 qs <- 0 for(i in min:max) { 
js <- if(i == 0) c(1:max) else c(min:max) 
for(j in js) { 
    arma <- corARMA(p = i, q = j) 
    browser() 
    assign(paste0(deparse(substitute(modIn)),'_',i,j), update(mod, .~., correlation = arma), envir = .GlobalEnv) 
    aic.arma <- c(aic.arma, AIC(get(paste0(deparse(substitute(mod)),'_',i,j)))) 
    ps <- c(ps, i) 
    qs <- c(qs, i) 
    } 
aic.arma.out <- data.frame(ps, qs, aic.arma) 
aic.arma.out 
} 
} 

надеюсь, что это то, что вы пытаетесь достичь.

Answer 2

Все еще не уверен, почему код работает в одиночку, но не в функции, но ясно, что deparse (substitute (mod)) по какой-то причине вытягивает mod отдельно ко всем его частям сначала в функции, vs. просто создавая имя самого объекта.

Вот мой новый код, который работает:

corARMA.chooser <- function(mod,p = 1,q = 0) { 
    #This function creates 1. object for each combo of ARMA {0:3} 2. AIC table comparing all of these models 
    mod.list <- NULL 
    nms <- NULL 
    aic.arma <- AIC(mod) 
    ps <- 0 
    qs <- 0 
    for(i in c(p)) { 
     js <- if(i == 0) c(q[q>0]) else c(q) 
     for(j in c(js)) { 
     arma <- corARMA(p = i, q = j) 
     mod.list <- c(mod.list, list(update(mod, .~., correlation = arma))) 
     names(mod.list) <- c(names(mod.list), paste0(deparse(substitute(mod)),'_',i,j)) 
     aic.arma <- c(aic.arma, AIC(eval(parse(text=(paste0('mod.list$',deparse(substitute(mod)),'_',i,j)))))) 
     ps <- c(ps, i) 
     qs <- c(qs, j) 
     } 
    } 
     assign(paste0(deparse(substitute(mod)),'_','ARMA'),mod.list, envir = .GlobalEnv)  
     aic.arma.out <- data.frame(p = ps, q = qs, AIC = aic.arma) 
     aic.arma.out 
    }