m <- 6 n <- 22 counts1 <- structure(c(54, 18, 32, 39, 32, 69, 54, 41, 156, 68, 83, 83, 88, 51, 67, 70, 19, 6, 5, 15, 44, 71, 8, 1, 8, 28, 7, 6, 9, 3, 21, 12, 13, 30, 8, 11, 14, 16, 7, 4, 2, 5, 6, 28, 84, 49, 21, 131, 34, 62, 28, 68, 63, 102, 49, 255, 41, 44, 110, 72, 7, 9, 6, 7, 63, 108, 1, 0, 0, 2, 0, 0, 0, 0, 0, 0, 1, 2, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 2, 1, 0, 2, 0, 0, 1, 2, 1, 2, 3, 5, 0, 3, 2, 3, 2, 0, 2, 0, 0, 1, 230, 158, 236, 440, 310, 242, 259, 394, 341, 279, 323, 692, 247, 156, 392, 233, 180, 87, 41, 134, 416, 304), .Dim = c(22, 6)) counts2 <- structure(c(54, 19, 29, 38, 32, 62, 52, 40, 157, 69, 81, 82, 91, 56, 65, 69, 17, 6, 5, 16, 45, 76, 8, 2, 9, 27, 6, 6, 8, 4, 21, 12, 12, 27, 11, 10, 13, 13, 6, 4, 2, 5, 5, 29, 70, 49, 27, 123, 33, 71, 33, 63, 62, 94, 50, 250, 36, 46, 112, 73, 4, 11, 6, 8, 66, 119, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 2, 0, 3, 0, 0, 0, 2, 3, 1, 4, 4, 0, 2, 2, 2, 1, 1, 1, 0, 1, 2, 240, 163, 224, 432, 318, 247, 254, 370, 336, 280, 299, 694, 239, 166, 394, 233, 179, 83, 39, 120, 427, 295), .Dim = c(22, 6)) x <- c(-0.508598071561408, -0.34350796037553, -0.333309343468294, 0.369757809574268, 0.266496813388506, 0.313665416584471, 0.382506080708313, 1.10788270823545, -0.439119993880865, 0.10778083776965, 0.712048889523366, 1.42531465947316, 0.329600755502027, 0.22761458642967, 1.10278339978184, 0.346173507976285, -1.45069530836731, -1.648293510945, -1.83696792372886, -1.9638132215126, 0.260122677821484, 1.57255719107138) logoffset <- c(5.63622008382947, 5.53182371175965, 5.60856770603342, 5.80107226615992, 5.22748275348865, 5.487632386555, 5.38581645358778, 5.23270672529503, 6.21468809522236, 5.87636228995183, 5.30014884836085, 5.57971825422116, 5.38634813881936, 5.40239752019213, 5.29597574960174, 5.559542088522, 5.54351741348397, 5.24691918666019, 4.72633258250038, 5.69651967794903, 5.43523163510263, 5.31025513055287)