R进行Logistic回归学习笔记

moonstone 发表于 2015-12-7 14:02
这是论坛对R语言编码显示的问题造成的，应该是infert$case,中间的“$”显示不全

我知道是怎么回事了。这里应该是infert$pred5吧，不是infert$pred6

xiaolovesfan 发表于 2015-12-7 14:25
我知道是怎么回事了。这里应该是infert$pred5吧，不是infert$pred6

谢谢，这个地方确实笔误了。^_^，已经修改

另外，对于$显示不正确的事情也已经修改，再看一下原始的编码。

具体解决方法是：同一行中不能出现两次美元号,不然就会被认为是特殊字符，显示为斜体。

也希望这个回复对论坛出现类似问题问题的坛友有点帮助吧

能问一下，Model5是什么意思嘛？是只对 education ="12+" 的子集建模吗？

moonstone 发表于 2015-12-7 14:32
谢谢，这个地方确实笔误了。^_^，已经修改

另外，对于$显示不正确的事情也已经修改，再看一下原始的编 ...

赞一个！看来是论坛里的常驻大神啊

xiaolovesfan 发表于 2015-12-7 14:36
能问一下，Model5是什么意思嘛？是只对 education ="12+" 的子集建模吗？

是，R没法像SAS里面的Proc logisitic过程中的by选项一样，根据by变量构建分层logisitc回归，只能选取亚集分别进行logisitic回归，这是SAS和R在完成logisitic回归上功能上最大的不一致。

moonstone 发表于 2015-12-7 14:48
是，R没法像SAS里面的Proc logisitic过程中的by选项一样，根据by变量构建分层logisitc回归，只能选取亚集 ...

我还有一个问题哈~这里只用education ="12+"，是只对education=12岁以上的建模是吗？education的其他组就没有建模？如果是的话，那么最后和model8的怎么有可比性呢？Model8是对所有education的子组都建模了啊。

xiaolovesfan 发表于 2015-12-7 19:51
我还有一个问题哈~这里只用education ="12+"，是只对education=12岁以上的建模是吗？education的其他组就 ...

非常感谢，示例只是一种举例，所以比较model5和model8确实没有特别明确的意义。

但是即使是示例，确实应该考虑实际意义。所以一开始的笔误应该是pred6而不是pred5，已经修改。

moonstone 发表于 2015-12-7 20:49
非常感谢，示例只是一种举例，所以比较model5和model8确实没有特别明确的意义。

但是即使是示例，确实 ...

运行差不多通过了。。。学R就应该这样

已经进一步完善“如何选取特定水平作为参考水平”的方法，也即采用stat程序包里面的relevel函数功能，这种方法比采用contrast函数更加实用，但是注意一点relevel函数只适用于factor变量，所以如果是数值型变量，需要先采用as.factor函数进行转换。具体函数已经在主函数中进行修改，这里再复制一下，强调一下relevel的功能：

1. infert$spontaneous <-
   
2.   as.factor(infert$spontaneous)
   
3. infert$induced <-
   
4.   as.factor(infert$induced)
   
5. 
   
6. infert$spontaneous <-
   
7.   relevel(infert$spontaneous, ref = "2")
   
8. infert$induced <-
   
9.   relevel(infert$induced, ref = "1")
   
10. 
    
11. model6 <- glm(case ~ spontaneous+induced,
    
12.               data = infert)
    
13. summary(model6)

复制代码

激动地哭了！！！今天有幸看到这个帖子！解决了困扰我多年的问题！！！谢谢楼主！！！