Function to compute B and S matrices for moment equations of Vardi's method (1996). It's not particularly efficient, but it works.
1  vardi.iteration(A, yBar, lambda, B, S)

A 
routing matrix (m x k) 
yBar 
numeric vector of mean link loads (length m) 
lambda 
value of lambda from last iteration 
B 
B matrix computed by

S 
S matrix computed by

numeric vector of length k with updated lambda
Y. Vardi. Network tomography: estimating sourcedestination traffic intensities from link data. Journal of the American Statistical Association, 91:365377, 1996.
Other vardi: vardi.algorithm
;
vardi.compute.BS
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.