参考McSpatial包里面的

1. qregbmat()

qregbmat(form,taumat=seq(.10,.90,.10),graphb=TRUE,graph.factor=FALSE,data=NULL)

2. qregsim1()

qregsim1(formall, formx, bmat, taumat, xvalues=NULL, ytarget=NULL,

xcolors=NULL,graphx=TRUE, graphy=TRUE, graphsim=TRUE, histogram=FALSE,

histfreq=FALSE, yname=NULL, xname=NULL, nsim=0, bwadjust=1,

legloc="topright",data=NULL)

还有更多的可以参考：http://cran.r-project.org/web/packages/McSpatial/McSpatial.pdf