错误: UndefVarError: y未定义

Question

我正试图在julia中运行一个混合效果模型(对于我的数据来说，R太慢了)，但是我一直收到这个错误。

我已经安装了DataArrays、DataFrames、MixedModels和RDatasets包，下面是教程--> http://dmbates.github.io/MixedModels.jl/latest/man/fitting/#Fitting-linear-mixed-effects-models-1。

以下是我的步骤：

using DataArrays, DataFrames , MixedModels, RDatasets 

我收到这些警告

警告:方法定义==(Base.Nullable{S}，Base.Nullable{T})位于nullable.jl:238的模块基础上，在nullable.jl:238中覆盖，在nullable.jl:238的模块NullableArrays中覆盖警告: /home/home/.julia/v0.6/DataFrames/src/statsmodels/formula.jl:352的DataFrames模块中的方法定义model_response(DataFrames.ModelFrame)覆盖在/home/home/.julia/v0.6/MixedModels/src/pls.jl:65.的MixedModels模块中警告:方法定义==(Base.Nullable{S}，Base.Nullable{T})位于nullable.jl:238的模块基础上，在nullable.jl:238中覆盖，在nullable.jl:238的模块NullableArrays中覆盖警告:方法定义==(Base.Nullable{S}，Base.Nullable{T})位于nullable.jl:238的模块基础上，在nullable.jl:238中覆盖，在nullable.jl:238的模块NullableArrays中覆盖警告: /home/home/.julia/v0.6/DataFrames/src/statsmodels/formula.jl:352的DataFrames模块中的方法定义model_response(DataFrames.ModelFrame)覆盖在/home/home/.julia/v0.6/MixedModels/src/pls.jl:65.的MixedModels模块中警告: /home/home/.julia/v0.6/DataFrames/src/statsmodels/formula.jl:352的DataFrames模块中的方法定义model_response(DataFrames.ModelFrame)覆盖在/home/home/.julia/v0.6/MixedModels/src/pls.jl:65.的MixedModels模块中

我从lme4包获得一个R数据集(在本教程中使用)

inst = dataset("lme4", "InstEval")

julia> head(inst)
6×7 DataFrames.DataFrame
│ Row │ S   │ D      │ Studage │ Lectage │ Service │ Dept │ Y │
├─────┼─────┼────────┼─────────┼─────────┼─────────┼──────┼───┤
│ 1   │ "1" │ "1002" │ "2"     │ "2"     │ "0"     │ "2"  │ 5 │
│ 2   │ "1" │ "1050" │ "2"     │ "1"     │ "1"     │ "6"  │ 2 │
│ 3   │ "1" │ "1582" │ "2"     │ "2"     │ "0"     │ "2"  │ 5 │
│ 4   │ "1" │ "2050" │ "2"     │ "2"     │ "1"     │ "3"  │ 3 │
│ 5   │ "2" │ "115"  │ "2"     │ "1"     │ "0"     │ "5"  │ 2 │
│ 6   │ "2" │ "756"  │ "2"     │ "1"     │ "0"     │ "5"  │ 4 │

如本教程所示，我运行该模型

m2 = fit!(lmm(y ~ 1 + dept*service + (1|s) + (1|d), inst))

并获取

ERROR: UndefVarError: y not defined
  Stacktrace:
  [1] macro expansion at ./REPL.jl:97 [inlined]
  [2] (::Base.REPL.##1#2{Base.REPL.REPLBackend})() at ./event.jl:73

当我使用DataFrames包中的"readtable“加载自己的数据时，也会发生同样的情况。

我正在运行Julia0.6.0，所有的软件包都是新安装的。我的系统是arch linux 4.11.7-1和所有最新的软件包。Julia安装时没有问题，但是有些包会发出警告(见上文)。

Answer 1

使用@formula宏：

julia> fit!(lmm(@formula(Y ~ (1 | Dept)), inst), true)
f_1: 250160.38873 [1.0]
f_2: 250175.99074 [1.75]
f_3: 250123.06531 [0.25]
f_4: 250602.3424 [0.0]
f_5: 250137.66303 [0.4375]
f_6: 250129.76244 [0.325]
f_7: 250125.94066 [0.280268]
f_8: 250121.15016 [0.23125]
f_9: 250119.12389 [0.2125]
f_10: 250114.7257 [0.175]
f_11: 250105.61264 [0.1]
f_12: 250602.3424 [0.0]
f_13: 250107.52714 [0.118027]
f_14: 250106.36924 [0.107778]
f_15: 250105.04638 [0.0925]
f_16: 250104.72722 [0.085]
f_17: 250104.93086 [0.0749222]
f_18: 250104.70046 [0.0831588]
f_19: 250104.70849 [0.0839088]
f_20: 250104.69659 [0.0824088]
f_21: 250104.69632 [0.0822501]
f_22: 250104.69625 [0.0821409]
f_23: 250104.69625 [0.0820659]
f_24: 250104.69624 [0.0821118]
f_25: 250104.69624 [0.0821193]
f_26: 250104.69624 [0.082111]
f_27: 250104.69624 [0.0821118]
f_28: 250104.69624 [0.0821117]
f_29: 250104.69624 [0.0821118]
f_30: 250104.69624 [0.0821118]
Linear mixed model fit by maximum likelihood
 Formula: Y ~ 1 | Dept
     logLik        -2 logLik          AIC             BIC       
 -1.25052348×10⁵  2.50104696×10⁵  2.50110696×10⁵  2.50138308×10⁵

Variance components:
              Column     Variance   Std.Dev.  
 Dept     (Intercept)  0.011897242 0.10907448
 Residual              1.764556375 1.32836605
 Number of obs: 73421; levels of grouping factors: 14

  Fixed-effects parameters:
             Estimate Std.Error z value P(>|z|)
(Intercept)   3.21373  0.029632 108.455  <1e-99

这些警告通常是"Julia (和包生态系统)仍在不断变化“的消息。但我不知道文档是否总是与代码同步。