2

我正在尝试在循环内使用包中的lme函数。这适用于包中的函数，但会为. 如何循环类似于下面示例中的函数的函数？nlmelapplylmerlme4lmelmelmer

library("nlme")
library("lme4")

set.seed(1)
dt <- data.frame(Resp1 = rnorm(100, 50, 23), Resp2 = rnorm(100, 80, 15), Pred = rnorm(100,10,2), group = factor(rep(LETTERS[1:10], each = 10)))

## Syntax:
lmer(Resp1 ~ Pred + (1 |group), data = dt)
lme(Resp1 ~ Pred, random = ~1 | group, data = dt)

## Works for lme4
lapply(c("Resp1", "Resp2"), function(k) {
  lmer(substitute(j ~ Pred + (1 | group), list(j = as.name(k))), data = dt)})

## Does not work for nlme
lapply(c("Resp1", "Resp2"), function(k) {
lme(substitute(j ~ Pred, list(j = as.name(k))), random = ~1 | group, data = dt)})

# Error in UseMethod("lme") : 
# no applicable method for 'lme' applied to an object of class "call"

PS。我知道存在此解决方案，但我想使用一种方法直接在模型函数中替换响应变量，而不是使用附加函数对数据进行子集化。

4

2 回答 2

3

您还可以执行以下操作，而不是摆弄substituteand ：eval

lapply(c("Resp1", "Resp2"), function(r) {
   f <- formula(paste(r, "Pred", sep = "~"))
   m <- lme(fixed = f, random = ~ 1 | group, data = dt)
   m$call$fixed <- f
   m})

如果您想为建模函数提供不同的数据集，您可以使用相同的技巧：

makeModel <- function(dat) {
   l <- lme(Resp1 ~ Pred, random = ~ 1 | group, data = dat)
   l$call$data <- as.symbol(deparse(substitute(dat)))
   l
}

当我想从一个函数中生成一个模型并想在之后更新它时，我会经常使用这个片段。

于 2014-10-14T14:14:33.427 回答
1

正如@CarlWitthoft 建议的那样，添加eval到函数中将解决问题：

lapply(c("Resp1", "Resp2"), function(k) {
lme(eval(substitute(j ~ Pred, list(j = as.name(k)))), random = ~1 | group, data = dt)})

另请参阅@thothal 的替代方案。

于 2014-10-14T14:23:59.947 回答