我optim()在 R 中应用函数来获得具有双变量随机效应的模型中固定效应和随机效应的最大似然估计。该代码运行良好，但在执行涉及该optim()函数的代码行时非常慢。我想知道 R 中是否有任何专门的优化函数/包来最大化随机效应模型中的可能性，而不是通用优化函数optim()。