Logistic模型【DOC精选】.docVIP

Logistic模型 例： 为研究高压电线对牲畜的影响，R.Norell研究小的电流对农场动物的影响。实验中选择了7头牛，6种电击强度（0、1、2、3、4、5mA）。每头牛被电击30下，每种强度5下，按随机次序进行。然后重复整个实验，每头牛总共被电击60下。对每次电击，响应变量——嘴巴运动或出现、或不出现。表d5.1的数据给出每种电击强度70次试验中的响应总次数，试分析电击对牛的影响。 这里，响应变量是分类的，它只有两个值：出现及未出现，这种情况下可用logistic回归： ，其中x是电流强度。概率p取0到1之间的值。 命令： ##在d5.1中选取A1：C7区域，然后拷贝 nk=read.table(clipboard,header=T) nk x n k 1 0 70 0 2 1 70 9 3 2 70 21 4 3 70 47 5 4 70 60 6 5 70 63 y-cbind(nk$k, nk$n-nk$k) glm.logit-glm(y~ nk$x,family=binomial) summary(glm.logit) Call: glm(formula = y ~ nk$x, family = binomial) Deviance Residuals: 1 2 3 4