Good day, It is implemented by the CRAN package multicon. The function is named popsd. But it does seem like something R should provide without creating a package dependency. -------------------------------------- Dario Strbenac University of Sydney Camperdown NSW 2050 Australia