如何从GLM(svyglm)模型对象中提取离散参数?
提取svyglm模型的离散参数
对于你的svyglm模型对象my_lm,要提取摘要中显示的离散参数并保存为变量,有两种直接的方法:
从模型摘要对象中提取:
disp_param <- summary(my_lm)$dispersion这个变量的值就是摘要里标注的
(Dispersion parameter for gaussian family taken to be xxx)对应的数值。针对高斯族模型,离散参数等于残差标准差的平方,也可以用
sigma()函数计算:disp_param <- sigma(my_lm)^2两种方法得到的结果完全一致。
提取完成后,直接调用disp_param就能查看或使用这个变量,无需仅通过打印摘要来获取数值。
内容的提问来源于stack exchange,提问作者Forklift17
相关产品推荐
相关产品推荐

