我想在一组360个观察中对~150个固定效应变量和7个随机效应变量进行模型选择。我决定使用Lasso程序进行混合模型，使用glmmLasso。我做了一些研究，找不到类似模型的例子但没有成功。以下是我的数据示例：
> str(RHI_12)
'data.frame': 350 obs. of 164 variables:
$ RHI_counts_12 : int 0 14 1 3 2 2 2 0 0 1 ...
$ Site : Factor w/ 6 levels "14_metzerlen",..: 1 1 1 1 1 1 1 1 1 1 ...
$ Location : Factor w/ 30 levels "1","2","3","4",..: 1 2 3 4 5 6 7 8 9 10 ...
$ Dist_roost : num 0.985 0.88 0.908 0.888 0.89 ...
$ Natural_light : num -0.194 -0.194 -0.194 -0.194 -0.194 ...
$ Mean_wind : num 0.836 0.836 0.836 0.836 0.836 ...
$ Mean_temp : num -0.427 -0.427 -0.427 -0.427 -0.427 ...
$ Day : num -0.993 -0.993 -0.993 -0.993 -0.993 ...
$ Artificial_light: num -0.2016 -0.2016 0.0772 -0.2016 -0.2016 ...
$ WBdi : num 1.14 1.14 1.14 1.14 1.14 ...
$ WCdi : num 1.49 1.49 1.49 1.49 1.47 ...
... (many more fixed-effect variables)
响应变量是计数（RHI_counts_12）。
我的问题是关于模型中随机效应变量的结构。 我有2个分类随机效应变量（＆＃34;网站＆＃34;和＆＃34;位置＆＃34 ;;＆＃34;位置＆＃34;嵌套在＆＃34;网站＆＃34;）和5个数字随机效应变量。我已经像这样构建了我的模型（仅使用固定效应变量的样本）：
lasso1<-glmmLasso(RHI_counts_12 ~ Artificial_light+WBdi+WCdi+BUdi+FOdi+TIdi, list(Site=~1,Location=~1+Dist_roost+Natural_light+Mean_wind+Mean_temp+Day),
lambda = 500,family = poisson(link = log), data = RHI_12)
如果我有这两种分类嵌套随机效应，我完全不相信构建随机效应的正确方法。我希望有一个模型，其中嵌套在Site中的位置，我不认为这是我得到的。以下是随机效果的输出（在此输出中，＆＃34; Loc＆＃34;代表位置，＆＃34; siteName＆＃34;代表网站）：
Random Effects:
StdDev:
[[1]]
siteName
siteName 1.180514
[[2]]
Loc Loc:Dist_roost Loc:Natural_light Loc:Mean_wind
Loc 1.15105859 -0.66317669 -0.35354821 -0.10805268
Loc:Dist_roost -0.66317669 1.42601945 0.46004662 -0.42795987
Loc:Natural_light -0.35354821 0.46004662 0.49532786 -0.15485395
Loc:Mean_wind -0.10805268 -0.42795987 -0.15485395 0.76175417
Loc:Mean_temp 0.02677276 0.03961902 -0.01431360 -0.03649499
Loc:Day 0.03756960 -0.02081360 0.02520654 -0.12082652
Loc:Mean_temp Loc:Day
Loc 0.02677276 0.03756960
Loc:Dist_roost 0.03961902 -0.02081360
Loc:Natural_light -0.01431360 0.02520654
Loc:Mean_wind -0.03649499 -0.12082652
Loc:Mean_temp 0.36923939 -0.08311209
Loc:Day -0.08311209 0.56876662
你认为这是对的吗？我无法使用＆＃34;位置＆＃34;来构建此模型。嵌套在＆＃34; Site＆＃34; （以及所有其他随机因素也将嵌套在＆＃34; Site＆＃34;。）我尝试了许多不同的方法但没有成功。
我已经非常感谢你读过我以及对glmmLasso中随机效果结构的任何建议！ ： - ）
托马斯