另一个rxGlm()问题...
已使用此功能安装了GLM,
MyGlm <- rxGlm(Response.Variable ~ 1 + F1 + F2 + F3,
family = gaussian,
data = my.df,
dropFirst = TRUE,
pweights = "myWeights"
)
我知道如何从MyGlm对象中提取某些信息,例如:
MyGlm$coefficients
MyGlm$coef.std.error
以此类推
我想要的是一个向量,它显示每个要素的每个层次的GLM拟合中使用的权重(以“系数”和“ coef.std.error”列出这些的顺序相同)。对于与该因子水平相对应的所有行,这只是原始数据中权重变量的总和。
因此,例如,如果我正在查看因子F1,并且我想知道级别6的权重,则希望这是原始数据中所有行的权重变量的总和(我。 df),其中F1取自因子级别6的值
如果这有意义...
不知道该怎么做。有人可以帮忙吗? 谢谢 艾伦