It appears that HSAURtable only works on two dimensional tables. Is this correct? For example, here is HairEyeColor: , , Sex = Male Eye Hair Brown Blue Hazel Green Black 32 11 10 3 Brown 38 50 25 15 Red 10 10 7 7 Blond 3 30 5 8 , , Sex = Female Eye Hair Brown Blue Hazel Green Black 36 9 5 2 Brown 81 34 29 14 Red 16 7 7 7 Blond 4 64 5 8 HSAURtable(HairEyeColor) only uses the first matrix above > HSAURtable(HairEyeColor) $xname [1] "HairEyeColor" $pkg NULL $varnames [1] "Hair" "Eye" "Sex" $data [,1] [,2] [,3] [,4] [,5] [1,] " " "Brown" "Blue" "Hazel" "Green" [2,] "Black" "32" "11" "10" "3" [3,] "Brown" "38" "50" "25" "15" [4,] "Red" "10" "10" "7" "7" [5,] "Blond" "3" "30" "5" "8" attr(,"class") [1] "tabtab"