我是R和统计建模的新手，我希望在r中使用 lmmlasso 库来拟合混合效果模型，只选择~300个可能变量中最好的固定效果 .
对于这个模型，我想包括固定拦截，随机效应和随机拦截 . 看一下CRAN的手册，我发现了以下几点：
x：维度为ntot x p的矩阵，包括固定效应协变量 . 截距必须包含在第一列中（1，...，1） . z：维度为ntot x q的随机效应矩阵 . 它必须是一个矩阵，即使q = 1 .
虽然很明显我需要为固定拦截做些什么，但我不太确定如何同时包含随机拦截和效果 . 它与固定矩阵完全相同，我在第一列中包含（1 ... 1）吗？
除此之外，我正在寻找验证我得到的另一个数据集的模型 . 对于 lmmlasso ， lme4 中是否存在类似于 predict 的函数，可用于根据我得到的输出计算新的预测？或者，使用 lmer 使用 lmmlasso 返回的非零系数变量构建新模型是否可行/正确，然后在新模型上使用 predict ？
提前致谢 .