# # data # ----------------------- # chapter 1 - example 2 # ----------------------- ex1_2.greece<-data.frame( game=1:6 ) ex1_2.greece$opponent<- c("Portugal", "Spain" , "Russia", "France", "Chech Rep.", "Portugal" ) ex1_2.greece$scored <-c(2, 1, 1, 1, 0, 1) ex1_2.greece$conceded<-c(1, 1, 2, 0, 0, 0) ex1_2.greece$total <- ex1_2.greece$scored + ex1_2.greece$conceded # ----------------------- # chapter 1 - example 3 # ----------------------- ex1_3<-c(1245,5) # ----------------------- # chapter 1 - example 4 # ----------------------- "ex1_4.kobe.total" <- structure(list(SEASON = structure(as.integer(c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11)), .Label = c(" 1996-97", " 1997-98", " 1998-99", " 1999-00", " 2000-01", " 2001-02", " 2002-03", " 2003-04", " 2004-05", " 2005-06", " 2006-07"), class = "factor"), TEAM = structure(as.integer(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)), .Label = "LAL ", class = "factor"), GAMES = c(71, 79, 50, 66, 68, 80, 82, 65, 66, 80, 42), TIME = c(1530, 2600, 3753, 3812, 4053, 3817, 4130, 3736, 4042, 4100, 3923 ), FGMEAN = c(2.5, 5, 7.2, 8.4, 10.3, 9.4, 10.6, 7.9, 8.7, 12.2, 9.5), FGATTEMP = c(5.9, 11.6, 15.6, 17.9, 22.2, 20, 23.5, 18.1, 20.1, 27.2, 20.1), FGPERC = c(41.7, 42.8, 46.5, 46.8, 46.4, 46.9, 45.1, 43.8, 43.3, 45, 47.2), TPTMEAN = c(0.7, 1, 0.5, 0.7, 0.9, 0.4, 1.5, 1.1, 2, 2.3, 1.7), TPTATTEM = c(1.9, 2.8, 2, 2.2, 2.9, 1.7, 4, 3.3, 5.9, 6.5, 4.7), TPTPERC = c(37.5, 34.1, 26.7, 31.9, 30.5, 25, 38.3, 32.7, 33.9, 34.8, 36.6), FTMEAN = c(1.9, 4.6, 4.9, 5, 7, 6.1, 7.3, 7, 8.2, 8.7, 8.1 ), FTATTEMP = c(2.3, 5.8, 5.8, 6.1, 8.2, 7.4, 8.7, 8.2, 10.1, 10.2, 9.4), FTPERC = c(81.9, 79.4, 83.9, 82.1, 85.3, 82.9, 84.3, 85.2, 81.6, 85, 85.4), REBOFF = c(0.7, 1, 1.1, 1.6, 1.5, 1.4, 1.3, 1.6, 1.4, 0.9, 0.8), REBDEF = c(1.2, 2.1, 4.2, 4.7, 4.3, 4.1, 5.6, 3.9, 4.5, 4.4, 4.6), REBTOT = c(1.9, 3.1, 5.3, 6.3, 5.9, 5.5, 6.9, 5.5, 5.9, 5.3, 5.5), AST = c(1.3, 2.5, 3.8, 4.9, 5, 5.5, 5.9, 5.1, 6, 4.5, 5.6), TO2 = c(1.6, 2, 3.1, 2.8, 3.2, 2.8, 3.5, 2.6, 4.1, 3.1, 3.5), STL = c(0.7, 0.9, 1.4, 1.6, 1.7, 1.5, 2.2, 1.7, 1.3, 1.8, 1.2), BLK = c(0.3, 0.5, 1, 0.9, 0.6, 0.4, 0.8, 0.4, 0.8, 0.4, 0.5), PF = c(1.4, 2.3, 3.1, 3.3, 3.3, 2.9, 2.7, 2.7, 2.6, 2.9, 3), PPG = c(7.6, 15.4, 19.9, 22.5, 28.5, 25.2, 30, 24, 27.6, 35.4, 28.8), FGTOTAL = c(177.5, 395, 360, 554, 701, 749, 868, 516, 573, 978, 399), FGTOTATT = c(418.9, 916.4, 780, 1183, 1510, 1597, 1924, 1178, 1324, 2173, 845)), .Names = c("SEASON", "TEAM", "GAMES", "TIME", "FGMEAN", "FGATTEMP", "FGPERC", "TPTMEAN", "TPTATTEM", "TPTPERC", "FTMEAN", "FTATTEMP", "FTPERC", "REBOFF", "REBDEF", "REBTOT", "AST", "TO2", "STL", "BLK", "PF", "PPG", "FGTOTAL", "FGTOTATT"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11"), class = "data.frame", variable.labels = structure(c("", "", "", "Time Played (min*100+seconds)", "Field Goals per game", "Field Goals attempts per game", "Percentage of Field Goals", "3 Point Goals per game", "3 Point Goal Attempts per game", "3 Point percentage", "Free Throws per Game", "Free Throws Attempts per game", "Free Throws percentage", "Offensive Rebounts per game", "Defensive Rebounts per game", "Total Rebounts per game", "Assist per game", "", "Steals per game", "Blocks per game", "", "Points per game", "", ""), .Names = c("SEASON", "TEAM", "GAMES", "TIME", "FGMEAN", "FGATTEMP", "FGPERC", "TPTMEAN", "TPTATTEM", "TPTPERC", "FTMEAN", "FTATTEMP", "FTPERC", "REBOFF", "REBDEF", "REBTOT", "AST", "TO2", "STL", "BLK", "PF", "PPG", "FGTOTAL", "FGTOTATT"))) "ex1_4.kobe.fg" <- structure(list(SEASON = structure(as.integer(c(4, 5, 6, 7, 8, 9, 10, 11)), .Label = c(" 1996-97", " 1997-98", " 1998-99", " 1999-00", " 2000-01", " 2001-02", " 2002-03", " 2003-04", " 2004-05", " 2005-06", " 2006-07"), class = "factor"), TEAM = structure(as.integer(c(1, 1, 1, 1, 1, 1, 1, 1)), .Label = "LAL ", class = "factor"), GAMES = c(66, 68, 80, 82, 65, 66, 80, 42), FGTOTAL = c(554, 701, 749, 868, 516, 573, 978, 399), FGTOTATT = c(1183, 1510, 1597, 1924, 1178, 1324, 2173, 845)), .Names = c("SEASON", "TEAM", "GAMES", "FGTOTAL", "FGTOTATT"), row.names = c("4", "5", "6", "7", "8", "9", "10", "11"), class = "data.frame") # ----------------------- # chapter 1 - example 5 # ----------------------- "ex1_5.bodytemp" <- structure(list(temperature = c(96.3, 96.7, 96.9, 97, 97.1, 97.1, 97.1, 97.2, 97.3, 97.4, 97.4, 97.4, 97.4, 97.5, 97.5, 97.6, 97.6, 97.6, 97.7, 97.8, 97.8, 97.8, 97.8, 97.9, 97.9, 98, 98, 98, 98, 98, 98, 98.1, 98.1, 98.2, 98.2, 98.2, 98.2, 98.3, 98.3, 98.4, 98.4, 98.4, 98.4, 98.5, 98.5, 98.6, 98.6, 98.6, 98.6, 98.6, 98.6, 98.7, 98.7, 98.8, 98.8, 98.8, 98.9, 99, 99, 99, 99.1, 99.2, 99.3, 99.4, 99.5, 96.4, 96.7, 96.8, 97.2, 97.2, 97.4, 97.6, 97.7, 97.7, 97.8, 97.8, 97.8, 97.9, 97.9, 97.9, 98, 98, 98, 98, 98, 98.1, 98.2, 98.2, 98.2, 98.2, 98.2, 98.2, 98.3, 98.3, 98.3, 98.4, 98.4, 98.4, 98.4, 98.4, 98.5, 98.6, 98.6, 98.6, 98.6, 98.7, 98.7, 98.7, 98.7, 98.7, 98.7, 98.8, 98.8, 98.8, 98.8, 98.8, 98.8, 98.8, 98.9, 99, 99, 99.1, 99.1, 99.2, 99.2, 99.3, 99.4, 99.9, 100, 100.8), gender = structure(as.integer(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2)), .Label = c("male", "female"), class = "factor"), heart.rate = as.integer(c(70, 71, 74, 80, 73, 75, 82, 64, 69, 70, 68, 72, 78, 70, 75, 74, 69, 73, 77, 58, 73, 65, 74, 76, 72, 78, 71, 74, 67, 64, 78, 73, 67, 66, 64, 71, 72, 86, 72, 68, 70, 82, 84, 68, 71, 77, 78, 83, 66, 70, 82, 73, 78, 78, 81, 78, 80, 75, 79, 81, 71, 83, 63, 70, 75, 69, 62, 75, 66, 68, 57, 61, 84, 61, 77, 62, 71, 68, 69, 79, 76, 87, 78, 73, 89, 81, 73, 64, 65, 73, 69, 57, 79, 78, 80, 79, 81, 73, 74, 84, 83, 82, 85, 86, 77, 72, 79, 59, 64, 65, 82, 64, 70, 83, 89, 69, 73, 84, 76, 79, 81, 80, 74, 77, 66, 68, 77, 79, 78, 77))), .Names = c("temperature", "gender", "heart.rate" ), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96", "97", "98", "99", "100", "101", "102", "103", "104", "105", "106", "107", "108", "109", "110", "111", "112", "113", "114", "115", "116", "117", "118", "119", "120", "121", "122", "123", "124", "125", "126", "127", "128", "129", "130"), class = "data.frame") # ----------------------- # chapter 2 - example 1 # ----------------------- y <- c( 25, 30) n <- c(300,900) ex2_1.risk <- data.frame(y=y,n=n) rm(y) rm(n) # ----------------------- # chapter 2 - example 2 # ----------------------- y<-399; N<-845 ex2_2.kobe2006 <- data.frame(y=y,N=N) rm(y) rm(N) # ----------------------- # chapter 2 - example 3 # ----------------------- "ex2_3.wais" <- structure(list(wais = as.integer(c(9, 13, 6, 8, 10, 4, 14, 8, 11, 7, 9, 7, 5, 14, 13, 16, 10, 12, 11, 14, 15, 18, 7, 16, 9, 9, 11, 13, 15, 13, 10, 11, 6, 17, 14, 19, 9, 11, 14, 10, 16, 10, 16, 14, 13, 13, 9, 15, 10, 11, 12, 4, 14, 20)), senility = as.integer(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0))), .Names = c("wais", "senility" ), class = "data.frame", row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54")) # ----------------------- # chapter 2 - example 4 # ----------------------- ex2_4 <- 'see ex1_5.bodytemp' # ----------------------- # chapter 2 - example 5 # ----------------------- ex2_5 <- 'see exex2_3.wais' # ----------------------- # chapter 3 - example 1 # ----------------------- ex3_1 <- c(-1.76, 0.38, 1.23, -0.67, -0.47, -1.36, 1.41,-0.07, -1.23, 2.35) # ----------------------- # chapter 4 - example 1 # ----------------------- ex4_1.kobes.fg <- list( YEARS=8, y=c(554,701,749,868,516,573,978,399), N=c(1183,1510,1597,1924,1178,1324,2173, 845) ) # ----------------------- # chapter 5 - example 1 # ----------------------- ex5_1.softdrinks.lst <- list( n=25, time = c(16.68, 11.5, 12.03, 14.88, 13.75, 18.11, 8, 17.83, 79.24, 21.5, 40.33, 21, 13.5, 19.75, 24, 29, 15.35, 19, 9.5, 35.1, 17.9, 52.32, 18.75, 19.83, 10.75), distance = c(560, 220, 340, 80, 150, 330, 110, 210, 1460, 605, 688, 215, 255, 462, 448, 776, 200, 132, 36, 770, 140, 810, 450, 635, 150), cases = c( 7, 3, 3, 4, 6, 7, 2, 7, 30, 5, 16, 10, 4, 6, 9, 10, 6, 7, 3, 17, 10, 26, 9, 8, 4) ) ex5_1.softdrinks <- data.frame( time= ex5_1.softdrinks.lst$time , distance=ex5_1.softdrinks.lst$distance , cases=ex5_1.softdrinks.lst$cases ) # ----------------------- # chapter 5 - example 2 # ----------------------- ex5_2.tutors.lst <-list( n=25, TUTORS=4, grade=c(84, 58, 100, 51, 28, 89, 97, 50, 76, 83, 45, 42, 83,64, 47, 83, 81, 83, 34, 61, 77, 69, 94, 80, 55, 79), class=c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4) ) ex5_2.tutors <- data.frame( grade=ex5_2.tutors.lst$grade, class=ex5_2.tutors.lst$class ) # ----------------------- # chapter 5 - example 3 # ----------------------- # data in tabular format (no missing values) y = structure(.Data=c( 9, 14, 18, 29, 25, 26, 22, 25, 23, 24, 12, 13),.Dim = c( 2,2,3 ) ) y2 <- array( dim=c(3,2,2) ) y2[,,1] <- t(y[1,,]) y2[,,2] <- t(y[2,,]) ex5_3.spq.tab.lst <- list( LA=3, LB=2, K=2, y=y2) rm(y) rm(y2) # data in individual data format (unbalanced dataset - with missing values) y<-c( 1, 1, 9, 1, 2, 25, 1, 3, 23, 2, 1, 18, 2, 2, 22, 2, 3, 12, 1, 1, 14, 1, 2, 26, 2, 1, 29, 2, 2, 25) y<-matrix(y, nrow=10,ncol=3,byrow=TRUE) ex5_3.spq.indiv <- data.frame( g=y[,1], e=y[,2], y=y[,3]) rm(y) # ----------------------- # chapter 6 - example 1 # ----------------------- y<-c(1, 1, 1, 9, 1, 2, 1, 25, 1, 3, 1, 23, 2, 1, 1, 18, 2, 2, 1, 22, 2, 3, 1, 12, 1, 1, 2, 14, 1, 2, 2, 26, 2, 1, 2, 29, 2, 2, 2, 25) y<-matrix(y, nrow=10,ncol=4,byrow=TRUE) ex6_1.spq.indiv <- data.frame( g=y[,1], e=y[,2], ci=y[,3], y=y[,4]) rm(y) # ----------------------- # chapter 6 - example 2 # ----------------------- ex6_2.bioassay.lst <- list( n=24, y=c(68.8, 67.6, 68.1, 67.6, 69.0, 67.9, 68.6, 68.3, 61.4, 59.8, 62.3, 60.6, 60.9, 60.3, 61.6, 61.8, 53.5, 51.9, 53.6, 52.2, 53.8, 54.9, 54.1, 54.2), dose=c(0.025, 0.025, 0.025, 0.025, 0.025, 0.025, 0.025, 0.025, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.050, 0.100, 0.100, 0.100, 0.100, 0.100, 0.100, 0.100, 0.100), drug=c(1,1,1,1,2,2,2,2, 1,1,1,1,2,2,2,2, 1,1,1,1,2,2,2,2) ) ex6_2.bioassay <- data.frame(y=ex6_2.bioassay.lst$y,dose=ex6_2.bioassay.lst$dose,drug=ex6_2.bioassay.lst$drug) # ----------------------- # chapter 6 - example 3 # ----------------------- ex6_3.oxygen.lst <-list(n=20, p=5, gamma=c(0,0,1,0,1), BOD = c(1125, 920, 835, 1000, 1150, 990, 840, 650, 640, 583, 570, 570, 510, 555, 460, 275, 510, 165, 244, 79), TKN = c(232, 268, 271, 237, 192, 202, 184, 200, 180, 165, 151, 171, 243, 147, 286, 198, 196, 210, 327, 334), TS = c(7160, 8804, 8108, 6370, 6441, 5154, 5896, 5336, 5041, 5012, 4825, 4391, 4320, 3709, 3969, 3558, 4361, 3301, 2964, 2777), TVS = c(85.9, 86.5, 85.2, 83.8, 82.1, 79.2, 81.2, 80.6, 78.4, 79.3, 78.7, 78, 72.3, 74.9, 74.4, 72.5, 57.7, 71.8, 72.5, 71.9), COD = c(8905, 7388, 5348, 8056, 6960, 5690, 6932, 5400, 3177, 4461, 3901, 5002, 4665, 4642, 4840, 4479, 4200, 3410, 3360, 2599), O2UP = c(36, 7.9, 5.6, 5.2, 2, 2.3, 1.3, 1.3, 0.6, 0.7, 1, 1, 0.8, 0.6, 0.4, 0.7, 0.6, 0.4, 0.3, 0.9)) ex6_3.oxygen <- data.frame( BOD=ex6_3.oxygen.lst[[4]], TKN=ex6_3.oxygen.lst[[5]], TS =ex6_3.oxygen.lst[[6]], TVS=ex6_3.oxygen.lst[[7]], COD=ex6_3.oxygen.lst[[8]], O2UP=ex6_3.oxygen.lst[[9]]) # ----------------------- # chapter 7 - example 1 # ----------------------- ex7_1.aircraft.lst <- list( damage = c(0, 1, 0, 0, 0, 0, 1, 0, 0, 2, 1, 1, 1, 1, 2, 3, 1, 1, 1, 2, 0, 1, 1, 2, 5, 1, 1, 5, 5, 7), type = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), bombload = c(4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 7, 7, 7, 10, 10, 10, 12, 12, 12, 8, 8, 8, 14, 14, 14), airexp = c(91.5, 84, 76.5, 69, 61.5, 80, 72.5, 65, 57.5, 50, 103, 95.5, 88, 80.5, 73, 116.1, 100.6, 85, 69.4, 53.9, 112.3, 96.7, 81.1, 65.6, 50, 120, 104.4, 88.9, 73.7, 57.8) ) ex7_1.aircraft<-as.data.frame(ex7_1.aircraft.lst) # ----------------------- # chapter 7 - example 2 # ----------------------- ex7_2_soccer<-structure(list(team1 = structure(c(16, 1, 7, 13, 14, 15, 19, 4, 11, 6, 18, 17, 2, 10, 3, 5, 8, 12, 9, 5, 8, 17, 18, 20, 10, 2, 3, 12, 7, 1, 4, 6, 13, 14, 16, 11, 19, 15, 5, 4, 7, 16, 18, 6, 3, 17, 19, 11, 9, 9, 1, 2, 8, 10, 12, 20, 15, 13, 14, 4, 5, 6, 7, 16, 19, 3, 11, 17, 18, 20, 1, 2, 9, 10, 12, 14, 15, 13, 8, 20, 5, 6, 7, 2, 11, 3, 12, 17, 15, 16, 1, 4, 8, 9, 14, 18, 13, 19, 10, 8, 4, 5, 9, 11, 18, 13, 19, 2, 17, 10, 6, 7, 12, 14, 16, 20, 3, 15, 1, 10, 1, 6, 7, 14, 15, 16, 12, 20, 3, 5, 2, 8, 9, 19, 4, 13, 17, 11, 18, 2, 8, 4, 9, 11, 1, 3, 14, 15, 16, 20, 12, 7, 10, 5, 17, 13, 19, 11, 3, 9, 12, 14, 17, 18, 4, 6, 16, 6, 20, 5, 1, 13, 15, 20, 2, 7, 10, 19, 8, 8, 1, 2, 9, 10, 12, 13, 14, 15, 20, 6, 17, 19, 3, 4, 7, 11, 16, 18, 5, 5, 3, 4, 6, 7, 11, 17, 19, 16, 9, 8, 10, 12, 14, 15, 20, 13, 1, 2, 18, 4, 5, 6, 11, 16, 19, 3, 7, 17, 9, 2, 8, 12, 13, 14, 15, 10, 20, 1, 18, 14, 16, 19, 15, 6, 4, 11, 13, 9, 2, 3, 5, 8, 10, 18, 20, 12, 17, 15, 6, 7, 11, 13, 16, 19, 14, 4, 1, 18, 7, 8, 5, 9, 12, 18, 20, 3, 17, 9, 1, 8, 10, 13, 16, 18, 14, 4, 19, 2, 10, 11, 6, 12, 15, 17, 20, 3, 5, 2, 7, 9, 4, 5, 8, 11, 13, 19, 18, 17, 2, 10, 7, 6, 1, 3, 12, 15, 16, 20, 14, 18, 2, 4, 8, 13, 5, 1, 10, 12, 14, 15, 16, 20, 7, 1, 11, 3, 9, 19, 17, 4, 5, 8, 9, 18, 19, 11, 13, 2, 6, 7, 3, 10, 12, 14, 16, 20, 1, 15, 10, 7, 8, 13, 15, 19, 20, 2, 1, 5, 6, 17, 3, 4, 6, 9, 11, 12, 14, 16, 17, 18), .Label = c("Arsenal ", "Aston Villa ", "Blackburn ", "Bolton ", "Charlton ", "Chelsea ", "Everton ", "Fulham ", "Liverpool ", "Man City ", "Man Utd ", "Middlesbrough", "Newcastle ", "Portsmouth ", "Reading ", "Sheff Utd ", "Tottenham ", "Watford ", "West Ham ", "Wigan "), class = "factor"), goals1 = c(1, 1, 2, 2, 3, 3, 3, 2, 5, 3, 1, 2, 2, 0, 1, 0, 1, 2, 2, 2, 1, 0, 1, 1, 1, 2, 0, 0, 3, 1, 1, 2, 1, 1, 0, 1, 1, 1, 0, 0, 2, 1, 0, 1, 4, 0, 0, 0, 2, 3, 3, 2, 0, 2, 0, 1, 1, 1, 0, 2, 1, 1, 1, 2, 0, 2, 2, 2, 3, 1, 3, 1, 1, 0, 2, 2, 0, 1, 2, 4, 0, 2, 2, 1, 2, 0, 1, 1, 0, 0, 1, 0, 0, 3, 3, 0, 0, 2, 1, 1, 0, 1, 2, 3, 2, 0, 1, 2, 2, 0, 4, 0, 1, 1, 2, 3, 0, 3, 3, 3, 1, 1, 1, 2, 2, 1, 0, 0, 1, 1, 1, 0, 1, 1, 3, 1, 3, 1, 0, 1, 2, 0, 0, 3, 3, 2, 2, 1, 2, 0, 1, 2, 0, 1, 2, 3, 0, 3, 1, 4, 1, 2, 5, 0, 4, 1, 2, 1, 0, 0, 2, 2, 1, 0, 0, 2, 1, 1, 2, 0, 6, 0, 2, 0, 2, 3, 3, 0, 2, 2, 2, 1, 1, 2, 0, 3, 0, 1, 2, 2, 2, 3, 2, 3, 3, 0, 0, 1, 3, 0, 2, 3, 1, 6, 0, 2, 4, 0, 0, 0, 1, 4, 3, 1, 3, 0, 1, 2, 2, 2, 1, 5, 2, 0, 3, 0, 0, 2, 2, 0, 2, 1, 3, 3, 1, 4, 3, 0, 1, 2, 0, 2, 0, 0, 1, 1, 0, 2, 3, 1, 2, 2, 2, 0, 2, 2, 2, 1, 1, 1, 4, 4, 2, 0, 1, 3, 4, 0, 2, 1, 0, 0, 1, 2, 0, 1, 3, 0, 0, 4, 3, 0, 0, 3, 0, 1, 2, 0, 1, 4, 1, 1, 1, 4, 0, 2, 0, 1, 1, 0, 4, 1, 0, 1, 4, 1, 1, 1, 2, 4, 1, 1, 1, 0, 0, 2, 0, 1, 2, 1, 3, 3, 2, 3, 2, 3, 2, 1, 2, 1, 1, 1, 2, 1, 1, 1, 0, 0, 2, 2, 4, 0, 2, 2, 1, 0, 3, 1, 0, 3, 1, 0, 0, 3, 0, 3, 1, 0, 0, 1, 3, 2, 1, 2, 0, 3, 0, 1, 2, 1), goals2 = c(1, 1, 1, 1, 0, 2, 1, 0, 1, 0, 1, 0, 1, 0, 1, 3, 1, 1, 1, 0, 0, 2, 2, 0, 0, 0, 2, 4, 0, 1, 0, 1, 2, 0, 0, 0, 1, 0, 1, 0, 2, 2, 0, 0, 2, 0, 2, 1, 0, 0, 0, 0, 2, 0, 1, 1, 1, 1, 1, 0, 2, 1, 1, 1, 1, 1, 0, 1, 3, 3, 0, 1, 1, 0, 1, 0, 1, 2, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 4, 2, 1, 4, 1, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 2, 2, 1, 1, 0, 1, 1, 0, 0, 1, 0, 2, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 3, 1, 1, 0, 0, 0, 0, 2, 0, 1, 4, 2, 0, 0, 0, 1, 2, 2, 1, 3, 0, 1, 0, 1, 0, 0, 1, 2, 0, 1, 3, 2, 1, 2, 1, 1, 3, 2, 0, 1, 0, 2, 3, 0, 2, 0, 1, 1, 2, 3, 2, 1, 2, 0, 1, 0, 1, 1, 2, 2, 1, 1, 2, 2, 0, 2, 1, 1, 0, 0, 0, 1, 1, 1, 0, 3, 2, 0, 0, 3, 0, 3, 0, 1, 1, 3, 2, 1, 3, 0, 0, 1, 1, 2, 1, 1, 3, 2, 1, 1, 0, 0, 2, 2, 0, 1, 0, 1, 0, 0, 1, 1, 1, 2, 1, 0, 1, 4, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 2, 2, 0, 0, 1, 3, 0, 0, 1, 1, 1, 1, 1, 0, 1, 2, 2, 2, 4, 1, 1, 1, 0, 2, 0, 1, 0, 2, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 2, 1, 2, 2, 3, 1, 2, 1, 1, 3, 0, 0, 1, 0, 3, 1, 0, 0, 3, 1, 1, 0, 1, 0, 4, 2, 3, 1, 1, 0, 1, 0, 1, 0, 0, 2, 4, 1, 2, 3, 1, 0, 3, 1, 0, 1, 0, 0, 2, 2, 1, 1, 0, 1, 2, 0, 1, 3, 2, 1, 2, 1, 1, 0, 2, 1, 1), team2 = structure(c(9, 2, 18, 20, 3, 12, 5, 17, 8, 10, 19, 16, 15, 14, 7, 11, 4, 6, 19, 4, 16, 7, 11, 15, 1, 13, 6, 14, 9, 12, 18, 5, 8, 20, 3, 17, 2, 10, 14, 12, 20, 15, 2, 9, 10, 8, 13, 1, 13, 17, 16, 5, 6, 19, 3, 18, 11, 7, 4, 9, 1, 2, 10, 12, 15, 20, 13, 14, 8, 11, 18, 17, 3, 16, 7, 19, 6, 4, 5, 10, 18, 14, 16, 8, 9, 4, 13, 19, 1, 6, 7, 11, 20, 2, 15, 17, 5, 3, 12, 7, 20, 10, 15, 14, 12, 16, 1, 3, 6, 13, 18, 2, 19, 8, 4, 5, 11, 17, 9, 8, 13, 19, 4, 18, 5, 11, 9, 2, 17, 7, 12, 15, 10, 16, 1, 14, 20, 6, 16, 10, 1, 6, 14, 7, 17, 8, 2, 4, 5, 9, 11, 19, 18, 3, 12, 15, 20, 10, 13, 8, 20, 7, 5, 15, 19, 1, 2, 13, 1, 9, 14, 18, 3, 16, 4, 6, 17, 11, 12, 19, 3, 11, 18, 4, 5, 17, 16, 7, 6, 15, 2, 14, 9, 13, 12, 20, 10, 1, 8, 2, 12, 14, 8, 13, 15, 9, 10, 1, 4, 18, 7, 16, 17, 19, 3, 11, 5, 6, 9, 10, 12, 20, 2, 14, 8, 1, 15, 13, 6, 18, 17, 4, 19, 5, 16, 3, 7, 11, 3, 12, 8, 9, 20, 3, 5, 18, 2, 7, 19, 16, 6, 13, 15, 4, 14, 1, 11, 2, 12, 3, 5, 9, 17, 18, 10, 8, 20, 20, 17, 11, 19, 16, 15, 7, 13, 14, 4, 11, 15, 2, 20, 12, 7, 5, 6, 3, 17, 1, 6, 4, 16, 10, 14, 18, 8, 19, 13, 9, 1, 1, 16, 20, 14, 3, 10, 12, 6, 15, 7, 5, 8, 17, 19, 2, 18, 9, 13, 4, 11, 14, 20, 7, 10, 1, 15, 4, 9, 2, 13, 8, 19, 17, 5, 10, 16, 18, 12, 6, 1, 15, 16, 3, 20, 10, 7, 12, 6, 14, 4, 11, 5, 2, 17, 9, 18, 19, 8, 13, 11, 14, 9, 3, 18, 4, 12, 16, 6, 17, 11, 3, 15, 2, 7, 5, 19, 8, 1, 20, 10, 13), .Label = c("Arsenal ", "Aston Villa ", "Blackburn ", "Bolton ", "Charlton ", "Chelsea ", "Everton ", "Fulham ", "Liverpool ", "Man City ", "Man Utd ", "Middlesbrough", "Newcastle ", "Portsmouth ", "Reading ", "Sheff Utd ", "Tottenham ", "Watford ", "West Ham ", "Wigan "), class = "factor"), ht = structure(c(16, 1, 7, 13, 14, 15, 19, 4, 11, 6, 18, 17, 2, 10, 3, 5, 8, 12, 9, 5, 8, 17, 18, 20, 10, 2, 3, 12, 7, 1, 4, 6, 13, 14, 16, 11, 19, 15, 5, 4, 7, 16, 18, 6, 3, 17, 19, 11, 9, 9, 1, 2, 8, 10, 12, 20, 15, 13, 14, 4, 5, 6, 7, 16, 19, 3, 11, 17, 18, 20, 1, 2, 9, 10, 12, 14, 15, 13, 8, 20, 5, 6, 7, 2, 11, 3, 12, 17, 15, 16, 1, 4, 8, 9, 14, 18, 13, 19, 10, 8, 4, 5, 9, 11, 18, 13, 19, 2, 17, 10, 6, 7, 12, 14, 16, 20, 3, 15, 1, 10, 1, 6, 7, 14, 15, 16, 12, 20, 3, 5, 2, 8, 9, 19, 4, 13, 17, 11, 18, 2, 8, 4, 9, 11, 1, 3, 14, 15, 16, 20, 12, 7, 10, 5, 17, 13, 19, 11, 3, 9, 12, 14, 17, 18, 4, 6, 16, 6, 20, 5, 1, 13, 15, 20, 2, 7, 10, 19, 8, 8, 1, 2, 9, 10, 12, 13, 14, 15, 20, 6, 17, 19, 3, 4, 7, 11, 16, 18, 5, 5, 3, 4, 6, 7, 11, 17, 19, 16, 9, 8, 10, 12, 14, 15, 20, 13, 1, 2, 18, 4, 5, 6, 11, 16, 19, 3, 7, 17, 9, 2, 8, 12, 13, 14, 15, 10, 20, 1, 18, 14, 16, 19, 15, 6, 4, 11, 13, 9, 2, 3, 5, 8, 10, 18, 20, 12, 17, 15, 6, 7, 11, 13, 16, 19, 14, 4, 1, 18, 7, 8, 5, 9, 12, 18, 20, 3, 17, 9, 1, 8, 10, 13, 16, 18, 14, 4, 19, 2, 10, 11, 6, 12, 15, 17, 20, 3, 5, 2, 7, 9, 4, 5, 8, 11, 13, 19, 18, 17, 2, 10, 7, 6, 1, 3, 12, 15, 16, 20, 14, 18, 2, 4, 8, 13, 5, 1, 10, 12, 14, 15, 16, 20, 7, 1, 11, 3, 9, 19, 17, 4, 5, 8, 9, 18, 19, 11, 13, 2, 6, 7, 3, 10, 12, 14, 16, 20, 1, 15, 10, 7, 8, 13, 15, 19, 20, 2, 1, 5, 6, 17, 3, 4, 6, 9, 11, 12, 14, 16, 17, 18), .Label = c("Arsenal", "Aston Villa", "Blackburn", "Bolton", "Charlton", "Chelsea", "Everton", "Fulham", "Liverpool", "Man City", "Man Utd", "Middlesbrough", "Newcastle", "Portsmouth", "Reading", "Sheff Utd", "Tottenham", "Watford", "West Ham", "Wigan"), class = "factor"), at = structure(c(9, 2, 18, 20, 3, 12, 5, 17, 8, 10, 19, 16, 15, 14, 7, 11, 4, 6, 19, 4, 16, 7, 11, 15, 1, 13, 6, 14, 9, 12, 18, 5, 8, 20, 3, 17, 2, 10, 14, 12, 20, 15, 2, 9, 10, 8, 13, 1, 13, 17, 16, 5, 6, 19, 3, 18, 11, 7, 4, 9, 1, 2, 10, 12, 15, 20, 13, 14, 8, 11, 18, 17, 3, 16, 7, 19, 6, 4, 5, 10, 18, 14, 16, 8, 9, 4, 13, 19, 1, 6, 7, 11, 20, 2, 15, 17, 5, 3, 12, 7, 20, 10, 15, 14, 12, 16, 1, 3, 6, 13, 18, 2, 19, 8, 4, 5, 11, 17, 9, 8, 13, 19, 4, 18, 5, 11, 9, 2, 17, 7, 12, 15, 10, 16, 1, 14, 20, 6, 16, 10, 1, 6, 14, 7, 17, 8, 2, 4, 5, 9, 11, 19, 18, 3, 12, 15, 20, 10, 13, 8, 20, 7, 5, 15, 19, 1, 2, 13, 1, 9, 14, 18, 3, 16, 4, 6, 17, 11, 12, 19, 3, 11, 18, 4, 5, 17, 16, 7, 6, 15, 2, 14, 9, 13, 12, 20, 10, 1, 8, 2, 12, 14, 8, 13, 15, 9, 10, 1, 4, 18, 7, 16, 17, 19, 3, 11, 5, 6, 9, 10, 12, 20, 2, 14, 8, 1, 15, 13, 6, 18, 17, 4, 19, 5, 16, 3, 7, 11, 3, 12, 8, 9, 20, 3, 5, 18, 2, 7, 19, 16, 6, 13, 15, 4, 14, 1, 11, 2, 12, 3, 5, 9, 17, 18, 10, 8, 20, 20, 17, 11, 19, 16, 15, 7, 13, 14, 4, 11, 15, 2, 20, 12, 7, 5, 6, 3, 17, 1, 6, 4, 16, 10, 14, 18, 8, 19, 13, 9, 1, 1, 16, 20, 14, 3, 10, 12, 6, 15, 7, 5, 8, 17, 19, 2, 18, 9, 13, 4, 11, 14, 20, 7, 10, 1, 15, 4, 9, 2, 13, 8, 19, 17, 5, 10, 16, 18, 12, 6, 1, 15, 16, 3, 20, 10, 7, 12, 6, 14, 4, 11, 5, 2, 17, 9, 18, 19, 8, 13, 11, 14, 9, 3, 18, 4, 12, 16, 6, 17, 11, 3, 15, 2, 7, 5, 19, 8, 1, 20, 10, 13), .Label = c("Arsenal", "Aston Villa", "Blackburn", "Bolton", "Charlton", "Chelsea", "Everton", "Fulham", "Liverpool", "Man City", "Man Utd", "Middlesbrough", "Newcastle", "Portsmouth", "Reading", "Sheff Utd", "Tottenham", "Watford", "West Ham", "Wigan"), class = "factor"), z = c(0, 0, 1, 1, 3, 1, 2, 2, 4, 3, 0, 2, 1, 0, 0, -3, 0, 1, 1, 2, 1, -2, -1, 1, 1, 2, -2, -4, 3, 0, 1, 1, -1, 1, 0, 1, 0, 1, -1, 0, 0, -1, 0, 1, 2, 0, -2, -1, 2, 3, 3, 2, -2, 2, -1, 0, 0, 0, -1, 2, -1, 0, 0, 1, -1, 1, 2, 1, 0, -2, 3, 0, 0, 0, 1, 2, -1, -1, 1, 4, 0, 1, 2, 0, 2, -1, 1, 1, -4, -2, 0, -4, -1, 2, 2, 0, 0, 1, 1, 1, -1, 1, 2, 3, 2, -1, 1, 2, 1, 0, 4, -1, 1, 0, 0, 1, -1, 2, 3, 2, 0, 1, 1, 1, 2, -1, 0, 0, 0, 0, 0, -1, 1, 1, 2, 1, 2, 0, -1, -2, 1, -1, 0, 3, 3, 2, 0, 1, 1, -4, -1, 2, 0, 1, 1, 1, -2, 2, -2, 4, 0, 2, 4, 0, 4, 0, 0, 1, -1, -3, 0, 1, -1, -1, -1, -1, -1, 1, 1, 0, 4, -3, 2, -2, 2, 2, 2, -2, -1, 0, 1, -1, 1, 1, 0, 2, -1, -1, 0, 1, 1, 1, 0, 3, 1, -1, -1, 1, 3, 0, 1, 2, 0, 6, -3, 0, 4, 0, -3, 0, -2, 4, 2, 0, 0, -2, 0, -1, 2, 2, 0, 4, 0, -1, 2, -3, -2, 1, 1, 0, 2, -1, 1, 3, 0, 4, 2, 0, 1, 1, -1, 1, -2, -1, 1, 0, -4, 2, 3, 1, 2, 1, 1, -1, 1, 1, 1, 0, -1, -1, 4, 4, 1, -3, 1, 3, 3, -1, 1, 0, -1, 0, 0, 0, -2, -1, -1, -1, -1, 3, 3, -2, 0, 2, 0, -1, 2, 0, 1, 3, 1, 1, 0, 3, -1, 2, -1, 1, 0, 0, 3, 1, -1, -1, 3, -1, -1, -2, 1, 2, 0, 0, -2, 0, 0, 1, 0, -2, 1, 1, 3, 0, 1, 2, 2, 2, 2, -3, 0, -2, 0, 0, 2, 0, 1, 0, 0, 0, 0, -2, 3, -2, -1, 1, 1, -3, 2, 1, -1, 3, 1, -2, -2, 2, -1, 3, 0, -2, 0, 0, 0, 0, 0, 0, -1, 2, 0, -1, 1, 0)), .Names = c("team1", "goals1", "goals2", "team2", "ht", "at", "z"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96", "97", "98", "99", "100", "101", "102", "103", "104", "105", "106", "107", "108", "109", "110", "111", "112", "113", "114", "115", "116", "117", "118", "119", "120", "121", "122", "123", "124", "125", "126", "127", "128", "129", "130", "131", "132", "133", "134", "135", "136", "137", "138", "139", "140", "141", "142", "143", "144", "145", "146", "147", "148", "149", "150", "151", "152", "153", "154", "155", "156", "157", "158", "159", "160", "161", "162", "163", "164", "165", "166", "167", "168", "169", "170", "171", "172", "173", "174", "175", "176", "177", "178", "179", "180", "181", "182", "183", "184", "185", "186", "187", "188", "189", "190", "191", "192", "193", "194", "195", "196", "197", "198", "199", "200", "201", "202", "203", "204", "205", "206", "207", "208", "209", "210", "211", "212", "213", "214", "215", "216", "217", "218", "219", "220", "221", "222", "223", "224", "225", "226", "227", "228", "229", "230", "231", "232", "233", "234", "235", "236", "237", "238", "239", "240", "241", "242", "243", "244", "245", "246", "247", "248", "249", "250", "251", "252", "253", "254", "255", "256", "257", "258", "259", "260", "261", "262", "263", "264", "265", "266", "267", "268", "269", "270", "271", "272", "273", "274", "275", "276", "277", "278", "279", "280", "281", "282", "283", "284", "285", "286", "287", "288", "289", "290", "291", "292", "293", "294", "295", "296", "297", "298", "299", "300", "301", "302", "303", "304", "305", "306", "307", "308", "309", "310", "311", "312", "313", "314", "315", "316", "317", "318", "319", "320", "321", "322", "323", "324", "325", "326", "327", "328", "329", "330", "331", "332", "333", "334", "335", "336", "337", "338", "339", "340", "341", "342", "343", "344", "345", "346", "347", "348", "349", "350", "351", "352", "353", "354", "355", "356", "357", "358", "359", "360", "361", "362", "363", "364", "365", "366", "367", "368", "369", "370", "371", "372", "373", "374", "375", "376", "377", "378", "379", "380"), class = "data.frame") ex7_2_soccer.lst<-list(n=380, K=20, ht=as.numeric(ex7_2_soccer$team1), at=as.numeric(ex7_2_soccer$team2), goals1=ex7_2_soccer$goals1, goals2=ex7_2_soccer$goals2) # ----------------------- # chapter 7 - example 3 # ----------------------- ex7_3.wais.lst <- list( n=54, wais = c(9, 13, 6, 8, 10, 4, 14, 8, 11, 7, 9, 7, 5, 14, 13, 16, 10, 12, 11, 14, 15, 18, 7, 16, 9, 9, 11, 13, 15, 13, 10, 11, 6, 17, 14, 19, 9, 11, 14, 10, 16, 10, 16, 14, 13, 13, 9, 15, 10, 11, 12, 4, 14, 20), senility = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) ) ex7_3.wais <-'see ex2_3.wais' # ----------------------- # chapter 8 - example 1 # ----------------------- ex8_1.igaussian.lst <-list( n=100, y = c(41.273, 72.628, 3.895, 56.98, 5.552, 5.744, 6.399, 10.324, 19.688, 14.923, 70.551, 36.876, 73.614, 0.63, 6.016, 14.83, 5.941, 8.226, 32.997, 5.993, 11.129, 6.265, 1.09, 34.514, 22.988, 13.241, 1.984, 1.587, 3.534, 125.5, 85.96, 15.998, 305.689, 68.492, 72.521, 7.445, 2.273, 29.444, 11.045, 3.131, 5.312, 33.617, 23.875, 41.723, 52.23, 5.472, 1.327, 3.533, 2.698, 34.683, 3.868, 20.195, 10.076, 8.767, 9.262, 28.838, 8.238, 8.719, 22.366, 15.731, 6.479, 5.372, 1150.242, 2.373, 5.696, 27.043, 10.092, 9.669, 21.733, 1.151, 11.928, 1.313, 3.558, 0.867, 8.658, 23.43, 0.361, 37.741, 6.074, 10.722, 5.246, 21.654, 4.344, 2.405, 12.9, 4.277, 4.003, 0.833, 3.549, 0.091, 12.15, 9.041, 2.137, 2.016, 8.747, 1.256, 25.758, 16.548, 9.812, 3.415), x1 = c(0.475, 1.958, -0.957, 0.442, 1.234, -0.407, 0.978, 1.733, 0.055, 0.863, 0.087, 0.836, -0.366, -1.478, 1.185, 1.058, 0.293, 0.665, 1.278, -1.578, -0.861, -0.889, -0.954, -0.248, 1.257, 2.029, -1.436, -0.484, 0.259, 1.477, 0.114, 0.137, 0.416, 0.672, 1.614, -0.246, -0.53, 0.704, -0.155, -0.738, -0.638, 1.591, -0.069, 1.692, 0.424, 0.214, -1.031, -0.393, -0.579, 0.546, -1.344, 0.105, 0.195, 1.313, -0.684, -0.003, -0.615, 0.466, -0.579, -0.927, -0.729, 0.136, 1.558, -0.136, -0.539, 0.637, -0.173, 1.673, -0.222, 0.14, 1.335, 0.884, 1.025, -1.228, 0.528, 0.678, -1.424, 1.997, -0.412, 0.406, 0.412, 2.916, -0.212, -0.711, -0.406, -0.524, 1.536, -1.537, -0.738, -2.417, 1.259, -0.548, -0.014, 0.957, 0.289, -1.214, 1.076, 0.056, -0.883, -0.069), x2 = c(0.404, -1.727, -0.556, 1.205, 0.077, 0.101, -0.334, -0.682, -1.702, 1.198, 0.304, 0.4, -1.397, 0.342, -1.669, 1.213, -3.665, -0.685, 1.116, -1.834, -1.285, -0.94, 1.201, -0.587, 0.348, 0.195, -0.986, 1.742, -1.063, 0.252, -0.706, -0.223, -0.228, -2.281, 1.148, -0.613, 0.249, 0.737, -0.966, -1.585, 0.373, 0.136, 0.115, 1.027, 0.357, 0.35, 0.588, 0.824, 1.338, -0.58, -0.788, -0.073, -2.042, -0.833, -1.021, 0.535, -0.171, 1.609, -0.446, -1.164, -0.2, 1.849, -0.203, 0.93, 0.586, -0.713, 0.744, -0.231, -1.904, -0.68, -1.232, 0.12, -0.265, 0.372, 0.233, 0.043, 0.961, -2.326, 0.506, 1.065, 0.838, -0.235, 0.28, 0.737, -0.327, -0.042, 0.512, 0.121, -0.035, 0.485, 1.534, 0.59, -0.067, -0.816, -1.2, 0.682, -0.727, -0.525, -0.186, 0.916), x3 = c(0.36, -0.975, 0.703, 1.386, 0.079, 0.389, 0.76, -2.64, -0.729, -0.973, 0.576, 2.475, -0.74, -0.347, 0.563, 1.4, -0.042, 0.748, -1.813, -2.182, 0.006, 0.149, 0.003, 0.135, -0.153, -1.331, -0.368, 0.046, 0.137, 1.868, -1.933, -0.284, 0.747, 0.678, -1.301, 0.117, -0.318, -0.58, -1.251, -1.549, -1.672, -1.433, 0.082, -0.261, -0.758, 0.039, -1.377, -1.238, -0.287, 0.171, 0.719, 0.184, 0.112, -2.449, 0.633, 0.667, -0.084, 1.321, 0.536, -0.932, -0.89, 0.125, 2.671, 0.265, -1.656, 0.291, -1.259, -0.082, -1.113, 0.645, -0.393, -0.166, -0.605, 0.228, 0.131, 0.55, 0.689, 1.1, 0.463, -0.712, -0.898, -0.737, 0.344, 2.71, 0.344, 1.332, 1.908, -1.562, 0.406, -0.505, -0.137, -0.8, -0.859, -0.629, 0.222, -0.718, 0.521, 0.197, -0.082, 0.921), x4 = c(-0.119, -1.501, -0.289, -0.002, -0.263, 1.287, -0.087, -0.743, -1.295, -0.341, -0.682, 0.798, -0.104, 0.03, 0.992, 1.023, 0.727, 2.839, -0.129, -0.063, 0.116, -0.425, 0.02, -0.276, -1.363, -0.33, -0.865, 0.679, -0.254, -0.837, 0.07, -0.13, 2.037, 0.099, 0.815, -0.549, -0.233, -0.163, 2.063, 0.88, 0.851, -1.447, -1.876, -0.44, 0.428, 1.344, -0.248, -0.466, 0.728, -0.823, -0.278, -2.204, -1.886, 0.027, 3.27, -0.868, 1.561, 0.364, 0.172, -0.185, -0.263, -1.786, -0.602, -1.981, 0.799, -1.218, -0.812, -0.22, -0.648, -1.874, -0.955, -1.107, -0.454, -0.615, 1.916, 1.688, 0.06, -0.006, -1.691, -0.22, -0.063, 0.285, -1.911, 0.173, 1.798, 0.921, -0.508, 0.809, -0.563, -0.265, -0.522, -0.411, 0.797, 0.886, 0.043, -0.639, -0.602, -0.389, 0.009, -0.109) ) ex8_1.igaussian <- data.frame( y = ex8_1.igaussian.lst$y, x1 = ex8_1.igaussian.lst$x1, x2 = ex8_1.igaussian.lst$x2, x3 = ex8_1.igaussian.lst$x3, x4 = ex8_1.igaussian.lst$x4 ) # ----------------------- # chapter 8 - example 2 # ----------------------- ex8_2.softdrinks <- 'see ex5_1.softdrinks & ex5_1.softdrinks.lst' # ----------------------- # chapter 8 - example 3 # ----------------------- # # original grouped data "ex8_3.gss1990.original" <- structure(list(y = as.integer(c(0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, 18, 20, 22, 23, 24, 25, 27, 30, 50, 60)), Male = as.integer(c(65, 11, 13, 14, 26, 13, 15, 7, 21, 2, 24, 6, 3, 0, 3, 3, 0, 0, 7, 0, 0, 1, 1, 0, 3, 1, 1)), Female = as.integer(c(128, 17, 23, 16, 19, 17, 17, 3, 15, 2, 13, 10, 3, 1, 10, 1, 1, 1, 6, 1, 1, 0, 3, 1, 1, 0, 0))), .Names = c("y", "Male", "Female" ), class = "data.frame", row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27" )) # # individual data in list format "ex8_3.gss1990.lst" <- list(n =550, y = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 12, 12, 12, 12, 12, 12, 13, 13, 13, 15, 15, 15, 16, 16, 16, 20, 20, 20, 20, 20, 20, 20, 24, 25, 30, 30, 30, 50, 60, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 16, 17, 18, 20, 20, 20, 20, 20, 20, 22, 23, 25, 25, 25, 27, 30), gender = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1) ) # # individual data in data.frame format ex8_3.gss1990 <- data.frame(y=ex8_3.gss1990.lst$y, gender=ex8_3.gss1990.lst$gender) # ----------------------- # chapter 8 - example 4 # ----------------------- "ex8_4.bivpois.lst" <- structure(list(y1 = c(2, 0, 4, 2, 16, 1, 0, 3, 20, 1, 4, 2, 2, 5, 1, 4, 7, 6, 2, 6, 8, 2, 0, 1, 1, 2, 1, 3, 2, 2, 0, 3, 3, 2, 1, 5, 3, 4, 1, 1, 2, 1, 1, 6, 10, 4, 1, 0, 4, 0, 1, 0, 1, 2, 2, 3, 3, 3, 1, 5, 3, 7, 0, 5, 2, 2, 4, 2, 4, 2, 3, 1, 3, 1, 1, 0, 3, 1, 4, 1, 6, 1, 2, 3, 1, 4, 3, 2, 3, 2, 5, 3, 4, 3, 1, 5, 5, 3, 3, 3, 2, 1, 2, 3, 1, 1, 2, 5, 4, 3, 2, 4, 2, 2, 3, 17, 4, 3, 3, 2, 4, 3, 0, 8, 1, 1, 0, 4, 5, 6, 3, 2, 5, 4, 2, 3, 4, 9, 7, 1, 3, 5, 0, 3, 0, 1, 1, 4, 5, 2, 9, 4, 3, 3, 2, 4, 7, 3, 6, 4, 1, 5, 3, 3, 3, 4, 12, 8, 2, 3, 2, 5, 9, 1, 3, 2, 1, 2, 5, 5, 3, 2, 6, 1, 1, 2, 2, 0, 1, 17, 3, 6, 2, 3, 4, 3, 1, 2, 1, 5), y2 = c(5, 0, 3, 2, 16, 10, 0, 10, 19, 1, 5, 4, 5, 1, 2, 3, 7, 5, 9, 4, 7, 0, 0, 0, 0, 4, 2, 7, 0, 1, 1, 2, 1, 4, 2, 5, 2, 4, 2, 4, 4, 5, 1, 11, 10, 7, 0, 3, 4, 2, 2, 0, 4, 1, 0, 3, 4, 4, 1, 5, 2, 6, 3, 8, 2, 7, 3, 1, 2, 1, 3, 1, 1, 2, 3, 2, 2, 1, 5, 5, 4, 1, 2, 4, 1, 4, 2, 5, 2, 2, 7, 3, 7, 1, 2, 4, 3, 4, 4, 2, 5, 4, 1, 2, 1, 1, 1, 10, 3, 2, 2, 0, 4, 1, 5, 16, 1, 4, 1, 5, 5, 5, 5, 8, 5, 3, 3, 5, 5, 5, 5, 2, 3, 4, 6, 5, 6, 6, 7, 0, 6, 8, 1, 4, 1, 1, 3, 3, 4, 14, 4, 4, 4, 2, 11, 2, 5, 1, 4, 4, 1, 5, 5, 2, 2, 7, 19, 4, 4, 3, 6, 3, 10, 6, 7, 4, 3, 3, 4, 3, 2, 5, 10, 2, 0, 1, 2, 12, 2, 27, 5, 4, 3, 1, 3, 1, 2, 4, 8, 6), x1 = c(1.23, -0.16, -0.33, 0.97, 0.08, -0.05, 0.27, 0.32, 0.05, 0.22, 0.09, 0.17, -0.44, -1.26, 0.11, -0.36, 0.65, 0.5, -0.42, 0.14, -0.82, -0.43, 0.09, 0.73, 0.17, 1.05, -1.49, -0.83, 1.92, -0.33, -1.44, 0.06, 0.09, 0.19, -0.05, 0.17, -0.25, -0.58, 0.59, -0.12, 1.1, 1, -1.93, -0.12, -2.06, -0.19, -1.45, 1.72, -0.74, 0.52, -0.01, -0.03, -0.03, 1.26, 0.25, 1.39, 0.62, -0.56, 0.13, -0.85, 0.16, -1.29, -1.08, 1.24, -0.51, 0.28, 0.17, -1.01, 0.15, 1.38, 0.1, -0.05, -0.64, -1.5, -1.64, -0.51, 1.55, 1, -0.43, -0.88, 1.54, 0.55, -0.1, 0.49, -1.2, -0.39, 0.91, -1.21, -0.16, -0.76, -1.31, 1.07, 0.21, 0.69, -0.23, -0.72, 1.63, -1.98, 1.21, -0.94, -2.89, 0.28, -0.07, 1.88, 1.52, 0.09, -0.47, 0.54, -0.01, 0.5, -0.53, 0.01, -0.3, 1.45, -0.25, -1.18, 1.59, -0.25, 0.13, 1.31, 1.84, -0.97, -0.56, 1.49, 0.26, -0.7, 1.38, -0.4, -0.87, 0.9, 0.18, 0.63, -1.05, 1.6, -1.74, 1.46, 0.54, 0.82, -1.2, -1.74, 0.5, 0.36, -1.15, 0.05, 2.14, -1.33, 0.65, 0.01, -0.24, 1.83, 2.39, 0.43, 0.22, -0.97, 0.01, 1.49, 0.16, -0.77, -0.89, 0.09, -0.11, 0.12, -0.03, -0.91, 0.81, -0.39, -2.67, 0.37, 0.06, -2.27, 1.5, 0.47, 0.11, 0.01, 1.28, -1.17, -0.79, 0.01, 0.59, -0.81, 1.23, 0.66, -0.06, 0.06, -0.15, 0.52, -0.99, 0.06, 1.06, -1.37, 0.19, -1.27, 1.37, -1.09, 2.17, -0.24, -0.71, 0.55, -1.78, -0.24 ), x2 = c(-0.32, -0.03, -0.1, -0.39, 0.2, -0.98, 0.96, 0.08, 0.06, 0.67, 0.73, -1.64, 0.07, -1.43, -0.65, 1.68, -1.03, -0.33, -0.94, -0.63, -0.06, 0.16, 1.34, -0.97, -0.37, 0.26, -0.74, -0.54, 1.51, -0.06, -0.12, 0.38, 0.02, -0.29, -2.53, 0.29, 0.27, -0.53, -0.8, -1.14, -2.27, -1.88, 0.73, -1.21, -0.05, -0.54, 1.28, -0.52, 0.93, -1.48, 0.08, -0.02, -0.34, 0.53, 0.54, 0.75, -0.93, -0.53, 0.15, -0.52, 0.13, -0.5, 0.56, -1.12, -2.33, -1.11, -0.06, 1.62, 0.29, 0.82, -1.49, 1.17, -1.41, 0.52, 1.65, 0.84, -1.1, -0.46, 0.41, -0.04, 1.45, -0.46, -1.66, -1.45, 0.25, 0.96, -0.64, -0.47, 1.15, -0.54, -1.98, -0.16, -0.23, -0.42, 0.92, -0.09, 0.51, -0.22, -0.47, 1.77, 1.16, 1.03, -0.02, -1, 0.62, 0.94, 1.47, -1.53, -0.91, -1.3, -0.51, 0.58, -1.46, 1.59, -0.44, 1.2, -0.56, -0.43, -0.29, -0.7, -1.5, 0.74, 0.54, -0.7, 0.36, 0.78, 1.28, -0.32, 0.57, 0.9, -0.55, -0.57, 1.1, -0.71, -1.13, -1.88, 0.28, 0.45, 0.15, 0.24, -0.31, 0.78, 0.34, 0.73, -0.95, 1.06, 1.51, -2.15, 2.46, -2.73, -1.28, 1.31, -0.73, 0.06, -0.28, -0.35, -0.16, 0.81, 1.04, 0.99, 2.3, 0.32, -0.18, 1.36, -0.62, 0.78, -1.54, -0.68, -0.48, -1.15, -2.15, -0.86, 1.38, -1.28, -1.39, -0.46, -0.68, -0.44, 2.07, -1.37, -1.01, -1.53, 0.56, 0.49, -0.68, 0.71, -2.25, -0.54, 0.13, -1.55, 0.56, 1.05, -0.53, 1.62, -1.45, 2.23, 0.69, -0.44, -0.58, 0.53), x3 = c(1.02, -0.86, -0.95, 0.06, 2.14, 0.66, -0.62, 1.12, 1.48, 0.53, 0.85, -0.11, -0.38, 0.44, 1.79, 1.79, -0.07, -0.32, 0.46, 0.09, -0.16, 0.7, 1.35, -0.87, -0.85, 1.94, -0.56, 0.03, -0.87, -0.75, -0.07, 1.58, -1.63, 1.16, 0.84, -0.82, 1.35, 1.26, -0.71, 2.37, -0.27, 1.05, -1.55, 0.5, -0.91, -0.22, 1.1, 0.09, -0.56, -1.13, -0.88, -0.3, 1.04, -0.5, 1.16, 0.53, 0.36, 0.45, -0.22, 0.44, 0.46, 1.01, -0.35, 1.26, 0.89, -0.18, 2.05, -1.55, 1.1, -0.13, 0.04, -1.65, -0.84, -0.76, 1.28, 0.12, -0.28, 0.05, -0.49, -2.08, 0.3, -0.92, -1.18, 1.05, -1.37, -0.89, 0.99, -1.18, -0.26, -0.52, 0.26, -0.59, 0.86, 1.82, -0.43, 1.17, 0.43, 1.39, -0.63, 0.11, -0.48, 1.7, -2.43, 0.47, 0.36, -0.87, -1.3, -0.56, -1.18, -0.34, -0.04, 1.12, -0.26, 0.39, -0.78, -0.42, 0.88, 1.21, 1.37, -1.36, 0.12, -0.59, 1.35, -1.98, -0.65, 0.41, 1.02, 0.85, -0.83, -0.49, 0.47, 0.63, 0.59, -0.27, 0.13, 1.29, -0.08, 1.03, -0.86, 0.72, 2.25, -0.4, 0.19, -1.15, -0.91, 0.48, 0.89, 0.27, 0.47, 1.04, -1.53, -1.13, 1.8, -0.87, -1.68, -0.43, 0.29, 0.53, -0.15, 0.76, -1.51, -0.43, 0.63, 0.09, 0.83, -0.01, -1.58, 1.3, -1.07, 0.11, -0.63, 0.39, -1.14, 0.2, 1.16, 0.42, -0.41, -0.08, 1.04, -0.16, 1.09, 1.16, -0.05, 1.02, -1.48, -0.53, 1.42, 0.8, 0.49, -1.12, 0.02, 0.97, 1.33, 0.13, 0.07, 0.16, 0.95, -0.1, -0.86, -1.01), x4 = c(-0.45, -0.26, 1.26, -2.08, 2.44, -0.72, 0.09, -0.19, 3, -2.4, 0.75, 0.43, -1.12, -0.76, -0.06, 1.26, 1.71, -0.39, 1.51, 0.15, 0.78, -0.04, 0.15, -0.76, -2.13, -0.65, -0.44, -3.27, -0.15, 0.65, -2.29, -0.85, -0.34, 1, 0.91, 0.37, 1.12, -0.6, -1.36, -1.26, -0.7, -0.08, 0.56, 1.62, 1.67, 0.86, 0.46, -0.18, 0.2, 0.66, -0.78, -0.77, 0.07, -1.13, -0.89, 0.03, 0.56, 1.07, 0.97, 1.5, -0.23, 0.64, -1.27, 0.62, -1.13, -0.48, 0.64, -1.29, 0.97, -1.14, 0.24, -0.44, 0.16, 0.29, 0.05, -0.3, -0.22, 0.74, 0.67, 0.69, 1.33, -0.13, -0.26, 0.12, 0.6, 0.66, -0.18, -0.29, 0.31, -0.19, 1.15, 0.5, 0.87, -0.61, -1.56, -1.34, 0.92, 0.77, -0.19, -0.37, 0.64, -0.65, -1.31, -0.9, -0.08, -0.59, -0.54, 0.62, 0.11, -0.53, 0.68, -2.56, -0.18, 0.01, 1.23, 2.86, -0.51, -0.23, -0.76, -0.09, -0.36, 1.16, -0.25, 1.35, 0.3, 0.14, -0.06, 0.69, 0.8, 0.55, -0.28, -0.77, 1.53, -0.53, -1.2, -0.85, 1.01, 1.17, 0.57, 0.35, -0.03, 1.52, 0.23, -0.95, -2.03, -0.41, -0.8, -1.37, 0.8, 0.31, 0.27, 0.65, -0.52, 0.44, 0.36, 0.48, 1.2, -0.04, 0.85, 1.62, -0.16, -0.15, 0.44, -0.42, 0.48, 0.38, 1.79, 1.08, -0.29, 1.24, 1.09, -0.68, 2.25, -0.34, -0.68, 0.3, -0.34, -0.5, 1.56, -1.44, -2.06, -0.47, 0.49, -1.72, -0.75, -1.02, 0.73, -1.14, -0.68, 3.13, 0.39, 1.57, 1.3, -0.73, 0.93, -1.11, 0.33, -0.62, -1.25, 0.56), x5 = c(-0.58, 0.96, -0.04, -0.78, 0.8, -1.39, 0.28, -1.53, 1.27, 0.25, -0.89, -0.12, -0.27, 0.19, -0.51, -0.29, 0.65, -1.2, -1.36, 0.8, 0.48, 0.15, 0.35, 1.69, 2.11, -0.76, -0.22, 0.31, 2.12, 0.31, 0.43, 0.8, 0.28, -1.64, 0.09, 0.57, 1.75, -0.22, 0.8, -0.41, 2.65, -0.93, -0.13, -0.64, -0.54, -0.49, 1.14, -1.1, -0.91, 1.58, 0.52, 0.06, -1.32, 0.13, 1.44, -0.37, -0.89, 0.52, -0.41, 0.59, -0.88, -0.41, -0.63, 1.01, 0.04, -2.15, 0.04, 1.27, 1.35, 0.38, 2.87, 0.36, 1.58, 0.81, -0.66, -1.41, -0.05, 1.94, -0.86, -3.16, 1.26, -0.77, 1.21, -0.07, 0.4, -1.36, 0.78, -0.68, 0.23, -0.63, -0.95, -0.96, -1.26, 1.91, 0.45, -1.08, -0.84, -0.56, -0.96, 0.41, 0.01, -0.69, -1.53, 2.5, -0.99, 0.07, -0.41, -0.19, 1.27, -0.9, 0.46, 0.08, -0.86, -0.2, 1.38, 0.04, 0.74, -0.74, 1.33, -1.79, -0.45, 0.85, 0.02, 0.35, -1.34, -0.47, -0.72, 0.18, -0.75, -0.29, -1.76, -0.04, 0.15, -0.26, -0.74, -0.42, -0.15, 0.5, 0.99, 0.88, -0.74, 0.22, 1.5, -1.56, 1.23, -1.53, -1.17, 0.47, 2.12, -1.26, -0.24, -1.76, -1.59, -0.9, -3.28, -0.86, -0.15, -0.01, 0.72, 1.09, 0.64, -1.37, -0.15, -1.7, 0.09, -1.4, -0.7, 0.77, -1.54, 1.26, -0.14, 0.24, -1.24, -1.02, -0.77, -0.47, 1.47, -0.31, 0.94, -0.94, 0.32, -0.56, -0.87, 0.64, 1.21, 0.52, 1.03, -2.44, 0.07, -1.13, -0.56, 0, 0.53, -0.1, 0.11, 1.01, 0.84, 0.74, -0.97, -1.21)), .Names = c("y1", "y2", "x1", "x2", "x3", "x4", "x5")) ex8_4.bivpois<-as.data.frame(ex8_4.bivpois.lst) # ----------------------- # chapter 8 - example 5 # ----------------------- ex8_5.pois.diff<- 'see ex8_4.bivpois & ex8_4.bivpois.lst' # ----------------------- # chapter 9 - example 1 # ----------------------- ex9_1.blood <- list(n=20, K=2, y=structure(.Data=c(108, 98, 91, 94, 93, 96, 104, 99, 99, 97, 95, 98, 93, 97, 99, 96, 90, 100, 92, 95, 101, 89, 97, 97, 97, 100, 96, 95, 106, 100, 100, 98, 90, 99, 88, 98, 92, 92, 100, 101), .Dim = c(2, 20) ) ) # ----------------------- # chapter 9 - example 2 # ----------------------- ex9_2.kobes <- 'see ex4_1.kobes.fg' # ----------------------- # chapter 9 - example 3 # ----------------------- ex9_3.gss1990 <-'see see ex8_3.gss1990.original, ex8_3.gss1990 & ex8_3.gss1990.lst ' # ----------------------- # chapter 9 - example 4 # ----------------------- y <-structure(.Data=c(49, 12, 29, 29958, 0, 4, 33, 89, 171, 70186, 118, 81), .Dim = c(3, 2, 2)) y2<-array(dim=c(2,2,3)) y2[,,1] <- t(y[1,,]) y2[,,2] <- t(y[2,,]) y2[,,3] <- t(y[3,,]) ex9_4.metaanalysis<-list( K1=7, K2=3, or=c(3.89, 3.97, 3.88, 17.47, 5.35, 9.1, 12.41), L=c(0.92, 2.2, 2.47, 14.24, 2.44, 5.57, 2.94), U=c(16.3, 7.16, 6.08, 21.43, 11.74, 14.86, 3.96), Y=y2) rm(y);rm(y2) # ----------------------- # chapter 9 - example 5 # ----------------------- ex9_5.crossover.lst<-list(N=188, n=94, p=3, patient = c(1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14, 15, 15, 16, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 22, 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27, 28, 28, 29, 29, 30, 30, 31, 31, 32, 32, 33, 33, 34, 34, 35, 35, 36, 36, 37, 37, 38, 38, 39, 39, 40, 40, 41, 41, 42, 42, 43, 43, 44, 44, 45, 45, 46, 46, 47, 47, 48, 48, 49, 49, 50, 50, 51, 51, 52, 52, 53, 53, 54, 54, 55, 55, 56, 56, 57, 57, 58, 58, 59, 59, 60, 60, 61, 61, 62, 62, 63, 63, 64, 64, 65, 65, 66, 66, 67, 67, 68, 68, 69, 69, 70, 70, 71, 71, 72, 72, 73, 73, 74, 74, 75, 75, 76, 76, 77, 77, 78, 78, 79, 79, 80, 80, 81, 81, 82, 82, 83, 83, 84, 84, 85, 85, 86, 86, 87, 87, 88, 88, 89, 89, 90, 90, 91, 91, 92, 92, 93, 93, 94, 94), treatment = c(2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1), period = c(1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2), oed = c(45, 45, 48, 48, 50, 52, 47, 47, 45, NA, 60, 59, 54, 54, 51, 52, 53, 54, 49, NA, 56, 56, 52, 51, 60.5, 60, 53, 53, 56, NA, 62, 61, 58, 58, 50, 48, 47, 48, 46, NA, 44.5, 42, 56, 58, 48, 48, 38, 39, 50, NA, 66, 66, 62, 62, 58, 58, 68, 66, 62, NA, 66, 66, 68, 68, 66, 66, 68, 66, 66, NA, 63, 62, 49.5, 49, 59, 57, 59, 60, 58, NA, 65, 66, 48, 49, 63, 63, 81, 80, 83, NA, 68, 68, 78, 77, 42, NA, 50, 49, 47, NA, 54, NA, 61, 61, 54, 51, 55, 55, 52, NA, 49, 49, 46, 44, 45, 43, 51, 52, 48, NA, 51, 51, 51.5, 51, 48, 50, 54, 56, 44, NA, 45, 45, 53, 53, 54, 55, 51, 50, 47, NA, 52, 52, 53, 53, 60, 60, 55, 55, 58, NA, 55, 54, 53, 53, 52, 52, 54, 53, 58, NA, 52, 52, 56, 56, 59, 59, 58, 58, 55, NA, 48, 48, 52, 50, 59, 57, 45, 45, 56, NA, 60, 59, 57, 57, 54, 54, 72, 70), dbp = c(55, 60, 60, 65, 70, 80, 60, 60, 60, NA, 94.6666666666667, 96.6666666666667, 86.6666666666667, 86.6666666666667, 81.3333333333333, 82, 94.6666666666667, 94.6666666666667, 87.3333333333333, NA, 92, 91.3333333333333, 89.3333333333333, 86.6666666666667, 93.3333333333333, 91.3333333333333, 94, 94.6666666666667, 88.6666666666667, NA, 97.3333333333333, 94, 94, 91.3333333333333, 65.3333333333333, 64.6666666666667, 69.3333333333333, 80.6666666666667, 71.3333333333333, NA, 90, 70, 60, 68.6666666666667, 80.6666666666667, 70, 67.3333333333333, 70, 84.6666666666667, NA, 80.6666666666667, 82, 89.3333333333333, 88, 74.6666666666667, 77.3333333333333, 90, 56.6666666666667, 87.3333333333333, NA, 90, 86.6666666666667, 90.6666666666667, 90, 90.6666666666667, 90.6666666666667, 89.3333333333333, 90.6666666666667, 86, NA, 84.6666666666667, 80.6666666666667, 81.3333333333333, 87.3333333333333, 87.3333333333333, 82.6666666666667, 100.666666666667, 94.3333333333333, 91.6666666666667, NA, 91.6666666666667, 95, 101.333333333333, 102, 90.6666666666667, 90.6666666666667, 84.6666666666667, 85, 92, NA, 90.6666666666667, 90, 82, 82.6666666666667, 85, NA, 70, 81.6666666666667, 66.6666666666667, NA, 76.6666666666667, NA, 83.3333333333333, 71.6666666666667, 80, 71.6666666666667, 80, 80, 76.6666666666667, NA, 81.6666666666667, 68.3333333333333, 85, 91.6666666666667, 90, 77.6666666666667, 83.3333333333333, 71.6666666666667, 62, NA, 73.3333333333333, 74, 77.3333333333333, 87, 86.6666666666667, 83.3333333333333, 73.3333333333333, 80, 92, NA, 93.6666666666667, 90.6666666666667, 98.6666666666667, 96, 93.3333333333333, 88.6666666666667, 99.3333333333333, 97.3333333333333, 96.6666666666667, NA, 85.3333333333333, 85.3333333333333, 85.3333333333333, 79.3333333333333, 85.3333333333333, 82.6666666666667, 88, 84.6666666666667, 90, NA, 82.3333333333333, 81.3333333333333, 94.3333333333333, 96, 80.6666666666667, 79.3333333333333, 83.3333333333333, 90, 90.6666666666667, NA, 88, 87.3333333333333, 96, 96.6666666666667, 90.6666666666667, 90.6666666666667, 94.6666666666667, 91.3333333333333, 92.6666666666667, NA, 79.3333333333333, 80, 80, 80, 86.6666666666667, 84.6666666666667, 84, 82.6666666666667, 85.3333333333333, NA, 80.6666666666667, 80.6666666666667, 76.6666666666667, 77.3333333333333, 80, 80.6666666666667, 85.3333333333333, 84.6666666666667) ) ex9_5.crossover <- data.frame( patient = ex9_5.crossover.lst$patient, treatment= ex9_5.crossover.lst$treatment, period = ex9_5.crossover.lst$period, oed = ex9_5.crossover.lst$oed, dbp = ex9_5.crossover.lst$dbp) # ----------------------- # chapter 9 - example 6 # ----------------------- ex9_6.matched.winbugs.lst <- list( n=60, F=c(0,7, 9,23,60), Y1=c(0,1,0,1), Y2=c(0,0,1,1) ) ex9_6.matched <- data.frame( f=c(7,2,14,37), Y1=c(0,1,0,1), Y2=c(0,0,1,1) ) ex9_6.matched.tab <- matrix( c(7,2,14,37), nrow=2, ncol=2) # ----------------------- # chapter 9 - example 7 # ----------------------- ex9_7.spq.items<- structure(c(0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, NA, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1), .Dim = c(74, 167), .Dimnames = list(c("SPQ18.1", "SPQ19.2", "SPQ20.3", "SPQ21.4", "SPQ22.5", "SPQ23.6", "SPQ24.7", "SPQ25.8", "SPQ26.9", "SPQ27.10", "SPQ28.11", "SPQ29.12", "SPQ30.13", "SPQ31.14", "SPQ32.15", "SPQ33.16", "SPQ34.17", "SPQ35.18", "SPQ36.19", "SPQ37.20", "SPQ38.21", "SPQ39.22", "SPQ40.23", "SPQ41.24", "SPQ42.25", "SPQ43.26", "SPQ44.27", "SPQ45.28", "SPQ46.29", "SPQ47.30", "SPQ48.31", "SPQ49.32", "SPQ50.33", "SPQ51.34", "SPQ52.35", "SPQ53.36", "SPQ54.37", "SPQ55.38", "SPQ56.39", "SPQ57.40", "SPQ58.41", "SPQ59.42", "SPQ60.43", "SPQ61.44", "SPQ62.45", "SPQ63.46", "SPQ64.47", "SPQ65.48", "SPQ66.49", "SPQ67.50", "SPQ68.51", "SPQ69.52", "SPQ70.53", "SPQ71.54", "SPQ72.55", "SPQ73.56", "SPQ74.57", "SPQ75.58", "SPQ76.59", "SPQ77.60", "SPQ78.61", "SPQ79.62", "SPQ80.63", "SPQ81.64", "SPQ82.65", "SPQ83.66", "SPQ84.67", "SPQ85.68", "SPQ86.69", "SPQ87.70", "SPQ88.71", "SPQ89.72", "SPQ90.73", "SPQ91.74"), c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96", "97", "98", "99", "100", "101", "102", "103", "104", "105", "106", "107", "108", "109", "110", "111", "112", "113", "114", "115", "116", "117", "118", "119", "120", "121", "122", "123", "124", "125", "126", "127", "128", "129", "130", "131", "132", "133", "134", "135", "136", "137", "138", "139", "140", "141", "142", "143", "144", "145", "146", "147", "148", "149", "150", "151", "152", "153", "154", "155", "156", "157", "158", "159", "160", "161", "162", "163", "164", "165", "166", "167"))) ex9_7.spq.items<-t(ex9_7.spq.items) # ----------------------- # chapter 9 - example 8 # ----------------------- ex9_8.water.polo<-list(n=60, teams=16, G1 = c(6, 9, 9, 8, 6, 10, 7, 7, 6, 8, 3, 1, 10, 9, 8, 4, 13, 5, 14, 5, 5, 7, 5, 4, 6, 13, 13, 7, 8, 4, 4, 3, 5, 12, 8, 5, 7, 6, 7, 8, 10, 10, 2, 5, 16, 5, 3, 8, 6, 5, 12, 9, 2, 9, 12, 3, 9, 5, 6, 4), G2 = c(7, 8, 10, 2, 1, 4, 3, 1, 3, 6, 3, 12, 11, 7, 5, 6, 5, 4, 6, 11, 6, 14, 6, 6, 6, 6, 3, 4, 9, 10, 3, 12, 5, 9, 7, 10, 10, 5, 8, 7, 6, 3, 4, 11, 6, 6, 2, 9, 7, 9, 11, 7, 4, 8, 6, 8, 10, 8, 7, 2), team1 = c(14, 6, 15, 4, 9, 12, 8, 5, 15, 4, 16, 10, 7, 13, 14, 6, 16, 7, 13, 14, 10, 6, 15, 4, 6, 2, 9, 13, 14, 15, 16, 12, 2, 3, 16, 12, 14, 15, 9, 13, 6, 2, 14, 15, 16, 12, 9, 13, 14, 11, 7, 8, 9, 16, 14, 1, 4, 7, 9, 5), team2 = c(16, 7, 13, 10, 2, 3, 11, 1, 3, 1, 9, 5, 8, 12, 2, 11, 2, 11, 3, 9, 1, 8, 12, 5, 3, 10, 11, 1, 8, 5, 7, 4, 6, 10, 8, 5, 11, 1, 7, 4, 10, 3, 7, 4, 11, 1, 8, 5, 1, 12, 4, 15, 5, 13, 11, 12, 15, 8, 13, 16), game = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60), phase = c(1, 2, 3, 4, 1, 3, 2, 4, 3, 4, 1, 4, 2, 3, 1, 2, 1, 2, 3, 1, 4, 2, 3, 4, 7, 7, 5, 6, 5, 6, 5, 6, 7, 7, 5, 6, 5, 6, 5, 6, 7, 7, 5, 6, 5, 6, 5, 6, 14, 14, 11, 11, 8, 8, 16, 15, 13, 12, 10, 9) ) "ex9_8.more.water.polo.data" <- structure(list(eng2000 = structure(list(GOALS1 = c(17, 8, 10, 5, 15, 8, 10, 7, 18, 8, 10, 6, 6, 15, 6, 13, 18, 8, 5, 7, 5, 7, 7, 5, 5, 6, 4, 7, 10, 11, 15, 9, 1, 8, 4, 9, 9, 11, 6, 10, 4, 4, 6, 4, 9, 2, 9, 14, 14, 12, 8, 7, 10, 6, 5, 6, 5, 7, 18, 10, 11, 10, 8, 13, 10, 8, 0, 15, 11, 12, 12, 6, 6, 8, 7, 2, 10, 10, 11, 7, 10, 3, 8, 12, 6, 9, 8, 5, 16, 13), GOAL2 = c(4, 7, 10, 10, 8, 2, 10, 11, 6, 11, 8, 14, 10, 7, 11, 11, 7, 11, 9, 6, 9, 4, 6, 9, 14, 13, 11, 4, 11, 3, 3, 4, 16, 9, 7, 11, 13, 10, 5, 6, 5, 6, 6, 6, 13, 11, 13, 4, 8, 2, 7, 8, 9, 5, 11, 12, 7, 12, 2, 21, 7, 13, 15, 10, 6, 7, 0, 6, 11, 9, 6, 10, 12, 5, 6, 13, 8, 10, 4, 9, 12, 17, 9, 11, 8, 9, 5, 13, 11, 8), TEAM1 = structure(as.integer(c(8, 4, 10, 1, 8, 4, 2, 2, 7, 9, 6, 1, 10, 4, 3, 7, 6, 9, 1, 8, 5, 4, 10, 5, 3, 9, 8, 2, 1, 7, 2, 3, 1, 2, 6, 3, 1, 8, 4, 1, 7, 2, 3, 7, 10, 5, 10, 2, 5, 8, 7, 6, 3, 7, 9, 1, 5, 3, 6, 9, 4, 3, 9, 10, 3, 7, 6, 2, 6, 10, 4, 8, 5, 2, 1, 9, 8, 6, 4, 9, 5, 10, 5, 7, 4, 10, 5, 9, 8, 6)), .Label = c("Birkenhead", "Bristol", "Cheltenham", "Lancaster", "Manchester", "Penguin", "Polytechnic", "Rotherham", "Sutton", "Tyldesley"), class = "factor"), TEAM2 = structure(as.integer(c(1, 9, 5, 6, 9, 3, 7, 6, 10, 5, 3, 2, 8, 5, 7, 8, 5, 10, 4, 6, 3, 2, 1, 7, 2, 6, 4, 10, 3, 9, 9, 10, 7, 5, 4, 8, 5, 5, 10, 9, 6, 8, 9, 4, 6, 4, 7, 1, 9, 10, 3, 2, 6, 2, 4, 8, 10, 4, 1, 8, 8, 5, 7, 2, 1, 5, 9, 3, 8, 9, 1, 7, 6, 4, 10, 3, 2, 10, 7, 1, 8, 4, 2, 1, 6, 3, 1, 2, 3, 7)), .Label = c("Birkenhead", "Bristol", "Cheltenham", "Lancaster", "Manchester", "Penguin", "Polytechnic", "Rotherham", "Sutton", "Tyldesley"), class = "factor")), .Names = c("GOALS1", "GOAL2", "TEAM1", "TEAM2"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", "90"), class = "data.frame"), engbc2000 = structure(list(GOALS1 = c(6, 7, 7, 7, 6, 8, 13, 11, 11, 6, 8, 9, 8, 9, 12, 6, 11, 8), GOALS2 = c(5, 6, 4, 6, 7, 5, 5, 10, 10, 15, 9, 11, 7, 12, 14, 15, 14, 5), TEAM1 = structure(as.integer(c(2, 1, 2, 5, 1, 2, 2, 4, 3, 8, 3, 4, 1, 5, 7, 6, 7, 5)), .Label = c("Bristol", "Lancaster", "Manchester", "Penguin", "Polytechnic", "Portobello", "Rotherham", "Sutton & Cheam"), class = "factor"), TEAM2 = structure(as.integer(c(4, 5, 1, 4, 4, 5, 3, 8, 4, 2, 8, 2, 7, 6, 5, 1, 6, 1)), .Label = c("Bristol", "Lancaster", "Manchester", "Penguin", "Polytechnic", "Portobello", "Rotherham", "Sutton & Cheam"), class = "factor"), PHASE = structure(as.integer(c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3)), .Label = c("Finals", "Group A", "Group B"), class = "factor")), .Names = c("GOALS1", "GOALS2", "TEAM1", "TEAM2", "PHASE"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18"), class = "data.frame"), world2000 = structure(list(GOALS1 = c(6, 9, 9, 8, 6, 10, 7, 7, 6, 8, 3, 1, 10, 9, 8, 4, 13, 5, 14, 5, 5, 7, 5, 4, 6, 13, 13, 7, 8, 4, 4, 3, 5, 12, 8, 5, 7, 6, 7, 8, 10, 10, 2, 5, 16, 5, 3, 8, 6, 5, 12, 9, 2, 9, 12, 3, 9, 5, 6, 4), GOALS2 = c(7, 8, 10, 2, 1, 4, 3, 1, 3, 6, 3, 12, 11, 7, 5, 6, 5, 4, 6, 11, 6, 14, 6, 6, 6, 6, 3, 4, 9, 10, 3, 12, 5, 9, 7, 10, 10, 5, 8, 7, 6, 3, 4, 11, 6, 6, 2, 9, 7, 9, 11, 7, 4, 8, 6, 8, 10, 8, 7, 2), TEAM1 = structure(as.integer(c(14, 6, 15, 4, 9, 12, 8, 5, 15, 4, 16, 10, 7, 13, 14, 6, 16, 7, 13, 14, 10, 6, 15, 4, 6, 2, 9, 13, 14, 15, 16, 12, 2, 3, 16, 12, 14, 15, 9, 13, 6, 2, 14, 15, 16, 12, 9, 13, 14, 11, 7, 8, 9, 16, 14, 1, 4, 7, 9, 5)), .Label = c("AUS", "BRA", "CAN", "CRO", "ESP", "GER", "GRE", "HUN", "ITA", "JPN", "KAZ", "NED", "RUS", "SVK", "USA", "YUG"), class = "factor"), TEAM2 = structure(as.integer(c(15, 7, 13, 10, 2, 3, 11, 1, 3, 1, 9, 5, 8, 12, 2, 11, 2, 11, 3, 9, 1, 8, 12, 5, 3, 10, 11, 1, 8, 5, 7, 4, 6, 10, 8, 5, 11, 1, 7, 4, 10, 3, 7, 4, 11, 1, 8, 5, 1, 12, 4, 14, 5, 13, 11, 12, 14, 8, 13, 15)), .Label = c("AUS", "BRA", "CAN", "CRO", "ESP", "GER", "GRE", "HUN", "ITA", "JPN", "KAZ", "NED", "RUS", "USA", "YUG"), class = "factor"), GAMES = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60 ), PHASE = structure(as.integer(c(2, 3, 4, 5, 2, 4, 3, 5, 4, 5, 2, 5, 3, 4, 2, 3, 2, 3, 4, 2, 5, 3, 4, 5, 8, 8, 6, 7, 6, 7, 6, 7, 8, 8, 6, 7, 6, 7, 6, 7, 8, 8, 6, 7, 6, 7, 6, 7, 16, 16, 15, 15, 14, 14, 9, 13, 12, 11, 10, 1)), .Label = c("FINAL", "Group A", "Group B", "Group C", "Group D", "Group E (Places 13-16)", "Group F", "Group G", "POS11-12", "POS3-4", "POS5-6", "POS7-8", "POS9-10", "SF (1-4)", "SF (5-8)", "SF (9-12)"), class = "factor")), .Names = c("GOALS1", "GOALS2", "TEAM1", "TEAM2", "GAMES", "PHASE"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60"), class = "data.frame"), euro99 = structure(list(GOALS1 = c(4, 9, 7, 7, 5, 7, 8, 3, 8, 3, 11, 7, 4, 4, 8, 3, 7, 7, 7, 4, 7, 7, 3, 5, 13, 6, 8, 9, 6, 9, 8, 8, 6, 9, 6, 15, 7, 5, 7, 7, 5, 14, 6, 12), GOAL2 = c(9, 5, 10, 9, 7, 7, 3, 5, 9, 11, 8, 6, 4, 6, 7, 3, 6, 6, 8, 9, 8, 10, 8, 6, 6, 11, 10, 6, 7, 7, 7, 9, 7, 7, 7, 4, 8, 8, 10, 5, 8, 12, 7, 15), TEAM1 = structure(as.integer(c(1, 2, 6, 9, 7, 4, 2, 7, 12, 10, 11, 5, 6, 8, 9, 1, 3, 4, 1, 10, 7, 12, 2, 9, 3, 6, 11, 5, 8, 4, 10, 9, 8, 5, 11, 3, 1, 12, 2, 3, 1, 8, 2, 5)), .Label = c("Germany", "Greece", "Hungary", "Italy", "Kroatia", "Netherlands", "Romania", "Russia", "Slovakia", "Slovenia", "Spain", "Yugoslavia"), class = "factor"), TEAM2 = structure(as.integer(c(11, 10, 8, 5, 12, 3, 9, 1, 8, 3, 6, 4, 7, 11, 10, 12, 5, 2, 6, 5, 8, 11, 3, 4, 9, 12, 7, 2, 1, 10, 6, 7, 4, 12, 2, 1, 8, 11, 5, 4, 12, 11, 4, 3)), .Label = c("Germany", "Greece", "Hungary", "Italy", "Kroatia", "Netherlands", "Romania", "Russia", "Slovakia", "Slovenia", "Spain", "Yugoslavia"), class = "factor"), GAME = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44), PHASE = structure(as.integer(c(3, 2, 3, 2, 3, 2, 2, 3, 3, 2, 3, 2, 3, 3, 2, 3, 2, 2, 3, 2, 3, 3, 2, 2, 2, 3, 3, 2, 3, 2, 4, 9, 11, 11, 11, 11, 7, 7, 10, 10, 8, 6, 5, 1)), .Label = c("Final", "Group A", "Group B", "Place 11", "Place 3", "Place 5", "Place 5/8", "Place 7", "Place 9", "Semifinal", "quarter finals"), class = "factor"), DATE = c(13155609600, 13155609600, 13155609600, 13155609600, 13155609600, 13155609600, 13155696000, 13155696000, 13155696000, 13155696000, 13155696000, 13155696000, 13155782400, 13155782400, 13155782400, 13155782400, 13155782400, 13155782400, 13155868800, 13155868800, 13155868800, 13155868800, 13155868800, 13155868800, 13155955200, 13155955200, 13155955200, 13155955200, 13155955200, 13155955200, 13156128000, 13156128000, 13156128000, 13156128000, 13156128000, 13156128000, 13156214400, 13156214400, 13156214400, 13156214400, 13156300800, 13156300800, 13156387200, 13156387200 )), .Names = c("GOALS1", "GOAL2", "TEAM1", "TEAM2", "GAME", "PHASE", "DATE"), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "43", "44", "45", "46", "47", "48", "49", "50"), class = "data.frame")), .Names = c("eng2000", "engbc2000", "world2000", "euro99")) ex9_8.water.polo.wc2000<-ex9_8.more.water.polo.data[[3]] # ----------------------- # chapter 10 - example 0 # ----------------------- x1<-c(1,3,NA) x2<-c(2,3,NA) ex10_0.spq.indiv <- rbind(ex5_3.spq.indiv,x1) ex10_0.spq.indiv <- rbind(ex10_0.spq.indiv,x2) rm(x1);rm(x2) ex10_0.spq.tab.lst <- ex5_3.spq.tab.lst ex10_0.spq.tab.lst$y[3,1,2] <- NA ex10_0.spq.tab.lst$y[3,2,2] <- NA # ----------------------- # chapter 10 - example 1 # ----------------------- ex10_1.insur.lst <- list( L=7, inf=c(1.000 , 1.204 , 1.439 , 1.666 , 1.906 , 2.142 , 2.356, 2.484 ,2.618, 2.736, 2.794, 2.875, 2.981), y = structure( .Data=c( 1546.6 , 647.5 , 382.2 , 246.8 , 211.8 , 63.6 , 146.3 , 2099.0 , 1001.8 , 487.2 , 293.0 , 318.8 , 269.9 , NA , 3422.2 , 1257.1 , 488.4 , 456.0 , 562.4 , NA , NA , 4948.8 , 1899.7 , 984.0 , 1253.3 , NA , NA , NA , 8161.3 , 2820.3 , 1304.8 , NA , NA , NA , NA , 10622.0 , 3897.7 , NA , NA , NA , NA , NA , 11744.9 , NA , NA , NA , NA , NA , NA ), .Dim = c( 7,7 ) ) ) ex10_1.insur.lst$y <- t(ex10_1.insur.lst$y) ex10_1.insur <- ex10_1.insur.lst$y # ----------------------- # chapter 10 - example 2 # ----------------------- x1<-list( nf=6, n.rep=19,y=c(0,1,2,3,4,5), fx=c(2, 3,4,6,3,1) ) # scored x2<-list( nf=6, n.rep=19, y=c(0,1,2,3,4,5), fx=c(8,10,1,0,0,0) ) # conceded ex10_2.manund.goals <- list( scored=x1, conceded=x2 ) rm(x1) rm(x2) # ----------------------- # chapter 10 - example 3 # ----------------------- x1<-list( n=19, y=c(0.51, 0.1, -2.53, 1, 0.65, -0.95, 2.76, 1.33, 0.25, 1.48, -0.25, -0.45, 2.11, -0.76, -1.51, -0.35, 0.18, 1.35, -0.18) ) x2<-list( n=20, y=c(0.51, 0.1, -2.53, 1, 0.65, -0.95, 2.76, 1.33, 0.25, 1.48, -0.25, -0.45, 2.11, -0.76, -1.51, -0.35, 0.18, 1.35, -0.18, 5) ) x3<-list( n=20, y=c(0.51, 0.1, -2.53, 1, 0.65, -0.95, 2.76, 1.33, 0.25, 1.48, -0.25, -0.45, 2.11, -0.76, -1.51, -0.35, 0.18, 1.35, -0.18, 10) ) x4<-list(label='from log(gamma(2,1)', n=19, y=c(1.43, 1.16, 1.64, 1.06, -1.76, -0.29, -1.37, -0.72, -0.94, 1.2, 1.41, 0.74, 0.77, 0.71, 1.83, 1.58, 0.18, 0.18, 1.43)) x5<-list(label='from log(gamma(2,1)', n=100, y=c(2.02, 0, 0.34, 0.19, 0.03, 0.02, 0.01, 0.21, 0.04, 0, 0, 0.79, 0, 0.04, 0.01, 0.04, 0, 0.22, 0.03, 0.04, 0, 0.11, 0, 0.1, 0.31, 0.05, 0.13, 0, 0.67, 0.06, 0.08, 0.26, 0, 0, 0, 0, 0.05, 0.98, 0.06, 0, 0, 0, 0, 1.07, 0.08, 0, 0, 2.11, 0, 0.01, 0, 0, 0.01, 0, 0.01, 0.48, 0.07, 0, 0, 0.22, 0.06, 0, 0, 0, 0, 0.04, 0.49, 1.36, 0.01, 0, 0, 0, 0, 0, 0, 0, 0, 0.06, 0.01, 0, 0.66, 0.24, 0.02, 0.02, 0.05, 0, 0.15, 0.07, 2.34, 0, 0.02, 0.12, 0.38, 0, 0.02, 0.01, 0, 0.79, 0.03, 0.21)) ex10_3.simulated <- list( dataset1=x1,dataset2=x2, dataset3=x3,dataset4=x4, dataset5=x5 ) rm(x1);rm(x2);rm(x3);rm(x4);rm(x5) # ----------------------- # chapter 10 - example 4 # ----------------------- ex10_4.softdrinks <- 'see ex5_1.softdrinks & ex5_1.softdrinks.lst' # ----------------------- # chapter 11 - example 1 # ----------------------- ex11_1.kobes <- 'see ex4_1.kobes.fg' # ----------------------- # chapter 11 - example 2 # ----------------------- ex11_2.softdrinks <- 'see ex5_1.softdrinks & ex5_1.softdrinks.lst' # ----------------------- # chapter 11 - example 3 # ----------------------- ex11_3.dellaportas_etal.lst<- list(n=50, p=15, prop.mean.beta=c(0.25690,-0.131, -0.2928, -0.00222, 1.384, 1.805, 0.4622, 0.07317, 0.5033, -0.5573, -0.4004, -0.4351, -1.234, -0.2499, 0.1184, -0.6131), prop.sd.beta=c(0.46600,0.4191, 0.4735, 0.4359, 0.3905, 0.3503, 0.4525, 0.4274, 0.5179, 0.4416, 0.5128, 0.4427, 0.6002, 0.486, 0.5677, 0.3765), y = c(-2.4386, -0.7768, -2.3082, -2.364, -5.415, 2.5454, 5.6997, -1.5727, 1.7478, 2.7069, 3.3031, -1.8023, 3.4703, 0.2567, 2.672, 5.1258, -0.291, 1.4621, 4.5879, -3.6333, 1.8433, -0.0457, 1.5518, -9.2514, -4.2667, 1.7334, -4.2685, 4.8231, -0.8476, 2.0168, 3.4958, 1.9169, 3.829, -0.5453, 4.9755, -2.4626, 1.534, -2.3165, -3.4395, 0.1536, 4.4862, 5.4609, -5.2323, 1.5587, 4.3393, 1.131, 1.8364, 3.7274, -2.2956, 1.5475), x=structure( .Data=c(0.2532, 1.8156, -0.483, -0.356, -1.8555, 0.9481, -1.3358, 0.1841, -0.4689, 0.6804, -0.7008, -0.2409, -0.8453, 0.9144, 0.1809, -1.2131, 0.7096, -0.0402, 0.1953, 0.1385, 0.0058, 1.2085, 0.2564, 1.4045, 0.9116, -0.6775, -1.675, 0.6242, -0.5213, 0.8702, -0.4847, -0.1016, -0.6435, 0.1701, -0.3215, -0.2552, -2.5718, 0.2579, -1.4463, 0.1115, 0.3261, -0.798, 0.0923, 0.5127, 1.4399, -0.8514, 1.5888, 1.3791, 0.1251, -0.8537, -1.3595, -0.1091, -0.5922, -0.9882, 0.4722, -0.0096, -0.715, 0.5835, 0.0069, 0.9861, -1.273, -0.6017, -0.7611, -0.8256, -0.7499, -0.5043, 0.9635, -1.0414, -0.4521, 1.0405, 1.1622, 0.7947, -2.0972, 0.2569, 0.7769, 0.0538, -1.6878, 0.1542, 0.0936, 1.5204, 0.5577, -1.0942, -0.7661, 0.1763, -1.0316, -1.313, 0.3309, 0.2222, -0.0018, 0.3009, -0.0371, 0.5701, 0.056, 0.7502, 1.9702, 0.0177, 0.7203, 0.183, -0.5191, 0.5726, 0.2503, -0.6135, -1.0759, -0.4996, 0.5746, 0.0343, -0.2322, -0.0152, 1.2714, -1.866, -0.1385, -1.6321, -0.2823, 0.9156, -0.8727, -0.3439, 0.3202, 0.4665, -1.3111, -0.9537, -0.1755, 0.7747, -1.4362, -0.7059, -0.1968, 0.8168, 0.3267, 1.2173, -1.5339, -0.0478, -0.0735, 0.3729, 1.1634, 0.7871, 0.638, 0.8531, 0.0236, 1.1302, 0.0606, 0.1544, 0.2727, -0.8014, 1.1449, -0.3672, 0.1691, 2.1495, -0.5756, 1.2238, -0.0029, -0.7291, -0.1384, 0.4197, -1.977, 0.1509, 0.2772, -0.8394, 0.58, -1.0992, -1.0541, 0.155, -0.233, -1.0515, -1.6162, 0.8783, 1.1748, 0.8032, -0.5469, -0.2908, -0.8209, 0.8454, -0.929, 0.608, -0.0652, -1.2, -0.2657, 0.4385, 0.604, 0.3884, -1.1925, -0.8979, -1.0286, 1.8212, -0.449, 1.0174, 1.2026, 0.3204, -1.2126, -0.734, 0.0475, -1.1168, 0.2648, 0.0406, -0.3929, -0.3909, 0.9785, 0.51, 1.0901, 1.0467, 0.2824, -0.7458, -1.337, -0.8706, -0.7959, 0.4132, -0.8158, -0.7637, 0.1194, -1.1479, 0.1261, 1.5405, -1.5259, 0.1946, -1.0717, 0.5397, 1.5064, -0.3946, -0.0142, 1.2302, -0.062, 1.0966, 1.1683, 0.3913, 0.7166, 0.3662, 0.8638, -0.3847, 0.8638, -0.0125, 1.5583, 0.2215, 0.5831, -0.7839, 0.6375, 1.9644, 1.2588, -1.4622, -0.3077, -1.1078, 1.6221, -2.3519, 1.4495, 1.4793, 0.7444, 1.7181, -2.28, -0.975, 1.494, 0.3758, 1.4219, 1.6146, 0.5361, -1.4991, 0.4427, -0.9266, -2.877, 0.715, 0.3277, -0.3406, 1.9136, -0.3782, 1.7287, -1.0935, -0.17, 1.0203, 0.9198, 0.8625, -0.188, -0.8651, 1.0775, -0.2109, -0.8634, -0.2092, -0.6246, 0.3993, 1.7807, 0.7549, 1.2155, -1.9109, -0.6407, 1.53, 0.0005, -0.1075, 2.0485, -1.936, -0.8069, 0.144, -0.2597, -0.9671, 0.6101, -1.228, -0.569, 0.7296, -2.1917, -1.2761, 1.4552, -2.0649, 1.1535, 0.316, 0.909, -0.1409, 0.6253, 0.7797, 0.4458, 0.7951, 0.2297, -0.9498, 0.0302, -0.7715, -0.7871, 0.6819, 0.5402, 0.327, 0.7586, 0.2633, -2.1068, -0.1977, -0.7025, -0.9499, 0.2206, -0.8347, -0.8338, 0.4483, -0.1515, -0.0043, -0.8642, -0.0132, 0.7021, -0.8466, -0.5904, -0.1068, 1.8521, -0.7045, -0.9314, 0.2761, 1.0288, -0.614, 1.7761, 1.6665, 1.4052, 0.0477, 2.1756, 0.3955, 1.2472, -0.1659, 0.8752, 1.8678, 1.2847, 0.4967, -0.3451, -2.194, 0.0322, -1.7619, 0.0413, 0.1795, 0.0882, 0.0741, 0.2199, -0.2527, 0.7575, -0.5475, -0.2652, -0.1864, -1.8443, -0.6213, -2.2861, 0.2397, 0.7165, -0.5088, 0.0723, 0.6554, 0.1944, 0.0627, -1.2333, -0.6668, -0.3643, 2.1336, -0.2696, 0.8352, -0.0905, 0.6053, 0.1278, -0.9659, -0.3242, 1.8595, -0.401, -0.5056, 0.9797, -0.3522, 1.5781, -0.6649, -0.0661, -0.0632, -1.0288, -0.3651, 1.0089, 1.0371, -1.2665, 0.062, -0.2491, 1.3932, -0.0286, 0.6036, 1.4613, 0.1986, -0.1378, 1.5776, 0.1716, 0.2851, 2.2335, 0.7554, -0.7355, -0.7447, 0.1274, 0.4182, -0.1174, 0.7349, -0.5518, -1.7765, -0.4166, -2.0655, 1.3805, 0.5255, -0.3012, -0.0047, -1.5563, 1.6014, -0.0095, 0.8014, -0.8466, 0.1172, -0.0097, -1.1792, -0.1674, 1.3232, 0.4011, 0.698, -0.0565, 1.2923, 0.6455, 0.8473, -0.4388, -0.0258, 0.6191, -0.4862, -0.3779, -2.221, -0.0139, -1.8425, 0.7279, 0.2261, -1.7503, -0.5339, -0.8615, 3.2967, 0.1566, -0.6211, 0.7122, 0.3329, 2.0307, 0.9708, 1.0827, 0.8803, 0.4843, 0.57, -1.0104, -1.6622, -0.9716, 0.3685, 1.1403, -0.4929, -1.8076, -1.0744, 1.4807, -1.2724, 0.358, 1.4126, 0.5296, 1.0928, -0.0743, 0.1657, 0.0948, -0.3091, -1.965, -0.3689, 1.6937, 0.4195, 0.4989, -1.6845, 0.1085, -0.3983, 1.0028, -0.8749, 0.3552, -0.8551, -0.7824, 0.0405, 1.7659, -1.1563, 1.7609, -0.1989, -0.3088, 0.9441, 0.2926, -0.549, -1.4652, 1.2621, 0.9841, 0.0035, -0.1114, -0.5988, 0.6284, -0.0279, 1.307, 0.8911, -0.0647, -0.4558, -0.2504, 0.8804, -0.1044, 0.3824, -0.0161, 0.0065, 0.5124, -1.5275, 0.6063, -0.0442, -0.3649, -0.1654, -1.1256, 2.0752, -2.0174, 1.9353, 2.0064, 1.6742, 0.1361, 0.4025, 0.0397, 0.1102, 0.3306, -0.1074, -1.3069, 1.9539, 0.9281, 0.2854, 0.1434, -1.9649, -1.1645, 1.7452, 0.6721, 0.8254, -0.0543, -0.3777, 0.2849, 1.8963, 0.0158, 2.1345, -0.4826, 1.7665, -0.1036, -0.3988, 0.2894, 0.9394, -0.2685, -0.068, 0.3147, -1.0013, 0.5726, 0.6401, -1.5506, -0.179, 2.5959, -1.2389, -1.0428, 0.0766, -0.7697, -0.3437, -0.0523, 0.4095, 2.4212, -0.5301, -1.0354, 0.8516, -1.484, 0.8148, 1.4214, -0.7104, 0.9131, -0.663, -0.2181, -0.2669, 1.6292, -0.1279, -0.5301, -1.0295, 0.0432, -0.6161, 0.8508, 0.2106, 1.2939, -0.5594, -0.5563, 1.0688, -1.3793, 1.3487, 0.7632, -1.9978, -1.8987, -0.2607, -0.5655, 0.1404, -0.1473, -1.1201, -0.3694, -0.7877, -1.2406, 0.0301, -1.4742, -0.0289, 1.5659, 1.1968, -0.5718, -1.1851, 0.7766, 0.1369, -0.3001, -0.2303, -0.1515, -0.1896, 0.942, 1.1239, -0.0113, -0.6285, -1.3864, 0.2737, -2.0212, -0.0423, -0.5297, 0.7101, -0.5613, -0.3932, 1.9499, -0.691, 0.7919, -0.5318, -0.5213, -0.6517, 0.5551, 0.6835, 0.6419, 0.0681, 1.4074, 0.3109, -0.3881, 0.8534, 0.3377, 0.6069, -1.0757, 1.1765, 0.9395, 1.1052, 1.2934, -0.3334, -1.201, 1.8051, -0.4944, 0.6722, -0.3811, 0.3006, -0.3241, 0.5145, 0.8625, -0.571, 0.8386, 0.14, 1.4898, -0.2371, -1.4977, 0.6917, -2.0052, -1.0114, 0.5536, 0.1943, -0.0383, -0.1806, -0.4167, -1.636, -0.6545, 0.7911, -0.4712, -1.7515, 0.4295, 0.2707, 1.8788, 0.4132, -0.431, -0.0985, -0.6662, 0.4337, 1.5325, -1.0914, 0.8049, -0.1756, -1.2483, 1.1944, -0.9778, 0.7469, 0.3923, -1.3149, -0.2081, -0.7817, 0.2014, 0.0257, 0.2161, -0.4244, 1.6736, -0.1795, -0.0329, -0.4735, 0.9362, -1.3239, -0.4339, 0.3462, 0.4204, 2.0448, -0.112, 0.9716, 0.4501, 0.2191, 0.5259, 1.5869, 1.4232, -0.6287, 0.9264, 0.0883, 0.5609, 0.3011, 0.0659, -1.4427, -1.8872, 1.4274, 2.1341, -0.2701, 0.5767, 0.0712, 1.5841, 2.1634, -0.1232, 0.4988, -0.2374, 0.5291), .Dim = c(15,50) ) ) ex11_3.dellaportas_etal.lst$x <- t(ex11_3.dellaportas_etal.lst$x) ex11_3.dellaportas_etal <- cbind(ex11_3.dellaportas_etal.lst$y, ex11_3.dellaportas_etal.lst$x) colnames(ex11_3.dellaportas_etal) <- c('y', paste('x',1:15,sep='')) # ----------------------- # chapter 11 - example 4 # ----------------------- ex11_4.softdrinks <- 'see ex5_1.softdrinks & ex5_1.softdrinks.lst' # ----------------------- # chapter 11 - example 5 # ----------------------- ex11_5.softdrinks <- 'see ex5_1.softdrinks & ex5_1.softdrinks.lst'