我正在寻求有关在 R 中编程新的随机效应方差-协方差矩阵/结构 (pdmat) 以用于 lme() 的建议？我已经检查了 lme 源代码（如果您想编写自己的 pdmat，则建议这样做），但仍然感到困惑。有人有好的参考吗？提前致谢！