[英]R : Pass argument to glm inside an R function
我想习惯于解决R中的问题。我想在函数内调用函数glm()
但它不起作用,显然是出于范围限制的原因,我没有设法用函数assign()
或eval()
进行修复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()
ao (y = Prop, x = Age, dataset = mydata, weights = Total)
这不起作用,但是以下起作用:
ao (y = mydata$Prop, x = mydata$Age, dataset = mydata, weights = mydata$Total)
有谁知道该怎么办?
任何帮助将不胜感激 !!!
顺便说一句,这是如何用我正在使用的数据集复制我的问题
library("MASS")
data(menarche)
mydata <- menarche
mydata$Prop <- mydata$Menarche / mydata$Total
替代解决方案(@DWin建议)。
function(y, x, dataset, weights){
f <- substitute(glm(y~x, data=dataset, weights=weights, family=binomial))
logLik(eval(f))
}
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
我建议使用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 。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.