R: передать аргумент glm внутри функции R

Question

Я пытаюсь привыкнуть к вопросам определения области в R. Я хотел бы вызвать функцию glm() внутри функции, но она не работает, по-видимому, по причинам, удалось установить с помощью функций assign() или eval().

Вот упрощенная версия:

ao <- function (y, x, phi = seq (0,1,0.1), dataset, weights) { 
    logLikvector <- rep(0,length(phi)) # vector of zeros to be replaced thereafter 
    for (i in 1:length(phi)) {   # loop to use glm() 
     fit <- glm (y ~ x, data = dataset, family = binomial, weights = weights)   
     logLikvector[i] <- logLik(fit)  # get log likelihood 
    } 
    logLikvector 
} 

Теперь я хочу использовать функцию ао() в моем наборе данных

ao (y = Prop, x = Age, dataset = mydata, weights = Total) 

Это не работает, но следующие работы:

ao (y = mydata$Prop, x = mydata$Age, dataset = mydata, weights = mydata$Total) 

Кто-нибудь знает, что делать?

Любая помощь будет принята с благодарностью!

Btw, вот как повторить мою проблему с набором данных я использую

library("MASS") 
data(menarche) 
mydata <- menarche 
mydata$Prop <- mydata$Menarche/mydata$Total 

Answer 1

Решение с заменой (предложение @DWin).

function(y, x, dataset, weights){ 
    f <- substitute(glm(y~x, data=dataset, weights=weights, family=binomial)) 
    logLik(eval(f)) 
} 

Answer 2

ao <- function (x, y, phi = seq (0,1,0.1), dataset, weights) { 
    logLikvector <- rep(0,length(phi)) 
    x <- dataset[,substitute(x)] 
    y <- dataset[,substitute(y)] 
    weights <- dataset[,substitute(weights)] 
     for (i in 1:length(phi)) {   # loop to use glm() 
     fit <- glm (y ~ x, data = dataset, family = binomial, weights = weights) 
     logLikvector[i] <- logLik(fit)  # get log likelihood 
    } 
    return(logLikvector) 
} 



library("MASS") 
data(menarche) 
mydata <- menarche 
mydata$Prop <- mydata$Menarche/mydata$Total 
ao(y = "Prop",x = "Age", dataset = mydata, weights = "Total") 


[1] -55.37763 -55.37763 -55.37763 -55.37763 -55.37763 -55.37763 
[7] -55.37763 -55.37763 -55.37763 -55.37763 -55.37763 

Answer 3

Я предлагаю создать формулу с paste и вызова функции с do.call.

ao <- function (y, x, phi = seq (0,1,0.1), dataset, weights) { 
    logLikvector <- rep(0,length(phi)) # vector of zeros to be replaced thereafter 
    for (i in 1:length(phi)) {   # loop to use glm() 
    f <- as.formula(paste(y, x, sep="~")) 
    fit <- do.call("glm", list(formula=f, data=as.name(dataset), 
        family="binomial", weights=as.name(weights))) 
    logLikvector[i] <- logLik(fit)  # get log likelihood 
    } 
    logLikvector 
} 

Затем вызовите его следующим образом:

ao("Prop", "Age", dataset="mydata", weights="Total") 

См https://stackoverflow.com/a/7668846/210673 для более подробной информации.