解释变量四列，sex,skin type,burn time,sun,都是分类变量。响应变量case，取值0或1.想做logstic 回归，看了一本书上的方法是先把数据分组，解释变量取值相同的归为一组，把case按分组计数整理好，形如
sex  skin  burn  sun  case  nocase
0       1      1       1       9       2

p<-case/(case+nocase)
把p嵌入数据中然后再做回归fit<-glm(p～sex,skin,burn,sun)

请问原始数据怎么分组计数？怎么做回归？

还有一个问题，the empirical logit is defined as epl=log{(case+0.5)/(nocase+0.5)}
why is the term"+0.5" included?

如果按照这个公式的话应该怎样写R代码做回归？