Dear friends,
  Both vcov(your.model) and summary(B)$cov.unscaled,summary(B)$cov.scaled
works, and vcov is the function that i'm looking for.
 Thanks very much!
-------------
with kind regards
zhijie

	[[alternative HTML version deleted]]

