稳健回归啊——
用robust 和 robustbase:
library(robustbase)
glmrob0<-glmrob(cbind(NumDeath, NumTested-NumDeath)~log(Conc),family=binomial,data=set6,method="Mqle")
library(robust)
glmrob4<-glmRob(cbind(NumDeath, NumTested-NumDeath)~log(Conc),family=binomial(),data=set6)