R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1 + ,1 + ,41 + ,38 + ,13 + ,12 + ,14 + ,1 + ,1 + ,39 + ,32 + ,16 + ,11 + ,18 + ,1 + ,1 + ,30 + ,35 + ,19 + ,15 + ,11 + ,1 + ,0 + ,31 + ,33 + ,15 + ,6 + ,12 + ,1 + ,1 + ,34 + ,37 + ,14 + ,13 + ,16 + ,1 + ,1 + ,35 + ,29 + ,13 + ,10 + ,18 + ,1 + ,1 + ,39 + ,31 + ,19 + ,12 + ,14 + ,1 + ,1 + ,34 + ,36 + ,15 + ,14 + ,14 + ,1 + ,1 + ,36 + ,35 + ,14 + ,12 + ,15 + ,1 + ,1 + ,37 + ,38 + ,15 + ,9 + ,15 + ,1 + ,0 + ,38 + ,31 + ,16 + ,10 + ,17 + ,1 + ,1 + ,36 + ,34 + ,16 + ,12 + ,19 + ,1 + ,0 + ,38 + ,35 + ,16 + ,12 + ,10 + ,1 + ,1 + ,39 + ,38 + ,16 + ,11 + ,16 + ,1 + ,1 + ,33 + ,37 + ,17 + ,15 + ,18 + ,1 + ,0 + ,32 + ,33 + ,15 + ,12 + ,14 + ,1 + ,0 + ,36 + ,32 + ,15 + ,10 + ,14 + ,1 + ,1 + ,38 + ,38 + ,20 + ,12 + ,17 + ,1 + ,0 + ,39 + ,38 + ,18 + ,11 + ,14 + ,1 + ,1 + ,32 + ,32 + ,16 + ,12 + ,16 + ,1 + ,0 + ,32 + ,33 + ,16 + ,11 + ,18 + ,1 + ,1 + ,31 + ,31 + ,16 + ,12 + ,11 + ,1 + ,1 + ,39 + ,38 + ,19 + ,13 + ,14 + ,1 + ,1 + ,37 + ,39 + ,16 + ,11 + ,12 + ,1 + ,0 + ,39 + ,32 + ,17 + ,12 + ,17 + ,1 + ,1 + ,41 + ,32 + ,17 + ,13 + ,9 + ,1 + ,0 + ,36 + ,35 + ,16 + ,10 + ,16 + ,1 + ,1 + ,33 + ,37 + ,15 + ,14 + ,14 + ,1 + ,1 + ,33 + ,33 + ,16 + ,12 + ,15 + ,1 + ,0 + ,34 + ,33 + ,14 + ,10 + ,11 + ,1 + ,1 + ,31 + ,31 + ,15 + ,12 + ,16 + ,1 + ,0 + ,27 + ,32 + ,12 + ,8 + ,13 + ,1 + ,1 + ,37 + ,31 + ,14 + ,10 + ,17 + ,1 + ,1 + ,34 + ,37 + ,16 + ,12 + ,15 + ,1 + ,0 + ,34 + ,30 + ,14 + ,12 + ,14 + ,1 + ,0 + ,32 + ,33 + ,10 + ,7 + ,16 + ,1 + ,0 + ,29 + ,31 + ,10 + ,9 + ,9 + ,1 + ,0 + ,36 + ,33 + ,14 + ,12 + ,15 + ,1 + ,1 + ,29 + ,31 + ,16 + ,10 + ,17 + ,1 + ,0 + ,35 + ,33 + ,16 + ,10 + ,13 + ,1 + ,0 + ,37 + ,32 + ,16 + ,10 + ,15 + ,1 + ,1 + ,34 + ,33 + ,14 + ,12 + ,16 + ,1 + ,0 + ,38 + ,32 + ,20 + ,15 + ,16 + ,1 + ,0 + ,35 + ,33 + ,14 + ,10 + ,12 + ,1 + ,1 + ,38 + ,28 + ,14 + ,10 + ,15 + ,1 + ,1 + ,37 + ,35 + ,11 + ,12 + ,11 + ,1 + ,1 + ,38 + ,39 + ,14 + ,13 + ,15 + ,1 + ,1 + ,33 + ,34 + ,15 + ,11 + ,15 + ,1 + ,1 + ,36 + ,38 + ,16 + ,11 + ,17 + ,1 + ,0 + ,38 + ,32 + ,14 + ,12 + ,13 + ,1 + ,1 + ,32 + ,38 + ,16 + ,14 + ,16 + ,1 + ,0 + ,32 + ,30 + ,14 + ,10 + ,14 + ,1 + ,0 + ,32 + ,33 + ,12 + ,12 + ,11 + ,1 + ,1 + ,34 + ,38 + ,16 + ,13 + ,12 + ,1 + ,0 + ,32 + ,32 + ,9 + ,5 + ,12 + ,1 + ,1 + ,37 + ,35 + ,14 + ,6 + ,15 + ,1 + ,1 + ,39 + ,34 + ,16 + ,12 + ,16 + ,1 + ,1 + ,29 + ,34 + ,16 + ,12 + ,15 + ,1 + ,0 + ,37 + ,36 + ,15 + ,11 + ,12 + ,1 + ,1 + ,35 + ,34 + ,16 + ,10 + ,12 + ,1 + ,0 + ,30 + ,28 + ,12 + ,7 + ,8 + ,1 + ,0 + ,38 + ,34 + ,16 + ,12 + ,13 + ,1 + ,1 + ,34 + ,35 + ,16 + ,14 + ,11 + ,1 + ,1 + ,31 + ,35 + ,14 + ,11 + ,14 + ,1 + ,1 + ,34 + ,31 + ,16 + ,12 + ,15 + ,1 + ,0 + ,35 + ,37 + ,17 + ,13 + ,10 + ,1 + ,1 + ,36 + ,35 + ,18 + ,14 + ,11 + ,1 + ,0 + ,30 + ,27 + ,18 + ,11 + ,12 + ,1 + ,1 + ,39 + ,40 + ,12 + ,12 + ,15 + ,1 + ,0 + ,35 + ,37 + ,16 + ,12 + ,15 + ,1 + ,0 + ,38 + ,36 + ,10 + ,8 + ,14 + ,1 + ,1 + ,31 + ,38 + ,14 + ,11 + ,16 + ,1 + ,1 + ,34 + ,39 + ,18 + ,14 + ,15 + ,1 + ,0 + ,38 + ,41 + ,18 + ,14 + ,15 + ,1 + ,0 + ,34 + ,27 + ,16 + ,12 + ,13 + ,1 + ,1 + ,39 + ,30 + ,17 + ,9 + ,12 + ,1 + ,1 + ,37 + ,37 + ,16 + ,13 + ,17 + ,1 + ,1 + ,34 + ,31 + ,16 + ,11 + ,13 + ,1 + ,0 + ,28 + ,31 + ,13 + ,12 + ,15 + ,1 + ,0 + ,37 + ,27 + ,16 + ,12 + ,13 + ,1 + ,0 + ,33 + ,36 + ,16 + ,12 + ,15 + ,1 + ,1 + ,35 + ,37 + ,16 + ,12 + ,15 + ,1 + ,0 + ,37 + ,33 + ,15 + ,12 + ,16 + ,1 + ,1 + ,32 + ,34 + ,15 + ,11 + ,15 + ,1 + ,1 + ,33 + ,31 + ,16 + ,10 + ,14 + ,1 + ,0 + ,38 + ,39 + ,14 + ,9 + ,15 + ,1 + ,1 + ,33 + ,34 + ,16 + ,12 + ,14 + ,1 + ,1 + ,29 + ,32 + ,16 + ,12 + ,13 + ,1 + ,1 + ,33 + ,33 + ,15 + ,12 + ,7 + ,1 + ,1 + ,31 + ,36 + ,12 + ,9 + ,17 + ,1 + ,1 + ,36 + ,32 + ,17 + ,15 + ,13 + ,1 + ,1 + ,35 + ,41 + ,16 + ,12 + ,15 + ,1 + ,1 + ,32 + ,28 + ,15 + ,12 + ,14 + ,1 + ,1 + ,29 + ,30 + ,13 + ,12 + ,13 + ,1 + ,1 + ,39 + ,36 + ,16 + ,10 + ,16 + ,1 + ,1 + ,37 + ,35 + ,16 + ,13 + ,12 + ,1 + ,1 + ,35 + ,31 + ,16 + ,9 + ,14 + ,1 + ,0 + ,37 + ,34 + ,16 + ,12 + ,17 + ,1 + ,0 + ,32 + ,36 + ,14 + ,10 + ,15 + ,1 + ,1 + ,38 + ,36 + ,16 + ,14 + ,17 + ,1 + ,0 + ,37 + ,35 + ,16 + ,11 + ,12 + ,1 + ,1 + ,36 + ,37 + ,20 + ,15 + ,16 + ,1 + ,0 + ,32 + ,28 + ,15 + ,11 + ,11 + ,1 + ,1 + ,33 + ,39 + ,16 + ,11 + ,15 + ,1 + ,0 + ,40 + ,32 + ,13 + ,12 + ,9 + ,1 + ,1 + ,38 + ,35 + ,17 + ,12 + ,16 + ,1 + ,0 + ,41 + ,39 + ,16 + ,12 + ,15 + ,1 + ,0 + ,36 + ,35 + ,16 + ,11 + ,10 + ,1 + ,1 + ,43 + ,42 + ,12 + ,7 + ,10 + ,1 + ,1 + ,30 + ,34 + ,16 + ,12 + ,15 + ,1 + ,1 + ,31 + ,33 + ,16 + ,14 + ,11 + ,1 + ,1 + ,32 + ,41 + ,17 + ,11 + ,13 + ,1 + ,1 + ,37 + ,34 + ,12 + ,10 + ,18 + ,1 + ,0 + ,37 + ,32 + ,18 + ,13 + ,16 + ,1 + ,1 + ,33 + ,40 + ,14 + ,13 + ,14 + ,1 + ,1 + ,34 + ,40 + ,14 + ,8 + ,14 + ,1 + ,1 + ,33 + ,35 + ,13 + ,11 + ,14 + ,1 + ,1 + ,38 + ,36 + ,16 + ,12 + ,14 + ,1 + ,0 + ,33 + ,37 + ,13 + ,11 + ,12 + ,1 + ,1 + ,31 + ,27 + ,16 + ,13 + ,14 + ,1 + ,1 + ,38 + ,39 + ,13 + ,12 + ,15 + ,1 + ,1 + ,37 + ,38 + ,16 + ,14 + ,15 + ,1 + ,1 + ,36 + ,31 + ,15 + ,13 + ,15 + ,1 + ,1 + ,31 + ,33 + ,16 + ,15 + ,13 + ,1 + ,0 + ,39 + ,32 + ,15 + ,10 + ,17 + ,1 + ,1 + ,44 + ,39 + ,17 + ,11 + ,17 + ,1 + ,1 + ,33 + ,36 + ,15 + ,9 + ,19 + ,1 + ,1 + ,35 + ,33 + ,12 + ,11 + ,15 + ,1 + ,0 + ,32 + ,33 + ,16 + ,10 + ,13 + ,1 + ,0 + ,28 + ,32 + ,10 + ,11 + ,9 + ,1 + ,1 + ,40 + ,37 + ,16 + ,8 + ,15 + ,1 + ,0 + ,27 + ,30 + ,12 + ,11 + ,15 + ,1 + ,0 + ,37 + ,38 + ,14 + ,12 + ,15 + ,1 + ,1 + ,32 + ,29 + ,15 + ,12 + ,16 + ,1 + ,0 + ,28 + ,22 + ,13 + ,9 + ,11 + ,1 + ,0 + ,34 + ,35 + ,15 + ,11 + ,14 + ,1 + ,1 + ,30 + ,35 + ,11 + ,10 + ,11 + ,1 + ,1 + ,35 + ,34 + ,12 + ,8 + ,15 + ,1 + ,0 + ,31 + ,35 + ,11 + ,9 + ,13 + ,1 + ,1 + ,32 + ,34 + ,16 + ,8 + ,15 + ,1 + ,0 + ,30 + ,37 + ,15 + ,9 + ,16 + ,1 + ,1 + ,30 + ,35 + ,17 + ,15 + ,14 + ,1 + ,0 + ,31 + ,23 + ,16 + ,11 + ,15 + ,1 + ,1 + ,40 + ,31 + ,10 + ,8 + ,16 + ,1 + ,1 + ,32 + ,27 + ,18 + ,13 + ,16 + ,1 + ,0 + ,36 + ,36 + ,13 + ,12 + ,11 + ,1 + ,0 + ,32 + ,31 + ,16 + ,12 + ,12 + ,1 + ,0 + ,35 + ,32 + ,13 + ,9 + ,9 + ,1 + ,1 + ,38 + ,39 + ,10 + ,7 + ,16 + ,1 + ,1 + ,42 + ,37 + ,15 + ,13 + ,13 + ,1 + ,0 + ,34 + ,38 + ,16 + ,9 + ,16 + ,1 + ,1 + ,35 + ,39 + ,16 + ,6 + ,12 + ,1 + ,1 + ,38 + ,34 + ,14 + ,8 + ,9 + ,1 + ,1 + ,33 + ,31 + ,10 + ,8 + ,13 + ,1 + ,1 + ,32 + ,37 + ,13 + ,6 + ,14 + ,1 + ,1 + ,33 + ,36 + ,15 + ,9 + ,19 + ,1 + ,1 + ,34 + ,32 + ,16 + ,11 + ,13 + ,1 + ,1 + ,32 + ,38 + ,12 + ,8 + ,12 + ,0 + ,0 + ,27 + ,26 + ,13 + ,10 + ,10 + ,0 + ,0 + ,31 + ,26 + ,12 + ,8 + ,14 + ,0 + ,0 + ,38 + ,33 + ,17 + ,14 + ,16 + ,0 + ,1 + ,34 + ,39 + ,15 + ,10 + ,10 + ,0 + ,0 + ,24 + ,30 + ,10 + ,8 + ,11 + ,0 + ,0 + ,30 + ,33 + ,14 + ,11 + ,14 + ,0 + ,1 + ,26 + ,25 + ,11 + ,12 + ,12 + ,0 + ,1 + ,34 + ,38 + ,13 + ,12 + ,9 + ,0 + ,0 + ,27 + ,37 + ,16 + ,12 + ,9 + ,0 + ,0 + ,37 + ,31 + ,12 + ,5 + ,11 + ,0 + ,1 + ,36 + ,37 + ,16 + ,12 + ,16 + ,0 + ,0 + ,41 + ,35 + ,12 + ,10 + ,9 + ,0 + ,1 + ,29 + ,25 + ,9 + ,7 + ,13 + ,0 + ,1 + ,36 + ,28 + ,12 + ,12 + ,16 + ,0 + ,0 + ,32 + ,35 + ,15 + ,11 + ,13 + ,0 + ,1 + ,37 + ,33 + ,12 + ,8 + ,9 + ,0 + ,0 + ,30 + ,30 + ,12 + ,9 + ,12 + ,0 + ,1 + ,31 + ,31 + ,14 + ,10 + ,16 + ,0 + ,1 + ,38 + ,37 + ,12 + ,9 + ,11 + ,0 + ,1 + ,36 + ,36 + ,16 + ,12 + ,14 + ,0 + ,0 + ,35 + ,30 + ,11 + ,6 + ,13 + ,0 + ,0 + ,31 + ,36 + ,19 + ,15 + ,15 + ,0 + ,0 + ,38 + ,32 + ,15 + ,12 + ,14 + ,0 + ,1 + ,22 + ,28 + ,8 + ,12 + ,16 + ,0 + ,1 + ,32 + ,36 + ,16 + ,12 + ,13 + ,0 + ,0 + ,36 + ,34 + ,17 + ,11 + ,14 + ,0 + ,1 + ,39 + ,31 + ,12 + ,7 + ,15 + ,0 + ,0 + ,28 + ,28 + ,11 + ,7 + ,13 + ,0 + ,0 + ,32 + ,36 + ,11 + ,5 + ,11 + ,0 + ,1 + ,32 + ,36 + ,14 + ,12 + ,11 + ,0 + ,1 + ,38 + ,40 + ,16 + ,12 + ,14 + ,0 + ,1 + ,32 + ,33 + ,12 + ,3 + ,15 + ,0 + ,1 + ,35 + ,37 + ,16 + ,11 + ,11 + ,0 + ,1 + ,32 + ,32 + ,13 + ,10 + ,15 + ,0 + ,0 + ,37 + ,38 + ,15 + ,12 + ,12 + ,0 + ,1 + ,34 + ,31 + ,16 + ,9 + ,14 + ,0 + ,1 + ,33 + ,37 + ,16 + ,12 + ,14 + ,0 + ,0 + ,33 + ,33 + ,14 + ,9 + ,8 + ,0 + ,0 + ,30 + ,30 + ,16 + ,12 + ,9 + ,0 + ,0 + ,24 + ,30 + ,14 + ,10 + ,15 + ,0 + ,0 + ,34 + ,31 + ,11 + ,9 + ,17 + ,0 + ,0 + ,34 + ,32 + ,12 + ,12 + ,13 + ,0 + ,1 + ,33 + ,34 + ,15 + ,8 + ,15 + ,0 + ,1 + ,34 + ,36 + ,15 + ,11 + ,15 + ,0 + ,1 + ,35 + ,37 + ,16 + ,11 + ,14 + ,0 + ,0 + ,35 + ,36 + ,16 + ,12 + ,16 + ,0 + ,0 + ,36 + ,33 + ,11 + ,10 + ,13 + ,0 + ,0 + ,34 + ,33 + ,15 + ,10 + ,16 + ,0 + ,1 + ,34 + ,33 + ,12 + ,12 + ,9 + ,0 + ,0 + ,41 + ,44 + ,12 + ,12 + ,16 + ,0 + ,0 + ,32 + ,39 + ,15 + ,11 + ,11 + ,0 + ,0 + ,30 + ,32 + ,15 + ,8 + ,10 + ,0 + ,1 + ,35 + ,35 + ,16 + ,12 + ,11 + ,0 + ,0 + ,28 + ,25 + ,14 + ,10 + ,15 + ,0 + ,1 + ,33 + ,35 + ,17 + ,11 + ,17 + ,0 + ,1 + ,39 + ,34 + ,14 + ,10 + ,14 + ,0 + ,0 + ,36 + ,35 + ,13 + ,8 + ,8 + ,0 + ,1 + ,36 + ,39 + ,15 + ,12 + ,15 + ,0 + ,0 + ,35 + ,33 + ,13 + ,12 + ,11 + ,0 + ,0 + ,38 + ,36 + ,14 + ,10 + ,16 + ,0 + ,1 + ,33 + ,32 + ,15 + ,12 + ,10 + ,0 + ,0 + ,31 + ,32 + ,12 + ,9 + ,15 + ,0 + ,1 + ,32 + ,36 + ,8 + ,6 + ,16 + ,0 + ,0 + ,31 + ,32 + ,14 + ,10 + ,19 + ,0 + ,0 + ,33 + ,34 + ,14 + ,9 + ,12 + ,0 + ,0 + ,34 + ,33 + ,11 + ,9 + ,8 + ,0 + ,0 + ,34 + ,35 + ,12 + ,9 + ,11 + ,0 + ,1 + ,34 + ,30 + ,13 + ,6 + ,14 + ,0 + ,0 + ,33 + ,38 + ,10 + ,10 + ,9 + ,0 + ,0 + ,32 + ,34 + ,16 + ,6 + ,15 + ,0 + ,1 + ,41 + ,33 + ,18 + ,14 + ,13 + ,0 + ,1 + ,34 + ,32 + ,13 + ,10 + ,16 + ,0 + ,0 + ,36 + ,31 + ,11 + ,10 + ,11 + ,0 + ,0 + ,37 + ,30 + ,4 + ,6 + ,12 + ,0 + ,0 + ,36 + ,27 + ,13 + ,12 + ,13 + ,0 + ,1 + ,29 + ,31 + ,16 + ,12 + ,10 + ,0 + ,0 + ,37 + ,30 + ,10 + ,7 + ,11 + ,0 + ,0 + ,27 + ,32 + ,12 + ,8 + ,12 + ,0 + ,0 + ,35 + ,35 + ,12 + ,11 + ,8 + ,0 + ,0 + ,28 + ,28 + ,10 + ,3 + ,12 + ,0 + ,0 + ,35 + ,33 + ,13 + ,6 + ,12 + ,0 + ,0 + ,29 + ,35 + ,12 + ,8 + ,11 + ,0 + ,0 + ,32 + ,35 + ,14 + ,9 + ,13 + ,0 + ,1 + ,36 + ,32 + ,10 + ,9 + ,14 + ,0 + ,1 + ,19 + ,21 + ,12 + ,8 + ,10 + ,0 + ,1 + ,21 + ,20 + ,12 + ,9 + ,12 + ,0 + ,0 + ,31 + ,34 + ,11 + ,7 + ,15 + ,0 + ,0 + ,33 + ,32 + ,10 + ,7 + ,13 + ,0 + ,1 + ,36 + ,34 + ,12 + ,6 + ,13 + ,0 + ,1 + ,33 + ,32 + ,16 + ,9 + ,13 + ,0 + ,0 + ,37 + ,33 + ,12 + ,10 + ,12 + ,0 + ,0 + ,34 + ,33 + ,14 + ,11 + ,12 + ,0 + ,0 + ,35 + ,37 + ,16 + ,12 + ,9 + ,0 + ,1 + ,31 + ,32 + ,14 + ,8 + ,9 + ,0 + ,1 + ,37 + ,34 + ,13 + ,11 + ,15 + ,0 + ,1 + ,35 + ,30 + ,4 + ,3 + ,10 + ,0 + ,1 + ,27 + ,30 + ,15 + ,11 + ,14 + ,0 + ,0 + ,34 + ,38 + ,11 + ,12 + ,15 + ,0 + ,0 + ,40 + ,36 + ,11 + ,7 + ,7 + ,0 + ,0 + ,29 + ,32 + ,14 + ,9 + ,14 + ,0 + ,0 + ,38 + ,34 + ,15 + ,12 + ,8 + ,0 + ,1 + ,34 + ,33 + ,14 + ,8 + ,10 + ,0 + ,0 + ,21 + ,27 + ,13 + ,11 + ,13 + ,0 + ,0 + ,36 + ,32 + ,11 + ,8 + ,13 + ,0 + ,1 + ,38 + ,34 + ,15 + ,10 + ,13 + ,0 + ,0 + ,30 + ,29 + ,11 + ,8 + ,8 + ,0 + ,0 + ,35 + ,35 + ,13 + ,7 + ,12 + ,0 + ,1 + ,30 + ,27 + ,13 + ,8 + ,13 + ,0 + ,1 + ,36 + ,33 + ,16 + ,10 + ,12 + ,0 + ,0 + ,34 + ,38 + ,13 + ,8 + ,10 + ,0 + ,1 + ,35 + ,36 + ,16 + ,12 + ,13 + ,0 + ,0 + ,34 + ,33 + ,16 + ,14 + ,12 + ,0 + ,0 + ,32 + ,39 + ,12 + ,7 + ,9 + ,0 + ,1 + ,33 + ,29 + ,7 + ,6 + ,15 + ,0 + ,0 + ,33 + ,32 + ,16 + ,11 + ,13 + ,0 + ,1 + ,26 + ,34 + ,5 + ,4 + ,13 + ,0 + ,0 + ,35 + ,38 + ,16 + ,9 + ,13 + ,0 + ,0 + ,21 + ,17 + ,4 + ,5 + ,15 + ,0 + ,0 + ,38 + ,35 + ,12 + ,9 + ,15 + ,0 + ,0 + ,35 + ,32 + ,15 + ,11 + ,14 + ,0 + ,1 + ,33 + ,34 + ,14 + ,12 + ,15 + ,0 + ,0 + ,37 + ,36 + ,11 + ,9 + ,11 + ,0 + ,0 + ,38 + ,31 + ,16 + ,12 + ,15 + ,0 + ,1 + ,34 + ,35 + ,15 + ,10 + ,14 + ,0 + ,0 + ,27 + ,29 + ,12 + ,9 + ,13 + ,0 + ,1 + ,16 + ,22 + ,6 + ,6 + ,12 + ,0 + ,0 + ,40 + ,41 + ,16 + ,10 + ,16 + ,0 + ,0 + ,36 + ,36 + ,10 + ,9 + ,16 + ,0 + ,1 + ,42 + ,42 + ,15 + ,13 + ,9 + ,0 + ,1 + ,30 + ,33 + ,14 + ,12 + ,14) + ,dim=c(7 + ,288) + ,dimnames=list(c('Pop' + ,'Gender' + ,'Connected' + ,'Separate' + ,'Learning' + ,'Software' + ,'Happiness') + ,1:288)) > y <- array(NA,dim=c(7,288),dimnames=list(c('Pop','Gender','Connected','Separate','Learning','Software','Happiness'),1:288)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par2 = 'No' > par1 = '1' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Dr. Ian E. Holliday > #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: > #Technical description: > library(rpart) > library(partykit) Loading required package: grid Loading required package: mvtnorm > par1 <- as.numeric(par1) > autoprune <- function ( tree, method='Minimum CV'){ + xerr <- tree$cptable[,'xerror'] + cpmin.id <- which.min(xerr) + if (method == 'Minimum CV Error plus 1 SD'){ + xstd <- tree$cptable[,'xstd'] + errt <- xerr[cpmin.id] + xstd[cpmin.id] + cpSE1.min <- which.min( errt < xerr ) + mycp <- (tree$cptable[,'CP'])[cpSE1.min] + } + if (method == 'Minimum CV') { + mycp <- (tree$cptable[,'CP'])[cpmin.id] + } + return (mycp) + } > conf.multi.mat <- function(true, new) + { + if ( all( is.na(match( levels(true),levels(new) ) )) ) + stop ( 'conflict of vector levels') + multi.t <- list() + for (mylev in levels(true) ) { + true.tmp <- true + new.tmp <- new + left.lev <- levels (true.tmp)[- match(mylev,levels(true) ) ] + levels(true.tmp) <- list ( mylev = mylev, all = left.lev ) + levels(new.tmp) <- list ( mylev = mylev, all = left.lev ) + curr.t <- conf.mat ( true.tmp , new.tmp ) + multi.t[[mylev]] <- curr.t + multi.t[[mylev]]$precision <- + round( curr.t$conf[1,1] / sum( curr.t$conf[1,] ), 2 ) + } + return (multi.t) + } > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > m <- rpart(as.data.frame(x1)) > par2 [1] "No" > if (par2 != 'No') { + mincp <- autoprune(m,method=par2) + print(mincp) + m <- prune(m,cp=mincp) + } > m$cptable CP nsplit rel error xerror xstd 1 0.10262902 0 1.0000000 1.0096809 0.01190953 2 0.03104743 1 0.8973710 0.9644950 0.03874415 3 0.02332768 2 0.8663235 0.9727186 0.04580101 4 0.01781639 5 0.7963405 0.9998951 0.05182210 5 0.01000000 7 0.7607077 1.0004916 0.05778712 > postscript(file="/var/www/html/rcomp/tmp/1cbmo1293206188.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(as.party(m),tp_args=list(id=FALSE)) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2cbmo1293206188.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plotcp(m) > dev.off() null device 1 > cbind(y=m$y,pred=predict(m),res=residuals(m)) y pred res 1 1 0.5185185 0.4814815 2 1 0.7647059 0.2352941 3 1 0.8787879 0.1212121 4 1 0.4509804 0.5490196 5 1 0.8787879 0.1212121 6 1 0.6000000 0.4000000 7 1 0.8750000 0.1250000 8 1 0.8787879 0.1212121 9 1 0.7647059 0.2352941 10 1 0.7647059 0.2352941 11 1 0.7647059 0.2352941 12 1 0.7647059 0.2352941 13 1 0.3333333 0.6666667 14 1 0.7647059 0.2352941 15 1 0.8787879 0.1212121 16 1 0.4509804 0.5490196 17 1 0.4509804 0.5490196 18 1 0.7647059 0.2352941 19 1 0.4509804 0.5490196 20 1 0.7647059 0.2352941 21 1 0.7647059 0.2352941 22 1 0.8750000 0.1250000 23 1 0.8787879 0.1212121 24 1 0.4509804 0.5490196 25 1 0.7647059 0.2352941 26 1 0.3333333 0.6666667 27 1 0.7647059 0.2352941 28 1 0.8787879 0.1212121 29 1 0.7647059 0.2352941 30 1 0.4509804 0.5490196 31 1 0.7647059 0.2352941 32 1 0.2153846 0.7846154 33 1 0.7647059 0.2352941 34 1 0.7647059 0.2352941 35 1 0.8750000 0.1250000 36 1 0.6000000 0.4000000 37 1 0.2153846 0.7846154 38 1 0.7647059 0.2352941 39 1 0.7647059 0.2352941 40 1 0.4509804 0.5490196 41 1 0.7647059 0.2352941 42 1 0.7647059 0.2352941 43 1 0.8787879 0.1212121 44 1 0.4509804 0.5490196 45 1 0.7647059 0.2352941 46 1 0.5185185 0.4814815 47 1 0.8787879 0.1212121 48 1 0.7647059 0.2352941 49 1 0.7647059 0.2352941 50 1 0.4509804 0.5490196 51 1 0.8787879 0.1212121 52 1 0.8750000 0.1250000 53 1 0.5185185 0.4814815 54 1 0.8787879 0.1212121 55 1 0.2153846 0.7846154 56 1 0.7647059 0.2352941 57 1 0.7647059 0.2352941 58 1 0.7647059 0.2352941 59 1 0.4509804 0.5490196 60 1 0.4509804 0.5490196 61 1 0.2153846 0.7846154 62 1 0.4509804 0.5490196 63 1 0.8787879 0.1212121 64 1 0.4509804 0.5490196 65 1 0.7647059 0.2352941 66 1 0.3333333 0.6666667 67 1 0.8787879 0.1212121 68 1 0.8750000 0.1250000 69 1 0.5185185 0.4814815 70 1 0.7647059 0.2352941 71 1 0.2153846 0.7846154 72 1 0.7647059 0.2352941 73 1 0.8787879 0.1212121 74 1 0.8787879 0.1212121 75 1 0.8750000 0.1250000 76 1 0.8750000 0.1250000 77 1 0.8787879 0.1212121 78 1 0.8750000 0.1250000 79 1 0.5185185 0.4814815 80 1 0.8750000 0.1250000 81 1 0.7647059 0.2352941 82 1 0.7647059 0.2352941 83 1 0.7647059 0.2352941 84 1 0.7647059 0.2352941 85 1 0.8750000 0.1250000 86 1 0.7647059 0.2352941 87 1 0.4509804 0.5490196 88 1 0.4509804 0.5490196 89 1 0.3333333 0.6666667 90 1 0.6000000 0.4000000 91 1 0.8787879 0.1212121 92 1 0.7647059 0.2352941 93 1 0.8750000 0.1250000 94 1 0.5185185 0.4814815 95 1 0.7647059 0.2352941 96 1 0.8787879 0.1212121 97 1 0.8750000 0.1250000 98 1 0.7647059 0.2352941 99 1 0.7647059 0.2352941 100 1 0.8787879 0.1212121 101 1 0.4509804 0.5490196 102 1 0.8787879 0.1212121 103 1 0.8750000 0.1250000 104 1 0.7647059 0.2352941 105 1 0.5185185 0.4814815 106 1 0.7647059 0.2352941 107 1 0.7647059 0.2352941 108 1 0.3333333 0.6666667 109 1 0.2153846 0.7846154 110 1 0.7647059 0.2352941 111 1 0.8787879 0.1212121 112 1 0.4509804 0.5490196 113 1 0.6000000 0.4000000 114 1 0.8787879 0.1212121 115 1 0.8787879 0.1212121 116 1 0.4509804 0.5490196 117 1 0.5185185 0.4814815 118 1 0.4509804 0.5490196 119 1 0.5185185 0.4814815 120 1 0.8787879 0.1212121 121 1 0.5185185 0.4814815 122 1 0.8787879 0.1212121 123 1 0.8787879 0.1212121 124 1 0.8787879 0.1212121 125 1 0.7647059 0.2352941 126 1 0.7647059 0.2352941 127 1 0.7647059 0.2352941 128 1 0.5185185 0.4814815 129 1 0.4509804 0.5490196 130 1 0.5185185 0.4814815 131 1 0.7647059 0.2352941 132 1 0.5185185 0.4814815 133 1 0.7647059 0.2352941 134 1 0.7647059 0.2352941 135 1 0.2153846 0.7846154 136 1 0.4509804 0.5490196 137 1 0.2153846 0.7846154 138 1 0.2153846 0.7846154 139 1 0.2153846 0.7846154 140 1 0.7647059 0.2352941 141 1 0.7647059 0.2352941 142 1 0.8787879 0.1212121 143 1 0.7647059 0.2352941 144 1 0.6000000 0.4000000 145 1 0.8787879 0.1212121 146 1 0.5185185 0.4814815 147 1 0.8750000 0.1250000 148 1 0.2153846 0.7846154 149 1 0.6000000 0.4000000 150 1 0.8787879 0.1212121 151 1 0.7647059 0.2352941 152 1 0.4509804 0.5490196 153 1 0.3333333 0.6666667 154 1 0.2153846 0.7846154 155 1 0.2153846 0.7846154 156 1 0.7647059 0.2352941 157 1 0.4509804 0.5490196 158 1 0.2153846 0.7846154 159 0 0.2153846 -0.2153846 160 0 0.2153846 -0.2153846 161 0 0.8787879 -0.8787879 162 0 0.3333333 -0.3333333 163 0 0.2153846 -0.2153846 164 0 0.4509804 -0.4509804 165 0 0.5185185 -0.5185185 166 0 0.5185185 -0.5185185 167 0 0.3333333 -0.3333333 168 0 0.2153846 -0.2153846 169 0 0.7647059 -0.7647059 170 0 0.2153846 -0.2153846 171 0 0.2153846 -0.2153846 172 0 0.5185185 -0.5185185 173 0 0.4509804 -0.4509804 174 0 0.2153846 -0.2153846 175 0 0.2153846 -0.2153846 176 0 0.7647059 -0.7647059 177 0 0.2153846 -0.2153846 178 0 0.4509804 -0.4509804 179 0 0.2153846 -0.2153846 180 0 0.8787879 -0.8787879 181 0 0.4509804 -0.4509804 182 0 0.5185185 -0.5185185 183 0 0.4509804 -0.4509804 184 0 0.4509804 -0.4509804 185 0 0.2153846 -0.2153846 186 0 0.2153846 -0.2153846 187 0 0.2153846 -0.2153846 188 0 0.4509804 -0.4509804 189 0 0.4509804 -0.4509804 190 0 0.2153846 -0.2153846 191 0 0.4509804 -0.4509804 192 0 0.2153846 -0.2153846 193 0 0.4509804 -0.4509804 194 0 0.8750000 -0.8750000 195 0 0.4509804 -0.4509804 196 0 0.3333333 -0.3333333 197 0 0.3333333 -0.3333333 198 0 0.7647059 -0.7647059 199 0 0.6000000 -0.6000000 200 0 0.5185185 -0.5185185 201 0 0.7647059 -0.7647059 202 0 0.7647059 -0.7647059 203 0 0.4509804 -0.4509804 204 0 0.7647059 -0.7647059 205 0 0.2153846 -0.2153846 206 0 0.7647059 -0.7647059 207 0 0.5185185 -0.5185185 208 0 0.5185185 -0.5185185 209 0 0.4509804 -0.4509804 210 0 0.3333333 -0.3333333 211 0 0.4509804 -0.4509804 212 0 0.7647059 -0.7647059 213 0 0.7647059 -0.7647059 214 0 0.4509804 -0.4509804 215 0 0.2153846 -0.2153846 216 0 0.7647059 -0.7647059 217 0 0.5185185 -0.5185185 218 0 0.7647059 -0.7647059 219 0 0.3333333 -0.3333333 220 0 0.2153846 -0.2153846 221 0 0.6000000 -0.6000000 222 0 0.7647059 -0.7647059 223 0 0.4509804 -0.4509804 224 0 0.2153846 -0.2153846 225 0 0.2153846 -0.2153846 226 0 0.2153846 -0.2153846 227 0 0.2153846 -0.2153846 228 0 0.7647059 -0.7647059 229 0 0.8787879 -0.8787879 230 0 0.6000000 -0.6000000 231 0 0.2153846 -0.2153846 232 0 0.2153846 -0.2153846 233 0 0.5185185 -0.5185185 234 0 0.3333333 -0.3333333 235 0 0.2153846 -0.2153846 236 0 0.2153846 -0.2153846 237 0 0.5185185 -0.5185185 238 0 0.2153846 -0.2153846 239 0 0.2153846 -0.2153846 240 0 0.2153846 -0.2153846 241 0 0.4509804 -0.4509804 242 0 0.2153846 -0.2153846 243 0 0.2153846 -0.2153846 244 0 0.2153846 -0.2153846 245 0 0.2153846 -0.2153846 246 0 0.2153846 -0.2153846 247 0 0.2153846 -0.2153846 248 0 0.4509804 -0.4509804 249 0 0.2153846 -0.2153846 250 0 0.4509804 -0.4509804 251 0 0.3333333 -0.3333333 252 0 0.3333333 -0.3333333 253 0 0.5185185 -0.5185185 254 0 0.2153846 -0.2153846 255 0 0.8750000 -0.8750000 256 0 0.5185185 -0.5185185 257 0 0.2153846 -0.2153846 258 0 0.4509804 -0.4509804 259 0 0.3333333 -0.3333333 260 0 0.3333333 -0.3333333 261 0 0.5185185 -0.5185185 262 0 0.2153846 -0.2153846 263 0 0.4509804 -0.4509804 264 0 0.2153846 -0.2153846 265 0 0.2153846 -0.2153846 266 0 0.2153846 -0.2153846 267 0 0.4509804 -0.4509804 268 0 0.2153846 -0.2153846 269 0 0.4509804 -0.4509804 270 0 0.8787879 -0.8787879 271 0 0.2153846 -0.2153846 272 0 0.2153846 -0.2153846 273 0 0.4509804 -0.4509804 274 0 0.2153846 -0.2153846 275 0 0.4509804 -0.4509804 276 0 0.2153846 -0.2153846 277 0 0.2153846 -0.2153846 278 0 0.4509804 -0.4509804 279 0 0.7647059 -0.7647059 280 0 0.2153846 -0.2153846 281 0 0.7647059 -0.7647059 282 0 0.4509804 -0.4509804 283 0 0.2153846 -0.2153846 284 0 0.2153846 -0.2153846 285 0 0.7647059 -0.7647059 286 0 0.6000000 -0.6000000 287 0 0.3333333 -0.3333333 288 0 0.4509804 -0.4509804 > myr <- residuals(m) > myp <- predict(m) > postscript(file="/var/www/html/rcomp/tmp/3424r1293206188.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > op <- par(mfrow=c(2,2)) > plot(myr,ylab='residuals') > plot(density(myr),main='Residual Kernel Density') > plot(myp,myr,xlab='predicted',ylab='residuals',main='Predicted vs Residuals') > plot(density(myp),main='Prediction Kernel Density') > par(op) > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Model Performance',6,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'#',header=TRUE) > a<-table.element(a,'Complexity',header=TRUE) > a<-table.element(a,'split',header=TRUE) > a<-table.element(a,'relative error',header=TRUE) > a<-table.element(a,'CV error',header=TRUE) > a<-table.element(a,'CV S.D.',header=TRUE) > a<-table.row.end(a) > for (i in 1:length(m$cptable[,1])) { + a<-table.row.start(a) + a<-table.element(a,i,header=TRUE) + a<-table.element(a,round(m$cptable[i,'CP'],3)) + a<-table.element(a,m$cptable[i,'nsplit']) + a<-table.element(a,round(m$cptable[i,'rel error'],3)) + a<-table.element(a,round(m$cptable[i,'xerror'],3)) + a<-table.element(a,round(m$cptable[i,'xstd'],3)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/48lkx1293206188.tab") > > try(system("convert tmp/1cbmo1293206188.ps tmp/1cbmo1293206188.png",intern=TRUE)) character(0) > try(system("convert tmp/2cbmo1293206188.ps tmp/2cbmo1293206188.png",intern=TRUE)) character(0) > try(system("convert tmp/3424r1293206188.ps tmp/3424r1293206188.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.182 0.613 5.104