当将数据名称作为另一个函数中的参数应用于gls函数时，将不会执行。

Question

使用包nlme，我有一个模型：

gls(temp.avg ~ year, data = df, method = 'ML')

这正是我所期望的。但是，当我在一个函数中创建带有gls函数的数据参数的模型时，我的模型不再使用数据'df‘，而是只打印"dat“。以下是功能：

function(dat) { gls(temp.avg ~ year, data = dat, method = 'ML') }

下面是当我(首先)在包装器函数之外创建模型时查看模型时的样子，然后(第二个)在函数内部看到"Data“行

Generalized least squares fit by maximum likelihood
  Model: temp.avg ~ year 
  Data: df 
  Log-likelihood: -3877.052

Coefficients:
   (Intercept)         (year) 
  15.135135363   -0.008796849 

Degrees of freedom: 1116 total; 1114 residual
Residual standard error: 7.807791

##########################

Generalized least squares fit by maximum likelihood
  Model: temp.avg ~ year
  Data: dat 
  Log-likelihood: -3877.052

Coefficients:
   (Intercept)         (year)  
  15.135135363   -0.008796849 

Degrees of freedom: 1116 total; 1114 residual
Residual standard error: 7.807791 

我真希望它不这么做。

如何才能让函数继承我指定的dat与"dat“本身？

Answer 1

library(nlme)
fun <- function(dat) { 
  res <- gls(follicles ~ sin(2*pi*Time) + cos(2*pi*Time), dat) 
  res$call$data <- substitute(dat)
  res
}
fun(Ovary)
#Generalized least squares fit by REML
#  Model: follicles ~ sin(2 * pi * Time) + cos(2 * pi * Time) 
#  Data: Ovary 
#  Log-restricted-likelihood: -898.434
#
#Coefficients:
#       (Intercept) sin(2 * pi * Time) cos(2 * pi * Time) 
#        12.2155822         -3.3396116         -0.8697358 
#
#Degrees of freedom: 308 total; 305 residual
#Residual standard error: 4.486121