我有两个模型的输出 - 第一个模型使用glm运行，第二个模型使用lmekin运行。两个模型的输入是相同的，除了lmekin模型拟合了亲属矩阵。
我想从glm模型中获取系数，并使用来自使用lmekin的模型运行的p值计算这些值的SE。
这可能吗？我知道你可以使用未缩放的协方差矩阵的对角元素的平方根来手动计算SE。因此，以下代码是否会成功实现我想要实现的目标？
beta <- glm_model$coefficients
nvar <- length(beta)
nfrail <- nrow(lmekin_model$var) - nvar
se <- sqrt(diag(lmekin_model$var)[nfrail + 1:nvar])