我希望能够绘制使用glm()R中的函数拟合的参数估计的轮廓偏差.在估计所有其他参数之后,轮廓偏差是所讨论的参数估计的不同值的偏差函数.我需要在拟合参数周围绘制几个值的偏差,以检查二次偏差函数的假设.
我的模型正在预测罪犯的重新定罪.该公式的形式为:
reconv ~ [other variables] + sex,其中reconv是二元是/否因子,sex是二元男/女因子.我想绘制性别=女性估计参数的轮廓偏差(性别=男性是参考水平).
该glm()函数将参数估计为-0.22,std误差为0.12.
[我问的是这个问题,因为我找不到答案,但是我把它解决了,并想发布一个对别人有用的解决方案.当然,欢迎提供额外的帮助.:-)]
见profileModel包由扬Kosmidis.他在R Journal(R News中会出现)中有一篇论文说明了这个包:
Ioannis Kosmidis.profilemodel R软件包:具有线性预测变量的模型的分析目标.R News,8(2):12-18,2008年10月.
PDF在这里(整个时事通讯).
请参阅包中的?profile.glm(和example("profile.glm"))MASS- 我认为它会执行您想要的所有操作(默认情况下不会加载,但是它被提及?profile,可能是您看到的第一个...)(注意,配置文件是通常在有符号的平方根刻度上绘制,以便真正的二次曲线将显示为直线.)
| 归档时间： |
|
| 查看次数： |
1130 次 |
| 最近记录： |