这是源代码：library(Zelig)library(MASS)
d<-data.frame(
  y=crabs$sp,
  x1=crabs$FL,
  x2=crabs$RW
)
r<-zelig(y~x1+x2,model="probit",data=d)
summary(r)
options<-par(mfrow=c(2,2),mar=c(4,4,3,2))
plot(r)
par(options)
错误在plot(r)之后： No simulated quantities of interest found