The ecodist package has a convenience function full() that converts a dist object to a symmetric matrix. After that subsetting works normally. lower() performs the reverse operation. Sarah -- Sarah Goslee http://www.functionaldiversity.org