没有人回答 啊？

就是做稳健性回归！stata里面就有这样的操作，不知R有没有啊！

try R packages: robust and robustbase. Below are my R codes:
library(robustbase)
# method="Mqle" fits a generalized linear model using Mallows or Huber type robust estimators, as described in Cantoni and Ronchetti (2001).
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)

实在感谢！早已经沉了很久的帖子了，又被翻出来，大家一起学习！谢谢！

glmrob, dlm 都是平稳模式的时间系列函数，dlm 还有电子书可以参考。