Thanks Rick, a very interesting article. How about using the DESIGN function to tabulate the data? . print (design(BP_Status)`* design(Sex)) [ rowname=(unique(BP_Status)) colname=(unique(Sex)) ];