你怎么还没解决?实在看不下去了。
PURPOSE: computes spatial durbin model estimates
(I-rho*W)y = a + X*B1 + W*X*B2 + e, using sparse algorithms
---------------------------------------------------
USAGE: results = sdm(y,x,W,info)
where: y = dependent variable vector
x = independent variables matrix
W = contiguity matrix (standardized)
info = an (optional) structure variable with input options:
info.rmin = (optional) minimum value of rho to use in search
info.rmax = (optional) maximum value of rho to use in search
info.convg = (optional) convergence criterion (default = 1e-8)
info.maxit = (optional) maximum # of iterations (default = 500)
info.lflag = 0 for full lndet computation (default = 1, fastest)
= 1 for MC lndet approximation (fast for very large problems)
= 2 for Spline lndet approximation (medium speed)
info.order = order to use with info.lflag = 1 option (default = 50)
info.iter = iterations to use with info.lflag = 1 option (default = 30)
info.lndet = a matrix returned by sar, sar_g, sarp_g, etc.
containing log-determinant information to save time