#################################################################################
# R functions used to produce MN county areal and edge plots			#
# Need to install R package "maps" first. Download it from www.r-project.org	#
# 			8/15/2005 Haijun Ma					#
#################################################################################

#########################################################################
# R function used to produce grayscale choropleth areal MN county maps  #
# figtitle: title for figure						#
# y: 		areal level observation					#
# n.col:	number of color levels					#
# type:	"q" produces equal quantile cutoffs				#
#		"e" produced equal interval cutoffs			#
#########################################################################

mn.Poly<-function(figtitle, y , n.col, type="q"){
# old.par <- par(no.readonly = TRUE) 
 

# my data order is a little bit different from R plotting order "orders$forPlot"
R_order<-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, 46, 43, 44, 45, 47, 48, 49, 
50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 
66, 67, 68, 72, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 80, 81, 
82, 83, 84, 85, 86, 87)

y<-y[R_order]

library(maps)
if(type=="q"){p <- seq(0,1, length=n.col+1)
              br <- round(quantile(y, probs=p),2)}
if(type=="e"){br <- round(seq(min(y), max(y), length=n.col+1),2)}

# 0: dark 1: light
shading<-gray((n.col-1):0/(n.col-1))
y.grp<-findInterval(y, vec=br, rightmost.closed = TRUE, all.inside = TRUE)
y.shad<-shading[y.grp]

 
map("county","minnesota", fill=T, plot=T, col=y.shad, interior=T)
leg.txt<-paste("[",br[n.col],",",br[n.col+1],"]",sep="")
for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste("[",br[i],",",br[i+1],")",sep=""),) }
   leg.txt<-rev(leg.txt)
#leg.pos<-matrix(c(-92,46.5,-91,45.5), byrow=T, nrow=2)
#legend(leg.pos,legend=leg.txt,fill=shading,cex=0.9,ncol=1,bty="n")
legend(-92.25,46.5,legend=leg.txt,fill=shading,cex=1,ncol=1,bty="n")

title(main=figtitle, cex.main=2)
#title(sub=figtitle,cex.sub=2)

#par(old.par)
}

#dev.print(device=postscript, "f:/Boundary/Latex/Plots/Esop_rSMR.eps",onefile=FALSE, horizontal=FALSE, paper="special")




#########################################################################
# R functions used to produce grayscale MN county edge maps  		#
# figtitle: title for figure						#
# y: 		edge level observation					#
# breaks:	a vector of given quantiles if qtle=T,			#
#		o.w. is a vector of cut-off points.			#
#########################################################################

# data should be probability of that the edge is part of the boundary, i.e. 1-wij/1-pij
# plotting order is absolutely alphabetic in county and in edges (R is not alphabetic in edge order)


mn.Edge.crisp<-function(figtitle,y,qtle=T,breaks){

library(maps)
map("county","minnesota", fill=F, plot=T, interior=F)
if(qtle==T) { br <- round(quantile(y, probs=breaks),2)
                n.col<- length(breaks)-1

    # 0: dark 1: light
    shading<-gray((n.col-1):0/n.col)
    y.grp<-findInterval(y, vec=br, rightmost.closed = TRUE, all.inside = TRUE)
    y.shad<-shading[y.grp]

   #leg.txt<-paste("[",br[n.col],",",br[n.col+1],"]",sep="")
   #for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste("[",br[i],",",br[i+1],,sep=""),) }
   leg.txt<-paste(breaks[n.col]*100,"%~",breaks[n.col+1]*100,"%",sep="")
   for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste( breaks[i]*100,"%~",breaks[i+1]*100,"%", sep=""),) }
   }

if(qtle==F) { br <- round(breaks,2)
             n.col<- length(breaks)-1

    # 0: dark 1: light
    shading<-gray((n.col-1):0/n.col)
    y.grp<-findInterval(y, vec=br, rightmost.closed = TRUE, all.inside = TRUE)
    y.shad<-shading[y.grp]

   #leg.txt<-paste("[",br[n.col],",",br[n.col+1],"]",sep="")
   #for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste("[",br[i],",",br[i+1],,sep=""),) }
   leg.txt<-paste("[",br[n.col],",",br[n.col+1],"]",sep="")
   for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste("[",br[i],",",br[i+1],")",sep=""),) }
   }

   leg.txt<-rev(leg.txt)
   leg.pos<-matrix(c(-92.25,46,-92,45.5), byrow=T, nrow=2)
   legend(leg.pos,legend=leg.txt,col=shading,lty="solid",lwd=5,cex=1,ncol=1,bty="n")

#title(sub=figtitle,cex.sub=2)
title(main=figtitle,cex.main=2)

edgelines<-structure(list(V1 = c(4, 5, 6, 7, 8, 24, 25, 26, 19, 20, 21, 
22, 23, 27, 28, 29, 12, 13, 14, 15, 16, 17, 18, 9, 10, 11, 1, 
2, 3, 36, 37, 38, 46, 47, 48, 49, 50, 51, 52, 53, 54, 33, 34, 
35, 42, 43, 44, 45, 30, 31, 32, 39, 40, 41, 58, 59, 60, 61, 62, 
69, 70, 71, 72, 73, 74, 75, 76, 66, 67, 68, 63, 64, 65, 55, 56, 
57, 77, 78, 79, 96, 97, 98, 99, 103, 104, 105, 106, 107, 108, 
109, 110, 111, 112, 113, 114, 115, 100, 101, 102, 93, 94, 95, 
88, 89, 90, 91, 92, 83, 84, 85, 86, 87, 119, 120, 121, 122, 123, 
116, 117, 118, 80, 81, 82, 124, 125, 126, 141, 142, 143, 144, 
145, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 
139, 140, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 
162, 163, 181, 182, 183, 184, 146, 147, 148, 149, 150, 178, 179, 
180, 188, 189, 190, 214, 215, 216, 207, 208, 209, 210, 217, 218, 
219, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 
203, 204, 205, 206, 211, 212, 213, 185, 186, 187, 239, 240, 241, 
242, 243, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 
231, 232, 233, 234, 235, 244, 245, 246, 247, 248, 249, 250, 251, 
252, 253, 236, 237, 238, 261, 262, 263, 254, 255, 256, 257, 289, 
290, 291, 292, 293, 283, 284, 285, 264, 265, 266, 267, 268, 269, 
270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 
286, 287, 288, 336, 337, 338, 339, 340, 341, 342, 375, 376, 377, 
378, 379, 380, 381, 382, 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, 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, 383, 384, 385, 407, 408, 409, 410, 411, 412, 413, 414, 415, 
416, 386, 387, 388, 389, 417, 418, 419, 390, 391, 392, 393, 394, 
395, 396, 397, 398, 399, 400, 401, 402, 403, 404, 405, 406, 435, 
436, 437, 438, 439, 440, 441, 442, 432, 433, 434, 467, 468, 469, 
447, 448, 449, 450, 451, 452, 453, 473, 474, 475, 476, 477, 478, 
479, 480, 470, 471, 472, 481, 482, 483, 484, 485, 534, 535, 536, 
537, 538, 542, 543, 544, 545, 546, 547, 548, 549, 550, 539, 540, 
541, 551, 552, 553, 554, 555, 556, 557, 558, 559, 560, 561, 562, 
563, 564, 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 611, 
612, 613, 576, 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 
587, 588, 565, 566, 567, 568, 569, 570, 571, 572, 573, 574, 575, 
589, 590, 591, 592, 593, 594, 595, 596, 597, 617, 618, 619, 624, 
625, 626, 620, 621, 622, 623, 614, 615, 616, 630, 631, 632, 633, 
634, 635, 636, 637, 627, 628, 629, 643, 644, 645, 638, 639, 640, 
641, 642, 658, 659, 660, 652, 653, 654, 655, 656, 657, 677, 678, 
679, 667, 668, 669, 670, 671, 672, 673, 674, 675, 676, 692, 693, 
694, 689, 690, 691, 686, 687, 688, 718, 719, 720, 695, 696, 697, 
712, 713, 714, 715, 716, 717, 732, 733, 734, 721, 722, 723, 724, 
725, 726, 727, 728, 729, 730, 731, 759, 760, 761, 762, 763, 735, 
736, 737, 738, 739, 740, 741, 742, 743, 744, 745, 746, 747, 748, 
749, 750, 751, 752, 753, 754, 755, 756, 757, 758, 770, 771, 772, 
788, 789, 790, 797, 798, 799, 794, 795, 796, 791, 792, 793, 800, 
801, 802, 803, 804, 805, 806, 807, 808, 809, 810, 811, 812, 813, 
826, 827, 828, 823, 824, 825, 829, 830, 831, 832, 833, 834, 835, 
836, 837, 838, 843, 844, 845, 846, 847, 856, 857, 858, 848, 849, 
850, 839, 840, 841, 842, 851, 852, 853, 854, 855, 882, 883, 884, 
877, 878, 879, 880, 881, 885, 886, 887, 908, 909, 910, 911, 912, 
913, 914, 915, 916, 917, 918, 919, 920, 921, 922, 923, 924, 925, 
942, 943, 944, 945, 946, 967, 968, 969, 970, 971, 972, 973, 974, 
975, 976, 977, 978, 979, 980, 981, 982, 983, 1022, 1023, 1024, 
1025, 1026, 1027, 1028, 1029, 1016, 1017, 1018, 1013, 1014, 1015, 
1030, 1031, 1032, 1033, 1034, 1035, 1036, 1019, 1020, 1021, 1046, 
1047, 1048, 1049, 1050, 1051, 1043, 1044, 1045, 1055, 1056, 1057, 
1052, 1053, 1054, 1058, 1059, 1060, 1061, 1062, 1079, 1080, 1081, 
1082, 1083, 1084, 1088, 1089, 1090, 1091, 1092, 1093, 1085, 1086, 
1087, 1110, 1111, 1112, 1066, 1067, 1068, 1069, 1063, 1064, 1065, 
1073, 1074, 1075, 1076, 1077, 1078, 1070, 1071, 1072, 1119, 1120, 
1121, 1122, 1123, 1124, 1125, 1126, 1116, 1117, 1118, 1127, 1128, 
1129, 1130, 1131, 1132, 1133, 1134, 1135, 1136, 1137, 1138, 1139, 
1140, 1141, 1142, 1143, 1144, 1145, 1146, 1159, 1160, 1161, 1162, 
1163, 1164, 1156, 1157, 1158, 1165, 1166, 1167, 1168, 1169, 1170, 
1174, 1175, 1176, 1193, 1194, 1195, 1196, 1197, 1198, 1199, 1200, 
1201, 1202, 1203, 1212, 1213, 1214, 1215, 1216, 1209, 1210, 1211, 
1204, 1205, 1206, 1207, 1208, 1217, 1218, 1219, 1224, 1225, 1226, 
1227, 1228, 1220, 1221, 1222, 1223, 1248, 1249, 1250, 1278, 1279, 
1280, 1281, 1282, 1283, 1284, 1285, 1286, 1287, 1296, 1297, 1298, 
1299, 1300, 1291, 1292, 1293, 1294, 1295, 1288, 1289, 1290, 1301, 
1302, 1303, 1304, 1305, 1310, 1311, 1312, 1313, 1314, 1315, 1316, 
1317, 1318, 1319, 1320, 1321, 1322, 1323, 1306, 1307, 1308, 1309, 
1324, 1325, 1326, 1327, 1328, 1329, 1330, 1331, 1332, 1333, 1334, 
1335, 1342, 1343, 1344, 1336, 1337, 1338, 1339, 1340, 1341, 1411, 
1412, 1413, 1414, 1415, 1416, 1417, 1418, 1419, 1440, 1441, 1442, 
1443, 1444, 1445, 1446, 1447, 1420, 1421, 1422, 1423, 1424, 1425, 
1426, 1427, 1428, 1429, 1430, 1431, 1432, 1433, 1434, 1435, 1436, 
1437, 1438, 1439, 1448, 1449, 1450, 1451, 1452, 1453, 1454, 1455, 
1456, 1457, 1458, 1459, 1460, 1461, 1462, 1463, 1464, 1465, 1466, 
1467, 1468, 1469, 1470, 1471, 1472, 1473, 1474, 1475, 1476, 1495, 
1496, 1497, 1498, 1499, 1500, 1501), V2 = c(-93.04834, -93.05408, 
-93.03689, -93.03689, NA, -93.77027, -93.77027, NA, -93.81038, 
-93.8161, -93.78174, -93.77027, NA, -93.03689, -93.77027, NA, 
-93.03689, -93.43223, NA, -93.43223, -93.43223, -93.81038, NA, 
-93.03689, -93.03689, NA, -93.04834, -93.03689, NA, -93.0197, 
-93.01397, NA, -93.23169, -93.29472, -93.289, -93.289, -93.32336, 
-93.3864, -93.48952, -93.51817, NA, -93.51817, -93.0197, NA, 
-93.01397, -93.23169, -93.23169, NA, -93.51817, -93.51817, NA, 
-93.01397, -93.01397, NA, -96.17669, -96.17669, -96.15378, -96.14805, 
NA, -95.16829, -95.54645, NA, -95.14537, -95.13965, -95.16256, 
-95.16829, NA, -96.05064, -95.54645, NA, -96.05064, -96.17669, 
NA, -95.14537, -96.14805, NA, -95.14537, -95.14537, NA, -94.38907, 
-94.3948, -94.66409, NA, -95.56936, -95.20267, -95.21413, -95.25423, 
-95.2657, -95.2657, -95.25423, -95.23132, -95.17975, -95.18548, 
-95.16829, -95.16829, NA, -94.66409, -95.16829, NA, -94.38907, 
-94.3948, NA, -94.3948, -94.40053, -94.3948, -94.3948, NA, -95.32299, 
-95.19694, -95.19694, -94.3948, NA, -95.58655, -95.58655, -95.57509, 
-95.57509, NA, -95.56936, -95.57509, NA, -95.58655, -95.32299, 
NA, -93.77027, -93.77027, NA, -93.77027, -94.35469, -94.32604, 
-94.28593, NA, -93.77027, -94.15988, NA, -94.15988, -94.1828, 
-94.20572, -94.19999, -94.22864, -94.2229, -94.2229, -94.25729, 
-94.26302, -94.28593, NA, -96.09075, -96.12513, -96.1595, -96.16524, 
-96.19389, -96.22826, -96.26836, -96.27983, -96.2741, -96.30274, 
-96.36004, -96.43452, NA, -96.10794, -96.234, -96.23972, NA, 
-96.09075, -96.08502, -96.10794, -96.10794, NA, -96.81268, -96.23972, 
NA, -94.37188, -94.37188, NA, -94.26302, -93.776, NA, -93.776, 
-93.776, -94.00519, NA, -94.37188, -94.26302, NA, -94.37188, 
-94.32604, -94.30312, -94.2401, -94.22864, -94.1828, -94.16561, 
-94.10832, -94.09113, -94.03383, -94.01665, -94.00519, -94.01092, 
-94.01092, -94.00519, NA, -93.776, -93.776, NA, -94.37188, -94.37188, 
NA, -94.86462, -95.09954, -95.09954, -95.11672, NA, -94.37188, 
-94.38907, -94.41199, -94.4349, -94.49792, -94.5094, -94.52657, 
-94.58387, -94.65263, -94.68127, -94.69846, -94.7042, -94.76149, 
-94.7844, -94.79014, NA, -94.85889, -94.87035, -95.11672, -95.11672, 
NA, -94.85889, -94.82452, -94.79587, -94.79014, NA, -94.37188, 
-94.86462, NA, -93.03689, -92.29205, NA, -93.04834, -92.29778, 
-92.3035, NA, -93.51817, -93.51817, -93.76454, -93.76454, NA, 
-94.01665, -94.01092, NA, -93.76454, -93.70152, -93.68433, -93.6614, 
-93.63276, -93.61558, -93.62703, -93.60985, -93.60985, -93.58692, 
-93.57547, -93.51817, NA, -93.76454, -93.76454, -93.87914, -93.88486, 
-94.01092, NA, -94.01665, -93.76454, NA, -94.34896, -94.32604, 
-94.31458, -94.33177, -94.33177, -93.77027, NA, -94.66409, -94.66982, 
-94.65836, -94.65836, -94.6469, -94.6469, -94.77868, NA, -94.38907, 
-94.36615, -94.3375, -94.31458, -94.29166, -94.28593, -94.26302, 
-94.25729, -94.24583, -94.24583, -94.2401, -94.22864, -94.21145, 
-94.19999, -94.17707, -94.15416, -94.10832, -94.07394, -94.04529, 
-94.0281, -93.99372, -93.99372, -94.01092, -94.02238, -94.02238, 
-94.01092, -93.95361, -93.91924, -93.89632, -93.88486, -93.86768, 
-93.83903, -93.81038, -93.78174, -93.77027, -93.776, -93.79319, 
-93.79892, -93.78746, -93.77027, -93.77027, NA, -94.34896, -94.35469, 
-94.35469, -94.38334, -94.40053, -94.41199, -94.4349, -94.45782, 
-94.4922, -94.5323, -94.55523, -94.57241, -94.6068, NA, -94.6068, 
-94.6297, -94.64117, -94.65836, -94.67554, -94.68701, -94.69273, 
-94.70992, -94.73284, -94.73857, -94.76149, -94.77295, NA, -94.77295, 
-94.76722, -94.77868, -94.77868, NA, -95.24277, -95.24277, NA, 
-96.03918, -95.97043, -95.9246, -95.89594, -95.83865, -95.81573, 
-95.80427, -95.76989, -95.74698, NA, -95.24277, -95.48341, -95.48914, 
NA, -96.03918, -95.24277, NA, -95.48914, -95.51207, -95.52925, 
-95.52352, -95.50061, -95.55218, -95.57509, -95.58082, -95.62093, 
-95.63811, -95.64384, -95.67822, -95.68395, -95.71833, -95.71833, 
-95.74698, NA, -93.14001, -93.14001, -93.01397, -93.0197, NA, 
-93.14001, -92.8478, NA, -92.76186, -93.01397, NA, -96.81841, 
-96.17669, NA, -96.2569, -96.15378, -96.14805, NA, -96.2569, 
-96.7783, NA, -95.16829, -95.17402, -95.16256, -95.16829, NA, 
-95.54645, -95.54645, NA, -95.56936, -95.56936, NA, -95.54645, 
-95.54645, -95.56936, -95.56936, NA, -91.00288, -91.00288, -91.02581, 
-91.01434, NA, -95.47768, -94.86462, NA, -95.47768, -95.46623, 
NA, -95.11672, -95.46623, NA, -94.86462, -94.86462, NA, -93.81038, 
-93.81038, NA, -93.81038, -94.40053, -94.41772, -94.41199, -94.38907, 
-94.3948, -94.3948, -94.3948, -94.37188, -94.34896, NA, -93.03689, 
-93.03689, -93.01397, -92.98532, -92.97386, -92.94521, -92.91656, 
-92.91656, -92.79625, -92.80197, -92.73322, -92.73895, NA, -93.1744, 
-93.19732, -93.24316, -93.28326, -93.3291, NA, -93.01397, -93.02543, 
-93.03689, -93.12283, -93.14001, -93.1744, NA, -93.28326, -93.2718, 
-93.03689, NA, -93.28326, -93.2718, -93.31763, -93.31763, -93.3291, 
-93.3291, NA, -92.80197, -92.84208, -92.89938, -92.91656, -93.01397, 
-93.01397, -93.00824, -93.01397, NA, -93.04834, -92.6702, NA, 
-93.05408, -92.69312, NA, -92.69312, -92.67592, -92.6702, NA, 
-93.05408, -93.04834, NA, -95.75844, -95.75844, -95.76989, -95.76989, 
NA, -95.14537, -95.76989, NA, -95.75844, -95.13965, NA, -95.13965, 
-95.13391, NA, -95.14537, -95.1511, -95.13391, -95.13391, NA, 
-93.64423, -93.64996, NA, -94.26302, -94.26302, NA, -93.776, 
-93.64996, NA, -91.72482, -91.71909, NA, -92.4582, -92.4582, 
NA, -92.4582, -92.45247, -92.08578, NA, -92.08578, -91.71909, 
NA, -93.05408, -93.0598, NA, -93.05408, -93.41504, NA, -93.64996, 
-93.41504, NA, -92.5556, -92.6702, NA, -93.04834, -93.03689, 
NA, -92.5556, -92.5556, -92.42383, -92.42955, -92.25194, NA, 
-95.76989, -96.25118, NA, -95.75844, -96.23972, NA, -96.25118, 
-96.2569, -96.23972, -96.23972, NA, -96.25118, -96.25118, NA, 
-93.23169, -93.20877, -93.20305, -93.1744, NA, -93.51817, -93.4838, 
-93.45515, -93.45515, -93.42077, -93.3864, -93.36347, -93.3291, 
NA, -93.76454, -93.76454, -93.75881, -93.72443, -93.70152, -93.68433, 
-93.64996, -93.64996, -93.6614, -93.64423, -93.60412, -93.58692, 
-93.53536, -93.51817, NA, -91.71909, -91.28937, NA, -95.14537, 
-94.77868, NA, -93.51817, -93.14001, NA, -93.51817, -93.51817, 
NA, -93.51817, -93.51817, NA, -94.3948, -93.75308, -93.75308, 
-93.0598, NA, -93.0598, -93.06554, -93.05408, -93.0598, -93.04834, 
-93.05408, -93.03689, -93.03689, NA, -94.86462, -94.86462, NA, 
-95.46623, -95.47768, NA, -93.51817, -93.51817, -93.43223, -93.43223, 
NA, -93.14001, -93.14001, -93.04834, -93.03689, NA, -94.76149, 
-94.75576, -94.75003, -94.75003, NA, -95.12818, -95.25996, NA, 
-95.24277, -94.75003, NA, -95.12818, -94.76149, -94.76149, NA, 
-95.25996, -95.25423, -95.24277, -95.24277, NA, -97.145, -96.34858, 
NA, -96.3715, -96.37724, -96.36004, -96.34858, NA, -94.40053, 
-94.3948, NA, -93.07127, -93.09418, -93.07127, -93.077, -93.06554, 
-93.0598, NA, -96.09075, -96.07356, -96.05637, -96.03918, NA, 
-96.43452, -95.86156, -95.86156, -95.74698, -95.74698, NA, -91.78784, 
-91.7821, -91.79357, -91.77065, NA, -95.32299, -95.32872, -95.08234, 
-95.07661, -95.09954, -95.12245, -95.14537, -95.17402, -95.2084, 
-95.2485, -95.27715, -95.28288, -95.29434, -95.28861, -95.31152, 
-95.29434, NA, -93.93643, -93.93643, -93.95935, -93.93643, -93.9479, 
-93.99372, -94.00519, NA, -93.5239, -93.5239, NA, -93.89632, 
-93.5239, NA, -93.89632, -93.89632, -93.89632, -93.90205, -93.92497, 
-93.93643, NA, -93.5239, -93.776, NA, -96.05064, -96.06783, NA, 
-96.05064, -96.4288, NA, -96.43452, -96.06783, NA, -95.59229, 
-96.05064, NA, -95.59229, -95.598, NA, -96.06783, -95.60947, 
-95.60947, -95.598, NA, -96.05064, -96.05064, NA, -95.54645, 
-96.05064, NA, -96.4689, -95.57509, NA, -96.4689, -97.145, NA, 
-95.58655, -96.34858, NA, -94.37188, -94.86462, NA, -94.25156, 
-94.50366, -94.49792, NA, -94.49792, -94.49792, NA, -94.49792, 
-94.4922, -94.25156, -94.25156, -94.01092, NA, -94.01665, -94.25156, 
NA, -94.49792, -94.75003, NA, -94.76149, -94.3776, -94.38334, 
-94.26302, NA, -94.25156, -94.26302, NA, -93.77027, -93.776, 
-93.81038, -93.81038, NA, -93.51817, -93.77027, NA, -94.63544, 
-94.28593, NA, -94.6068, -94.61825, -94.65263, -94.63544, NA, 
-92.69312, -92.68738, -92.4582, NA, -96.05064, -95.47768, NA, 
-96.05064, -96.05064, NA, -95.59229, -95.46623, NA, -94.79014, 
-94.61825, NA, -94.61825, -93.93643, NA, -96.04491, -96.05064, 
NA, -96.05064, -96.8356, NA, -92.08005, -92.32069, -92.32069, 
-92.5556, NA, -92.08005, -92.08578, NA, -95.14537, -95.13965, 
-95.1511, -95.14537, NA, -95.14537, -95.14537, NA, -96.2569, 
-96.2741, -96.2569, -96.25118, NA, -95.56936, -95.69541, NA, 
-96.4689, -96.46317, -96.45745, -96.45745, NA, -95.69541, -95.69541, 
-96.45745, NA, -96.4288, -96.05064, NA, -95.69541, -95.70114, 
-95.82146, -95.82146, -96.0793, -96.3314, -96.3314, -96.45745, 
-96.45745, NA, -95.12818, -95.12818, -95.13391, -95.13965, NA, 
-95.75844, -95.75844, -95.74698, -95.74698, NA, -95.25996, -95.74698, 
NA, -93.01397, -92.9796, -92.9796, -93.01397, NA, -94.85889, 
-94.899, -95.0537, -95.07661, -95.08807, -95.08234, -95.13965, 
-95.14537, -95.18548, -95.23132, -95.28288, -95.31152, -95.35163, 
NA, -95.35163, -95.35736, -95.598, NA, -94.49792, -94.62398, 
-94.6297, -94.61252, -94.61825, NA, -95.35163, -95.36883, -95.40894, 
-95.44904, -95.48914, NA, -93.28326, -93.5239, NA, -93.04834, 
-93.40931, NA, -93.40931, -93.5239, NA, -93.76454, -93.76454, 
-93.78746, -93.8333, -93.89632, -93.8906, -93.90205, -93.89632, 
NA, -94.15988, -94.14843, -94.14843, -94.11977, -94.1026, -94.05103, 
-94.04529, NA, -94.04529, -94.03383, -93.988, -93.96508, -93.93643, 
-93.91924, -93.88486, -93.82757, -93.8161, -93.73016, -93.70152, 
-93.67287, -93.63849, -93.60985, -93.57547, -93.56401, -93.55255, 
-93.53536, -93.51817, NA, -95.13391, -94.63544, NA, -94.26302, 
-94.25729, -94.2229, -94.20572, -94.1828, -94.15988, -94.16561, 
-94.15416, -94.11977, -94.1026, -94.0854, -94.06821, -94.04529, 
NA, -93.40931, -93.41504, NA, -95.74698, -96.10794, NA, -96.23972, 
-96.23972, NA, -94.77295, -95.14537, NA, -96.25118, -96.54912, 
NA, -92.08005, -92.08005, -91.83368, NA), V3 = c(46.76482, 46.5872, 
46.57574, 46.42677, NA, 46.79919, 47.01691, NA, 46.25488, 46.59293, 
46.5872, 46.79919, NA, 47.02264, 47.01691, NA, 46.16894, 46.16321, 
NA, 46.16321, 46.26061, 46.25488, NA, 46.42677, 46.16894, NA, 
46.76482, 47.02264, NA, 45.41837, 45.30378, NA, 45.05167, 45.0574, 
45.09751, 45.13189, 45.1548, 45.20637, 45.23502, 45.24648, NA, 
45.42982, 45.41837, NA, 45.13189, 45.13189, 45.05167, NA, 45.42982, 
45.24648, NA, 45.13189, 45.30378, NA, 47.15443, 46.98254, 46.98254, 
46.71325, NA, 47.1487, 47.15443, NA, 46.80492, 46.98254, 46.98254, 
47.1487, NA, 47.15443, 47.15443, NA, 47.15443, 47.15443, NA, 
46.71898, 46.71325, NA, 46.71898, 46.80492, NA, 47.42372, 47.40653, 
47.40653, NA, 48.02532, 48.02532, 48.00813, 47.97375, 47.93938, 
47.91073, 47.88781, 47.87635, 47.87635, 47.67582, 47.67582, 47.4008, 
NA, 47.40653, 47.4008, NA, 47.42372, 47.84771, NA, 48.35764, 
48.02532, 48.01959, 47.84771, NA, 48.53526, 48.53526, 48.3691, 
48.35764, NA, 48.53526, 48.35764, 48.35191, 48.17429, NA, 48.02532, 
48.17429, NA, 48.53526, 48.53526, NA, 45.83662, 45.57306, NA, 
45.83662, 45.83662, 45.80798, 45.78506, NA, 45.57306, 45.57306, 
NA, 45.57306, 45.60171, 45.63609, 45.68765, 45.72776, 45.73922, 
45.74495, 45.74495, 45.7736, 45.78506, NA, 45.17772, 45.18918, 
45.21783, 45.2121, 45.21783, 45.22356, 45.22356, 45.22929, 45.24648, 
45.24648, 45.25794, 45.26367, NA, 45.41837, 45.41837, 45.59025, 
NA, 45.17772, 45.24648, 45.24075, 45.41837, NA, 45.58452, 45.59025, 
NA, 44.28391, 44.11775, NA, 43.85992, 43.85992, NA, 44.20369, 
44.25526, 44.25526, NA, 43.85992, 43.85992, NA, 44.28391, 44.27245, 
44.26099, 44.2438, 44.22661, 44.21515, 44.20369, 44.20369, 44.18077, 
44.17505, 44.18077, 44.20369, 44.22661, 44.23807, 44.25526, NA, 
43.85992, 44.20369, NA, 43.85992, 44.11775, NA, 44.11775, 44.12348, 
44.20369, 44.20369, NA, 44.28391, 44.28391, 44.30109, 44.3412, 
44.36985, 44.38131, 44.38131, 44.38704, 44.42142, 44.42142, 44.42715, 
44.43861, 44.46152, 44.46152, 44.47298, NA, 44.50736, 44.29537, 
44.29537, 44.20369, NA, 44.50736, 44.4959, 44.49017, 44.47298, 
NA, 44.11775, 44.11775, NA, 46.42677, 46.42104, NA, 46.76482, 
46.76482, 46.65595, NA, 44.82249, 44.89697, 44.89697, 44.98865, 
NA, 44.98291, 44.72509, NA, 44.66206, 44.67925, 44.70217, 44.7079, 
44.71363, 44.72509, 44.748, 44.76519, 44.77665, 44.78238, 44.8053, 
44.82249, NA, 44.66206, 44.68498, 44.68498, 44.72509, 44.72509, 
NA, 44.98291, 44.98865, NA, 46.28353, 46.30645, 46.6273, 46.6273, 
46.80492, 46.79919, NA, 47.40653, 47.33204, 47.32631, 46.97681, 
46.97681, 46.80492, 46.80492, NA, 47.42372, 47.41799, 47.42372, 
47.41799, 47.41799, 47.42944, 47.44091, 47.45236, 47.4581, 47.46955, 
47.49247, 47.50966, 47.51539, 47.50393, 47.51539, 47.51539, 47.4982, 
47.48674, 47.46383, 47.43518, 47.41799, 47.40653, 47.4008, 47.38934, 
47.36642, 47.35496, 47.33204, 47.29767, 47.29767, 47.30912, 47.32059, 
47.31485, 47.29193, 47.29193, 47.28048, 47.26329, 47.2461, 47.22891, 
47.21172, 47.20599, 47.01691, NA, 46.28353, 46.30645, 46.31218, 
46.31791, 46.31791, 46.3351, 46.34082, 46.34082, 46.32937, 46.32364, 
46.31791, 46.32937, 46.3351, NA, 46.3351, 46.34082, 46.36374, 
46.36948, 46.36374, 46.35801, 46.35229, 46.34082, 46.34082, 46.36948, 
46.37521, 46.37521, NA, 46.37521, 46.62158, 46.6273, 46.80492, 
NA, 44.89697, 45.1548, NA, 45.1548, 45.1147, 45.07459, 45.04594, 
45.00583, 44.97146, 44.95427, 44.94281, 44.93135, NA, 44.89697, 
44.89697, 44.77665, NA, 45.1548, 45.1548, NA, 44.77665, 44.78238, 
44.78238, 44.79957, 44.8053, 44.83968, 44.83968, 44.85114, 44.87406, 
44.87978, 44.89124, 44.89697, 44.89697, 44.90843, 44.92562, 44.93135, 
NA, 45.74495, 45.57306, 45.57306, 45.41837, NA, 45.74495, 45.73349, 
NA, 45.29804, 45.30378, NA, 47.1487, 47.15443, NA, 46.63303, 
46.63303, 46.71325, NA, 46.63303, 46.63303, NA, 47.4008, 47.32631, 
47.32059, 47.1487, NA, 47.49247, 47.15443, NA, 48.02532, 47.93365, 
NA, 47.49247, 47.67009, 47.67009, 47.93365, NA, 47.45236, 48.06543, 
48.0597, 48.19721, NA, 43.85992, 43.85992, NA, 43.85992, 44.21515, 
NA, 44.20369, 44.21515, NA, 44.11775, 43.85992, NA, 46.16321, 
46.25488, NA, 46.16321, 46.16321, 46.16894, 46.19186, 46.20332, 
46.22623, 46.24342, 46.26061, 46.26634, 46.28353, NA, 44.47871, 
44.52455, 44.52455, 44.52455, 44.51309, 44.51309, 44.52455, 44.5532, 
44.5532, 44.64487, 44.64487, 44.71363, NA, 44.89697, 44.87406, 
44.83968, 44.81103, 44.81103, NA, 44.89697, 44.92562, 44.93708, 
44.93708, 44.91416, 44.89697, NA, 44.5532, 44.48444, 44.47871, 
NA, 44.5532, 44.63914, 44.63914, 44.72509, 44.72509, 44.81103, 
NA, 44.75373, 44.75946, 44.78238, 44.77665, 44.78238, 44.8053, 
44.86259, 44.89697, NA, 44.20369, 44.20369, NA, 43.85419, 43.85992, 
NA, 43.85992, 43.85992, 44.20369, NA, 43.85419, 44.20369, NA, 
45.76787, 45.95694, 45.95694, 46.11164, NA, 46.11164, 46.11164, 
NA, 45.76787, 45.76787, NA, 45.76787, 45.78506, NA, 46.11164, 
45.9283, 45.9283, 45.78506, NA, 43.51041, 43.85992, NA, 43.85992, 
43.50468, NA, 43.85992, 43.85992, NA, 43.51041, 43.85992, NA, 
43.50468, 43.84273, NA, 43.84273, 43.85992, 43.85992, NA, 43.85992, 
43.85992, NA, 43.85419, 43.51041, NA, 43.85419, 43.85992, NA, 
43.85992, 43.85992, NA, 44.20369, 44.20369, NA, 44.20369, 44.47871, 
NA, 44.20369, 44.38131, 44.37558, 44.46725, 44.4558, NA, 46.11164, 
46.11164, NA, 45.76787, 45.76787, NA, 46.03143, 45.93976, 45.93402, 
45.76787, NA, 46.11164, 46.03143, NA, 45.05167, 45.05167, 44.93708, 
44.89697, NA, 44.82249, 44.81676, 44.81676, 44.82249, 44.82249, 
44.81103, 44.79957, 44.81103, NA, 44.98865, 45.06886, 45.09178, 
45.10897, 45.1548, 45.16053, 45.17199, 45.18345, 45.20637, 45.21783, 
45.21783, 45.23502, 45.22929, 45.24648, NA, 43.85992, 43.85419, 
NA, 46.80492, 46.80492, NA, 45.74495, 45.74495, NA, 45.57306, 
45.74495, NA, 45.42982, 45.57306, NA, 47.84771, 47.84771, 47.89354, 
47.88208, NA, 47.88208, 47.71592, 47.71592, 47.36642, 47.36642, 
47.10859, 47.10859, 47.02264, NA, 43.50468, 43.85992, NA, 43.50468, 
43.85992, NA, 45.74495, 45.99705, 45.99705, 46.16321, NA, 45.74495, 
45.99705, 45.99705, 46.16894, NA, 45.33242, 45.24075, 45.23502, 
44.89697, NA, 45.41837, 45.41837, NA, 44.89697, 44.89697, NA, 
45.41837, 45.41837, 45.33242, NA, 45.41837, 45.24648, 45.24648, 
45.1548, NA, 48.54671, 48.54098, NA, 48.99935, 48.72433, 48.7186, 
48.54098, NA, 48.70141, 48.35764, NA, 48.62693, 48.06543, 48.0597, 
47.97375, 47.97375, 47.88208, NA, 45.17772, 45.16626, 45.16053, 
45.1548, NA, 44.81676, 44.81103, 44.89697, 44.89697, 44.93135, 
NA, 46.94816, 47.49247, 47.72739, 48.20294, NA, 48.53526, 48.71287, 
48.70714, 48.9306, 48.90768, 48.89622, 48.89049, 48.88476, 48.88476, 
48.89049, 48.89049, 48.91341, 48.92487, 48.94778, 48.9707, 48.99935, 
NA, 44.46725, 44.4558, 44.41569, 44.35839, 44.32975, 44.27818, 
44.25526, NA, 44.20369, 44.5532, NA, 44.5532, 44.5532, NA, 44.5532, 
44.54174, 44.53028, 44.4959, 44.48444, 44.46725, NA, 44.20369, 
44.20369, NA, 44.21515, 44.63914, NA, 44.21515, 44.20369, NA, 
44.63914, 44.63914, NA, 44.21515, 44.21515, NA, 44.21515, 44.5532, 
NA, 44.63914, 44.64487, 44.5532, 44.5532, NA, 47.50393, 47.15443, 
NA, 47.49247, 47.50393, NA, 48.16856, 48.17429, NA, 48.16856, 
48.16856, NA, 48.53526, 48.54098, NA, 43.85992, 43.85992, NA, 
44.98291, 44.98291, 44.89697, NA, 44.72509, 44.89697, NA, 44.72509, 
44.63914, 44.63914, 44.72509, 44.72509, NA, 44.98291, 44.98291, 
NA, 44.89697, 44.89697, NA, 45.33242, 45.33242, 45.29231, 45.29231, 
NA, 44.98291, 45.29231, NA, 45.83662, 45.99132, 45.99705, 46.16321, 
NA, 45.57306, 45.57306, NA, 45.78506, 45.78506, NA, 46.3351, 
46.3351, 46.3351, 45.78506, NA, 43.85992, 43.84273, 43.84273, 
NA, 43.85992, 43.85992, NA, 44.21515, 43.85992, NA, 44.21515, 
44.21515, NA, 44.47298, 44.47298, NA, 44.47298, 44.46725, NA, 
43.49895, 43.85992, NA, 47.50393, 47.4982, NA, 44.11775, 44.12348, 
44.20369, 44.20369, NA, 44.11775, 43.85992, NA, 46.11164, 46.28353, 
46.28353, 46.37521, NA, 46.71898, 46.37521, NA, 46.63303, 46.28353, 
46.28353, 46.11164, NA, 47.93365, 47.93365, NA, 48.16856, 48.03105, 
48.03105, 47.96803, NA, 47.93365, 47.9623, 47.96803, NA, 43.85419, 
43.85992, NA, 47.93365, 47.84771, 47.84771, 47.76176, 47.76176, 
47.76176, 47.85344, 47.85917, 47.96803, NA, 45.41837, 45.59025, 
45.59025, 45.76787, NA, 45.76787, 45.60171, 45.60171, 45.41837, 
NA, 45.41837, 45.41837, NA, 45.13189, 45.13189, 44.89697, 44.89697, 
NA, 44.50736, 44.53028, 44.58757, 44.58757, 44.5933, 44.59903, 
44.61049, 44.63341, 44.63914, 44.67925, 44.67925, 44.7079, 44.71363, 
NA, 44.71363, 44.5532, 44.5532, NA, 44.72509, 44.72509, 44.5532, 
44.5532, 44.47298, NA, 44.71363, 44.71936, 44.75373, 44.75946, 
44.77665, NA, 44.5532, 44.5532, NA, 44.20369, 44.20369, NA, 44.20369, 
44.20369, NA, 44.66206, 44.63341, 44.63914, 44.62769, 44.61049, 
44.58757, 44.57038, 44.5532, NA, 45.57306, 45.55587, 45.52723, 
45.49285, 45.45274, 45.44701, 45.41837, NA, 45.41837, 45.41263, 
45.40691, 45.38972, 45.39545, 45.38972, 45.36107, 45.34961, 45.34388, 
45.30378, 45.31523, 45.30378, 45.31523, 45.30378, 45.30378, 45.30378, 
45.27512, 45.25794, 45.24648, NA, 45.78506, 45.78506, NA, 45.29231, 
45.30378, 45.32096, 45.32669, 45.30378, 45.30378, 45.32669, 45.33815, 
45.33815, 45.34961, 45.38972, 45.41263, 45.41837, NA, 44.20369, 
43.85992, NA, 45.41837, 45.41837, NA, 45.59025, 45.76787, NA, 
46.37521, 46.37521, NA, 46.03143, 46.03143, NA, 44.11775, 44.20369, 
44.19223, NA), V4 = 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, 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, 3, 
3, 3, 3, 3, 3, 3, 3, 3, 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, 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, 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, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 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, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 
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, 10, 10, 10, 10, 10, 10, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 
12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 
12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 
13, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 
14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 
15, 15, 15, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 17, 17, 
17, 17, 17, 17, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 
18, 18, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 
19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 
19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 
20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 21, 21, 21, 
21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 
22, 22, 22, 22, 22, 22, 22, 22, 22, 23, 23, 23, 23, 23, 23, 23, 
23, 23, 23, 23, 23, 23, 24, 24, 24, 24, 24, 24, 24, 24, 24, 25, 
25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 26, 26, 26, 26, 26, 
26, 26, 26, 26, 26, 26, 26, 26, 26, 27, 27, 27, 27, 27, 27, 27, 
27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 
27, 27, 27, 27, 27, 27, 28, 28, 28, 29, 29, 29, 30, 30, 30, 30, 
30, 30, 30, 30, 30, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 
31, 31, 31, 32, 32, 32, 32, 32, 32, 33, 33, 33, 33, 33, 33, 33, 
33, 33, 33, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 
34, 34, 34, 34, 34, 34, 34, 35, 35, 35, 35, 35, 35, 35, 35, 36, 
36, 36, 36, 36, 36, 36, 36, 36, 36, 37, 37, 37, 37, 37, 37, 37, 
37, 37, 37, 37, 38, 38, 38, 38, 38, 39, 39, 39, 39, 39, 39, 39, 
39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 40, 40, 40, 40, 40, 40, 
40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 
40, 40, 41, 41, 41, 41, 41, 41, 41, 41, 41, 42, 42, 42, 42, 42, 
42, 42, 42, 42, 42, 42, 43, 43, 43, 43, 43, 43, 44, 44, 44, 44, 
44, 44, 44, 44, 44, 45, 45, 45, 46, 46, 46, 46, 46, 46, 46, 46, 
46, 46, 46, 46, 46, 46, 46, 46, 47, 47, 47, 47, 47, 47, 47, 47, 
47, 47, 47, 48, 48, 48, 48, 48, 48, 48, 48, 49, 49, 49, 49, 49, 
49, 49, 49, 50, 50, 50, 50, 51, 51, 51, 51, 51, 51, 51, 51, 51, 
52, 52, 52, 52, 52, 52, 53, 53, 53, 54, 54, 54, 55, 55, 55, 55, 
55, 55, 55, 55, 56, 56, 56, 56, 56, 56, 56, 56, 56, 56, 56, 56, 
56, 57, 57, 57, 57, 57, 57, 57, 57, 57, 57, 57, 57, 59, 59, 59, 
60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 61, 61, 61, 61, 61, 61, 
61, 61, 61, 61, 61, 61, 61, 62, 62, 62, 62, 62, 64, 64, 64, 64, 
64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 65, 65, 
65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 66, 66, 66, 66, 66, 66, 
66, 66, 66, 69, 69, 69, 69, 69, 69, 69, 69, 69, 70, 70, 70, 70, 
70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 
70, 70, 70, 70, 70, 70, 70, 70, 73, 73, 73, 73, 73, 73, 73, 73, 
73, 73, 73, 73, 73, 73, 73, 73, 73, 74, 74, 74, 75, 75, 75, 75, 
75, 75, 77, 77, 77, 78, 78, 78, 79, 79, 79, 79), V5 = c(9, 9, 
9, 9, 9, 11, 11, 11, 18, 18, 18, 18, 18, 31, 31, 31, 33, 33, 
33, 48, 48, 48, 48, 58, 58, 58, 72, 72, 72, 13, 13, 13, 27, 27, 
27, 27, 27, 27, 27, 27, 27, 30, 30, 30, 62, 62, 62, 62, 70, 70, 
70, 82, 82, 82, 14, 14, 14, 14, 14, 15, 15, 15, 29, 29, 29, 29, 
29, 43, 43, 43, 54, 54, 54, 56, 56, 56, 80, 80, 80, 11, 11, 11, 
11, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 29, 29, 
29, 31, 31, 31, 36, 36, 36, 36, 36, 39, 39, 39, 39, 39, 44, 44, 
44, 44, 44, 57, 57, 57, 68, 68, 68, 48, 48, 48, 49, 49, 49, 49, 
49, 70, 70, 70, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 37, 
37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 75, 75, 75, 75, 
76, 76, 76, 76, 76, 78, 78, 78, 8, 8, 8, 22, 22, 22, 40, 40, 
40, 40, 45, 45, 45, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 
52, 52, 52, 52, 52, 81, 81, 81, 83, 83, 83, 17, 17, 17, 17, 17, 
52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 
64, 64, 64, 64, 64, 65, 65, 65, 65, 65, 83, 83, 83, 58, 58, 58, 
72, 72, 72, 72, 27, 27, 27, 27, 27, 46, 46, 46, 69, 69, 69, 69, 
69, 69, 69, 69, 69, 69, 69, 69, 69, 71, 71, 71, 71, 71, 71, 86, 
86, 86, 18, 18, 18, 18, 18, 18, 18, 29, 29, 29, 29, 29, 29, 29, 
29, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 
31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 
31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 49, 49, 49, 49, 49, 
49, 49, 49, 49, 49, 49, 49, 49, 49, 77, 77, 77, 77, 77, 77, 77, 
77, 77, 77, 77, 77, 77, 80, 80, 80, 80, 80, 34, 34, 34, 37, 37, 
37, 37, 37, 37, 37, 37, 37, 37, 65, 65, 65, 65, 76, 76, 76, 87, 
87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 
30, 30, 30, 30, 30, 58, 58, 58, 82, 82, 82, 54, 54, 54, 56, 56, 
56, 56, 84, 84, 84, 29, 29, 29, 29, 29, 43, 43, 43, 57, 57, 57, 
60, 60, 60, 60, 60, 38, 38, 38, 38, 38, 32, 32, 32, 51, 51, 51, 
64, 64, 64, 83, 83, 83, 48, 48, 48, 49, 49, 49, 49, 49, 49, 49, 
49, 49, 49, 49, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 
25, 27, 27, 27, 27, 27, 27, 62, 62, 62, 62, 62, 62, 62, 66, 66, 
66, 66, 69, 69, 69, 69, 69, 69, 69, 82, 82, 82, 82, 82, 82, 82, 
82, 82, 25, 25, 25, 50, 50, 50, 55, 55, 55, 55, 74, 74, 74, 26, 
26, 26, 26, 26, 56, 56, 56, 61, 61, 61, 73, 73, 73, 77, 77, 77, 
77, 77, 24, 24, 24, 45, 45, 45, 81, 81, 81, 28, 28, 28, 50, 50, 
50, 55, 55, 55, 55, 85, 85, 85, 50, 50, 50, 74, 74, 74, 81, 81, 
81, 55, 55, 55, 66, 66, 66, 79, 79, 79, 79, 79, 79, 56, 56, 56, 
75, 75, 75, 78, 78, 78, 78, 78, 84, 84, 84, 62, 62, 62, 62, 62, 
69, 69, 69, 69, 69, 69, 69, 69, 69, 86, 86, 86, 86, 86, 86, 86, 
86, 86, 86, 86, 86, 86, 86, 86, 85, 85, 85, 80, 80, 80, 33, 33, 
33, 48, 48, 48, 70, 70, 70, 36, 36, 36, 36, 36, 72, 72, 72, 72, 
72, 72, 72, 72, 72, 45, 45, 45, 53, 53, 53, 48, 48, 48, 48, 48, 
58, 58, 58, 58, 58, 47, 47, 47, 47, 47, 61, 61, 61, 65, 65, 65, 
73, 73, 73, 73, 76, 76, 76, 76, 76, 44, 44, 44, 68, 68, 68, 68, 
68, 39, 39, 39, 72, 72, 72, 72, 72, 72, 72, 76, 76, 76, 76, 76, 
87, 87, 87, 87, 87, 87, 72, 72, 72, 72, 72, 68, 68, 68, 68, 68, 
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 52, 52, 52, 52, 
52, 52, 52, 52, 66, 66, 66, 69, 69, 69, 71, 71, 71, 71, 71, 71, 
71, 81, 81, 81, 42, 42, 42, 59, 59, 59, 87, 87, 87, 51, 51, 51, 
64, 64, 64, 87, 87, 87, 87, 87, 54, 54, 54, 60, 60, 60, 57, 57, 
57, 60, 60, 60, 68, 68, 68, 83, 83, 83, 47, 47, 47, 47, 65, 65, 
65, 71, 71, 71, 71, 71, 71, 86, 86, 86, 65, 65, 65, 73, 73, 73, 
73, 73, 86, 86, 86, 49, 49, 49, 49, 49, 70, 70, 70, 73, 73, 73, 
77, 77, 77, 77, 77, 55, 55, 55, 55, 53, 53, 53, 59, 59, 59, 64, 
64, 64, 65, 65, 65, 71, 71, 71, 67, 67, 67, 60, 60, 60, 79, 79, 
79, 79, 79, 85, 85, 85, 77, 77, 77, 77, 77, 80, 80, 80, 84, 84, 
84, 84, 84, 60, 60, 60, 60, 60, 60, 60, 60, 63, 63, 63, 63, 67, 
67, 67, 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, 73, 73, 73, 73, 
73, 75, 75, 75, 75, 75, 76, 76, 76, 82, 82, 82, 82, 82, 65, 65, 
65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 87, 87, 87, 87, 
71, 71, 71, 71, 71, 71, 87, 87, 87, 87, 87, 87, 69, 69, 69, 74, 
74, 74, 81, 81, 81, 71, 71, 71, 71, 71, 71, 71, 71, 71, 73, 73, 
73, 73, 73, 73, 73, 73, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 
86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 77, 77, 77, 86, 86, 86, 
86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 81, 81, 81, 76, 76, 
76, 78, 78, 78, 80, 80, 80, 84, 84, 84, 85, 85, 85, 85)), .Names = c("V1", 
"V2", "V3", "V4", "V5"), 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", "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", "381", "382", "383", "384", "385", "386", "387", "388", 
"389", "390", "391", "392", "393", "394", "395", "396", "397", 
"398", "399", "400", "401", "402", "403", "404", "405", "406", 
"407", "408", "409", "410", "411", "412", "413", "414", "415", 
"416", "417", "418", "419", "420", "421", "422", "423", "424", 
"425", "426", "427", "428", "429", "430", "431", "432", "433", 
"434", "435", "436", "437", "438", "439", "440", "441", "442", 
"443", "444", "445", "446", "447", "448", "449", "450", "451", 
"452", "453", "454", "455", "456", "457", "458", "459", "460", 
"461", "462", "463", "464", "465", "466", "467", "468", "469", 
"470", "471", "472", "473", "474", "475", "476", "477", "478", 
"479", "480", "481", "482", "483", "484", "485", "486", "487", 
"488", "489", "490", "491", "492", "493", "494", "495", "496", 
"497", "498", "499", "500", "501", "502", "503", "504", "505", 
"506", "507", "508", "509", "510", "511", "512", "513", "514", 
"515", "516", "517", "518", "519", "520", "521", "522", "523", 
"524", "525", "526", "527", "528", "529", "530", "531", "532", 
"533", "534", "535", "536", "537", "538", "539", "540", "541", 
"542", "543", "544", "545", "546", "547", "548", "549", "550", 
"551", "552", "553", "554", "555", "556", "557", "558", "559", 
"560", "561", "562", "563", "564", "565", "566", "567", "568", 
"569", "570", "571", "572", "573", "574", "575", "576", "577", 
"578", "579", "580", "581", "582", "583", "584", "585", "586", 
"587", "588", "589", "590", "591", "592", "593", "594", "595", 
"596", "597", "598", "599", "600", "601", "602", "603", "604", 
"605", "606", "607", "608", "609", "610", "611", "612", "613", 
"614", "615", "616", "617", "618", "619", "620", "621", "622", 
"623", "624", "625", "626", "627", "628", "629", "630", "631", 
"632", "633", "634", "635", "636", "637", "638", "639", "640", 
"641", "642", "643", "644", "645", "646", "647", "648", "649", 
"650", "651", "652", "653", "654", "655", "656", "657", "658", 
"659", "660", "661", "662", "663", "664", "665", "666", "667", 
"668", "669", "670", "671", "672", "673", "674", "675", "676", 
"677", "678", "679", "680", "681", "682", "683", "684", "685", 
"686", "687", "688", "689", "690", "691", "692", "693", "694", 
"695", "696", "697", "698", "699", "700", "701", "702", "703", 
"704", "705", "706", "707", "708", "709", "710", "711", "712", 
"713", "714", "715", "716", "717", "718", "719", "720", "721", 
"722", "723", "724", "725", "726", "727", "728", "729", "730", 
"731", "732", "733", "734", "735", "736", "737", "738", "739", 
"740", "741", "742", "743", "744", "745", "746", "747", "748", 
"749", "750", "751", "752", "753", "754", "755", "756", "757", 
"758", "759", "760", "761", "762", "763", "764", "765", "766", 
"767", "768", "769", "770", "771", "772", "773", "774", "775", 
"776", "777", "778", "779", "780", "781", "782", "783", "784", 
"785", "786", "787", "788", "789", "790", "791", "792", "793", 
"794", "795", "796", "797", "798", "799", "800", "801", "802", 
"803", "804", "805", "806", "807", "808", "809", "810", "811", 
"812", "813", "814", "815", "816", "817", "818", "819", "820", 
"821", "822", "823", "824", "825", "826", "827", "828", "829", 
"830", "831", "832", "833", "834", "835", "836", "837", "838", 
"839", "840", "841", "842", "843", "844", "845", "846", "847", 
"848", "849", "850", "851", "852", "853", "854", "855", "856", 
"857", "858", "859", "860", "861", "862", "863", "864", "865", 
"866", "867", "868", "869", "870", "871", "872", "873", "874", 
"875", "876", "877", "878", "879", "880", "881", "882", "883", 
"884", "885", "886", "887", "888", "889", "890", "891", "892", 
"893", "894", "895", "896", "897", "898", "899", "900", "901", 
"902", "903", "904", "905", "906", "907", "908", "909", "910", 
"911", "912", "913", "914", "915", "916", "917", "918", "919", 
"920", "921", "922", "923", "924", "925", "926", "927", "928", 
"929", "930", "931", "932", "933", "934", "935", "936", "937", 
"938", "939", "940", "941", "942", "943", "944", "945", "946", 
"947", "948", "949", "950", "951", "952", "953", "954", "955", 
"956", "957", "958", "959", "960", "961", "962", "963", "964", 
"965", "966", "967", "968", "969", "970", "971", "972", "973", 
"974", "975", "976", "977", "978", "979", "980", "981", "982", 
"983", "984", "985", "986", "987", "988", "989", "990", "991", 
"992", "993", "994", "995", "996", "997", "998", "999", "1000", 
"1001", "1002", "1003", "1004", "1005", "1006", "1007", "1008", 
"1009", "1010", "1011", "1012", "1013", "1014", "1015", "1016", 
"1017", "1018", "1019", "1020", "1021", "1022", "1023", "1024", 
"1025", "1026", "1027", "1028", "1029", "1030", "1031", "1032", 
"1033", "1034", "1035", "1036", "1037", "1038", "1039", "1040", 
"1041", "1042", "1043", "1044", "1045", "1046", "1047", "1048", 
"1049", "1050", "1051", "1052", "1053", "1054", "1055", "1056"
))


lines(edgelines[1:5,2:3],col=y.shad[1],lwd=5)
lines(edgelines[6:8,2:3],col=y.shad[2],lwd=5)
lines(edgelines[9:13,2:3],col=y.shad[3],lwd=5)
lines(edgelines[14:16,2:3],col=y.shad[4],lwd=5)
lines(edgelines[17:19,2:3],col=y.shad[5],lwd=5)
lines(edgelines[20:23,2:3],col=y.shad[6],lwd=5)
lines(edgelines[24:26,2:3],col=y.shad[7],lwd=5)
lines(edgelines[27:29,2:3],col=y.shad[8],lwd=5)
lines(edgelines[30:32,2:3],col=y.shad[9],lwd=5)
lines(edgelines[33:41,2:3],col=y.shad[10],lwd=5)
lines(edgelines[42:44,2:3],col=y.shad[11],lwd=5)
lines(edgelines[45:48,2:3],col=y.shad[12],lwd=5)
lines(edgelines[49:51,2:3],col=y.shad[13],lwd=5)
lines(edgelines[52:54,2:3],col=y.shad[14],lwd=5)
lines(edgelines[55:59,2:3],col=y.shad[15],lwd=5)
lines(edgelines[60:62,2:3],col=y.shad[16],lwd=5)
lines(edgelines[63:67,2:3],col=y.shad[17],lwd=5)
lines(edgelines[68:70,2:3],col=y.shad[18],lwd=5)
lines(edgelines[71:73,2:3],col=y.shad[19],lwd=5)
lines(edgelines[74:76,2:3],col=y.shad[20],lwd=5)
lines(edgelines[77:79,2:3],col=y.shad[21],lwd=5)
lines(edgelines[80:83,2:3],col=y.shad[22],lwd=5)
lines(edgelines[84:96,2:3],col=y.shad[23],lwd=5)
lines(edgelines[97:99,2:3],col=y.shad[24],lwd=5)
lines(edgelines[100:102,2:3],col=y.shad[25],lwd=5)
lines(edgelines[103:107,2:3],col=y.shad[26],lwd=5)
lines(edgelines[108:112,2:3],col=y.shad[27],lwd=5)
lines(edgelines[113:117,2:3],col=y.shad[28],lwd=5)
lines(edgelines[118:120,2:3],col=y.shad[29],lwd=5)
lines(edgelines[121:123,2:3],col=y.shad[30],lwd=5)
lines(edgelines[124:126,2:3],col=y.shad[31],lwd=5)
lines(edgelines[127:131,2:3],col=y.shad[32],lwd=5)
lines(edgelines[132:134,2:3],col=y.shad[33],lwd=5)
lines(edgelines[135:145,2:3],col=y.shad[34],lwd=5)
lines(edgelines[146:158,2:3],col=y.shad[35],lwd=5)
lines(edgelines[159:162,2:3],col=y.shad[36],lwd=5)
lines(edgelines[163:167,2:3],col=y.shad[37],lwd=5)
lines(edgelines[168:170,2:3],col=y.shad[38],lwd=5)
lines(edgelines[171:173,2:3],col=y.shad[39],lwd=5)
lines(edgelines[174:176,2:3],col=y.shad[40],lwd=5)
lines(edgelines[177:180,2:3],col=y.shad[41],lwd=5)
lines(edgelines[181:183,2:3],col=y.shad[42],lwd=5)
lines(edgelines[184:199,2:3],col=y.shad[43],lwd=5)
lines(edgelines[200:202,2:3],col=y.shad[44],lwd=5)
lines(edgelines[203:205,2:3],col=y.shad[45],lwd=5)
lines(edgelines[206:210,2:3],col=y.shad[46],lwd=5)
lines(edgelines[211:226,2:3],col=y.shad[47],lwd=5)
lines(edgelines[227:231,2:3],col=y.shad[48],lwd=5)
lines(edgelines[232:236,2:3],col=y.shad[49],lwd=5)
lines(edgelines[237:239,2:3],col=y.shad[50],lwd=5)
lines(edgelines[240:242,2:3],col=y.shad[51],lwd=5)
lines(edgelines[243:246,2:3],col=y.shad[52],lwd=5)
lines(edgelines[247:251,2:3],col=y.shad[53],lwd=5)
lines(edgelines[252:254,2:3],col=y.shad[54],lwd=5)
lines(edgelines[255:267,2:3],col=y.shad[55],lwd=5)
lines(edgelines[268:273,2:3],col=y.shad[56],lwd=5)
lines(edgelines[274:276,2:3],col=y.shad[57],lwd=5)
lines(edgelines[277:283,2:3],col=y.shad[58],lwd=5)
lines(edgelines[284:291,2:3],col=y.shad[59],lwd=5)
lines(edgelines[292:333,2:3],col=y.shad[60],lwd=5)
lines(edgelines[334:347,2:3],col=y.shad[61],lwd=5)
lines(edgelines[348:360,2:3],col=y.shad[62],lwd=5)
lines(edgelines[361:365,2:3],col=y.shad[63],lwd=5)
lines(edgelines[366:368,2:3],col=y.shad[64],lwd=5)
lines(edgelines[369:378,2:3],col=y.shad[65],lwd=5)
lines(edgelines[379:382,2:3],col=y.shad[66],lwd=5)
lines(edgelines[383:385,2:3],col=y.shad[67],lwd=5)
lines(edgelines[386:402,2:3],col=y.shad[68],lwd=5)
lines(edgelines[403:407,2:3],col=y.shad[69],lwd=5)
lines(edgelines[408:410,2:3],col=y.shad[70],lwd=5)
lines(edgelines[411:413,2:3],col=y.shad[71],lwd=5)
lines(edgelines[414:416,2:3],col=y.shad[72],lwd=5)
lines(edgelines[417:420,2:3],col=y.shad[73],lwd=5)
lines(edgelines[421:423,2:3],col=y.shad[74],lwd=5)
lines(edgelines[424:428,2:3],col=y.shad[75],lwd=5)
lines(edgelines[429:431,2:3],col=y.shad[76],lwd=5)
lines(edgelines[432:434,2:3],col=y.shad[77],lwd=5)
lines(edgelines[435:439,2:3],col=y.shad[78],lwd=5)
lines(edgelines[440:444,2:3],col=y.shad[79],lwd=5)
lines(edgelines[445:447,2:3],col=y.shad[80],lwd=5)
lines(edgelines[448:450,2:3],col=y.shad[81],lwd=5)
lines(edgelines[451:453,2:3],col=y.shad[82],lwd=5)
lines(edgelines[454:456,2:3],col=y.shad[83],lwd=5)
lines(edgelines[457:459,2:3],col=y.shad[84],lwd=5)
lines(edgelines[460:470,2:3],col=y.shad[85],lwd=5)
lines(edgelines[471:483,2:3],col=y.shad[86],lwd=5)
lines(edgelines[484:489,2:3],col=y.shad[87],lwd=5)
lines(edgelines[490:496,2:3],col=y.shad[88],lwd=5)
lines(edgelines[497:500,2:3],col=y.shad[89],lwd=5)
lines(edgelines[501:507,2:3],col=y.shad[90],lwd=5)
lines(edgelines[508:516,2:3],col=y.shad[91],lwd=5)
lines(edgelines[517:519,2:3],col=y.shad[92],lwd=5)
lines(edgelines[520:522,2:3],col=y.shad[93],lwd=5)
lines(edgelines[523:526,2:3],col=y.shad[94],lwd=5)
lines(edgelines[527:529,2:3],col=y.shad[95],lwd=5)
lines(edgelines[530:534,2:3],col=y.shad[96],lwd=5)
lines(edgelines[535:537,2:3],col=y.shad[97],lwd=5)
lines(edgelines[538:540,2:3],col=y.shad[98],lwd=5)
lines(edgelines[541:543,2:3],col=y.shad[99],lwd=5)
lines(edgelines[544:548,2:3],col=y.shad[100],lwd=5)
lines(edgelines[549:551,2:3],col=y.shad[101],lwd=5)
lines(edgelines[552:554,2:3],col=y.shad[102],lwd=5)
lines(edgelines[555:557,2:3],col=y.shad[103],lwd=5)
lines(edgelines[558:560,2:3],col=y.shad[104],lwd=5)
lines(edgelines[561:563,2:3],col=y.shad[105],lwd=5)
lines(edgelines[564:567,2:3],col=y.shad[106],lwd=5)
lines(edgelines[568:570,2:3],col=y.shad[107],lwd=5)
lines(edgelines[571:573,2:3],col=y.shad[108],lwd=5)
lines(edgelines[574:576,2:3],col=y.shad[109],lwd=5)
lines(edgelines[577:579,2:3],col=y.shad[110],lwd=5)
lines(edgelines[580:582,2:3],col=y.shad[111],lwd=5)
lines(edgelines[583:585,2:3],col=y.shad[112],lwd=5)
lines(edgelines[586:591,2:3],col=y.shad[113],lwd=5)
lines(edgelines[592:594,2:3],col=y.shad[114],lwd=5)
lines(edgelines[595:597,2:3],col=y.shad[115],lwd=5)
lines(edgelines[598:602,2:3],col=y.shad[116],lwd=5)
lines(edgelines[603:605,2:3],col=y.shad[117],lwd=5)
lines(edgelines[606:610,2:3],col=y.shad[118],lwd=5)
lines(edgelines[611:619,2:3],col=y.shad[119],lwd=5)
lines(edgelines[620:634,2:3],col=y.shad[120],lwd=5)
lines(edgelines[635:637,2:3],col=y.shad[121],lwd=5)
lines(edgelines[638:640,2:3],col=y.shad[122],lwd=5)
lines(edgelines[641:643,2:3],col=y.shad[123],lwd=5)
lines(edgelines[644:646,2:3],col=y.shad[124],lwd=5)
lines(edgelines[647:649,2:3],col=y.shad[125],lwd=5)
lines(edgelines[650:654,2:3],col=y.shad[126],lwd=5)
lines(edgelines[655:663,2:3],col=y.shad[127],lwd=5)
lines(edgelines[664:666,2:3],col=y.shad[128],lwd=5)
lines(edgelines[667:669,2:3],col=y.shad[129],lwd=5)
lines(edgelines[670:674,2:3],col=y.shad[130],lwd=5)
lines(edgelines[675:679,2:3],col=y.shad[131],lwd=5)
lines(edgelines[680:684,2:3],col=y.shad[132],lwd=5)
lines(edgelines[685:687,2:3],col=y.shad[133],lwd=5)
lines(edgelines[688:690,2:3],col=y.shad[134],lwd=5)
lines(edgelines[691:694,2:3],col=y.shad[135],lwd=5)
lines(edgelines[695:699,2:3],col=y.shad[136],lwd=5)
lines(edgelines[700:702,2:3],col=y.shad[137],lwd=5)
lines(edgelines[703:707,2:3],col=y.shad[138],lwd=5)
lines(edgelines[708:710,2:3],col=y.shad[139],lwd=5)
lines(edgelines[711:717,2:3],col=y.shad[140],lwd=5)
lines(edgelines[718:722,2:3],col=y.shad[141],lwd=5)
lines(edgelines[723:728,2:3],col=y.shad[142],lwd=5)
lines(edgelines[729:733,2:3],col=y.shad[143],lwd=5)
lines(edgelines[734:750,2:3],col=y.shad[144],lwd=5)
lines(edgelines[751:758,2:3],col=y.shad[145],lwd=5)
lines(edgelines[759:761,2:3],col=y.shad[146],lwd=5)
lines(edgelines[762:764,2:3],col=y.shad[147],lwd=5)
lines(edgelines[765:771,2:3],col=y.shad[148],lwd=5)
lines(edgelines[772:774,2:3],col=y.shad[149],lwd=5)
lines(edgelines[775:777,2:3],col=y.shad[150],lwd=5)
lines(edgelines[778:780,2:3],col=y.shad[151],lwd=5)
lines(edgelines[781:783,2:3],col=y.shad[152],lwd=5)
lines(edgelines[784:786,2:3],col=y.shad[153],lwd=5)
lines(edgelines[787:789,2:3],col=y.shad[154],lwd=5)
lines(edgelines[790:794,2:3],col=y.shad[155],lwd=5)
lines(edgelines[795:797,2:3],col=y.shad[156],lwd=5)
lines(edgelines[798:800,2:3],col=y.shad[157],lwd=5)
lines(edgelines[801:803,2:3],col=y.shad[158],lwd=5)
lines(edgelines[804:806,2:3],col=y.shad[159],lwd=5)
lines(edgelines[807:809,2:3],col=y.shad[160],lwd=5)
lines(edgelines[810:812,2:3],col=y.shad[161],lwd=5)
lines(edgelines[813:816,2:3],col=y.shad[162],lwd=5)
lines(edgelines[817:819,2:3],col=y.shad[163],lwd=5)
lines(edgelines[820:825,2:3],col=y.shad[164],lwd=5)
lines(edgelines[826:828,2:3],col=y.shad[165],lwd=5)
lines(edgelines[829:831,2:3],col=y.shad[166],lwd=5)
lines(edgelines[832:836,2:3],col=y.shad[167],lwd=5)
lines(edgelines[837:839,2:3],col=y.shad[168],lwd=5)
lines(edgelines[840:844,2:3],col=y.shad[169],lwd=5)
lines(edgelines[845:847,2:3],col=y.shad[170],lwd=5)
lines(edgelines[848:850,2:3],col=y.shad[171],lwd=5)
lines(edgelines[851:855,2:3],col=y.shad[172],lwd=5)
lines(edgelines[856:859,2:3],col=y.shad[173],lwd=5)
lines(edgelines[860:862,2:3],col=y.shad[174],lwd=5)
lines(edgelines[863:865,2:3],col=y.shad[175],lwd=5)
lines(edgelines[866:868,2:3],col=y.shad[176],lwd=5)
lines(edgelines[869:871,2:3],col=y.shad[177],lwd=5)
lines(edgelines[872:874,2:3],col=y.shad[178],lwd=5)
lines(edgelines[875:877,2:3],col=y.shad[179],lwd=5)
lines(edgelines[878:880,2:3],col=y.shad[180],lwd=5)
lines(edgelines[881:885,2:3],col=y.shad[181],lwd=5)
lines(edgelines[886:888,2:3],col=y.shad[182],lwd=5)
lines(edgelines[889:893,2:3],col=y.shad[183],lwd=5)
lines(edgelines[894:896,2:3],col=y.shad[184],lwd=5)
lines(edgelines[897:901,2:3],col=y.shad[185],lwd=5)
lines(edgelines[902:909,2:3],col=y.shad[186],lwd=5)
lines(edgelines[910:913,2:3],col=y.shad[187],lwd=5)
lines(edgelines[914:916,2:3],col=y.shad[188],lwd=5)
lines(edgelines[917:926,2:3],col=y.shad[189],lwd=5)
lines(edgelines[927:931,2:3],col=y.shad[190],lwd=5)
lines(edgelines[932:936,2:3],col=y.shad[191],lwd=5)
lines(edgelines[937:939,2:3],col=y.shad[192],lwd=5)
lines(edgelines[940:944,2:3],col=y.shad[193],lwd=5)
lines(edgelines[945:958,2:3],col=y.shad[194],lwd=5)
lines(edgelines[959:962,2:3],col=y.shad[195],lwd=5)
lines(edgelines[963:968,2:3],col=y.shad[196],lwd=5)
lines(edgelines[969:974,2:3],col=y.shad[197],lwd=5)
lines(edgelines[975:977,2:3],col=y.shad[198],lwd=5)
lines(edgelines[978:980,2:3],col=y.shad[199],lwd=5)
lines(edgelines[981:983,2:3],col=y.shad[200],lwd=5)
lines(edgelines[984:992,2:3],col=y.shad[201],lwd=5)
lines(edgelines[993:1000,2:3],col=y.shad[202],lwd=5)
lines(edgelines[1001:1020,2:3],col=y.shad[203],lwd=5)
lines(edgelines[1021:1023,2:3],col=y.shad[204],lwd=5)
lines(edgelines[1024:1037,2:3],col=y.shad[205],lwd=5)
lines(edgelines[1038:1040,2:3],col=y.shad[206],lwd=5)
lines(edgelines[1041:1043,2:3],col=y.shad[207],lwd=5)
lines(edgelines[1044:1046,2:3],col=y.shad[208],lwd=5)
lines(edgelines[1047:1049,2:3],col=y.shad[209],lwd=5)
lines(edgelines[1050:1052,2:3],col=y.shad[210],lwd=5)
lines(edgelines[1053:1056,2:3],col=y.shad[211],lwd=5)
#par(old.par)
}





#########################################################################
##%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%##
##%%                                                                 %%##
##%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%##
#########################################################################

# data should be prob(boundary), i.e. 1-wij/1-pij
# plotting order is absolutely alphabetic in county and in edges (R is not alphabetic in edge order)


mn.Edge.crisp.rev<-function(figtitle,y,qtle=T,breaks){
#old.par <- par(no.readonly = TRUE)
library(maps)
map("county","minnesota", fill=F, plot=T, interior=F)
if(qtle==T) { br <- round(quantile(y, probs=breaks),2)
                n.col<- length(breaks)-1

    # 0: dark 1: light
    shading<-gray(0:(n.col-1)/n.col)
    y.grp<-findInterval(y, vec=br, rightmost.closed = TRUE, all.inside = TRUE)
    y.shad<-shading[y.grp]

   #leg.txt<-paste("[",br[n.col],",",br[n.col+1],"]",sep="")
   #for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste("[",br[i],",",br[i+1],,sep=""),) }
   leg.txt<-paste(breaks[n.col]*100,"%~",breaks[n.col+1]*100,"%",sep="")
   for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste( breaks[i]*100,"%~",breaks[i+1]*100,"%", sep=""),) }
   }

if(qtle==F) { br <- round(breaks,2)
             n.col<- length(breaks)-1

    # 0: dark 1: light
    shading<-gray(0:(n.col-1)/n.col)
    y.grp<-findInterval(y, vec=br, rightmost.closed = TRUE, all.inside = TRUE)
    y.shad<-shading[y.grp]

   #leg.txt<-paste("[",br[n.col],",",br[n.col+1],"]",sep="")
   #for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste("[",br[i],",",br[i+1],,sep=""),) }
   leg.txt<-paste("[",br[n.col],",",br[n.col+1],"]",sep="")
   for(i in (n.col-1):1){    leg.txt<-append(leg.txt,paste("[",br[i],",",br[i+1],")",sep=""),) }
   }

   leg.txt<-rev(leg.txt)
   leg.pos<-matrix(c(-92.25,46,-92,45.5), byrow=T, nrow=2)
   legend(leg.pos,legend=leg.txt,col=shading,lty="solid",lwd=5,cex=1,ncol=1,bty="n")

#title(sub=figtitle,cex.sub=2)
title(main=figtitle,cex.main=2)

edgelines<-structure(list(V1 = c(4, 5, 6, 7, 8, 24, 25, 26, 19, 20, 21, 
22, 23, 27, 28, 29, 12, 13, 14, 15, 16, 17, 18, 9, 10, 11, 1, 
2, 3, 36, 37, 38, 46, 47, 48, 49, 50, 51, 52, 53, 54, 33, 34, 
35, 42, 43, 44, 45, 30, 31, 32, 39, 40, 41, 58, 59, 60, 61, 62, 
69, 70, 71, 72, 73, 74, 75, 76, 66, 67, 68, 63, 64, 65, 55, 56, 
57, 77, 78, 79, 96, 97, 98, 99, 103, 104, 105, 106, 107, 108, 
109, 110, 111, 112, 113, 114, 115, 100, 101, 102, 93, 94, 95, 
88, 89, 90, 91, 92, 83, 84, 85, 86, 87, 119, 120, 121, 122, 123, 
116, 117, 118, 80, 81, 82, 124, 125, 126, 141, 142, 143, 144, 
145, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 
139, 140, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 
162, 163, 181, 182, 183, 184, 146, 147, 148, 149, 150, 178, 179, 
180, 188, 189, 190, 214, 215, 216, 207, 208, 209, 210, 217, 218, 
219, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 
203, 204, 205, 206, 211, 212, 213, 185, 186, 187, 239, 240, 241, 
242, 243, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 
231, 232, 233, 234, 235, 244, 245, 246, 247, 248, 249, 250, 251, 
252, 253, 236, 237, 238, 261, 262, 263, 254, 255, 256, 257, 289, 
290, 291, 292, 293, 283, 284, 285, 264, 265, 266, 267, 268, 269, 
270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 
286, 287, 288, 336, 337, 338, 339, 340, 341, 342, 375, 376, 377, 
378, 379, 380, 381, 382, 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, 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, 383, 384, 385, 407, 408, 409, 410, 411, 412, 413, 414, 415, 
416, 386, 387, 388, 389, 417, 418, 419, 390, 391, 392, 393, 394, 
395, 396, 397, 398, 399, 400, 401, 402, 403, 404, 405, 406, 435, 
436, 437, 438, 439, 440, 441, 442, 432, 433, 434, 467, 468, 469, 
447, 448, 449, 450, 451, 452, 453, 473, 474, 475, 476, 477, 478, 
479, 480, 470, 471, 472, 481, 482, 483, 484, 485, 534, 535, 536, 
537, 538, 542, 543, 544, 545, 546, 547, 548, 549, 550, 539, 540, 
541, 551, 552, 553, 554, 555, 556, 557, 558, 559, 560, 561, 562, 
563, 564, 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 611, 
612, 613, 576, 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 
587, 588, 565, 566, 567, 568, 569, 570, 571, 572, 573, 574, 575, 
589, 590, 591, 592, 593, 594, 595, 596, 597, 617, 618, 619, 624, 
625, 626, 620, 621, 622, 623, 614, 615, 616, 630, 631, 632, 633, 
634, 635, 636, 637, 627, 628, 629, 643, 644, 645, 638, 639, 640, 
641, 642, 658, 659, 660, 652, 653, 654, 655, 656, 657, 677, 678, 
679, 667, 668, 669, 670, 671, 672, 673, 674, 675, 676, 692, 693, 
694, 689, 690, 691, 686, 687, 688, 718, 719, 720, 695, 696, 697, 
712, 713, 714, 715, 716, 717, 732, 733, 734, 721, 722, 723, 724, 
725, 726, 727, 728, 729, 730, 731, 759, 760, 761, 762, 763, 735, 
736, 737, 738, 739, 740, 741, 742, 743, 744, 745, 746, 747, 748, 
749, 750, 751, 752, 753, 754, 755, 756, 757, 758, 770, 771, 772, 
788, 789, 790, 797, 798, 799, 794, 795, 796, 791, 792, 793, 800, 
801, 802, 803, 804, 805, 806, 807, 808, 809, 810, 811, 812, 813, 
826, 827, 828, 823, 824, 825, 829, 830, 831, 832, 833, 834, 835, 
836, 837, 838, 843, 844, 845, 846, 847, 856, 857, 858, 848, 849, 
850, 839, 840, 841, 842, 851, 852, 853, 854, 855, 882, 883, 884, 
877, 878, 879, 880, 881, 885, 886, 887, 908, 909, 910, 911, 912, 
913, 914, 915, 916, 917, 918, 919, 920, 921, 922, 923, 924, 925, 
942, 943, 944, 945, 946, 967, 968, 969, 970, 971, 972, 973, 974, 
975, 976, 977, 978, 979, 980, 981, 982, 983, 1022, 1023, 1024, 
1025, 1026, 1027, 1028, 1029, 1016, 1017, 1018, 1013, 1014, 1015, 
1030, 1031, 1032, 1033, 1034, 1035, 1036, 1019, 1020, 1021, 1046, 
1047, 1048, 1049, 1050, 1051, 1043, 1044, 1045, 1055, 1056, 1057, 
1052, 1053, 1054, 1058, 1059, 1060, 1061, 1062, 1079, 1080, 1081, 
1082, 1083, 1084, 1088, 1089, 1090, 1091, 1092, 1093, 1085, 1086, 
1087, 1110, 1111, 1112, 1066, 1067, 1068, 1069, 1063, 1064, 1065, 
1073, 1074, 1075, 1076, 1077, 1078, 1070, 1071, 1072, 1119, 1120, 
1121, 1122, 1123, 1124, 1125, 1126, 1116, 1117, 1118, 1127, 1128, 
1129, 1130, 1131, 1132, 1133, 1134, 1135, 1136, 1137, 1138, 1139, 
1140, 1141, 1142, 1143, 1144, 1145, 1146, 1159, 1160, 1161, 1162, 
1163, 1164, 1156, 1157, 1158, 1165, 1166, 1167, 1168, 1169, 1170, 
1174, 1175, 1176, 1193, 1194, 1195, 1196, 1197, 1198, 1199, 1200, 
1201, 1202, 1203, 1212, 1213, 1214, 1215, 1216, 1209, 1210, 1211, 
1204, 1205, 1206, 1207, 1208, 1217, 1218, 1219, 1224, 1225, 1226, 
1227, 1228, 1220, 1221, 1222, 1223, 1248, 1249, 1250, 1278, 1279, 
1280, 1281, 1282, 1283, 1284, 1285, 1286, 1287, 1296, 1297, 1298, 
1299, 1300, 1291, 1292, 1293, 1294, 1295, 1288, 1289, 1290, 1301, 
1302, 1303, 1304, 1305, 1310, 1311, 1312, 1313, 1314, 1315, 1316, 
1317, 1318, 1319, 1320, 1321, 1322, 1323, 1306, 1307, 1308, 1309, 
1324, 1325, 1326, 1327, 1328, 1329, 1330, 1331, 1332, 1333, 1334, 
1335, 1342, 1343, 1344, 1336, 1337, 1338, 1339, 1340, 1341, 1411, 
1412, 1413, 1414, 1415, 1416, 1417, 1418, 1419, 1440, 1441, 1442, 
1443, 1444, 1445, 1446, 1447, 1420, 1421, 1422, 1423, 1424, 1425, 
1426, 1427, 1428, 1429, 1430, 1431, 1432, 1433, 1434, 1435, 1436, 
1437, 1438, 1439, 1448, 1449, 1450, 1451, 1452, 1453, 1454, 1455, 
1456, 1457, 1458, 1459, 1460, 1461, 1462, 1463, 1464, 1465, 1466, 
1467, 1468, 1469, 1470, 1471, 1472, 1473, 1474, 1475, 1476, 1495, 
1496, 1497, 1498, 1499, 1500, 1501), V2 = c(-93.04834, -93.05408, 
-93.03689, -93.03689, NA, -93.77027, -93.77027, NA, -93.81038, 
-93.8161, -93.78174, -93.77027, NA, -93.03689, -93.77027, NA, 
-93.03689, -93.43223, NA, -93.43223, -93.43223, -93.81038, NA, 
-93.03689, -93.03689, NA, -93.04834, -93.03689, NA, -93.0197, 
-93.01397, NA, -93.23169, -93.29472, -93.289, -93.289, -93.32336, 
-93.3864, -93.48952, -93.51817, NA, -93.51817, -93.0197, NA, 
-93.01397, -93.23169, -93.23169, NA, -93.51817, -93.51817, NA, 
-93.01397, -93.01397, NA, -96.17669, -96.17669, -96.15378, -96.14805, 
NA, -95.16829, -95.54645, NA, -95.14537, -95.13965, -95.16256, 
-95.16829, NA, -96.05064, -95.54645, NA, -96.05064, -96.17669, 
NA, -95.14537, -96.14805, NA, -95.14537, -95.14537, NA, -94.38907, 
-94.3948, -94.66409, NA, -95.56936, -95.20267, -95.21413, -95.25423, 
-95.2657, -95.2657, -95.25423, -95.23132, -95.17975, -95.18548, 
-95.16829, -95.16829, NA, -94.66409, -95.16829, NA, -94.38907, 
-94.3948, NA, -94.3948, -94.40053, -94.3948, -94.3948, NA, -95.32299, 
-95.19694, -95.19694, -94.3948, NA, -95.58655, -95.58655, -95.57509, 
-95.57509, NA, -95.56936, -95.57509, NA, -95.58655, -95.32299, 
NA, -93.77027, -93.77027, NA, -93.77027, -94.35469, -94.32604, 
-94.28593, NA, -93.77027, -94.15988, NA, -94.15988, -94.1828, 
-94.20572, -94.19999, -94.22864, -94.2229, -94.2229, -94.25729, 
-94.26302, -94.28593, NA, -96.09075, -96.12513, -96.1595, -96.16524, 
-96.19389, -96.22826, -96.26836, -96.27983, -96.2741, -96.30274, 
-96.36004, -96.43452, NA, -96.10794, -96.234, -96.23972, NA, 
-96.09075, -96.08502, -96.10794, -96.10794, NA, -96.81268, -96.23972, 
NA, -94.37188, -94.37188, NA, -94.26302, -93.776, NA, -93.776, 
-93.776, -94.00519, NA, -94.37188, -94.26302, NA, -94.37188, 
-94.32604, -94.30312, -94.2401, -94.22864, -94.1828, -94.16561, 
-94.10832, -94.09113, -94.03383, -94.01665, -94.00519, -94.01092, 
-94.01092, -94.00519, NA, -93.776, -93.776, NA, -94.37188, -94.37188, 
NA, -94.86462, -95.09954, -95.09954, -95.11672, NA, -94.37188, 
-94.38907, -94.41199, -94.4349, -94.49792, -94.5094, -94.52657, 
-94.58387, -94.65263, -94.68127, -94.69846, -94.7042, -94.76149, 
-94.7844, -94.79014, NA, -94.85889, -94.87035, -95.11672, -95.11672, 
NA, -94.85889, -94.82452, -94.79587, -94.79014, NA, -94.37188, 
-94.86462, NA, -93.03689, -92.29205, NA, -93.04834, -92.29778, 
-92.3035, NA, -93.51817, -93.51817, -93.76454, -93.76454, NA, 
-94.01665, -94.01092, NA, -93.76454, -93.70152, -93.68433, -93.6614, 
-93.63276, -93.61558, -93.62703, -93.60985, -93.60985, -93.58692, 
-93.57547, -93.51817, NA, -93.76454, -93.76454, -93.87914, -93.88486, 
-94.01092, NA, -94.01665, -93.76454, NA, -94.34896, -94.32604, 
-94.31458, -94.33177, -94.33177, -93.77027, NA, -94.66409, -94.66982, 
-94.65836, -94.65836, -94.6469, -94.6469, -94.77868, NA, -94.38907, 
-94.36615, -94.3375, -94.31458, -94.29166, -94.28593, -94.26302, 
-94.25729, -94.24583, -94.24583, -94.2401, -94.22864, -94.21145, 
-94.19999, -94.17707, -94.15416, -94.10832, -94.07394, -94.04529, 
-94.0281, -93.99372, -93.99372, -94.01092, -94.02238, -94.02238, 
-94.01092, -93.95361, -93.91924, -93.89632, -93.88486, -93.86768, 
-93.83903, -93.81038, -93.78174, -93.77027, -93.776, -93.79319, 
-93.79892, -93.78746, -93.77027, -93.77027, NA, -94.34896, -94.35469, 
-94.35469, -94.38334, -94.40053, -94.41199, -94.4349, -94.45782, 
-94.4922, -94.5323, -94.55523, -94.57241, -94.6068, NA, -94.6068, 
-94.6297, -94.64117, -94.65836, -94.67554, -94.68701, -94.69273, 
-94.70992, -94.73284, -94.73857, -94.76149, -94.77295, NA, -94.77295, 
-94.76722, -94.77868, -94.77868, NA, -95.24277, -95.24277, NA, 
-96.03918, -95.97043, -95.9246, -95.89594, -95.83865, -95.81573, 
-95.80427, -95.76989, -95.74698, NA, -95.24277, -95.48341, -95.48914, 
NA, -96.03918, -95.24277, NA, -95.48914, -95.51207, -95.52925, 
-95.52352, -95.50061, -95.55218, -95.57509, -95.58082, -95.62093, 
-95.63811, -95.64384, -95.67822, -95.68395, -95.71833, -95.71833, 
-95.74698, NA, -93.14001, -93.14001, -93.01397, -93.0197, NA, 
-93.14001, -92.8478, NA, -92.76186, -93.01397, NA, -96.81841, 
-96.17669, NA, -96.2569, -96.15378, -96.14805, NA, -96.2569, 
-96.7783, NA, -95.16829, -95.17402, -95.16256, -95.16829, NA, 
-95.54645, -95.54645, NA, -95.56936, -95.56936, NA, -95.54645, 
-95.54645, -95.56936, -95.56936, NA, -91.00288, -91.00288, -91.02581, 
-91.01434, NA, -95.47768, -94.86462, NA, -95.47768, -95.46623, 
NA, -95.11672, -95.46623, NA, -94.86462, -94.86462, NA, -93.81038, 
-93.81038, NA, -93.81038, -94.40053, -94.41772, -94.41199, -94.38907, 
-94.3948, -94.3948, -94.3948, -94.37188, -94.34896, NA, -93.03689, 
-93.03689, -93.01397, -92.98532, -92.97386, -92.94521, -92.91656, 
-92.91656, -92.79625, -92.80197, -92.73322, -92.73895, NA, -93.1744, 
-93.19732, -93.24316, -93.28326, -93.3291, NA, -93.01397, -93.02543, 
-93.03689, -93.12283, -93.14001, -93.1744, NA, -93.28326, -93.2718, 
-93.03689, NA, -93.28326, -93.2718, -93.31763, -93.31763, -93.3291, 
-93.3291, NA, -92.80197, -92.84208, -92.89938, -92.91656, -93.01397, 
-93.01397, -93.00824, -93.01397, NA, -93.04834, -92.6702, NA, 
-93.05408, -92.69312, NA, -92.69312, -92.67592, -92.6702, NA, 
-93.05408, -93.04834, NA, -95.75844, -95.75844, -95.76989, -95.76989, 
NA, -95.14537, -95.76989, NA, -95.75844, -95.13965, NA, -95.13965, 
-95.13391, NA, -95.14537, -95.1511, -95.13391, -95.13391, NA, 
-93.64423, -93.64996, NA, -94.26302, -94.26302, NA, -93.776, 
-93.64996, NA, -91.72482, -91.71909, NA, -92.4582, -92.4582, 
NA, -92.4582, -92.45247, -92.08578, NA, -92.08578, -91.71909, 
NA, -93.05408, -93.0598, NA, -93.05408, -93.41504, NA, -93.64996, 
-93.41504, NA, -92.5556, -92.6702, NA, -93.04834, -93.03689, 
NA, -92.5556, -92.5556, -92.42383, -92.42955, -92.25194, NA, 
-95.76989, -96.25118, NA, -95.75844, -96.23972, NA, -96.25118, 
-96.2569, -96.23972, -96.23972, NA, -96.25118, -96.25118, NA, 
-93.23169, -93.20877, -93.20305, -93.1744, NA, -93.51817, -93.4838, 
-93.45515, -93.45515, -93.42077, -93.3864, -93.36347, -93.3291, 
NA, -93.76454, -93.76454, -93.75881, -93.72443, -93.70152, -93.68433, 
-93.64996, -93.64996, -93.6614, -93.64423, -93.60412, -93.58692, 
-93.53536, -93.51817, NA, -91.71909, -91.28937, NA, -95.14537, 
-94.77868, NA, -93.51817, -93.14001, NA, -93.51817, -93.51817, 
NA, -93.51817, -93.51817, NA, -94.3948, -93.75308, -93.75308, 
-93.0598, NA, -93.0598, -93.06554, -93.05408, -93.0598, -93.04834, 
-93.05408, -93.03689, -93.03689, NA, -94.86462, -94.86462, NA, 
-95.46623, -95.47768, NA, -93.51817, -93.51817, -93.43223, -93.43223, 
NA, -93.14001, -93.14001, -93.04834, -93.03689, NA, -94.76149, 
-94.75576, -94.75003, -94.75003, NA, -95.12818, -95.25996, NA, 
-95.24277, -94.75003, NA, -95.12818, -94.76149, -94.76149, NA, 
-95.25996, -95.25423, -95.24277, -95.24277, NA, -97.145, -96.34858, 
NA, -96.3715, -96.37724, -96.36004, -96.34858, NA, -94.40053, 
-94.3948, NA, -93.07127, -93.09418, -93.07127, -93.077, -93.06554, 
-93.0598, NA, -96.09075, -96.07356, -96.05637, -96.03918, NA, 
-96.43452, -95.86156, -95.86156, -95.74698, -95.74698, NA, -91.78784, 
-91.7821, -91.79357, -91.77065, NA, -95.32299, -95.32872, -95.08234, 
-95.07661, -95.09954, -95.12245, -95.14537, -95.17402, -95.2084, 
-95.2485, -95.27715, -95.28288, -95.29434, -95.28861, -95.31152, 
-95.29434, NA, -93.93643, -93.93643, -93.95935, -93.93643, -93.9479, 
-93.99372, -94.00519, NA, -93.5239, -93.5239, NA, -93.89632, 
-93.5239, NA, -93.89632, -93.89632, -93.89632, -93.90205, -93.92497, 
-93.93643, NA, -93.5239, -93.776, NA, -96.05064, -96.06783, NA, 
-96.05064, -96.4288, NA, -96.43452, -96.06783, NA, -95.59229, 
-96.05064, NA, -95.59229, -95.598, NA, -96.06783, -95.60947, 
-95.60947, -95.598, NA, -96.05064, -96.05064, NA, -95.54645, 
-96.05064, NA, -96.4689, -95.57509, NA, -96.4689, -97.145, NA, 
-95.58655, -96.34858, NA, -94.37188, -94.86462, NA, -94.25156, 
-94.50366, -94.49792, NA, -94.49792, -94.49792, NA, -94.49792, 
-94.4922, -94.25156, -94.25156, -94.01092, NA, -94.01665, -94.25156, 
NA, -94.49792, -94.75003, NA, -94.76149, -94.3776, -94.38334, 
-94.26302, NA, -94.25156, -94.26302, NA, -93.77027, -93.776, 
-93.81038, -93.81038, NA, -93.51817, -93.77027, NA, -94.63544, 
-94.28593, NA, -94.6068, -94.61825, -94.65263, -94.63544, NA, 
-92.69312, -92.68738, -92.4582, NA, -96.05064, -95.47768, NA, 
-96.05064, -96.05064, NA, -95.59229, -95.46623, NA, -94.79014, 
-94.61825, NA, -94.61825, -93.93643, NA, -96.04491, -96.05064, 
NA, -96.05064, -96.8356, NA, -92.08005, -92.32069, -92.32069, 
-92.5556, NA, -92.08005, -92.08578, NA, -95.14537, -95.13965, 
-95.1511, -95.14537, NA, -95.14537, -95.14537, NA, -96.2569, 
-96.2741, -96.2569, -96.25118, NA, -95.56936, -95.69541, NA, 
-96.4689, -96.46317, -96.45745, -96.45745, NA, -95.69541, -95.69541, 
-96.45745, NA, -96.4288, -96.05064, NA, -95.69541, -95.70114, 
-95.82146, -95.82146, -96.0793, -96.3314, -96.3314, -96.45745, 
-96.45745, NA, -95.12818, -95.12818, -95.13391, -95.13965, NA, 
-95.75844, -95.75844, -95.74698, -95.74698, NA, -95.25996, -95.74698, 
NA, -93.01397, -92.9796, -92.9796, -93.01397, NA, -94.85889, 
-94.899, -95.0537, -95.07661, -95.08807, -95.08234, -95.13965, 
-95.14537, -95.18548, -95.23132, -95.28288, -95.31152, -95.35163, 
NA, -95.35163, -95.35736, -95.598, NA, -94.49792, -94.62398, 
-94.6297, -94.61252, -94.61825, NA, -95.35163, -95.36883, -95.40894, 
-95.44904, -95.48914, NA, -93.28326, -93.5239, NA, -93.04834, 
-93.40931, NA, -93.40931, -93.5239, NA, -93.76454, -93.76454, 
-93.78746, -93.8333, -93.89632, -93.8906, -93.90205, -93.89632, 
NA, -94.15988, -94.14843, -94.14843, -94.11977, -94.1026, -94.05103, 
-94.04529, NA, -94.04529, -94.03383, -93.988, -93.96508, -93.93643, 
-93.91924, -93.88486, -93.82757, -93.8161, -93.73016, -93.70152, 
-93.67287, -93.63849, -93.60985, -93.57547, -93.56401, -93.55255, 
-93.53536, -93.51817, NA, -95.13391, -94.63544, NA, -94.26302, 
-94.25729, -94.2229, -94.20572, -94.1828, -94.15988, -94.16561, 
-94.15416, -94.11977, -94.1026, -94.0854, -94.06821, -94.04529, 
NA, -93.40931, -93.41504, NA, -95.74698, -96.10794, NA, -96.23972, 
-96.23972, NA, -94.77295, -95.14537, NA, -96.25118, -96.54912, 
NA, -92.08005, -92.08005, -91.83368, NA), V3 = c(46.76482, 46.5872, 
46.57574, 46.42677, NA, 46.79919, 47.01691, NA, 46.25488, 46.59293, 
46.5872, 46.79919, NA, 47.02264, 47.01691, NA, 46.16894, 46.16321, 
NA, 46.16321, 46.26061, 46.25488, NA, 46.42677, 46.16894, NA, 
46.76482, 47.02264, NA, 45.41837, 45.30378, NA, 45.05167, 45.0574, 
45.09751, 45.13189, 45.1548, 45.20637, 45.23502, 45.24648, NA, 
45.42982, 45.41837, NA, 45.13189, 45.13189, 45.05167, NA, 45.42982, 
45.24648, NA, 45.13189, 45.30378, NA, 47.15443, 46.98254, 46.98254, 
46.71325, NA, 47.1487, 47.15443, NA, 46.80492, 46.98254, 46.98254, 
47.1487, NA, 47.15443, 47.15443, NA, 47.15443, 47.15443, NA, 
46.71898, 46.71325, NA, 46.71898, 46.80492, NA, 47.42372, 47.40653, 
47.40653, NA, 48.02532, 48.02532, 48.00813, 47.97375, 47.93938, 
47.91073, 47.88781, 47.87635, 47.87635, 47.67582, 47.67582, 47.4008, 
NA, 47.40653, 47.4008, NA, 47.42372, 47.84771, NA, 48.35764, 
48.02532, 48.01959, 47.84771, NA, 48.53526, 48.53526, 48.3691, 
48.35764, NA, 48.53526, 48.35764, 48.35191, 48.17429, NA, 48.02532, 
48.17429, NA, 48.53526, 48.53526, NA, 45.83662, 45.57306, NA, 
45.83662, 45.83662, 45.80798, 45.78506, NA, 45.57306, 45.57306, 
NA, 45.57306, 45.60171, 45.63609, 45.68765, 45.72776, 45.73922, 
45.74495, 45.74495, 45.7736, 45.78506, NA, 45.17772, 45.18918, 
45.21783, 45.2121, 45.21783, 45.22356, 45.22356, 45.22929, 45.24648, 
45.24648, 45.25794, 45.26367, NA, 45.41837, 45.41837, 45.59025, 
NA, 45.17772, 45.24648, 45.24075, 45.41837, NA, 45.58452, 45.59025, 
NA, 44.28391, 44.11775, NA, 43.85992, 43.85992, NA, 44.20369, 
44.25526, 44.25526, NA, 43.85992, 43.85992, NA, 44.28391, 44.27245, 
44.26099, 44.2438, 44.22661, 44.21515, 44.20369, 44.20369, 44.18077, 
44.17505, 44.18077, 44.20369, 44.22661, 44.23807, 44.25526, NA, 
43.85992, 44.20369, NA, 43.85992, 44.11775, NA, 44.11775, 44.12348, 
44.20369, 44.20369, NA, 44.28391, 44.28391, 44.30109, 44.3412, 
44.36985, 44.38131, 44.38131, 44.38704, 44.42142, 44.42142, 44.42715, 
44.43861, 44.46152, 44.46152, 44.47298, NA, 44.50736, 44.29537, 
44.29537, 44.20369, NA, 44.50736, 44.4959, 44.49017, 44.47298, 
NA, 44.11775, 44.11775, NA, 46.42677, 46.42104, NA, 46.76482, 
46.76482, 46.65595, NA, 44.82249, 44.89697, 44.89697, 44.98865, 
NA, 44.98291, 44.72509, NA, 44.66206, 44.67925, 44.70217, 44.7079, 
44.71363, 44.72509, 44.748, 44.76519, 44.77665, 44.78238, 44.8053, 
44.82249, NA, 44.66206, 44.68498, 44.68498, 44.72509, 44.72509, 
NA, 44.98291, 44.98865, NA, 46.28353, 46.30645, 46.6273, 46.6273, 
46.80492, 46.79919, NA, 47.40653, 47.33204, 47.32631, 46.97681, 
46.97681, 46.80492, 46.80492, NA, 47.42372, 47.41799, 47.42372, 
47.41799, 47.41799, 47.42944, 47.44091, 47.45236, 47.4581, 47.46955, 
47.49247, 47.50966, 47.51539, 47.50393, 47.51539, 47.51539, 47.4982, 
47.48674, 47.46383, 47.43518, 47.41799, 47.40653, 47.4008, 47.38934, 
47.36642, 47.35496, 47.33204, 47.29767, 47.29767, 47.30912, 47.32059, 
47.31485, 47.29193, 47.29193, 47.28048, 47.26329, 47.2461, 47.22891, 
47.21172, 47.20599, 47.01691, NA, 46.28353, 46.30645, 46.31218, 
46.31791, 46.31791, 46.3351, 46.34082, 46.34082, 46.32937, 46.32364, 
46.31791, 46.32937, 46.3351, NA, 46.3351, 46.34082, 46.36374, 
46.36948, 46.36374, 46.35801, 46.35229, 46.34082, 46.34082, 46.36948, 
46.37521, 46.37521, NA, 46.37521, 46.62158, 46.6273, 46.80492, 
NA, 44.89697, 45.1548, NA, 45.1548, 45.1147, 45.07459, 45.04594, 
45.00583, 44.97146, 44.95427, 44.94281, 44.93135, NA, 44.89697, 
44.89697, 44.77665, NA, 45.1548, 45.1548, NA, 44.77665, 44.78238, 
44.78238, 44.79957, 44.8053, 44.83968, 44.83968, 44.85114, 44.87406, 
44.87978, 44.89124, 44.89697, 44.89697, 44.90843, 44.92562, 44.93135, 
NA, 45.74495, 45.57306, 45.57306, 45.41837, NA, 45.74495, 45.73349, 
NA, 45.29804, 45.30378, NA, 47.1487, 47.15443, NA, 46.63303, 
46.63303, 46.71325, NA, 46.63303, 46.63303, NA, 47.4008, 47.32631, 
47.32059, 47.1487, NA, 47.49247, 47.15443, NA, 48.02532, 47.93365, 
NA, 47.49247, 47.67009, 47.67009, 47.93365, NA, 47.45236, 48.06543, 
48.0597, 48.19721, NA, 43.85992, 43.85992, NA, 43.85992, 44.21515, 
NA, 44.20369, 44.21515, NA, 44.11775, 43.85992, NA, 46.16321, 
46.25488, NA, 46.16321, 46.16321, 46.16894, 46.19186, 46.20332, 
46.22623, 46.24342, 46.26061, 46.26634, 46.28353, NA, 44.47871, 
44.52455, 44.52455, 44.52455, 44.51309, 44.51309, 44.52455, 44.5532, 
44.5532, 44.64487, 44.64487, 44.71363, NA, 44.89697, 44.87406, 
44.83968, 44.81103, 44.81103, NA, 44.89697, 44.92562, 44.93708, 
44.93708, 44.91416, 44.89697, NA, 44.5532, 44.48444, 44.47871, 
NA, 44.5532, 44.63914, 44.63914, 44.72509, 44.72509, 44.81103, 
NA, 44.75373, 44.75946, 44.78238, 44.77665, 44.78238, 44.8053, 
44.86259, 44.89697, NA, 44.20369, 44.20369, NA, 43.85419, 43.85992, 
NA, 43.85992, 43.85992, 44.20369, NA, 43.85419, 44.20369, NA, 
45.76787, 45.95694, 45.95694, 46.11164, NA, 46.11164, 46.11164, 
NA, 45.76787, 45.76787, NA, 45.76787, 45.78506, NA, 46.11164, 
45.9283, 45.9283, 45.78506, NA, 43.51041, 43.85992, NA, 43.85992, 
43.50468, NA, 43.85992, 43.85992, NA, 43.51041, 43.85992, NA, 
43.50468, 43.84273, NA, 43.84273, 43.85992, 43.85992, NA, 43.85992, 
43.85992, NA, 43.85419, 43.51041, NA, 43.85419, 43.85992, NA, 
43.85992, 43.85992, NA, 44.20369, 44.20369, NA, 44.20369, 44.47871, 
NA, 44.20369, 44.38131, 44.37558, 44.46725, 44.4558, NA, 46.11164, 
46.11164, NA, 45.76787, 45.76787, NA, 46.03143, 45.93976, 45.93402, 
45.76787, NA, 46.11164, 46.03143, NA, 45.05167, 45.05167, 44.93708, 
44.89697, NA, 44.82249, 44.81676, 44.81676, 44.82249, 44.82249, 
44.81103, 44.79957, 44.81103, NA, 44.98865, 45.06886, 45.09178, 
45.10897, 45.1548, 45.16053, 45.17199, 45.18345, 45.20637, 45.21783, 
45.21783, 45.23502, 45.22929, 45.24648, NA, 43.85992, 43.85419, 
NA, 46.80492, 46.80492, NA, 45.74495, 45.74495, NA, 45.57306, 
45.74495, NA, 45.42982, 45.57306, NA, 47.84771, 47.84771, 47.89354, 
47.88208, NA, 47.88208, 47.71592, 47.71592, 47.36642, 47.36642, 
47.10859, 47.10859, 47.02264, NA, 43.50468, 43.85992, NA, 43.50468, 
43.85992, NA, 45.74495, 45.99705, 45.99705, 46.16321, NA, 45.74495, 
45.99705, 45.99705, 46.16894, NA, 45.33242, 45.24075, 45.23502, 
44.89697, NA, 45.41837, 45.41837, NA, 44.89697, 44.89697, NA, 
45.41837, 45.41837, 45.33242, NA, 45.41837, 45.24648, 45.24648, 
45.1548, NA, 48.54671, 48.54098, NA, 48.99935, 48.72433, 48.7186, 
48.54098, NA, 48.70141, 48.35764, NA, 48.62693, 48.06543, 48.0597, 
47.97375, 47.97375, 47.88208, NA, 45.17772, 45.16626, 45.16053, 
45.1548, NA, 44.81676, 44.81103, 44.89697, 44.89697, 44.93135, 
NA, 46.94816, 47.49247, 47.72739, 48.20294, NA, 48.53526, 48.71287, 
48.70714, 48.9306, 48.90768, 48.89622, 48.89049, 48.88476, 48.88476, 
48.89049, 48.89049, 48.91341, 48.92487, 48.94778, 48.9707, 48.99935, 
NA, 44.46725, 44.4558, 44.41569, 44.35839, 44.32975, 44.27818, 
44.25526, NA, 44.20369, 44.5532, NA, 44.5532, 44.5532, NA, 44.5532, 
44.54174, 44.53028, 44.4959, 44.48444, 44.46725, NA, 44.20369, 
44.20369, NA, 44.21515, 44.63914, NA, 44.21515, 44.20369, NA, 
44.63914, 44.63914, NA, 44.21515, 44.21515, NA, 44.21515, 44.5532, 
NA, 44.63914, 44.64487, 44.5532, 44.5532, NA, 47.50393, 47.15443, 
NA, 47.49247, 47.50393, NA, 48.16856, 48.17429, NA, 48.16856, 
48.16856, NA, 48.53526, 48.54098, NA, 43.85992, 43.85992, NA, 
44.98291, 44.98291, 44.89697, NA, 44.72509, 44.89697, NA, 44.72509, 
44.63914, 44.63914, 44.72509, 44.72509, NA, 44.98291, 44.98291, 
NA, 44.89697, 44.89697, NA, 45.33242, 45.33242, 45.29231, 45.29231, 
NA, 44.98291, 45.29231, NA, 45.83662, 45.99132, 45.99705, 46.16321, 
NA, 45.57306, 45.57306, NA, 45.78506, 45.78506, NA, 46.3351, 
46.3351, 46.3351, 45.78506, NA, 43.85992, 43.84273, 43.84273, 
NA, 43.85992, 43.85992, NA, 44.21515, 43.85992, NA, 44.21515, 
44.21515, NA, 44.47298, 44.47298, NA, 44.47298, 44.46725, NA, 
43.49895, 43.85992, NA, 47.50393, 47.4982, NA, 44.11775, 44.12348, 
44.20369, 44.20369, NA, 44.11775, 43.85992, NA, 46.11164, 46.28353, 
46.28353, 46.37521, NA, 46.71898, 46.37521, NA, 46.63303, 46.28353, 
46.28353, 46.11164, NA, 47.93365, 47.93365, NA, 48.16856, 48.03105, 
48.03105, 47.96803, NA, 47.93365, 47.9623, 47.96803, NA, 43.85419, 
43.85992, NA, 47.93365, 47.84771, 47.84771, 47.76176, 47.76176, 
47.76176, 47.85344, 47.85917, 47.96803, NA, 45.41837, 45.59025, 
45.59025, 45.76787, NA, 45.76787, 45.60171, 45.60171, 45.41837, 
NA, 45.41837, 45.41837, NA, 45.13189, 45.13189, 44.89697, 44.89697, 
NA, 44.50736, 44.53028, 44.58757, 44.58757, 44.5933, 44.59903, 
44.61049, 44.63341, 44.63914, 44.67925, 44.67925, 44.7079, 44.71363, 
NA, 44.71363, 44.5532, 44.5532, NA, 44.72509, 44.72509, 44.5532, 
44.5532, 44.47298, NA, 44.71363, 44.71936, 44.75373, 44.75946, 
44.77665, NA, 44.5532, 44.5532, NA, 44.20369, 44.20369, NA, 44.20369, 
44.20369, NA, 44.66206, 44.63341, 44.63914, 44.62769, 44.61049, 
44.58757, 44.57038, 44.5532, NA, 45.57306, 45.55587, 45.52723, 
45.49285, 45.45274, 45.44701, 45.41837, NA, 45.41837, 45.41263, 
45.40691, 45.38972, 45.39545, 45.38972, 45.36107, 45.34961, 45.34388, 
45.30378, 45.31523, 45.30378, 45.31523, 45.30378, 45.30378, 45.30378, 
45.27512, 45.25794, 45.24648, NA, 45.78506, 45.78506, NA, 45.29231, 
45.30378, 45.32096, 45.32669, 45.30378, 45.30378, 45.32669, 45.33815, 
45.33815, 45.34961, 45.38972, 45.41263, 45.41837, NA, 44.20369, 
43.85992, NA, 45.41837, 45.41837, NA, 45.59025, 45.76787, NA, 
46.37521, 46.37521, NA, 46.03143, 46.03143, NA, 44.11775, 44.20369, 
44.19223, NA), V4 = 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, 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, 3, 
3, 3, 3, 3, 3, 3, 3, 3, 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, 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, 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, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 
7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 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, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 
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, 10, 10, 10, 10, 10, 10, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 
11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 
12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 
12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 
13, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 
14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 
15, 15, 15, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 17, 17, 
17, 17, 17, 17, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 
18, 18, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 
19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 
19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 
20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 21, 21, 21, 
21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 
22, 22, 22, 22, 22, 22, 22, 22, 22, 23, 23, 23, 23, 23, 23, 23, 
23, 23, 23, 23, 23, 23, 24, 24, 24, 24, 24, 24, 24, 24, 24, 25, 
25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 26, 26, 26, 26, 26, 
26, 26, 26, 26, 26, 26, 26, 26, 26, 27, 27, 27, 27, 27, 27, 27, 
27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 27, 
27, 27, 27, 27, 27, 27, 28, 28, 28, 29, 29, 29, 30, 30, 30, 30, 
30, 30, 30, 30, 30, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 
31, 31, 31, 32, 32, 32, 32, 32, 32, 33, 33, 33, 33, 33, 33, 33, 
33, 33, 33, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 
34, 34, 34, 34, 34, 34, 34, 35, 35, 35, 35, 35, 35, 35, 35, 36, 
36, 36, 36, 36, 36, 36, 36, 36, 36, 37, 37, 37, 37, 37, 37, 37, 
37, 37, 37, 37, 38, 38, 38, 38, 38, 39, 39, 39, 39, 39, 39, 39, 
39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 40, 40, 40, 40, 40, 40, 
40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 
40, 40, 41, 41, 41, 41, 41, 41, 41, 41, 41, 42, 42, 42, 42, 42, 
42, 42, 42, 42, 42, 42, 43, 43, 43, 43, 43, 43, 44, 44, 44, 44, 
44, 44, 44, 44, 44, 45, 45, 45, 46, 46, 46, 46, 46, 46, 46, 46, 
46, 46, 46, 46, 46, 46, 46, 46, 47, 47, 47, 47, 47, 47, 47, 47, 
47, 47, 47, 48, 48, 48, 48, 48, 48, 48, 48, 49, 49, 49, 49, 49, 
49, 49, 49, 50, 50, 50, 50, 51, 51, 51, 51, 51, 51, 51, 51, 51, 
52, 52, 52, 52, 52, 52, 53, 53, 53, 54, 54, 54, 55, 55, 55, 55, 
55, 55, 55, 55, 56, 56, 56, 56, 56, 56, 56, 56, 56, 56, 56, 56, 
56, 57, 57, 57, 57, 57, 57, 57, 57, 57, 57, 57, 57, 59, 59, 59, 
60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 61, 61, 61, 61, 61, 61, 
61, 61, 61, 61, 61, 61, 61, 62, 62, 62, 62, 62, 64, 64, 64, 64, 
64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 65, 65, 
65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 66, 66, 66, 66, 66, 66, 
66, 66, 66, 69, 69, 69, 69, 69, 69, 69, 69, 69, 70, 70, 70, 70, 
70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 70, 
70, 70, 70, 70, 70, 70, 70, 70, 73, 73, 73, 73, 73, 73, 73, 73, 
73, 73, 73, 73, 73, 73, 73, 73, 73, 74, 74, 74, 75, 75, 75, 75, 
75, 75, 77, 77, 77, 78, 78, 78, 79, 79, 79, 79), V5 = c(9, 9, 
9, 9, 9, 11, 11, 11, 18, 18, 18, 18, 18, 31, 31, 31, 33, 33, 
33, 48, 48, 48, 48, 58, 58, 58, 72, 72, 72, 13, 13, 13, 27, 27, 
27, 27, 27, 27, 27, 27, 27, 30, 30, 30, 62, 62, 62, 62, 70, 70, 
70, 82, 82, 82, 14, 14, 14, 14, 14, 15, 15, 15, 29, 29, 29, 29, 
29, 43, 43, 43, 54, 54, 54, 56, 56, 56, 80, 80, 80, 11, 11, 11, 
11, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 29, 29, 
29, 31, 31, 31, 36, 36, 36, 36, 36, 39, 39, 39, 39, 39, 44, 44, 
44, 44, 44, 57, 57, 57, 68, 68, 68, 48, 48, 48, 49, 49, 49, 49, 
49, 70, 70, 70, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 37, 
37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 75, 75, 75, 75, 
76, 76, 76, 76, 76, 78, 78, 78, 8, 8, 8, 22, 22, 22, 40, 40, 
40, 40, 45, 45, 45, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 
52, 52, 52, 52, 52, 81, 81, 81, 83, 83, 83, 17, 17, 17, 17, 17, 
52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 
64, 64, 64, 64, 64, 65, 65, 65, 65, 65, 83, 83, 83, 58, 58, 58, 
72, 72, 72, 72, 27, 27, 27, 27, 27, 46, 46, 46, 69, 69, 69, 69, 
69, 69, 69, 69, 69, 69, 69, 69, 69, 71, 71, 71, 71, 71, 71, 86, 
86, 86, 18, 18, 18, 18, 18, 18, 18, 29, 29, 29, 29, 29, 29, 29, 
29, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 
31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 
31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 49, 49, 49, 49, 49, 
49, 49, 49, 49, 49, 49, 49, 49, 49, 77, 77, 77, 77, 77, 77, 77, 
77, 77, 77, 77, 77, 77, 80, 80, 80, 80, 80, 34, 34, 34, 37, 37, 
37, 37, 37, 37, 37, 37, 37, 37, 65, 65, 65, 65, 76, 76, 76, 87, 
87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 
30, 30, 30, 30, 30, 58, 58, 58, 82, 82, 82, 54, 54, 54, 56, 56, 
56, 56, 84, 84, 84, 29, 29, 29, 29, 29, 43, 43, 43, 57, 57, 57, 
60, 60, 60, 60, 60, 38, 38, 38, 38, 38, 32, 32, 32, 51, 51, 51, 
64, 64, 64, 83, 83, 83, 48, 48, 48, 49, 49, 49, 49, 49, 49, 49, 
49, 49, 49, 49, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 
25, 27, 27, 27, 27, 27, 27, 62, 62, 62, 62, 62, 62, 62, 66, 66, 
66, 66, 69, 69, 69, 69, 69, 69, 69, 82, 82, 82, 82, 82, 82, 82, 
82, 82, 25, 25, 25, 50, 50, 50, 55, 55, 55, 55, 74, 74, 74, 26, 
26, 26, 26, 26, 56, 56, 56, 61, 61, 61, 73, 73, 73, 77, 77, 77, 
77, 77, 24, 24, 24, 45, 45, 45, 81, 81, 81, 28, 28, 28, 50, 50, 
50, 55, 55, 55, 55, 85, 85, 85, 50, 50, 50, 74, 74, 74, 81, 81, 
81, 55, 55, 55, 66, 66, 66, 79, 79, 79, 79, 79, 79, 56, 56, 56, 
75, 75, 75, 78, 78, 78, 78, 78, 84, 84, 84, 62, 62, 62, 62, 62, 
69, 69, 69, 69, 69, 69, 69, 69, 69, 86, 86, 86, 86, 86, 86, 86, 
86, 86, 86, 86, 86, 86, 86, 86, 85, 85, 85, 80, 80, 80, 33, 33, 
33, 48, 48, 48, 70, 70, 70, 36, 36, 36, 36, 36, 72, 72, 72, 72, 
72, 72, 72, 72, 72, 45, 45, 45, 53, 53, 53, 48, 48, 48, 48, 48, 
58, 58, 58, 58, 58, 47, 47, 47, 47, 47, 61, 61, 61, 65, 65, 65, 
73, 73, 73, 73, 76, 76, 76, 76, 76, 44, 44, 44, 68, 68, 68, 68, 
68, 39, 39, 39, 72, 72, 72, 72, 72, 72, 72, 76, 76, 76, 76, 76, 
87, 87, 87, 87, 87, 87, 72, 72, 72, 72, 72, 68, 68, 68, 68, 68, 
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 52, 52, 52, 52, 
52, 52, 52, 52, 66, 66, 66, 69, 69, 69, 71, 71, 71, 71, 71, 71, 
71, 81, 81, 81, 42, 42, 42, 59, 59, 59, 87, 87, 87, 51, 51, 51, 
64, 64, 64, 87, 87, 87, 87, 87, 54, 54, 54, 60, 60, 60, 57, 57, 
57, 60, 60, 60, 68, 68, 68, 83, 83, 83, 47, 47, 47, 47, 65, 65, 
65, 71, 71, 71, 71, 71, 71, 86, 86, 86, 65, 65, 65, 73, 73, 73, 
73, 73, 86, 86, 86, 49, 49, 49, 49, 49, 70, 70, 70, 73, 73, 73, 
77, 77, 77, 77, 77, 55, 55, 55, 55, 53, 53, 53, 59, 59, 59, 64, 
64, 64, 65, 65, 65, 71, 71, 71, 67, 67, 67, 60, 60, 60, 79, 79, 
79, 79, 79, 85, 85, 85, 77, 77, 77, 77, 77, 80, 80, 80, 84, 84, 
84, 84, 84, 60, 60, 60, 60, 60, 60, 60, 60, 63, 63, 63, 63, 67, 
67, 67, 63, 63, 63, 63, 63, 63, 63, 63, 63, 63, 73, 73, 73, 73, 
73, 75, 75, 75, 75, 75, 76, 76, 76, 82, 82, 82, 82, 82, 65, 65, 
65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 65, 87, 87, 87, 87, 
71, 71, 71, 71, 71, 71, 87, 87, 87, 87, 87, 87, 69, 69, 69, 74, 
74, 74, 81, 81, 81, 71, 71, 71, 71, 71, 71, 71, 71, 71, 73, 73, 
73, 73, 73, 73, 73, 73, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 
86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 77, 77, 77, 86, 86, 86, 
86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 81, 81, 81, 76, 76, 
76, 78, 78, 78, 80, 80, 80, 84, 84, 84, 85, 85, 85, 85)), .Names = c("V1", 
"V2", "V3", "V4", "V5"), 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", "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", "381", "382", "383", "384", "385", "386", "387", "388", 
"389", "390", "391", "392", "393", "394", "395", "396", "397", 
"398", "399", "400", "401", "402", "403", "404", "405", "406", 
"407", "408", "409", "410", "411", "412", "413", "414", "415", 
"416", "417", "418", "419", "420", "421", "422", "423", "424", 
"425", "426", "427", "428", "429", "430", "431", "432", "433", 
"434", "435", "436", "437", "438", "439", "440", "441", "442", 
"443", "444", "445", "446", "447", "448", "449", "450", "451", 
"452", "453", "454", "455", "456", "457", "458", "459", "460", 
"461", "462", "463", "464", "465", "466", "467", "468", "469", 
"470", "471", "472", "473", "474", "475", "476", "477", "478", 
"479", "480", "481", "482", "483", "484", "485", "486", "487", 
"488", "489", "490", "491", "492", "493", "494", "495", "496", 
"497", "498", "499", "500", "501", "502", "503", "504", "505", 
"506", "507", "508", "509", "510", "511", "512", "513", "514", 
"515", "516", "517", "518", "519", "520", "521", "522", "523", 
"524", "525", "526", "527", "528", "529", "530", "531", "532", 
"533", "534", "535", "536", "537", "538", "539", "540", "541", 
"542", "543", "544", "545", "546", "547", "548", "549", "550", 
"551", "552", "553", "554", "555", "556", "557", "558", "559", 
"560", "561", "562", "563", "564", "565", "566", "567", "568", 
"569", "570", "571", "572", "573", "574", "575", "576", "577", 
"578", "579", "580", "581", "582", "583", "584", "585", "586", 
"587", "588", "589", "590", "591", "592", "593", "594", "595", 
"596", "597", "598", "599", "600", "601", "602", "603", "604", 
"605", "606", "607", "608", "609", "610", "611", "612", "613", 
"614", "615", "616", "617", "618", "619", "620", "621", "622", 
"623", "624", "625", "626", "627", "628", "629", "630", "631", 
"632", "633", "634", "635", "636", "637", "638", "639", "640", 
"641", "642", "643", "644", "645", "646", "647", "648", "649", 
"650", "651", "652", "653", "654", "655", "656", "657", "658", 
"659", "660", "661", "662", "663", "664", "665", "666", "667", 
"668", "669", "670", "671", "672", "673", "674", "675", "676", 
"677", "678", "679", "680", "681", "682", "683", "684", "685", 
"686", "687", "688", "689", "690", "691", "692", "693", "694", 
"695", "696", "697", "698", "699", "700", "701", "702", "703", 
"704", "705", "706", "707", "708", "709", "710", "711", "712", 
"713", "714", "715", "716", "717", "718", "719", "720", "721", 
"722", "723", "724", "725", "726", "727", "728", "729", "730", 
"731", "732", "733", "734", "735", "736", "737", "738", "739", 
"740", "741", "742", "743", "744", "745", "746", "747", "748", 
"749", "750", "751", "752", "753", "754", "755", "756", "757", 
"758", "759", "760", "761", "762", "763", "764", "765", "766", 
"767", "768", "769", "770", "771", "772", "773", "774", "775", 
"776", "777", "778", "779", "780", "781", "782", "783", "784", 
"785", "786", "787", "788", "789", "790", "791", "792", "793", 
"794", "795", "796", "797", "798", "799", "800", "801", "802", 
"803", "804", "805", "806", "807", "808", "809", "810", "811", 
"812", "813", "814", "815", "816", "817", "818", "819", "820", 
"821", "822", "823", "824", "825", "826", "827", "828", "829", 
"830", "831", "832", "833", "834", "835", "836", "837", "838", 
"839", "840", "841", "842", "843", "844", "845", "846", "847", 
"848", "849", "850", "851", "852", "853", "854", "855", "856", 
"857", "858", "859", "860", "861", "862", "863", "864", "865", 
"866", "867", "868", "869", "870", "871", "872", "873", "874", 
"875", "876", "877", "878", "879", "880", "881", "882", "883", 
"884", "885", "886", "887", "888", "889", "890", "891", "892", 
"893", "894", "895", "896", "897", "898", "899", "900", "901", 
"902", "903", "904", "905", "906", "907", "908", "909", "910", 
"911", "912", "913", "914", "915", "916", "917", "918", "919", 
"920", "921", "922", "923", "924", "925", "926", "927", "928", 
"929", "930", "931", "932", "933", "934", "935", "936", "937", 
"938", "939", "940", "941", "942", "943", "944", "945", "946", 
"947", "948", "949", "950", "951", "952", "953", "954", "955", 
"956", "957", "958", "959", "960", "961", "962", "963", "964", 
"965", "966", "967", "968", "969", "970", "971", "972", "973", 
"974", "975", "976", "977", "978", "979", "980", "981", "982", 
"983", "984", "985", "986", "987", "988", "989", "990", "991", 
"992", "993", "994", "995", "996", "997", "998", "999", "1000", 
"1001", "1002", "1003", "1004", "1005", "1006", "1007", "1008", 
"1009", "1010", "1011", "1012", "1013", "1014", "1015", "1016", 
"1017", "1018", "1019", "1020", "1021", "1022", "1023", "1024", 
"1025", "1026", "1027", "1028", "1029", "1030", "1031", "1032", 
"1033", "1034", "1035", "1036", "1037", "1038", "1039", "1040", 
"1041", "1042", "1043", "1044", "1045", "1046", "1047", "1048", 
"1049", "1050", "1051", "1052", "1053", "1054", "1055", "1056"
))


lines(edgelines[1:5,2:3],col=y.shad[1],lwd=5)
lines(edgelines[6:8,2:3],col=y.shad[2],lwd=5)
lines(edgelines[9:13,2:3],col=y.shad[3],lwd=5)
lines(edgelines[14:16,2:3],col=y.shad[4],lwd=5)
lines(edgelines[17:19,2:3],col=y.shad[5],lwd=5)
lines(edgelines[20:23,2:3],col=y.shad[6],lwd=5)
lines(edgelines[24:26,2:3],col=y.shad[7],lwd=5)
lines(edgelines[27:29,2:3],col=y.shad[8],lwd=5)
lines(edgelines[30:32,2:3],col=y.shad[9],lwd=5)
lines(edgelines[33:41,2:3],col=y.shad[10],lwd=5)
lines(edgelines[42:44,2:3],col=y.shad[11],lwd=5)
lines(edgelines[45:48,2:3],col=y.shad[12],lwd=5)
lines(edgelines[49:51,2:3],col=y.shad[13],lwd=5)
lines(edgelines[52:54,2:3],col=y.shad[14],lwd=5)
lines(edgelines[55:59,2:3],col=y.shad[15],lwd=5)
lines(edgelines[60:62,2:3],col=y.shad[16],lwd=5)
lines(edgelines[63:67,2:3],col=y.shad[17],lwd=5)
lines(edgelines[68:70,2:3],col=y.shad[18],lwd=5)
lines(edgelines[71:73,2:3],col=y.shad[19],lwd=5)
lines(edgelines[74:76,2:3],col=y.shad[20],lwd=5)
lines(edgelines[77:79,2:3],col=y.shad[21],lwd=5)
lines(edgelines[80:83,2:3],col=y.shad[22],lwd=5)
lines(edgelines[84:96,2:3],col=y.shad[23],lwd=5)
lines(edgelines[97:99,2:3],col=y.shad[24],lwd=5)
lines(edgelines[100:102,2:3],col=y.shad[25],lwd=5)
lines(edgelines[103:107,2:3],col=y.shad[26],lwd=5)
lines(edgelines[108:112,2:3],col=y.shad[27],lwd=5)
lines(edgelines[113:117,2:3],col=y.shad[28],lwd=5)
lines(edgelines[118:120,2:3],col=y.shad[29],lwd=5)
lines(edgelines[121:123,2:3],col=y.shad[30],lwd=5)
lines(edgelines[124:126,2:3],col=y.shad[31],lwd=5)
lines(edgelines[127:131,2:3],col=y.shad[32],lwd=5)
lines(edgelines[132:134,2:3],col=y.shad[33],lwd=5)
lines(edgelines[135:145,2:3],col=y.shad[34],lwd=5)
lines(edgelines[146:158,2:3],col=y.shad[35],lwd=5)
lines(edgelines[159:162,2:3],col=y.shad[36],lwd=5)
lines(edgelines[163:167,2:3],col=y.shad[37],lwd=5)
lines(edgelines[168:170,2:3],col=y.shad[38],lwd=5)
lines(edgelines[171:173,2:3],col=y.shad[39],lwd=5)
lines(edgelines[174:176,2:3],col=y.shad[40],lwd=5)
lines(edgelines[177:180,2:3],col=y.shad[41],lwd=5)
lines(edgelines[181:183,2:3],col=y.shad[42],lwd=5)
lines(edgelines[184:199,2:3],col=y.shad[43],lwd=5)
lines(edgelines[200:202,2:3],col=y.shad[44],lwd=5)
lines(edgelines[203:205,2:3],col=y.shad[45],lwd=5)
lines(edgelines[206:210,2:3],col=y.shad[46],lwd=5)
lines(edgelines[211:226,2:3],col=y.shad[47],lwd=5)
lines(edgelines[227:231,2:3],col=y.shad[48],lwd=5)
lines(edgelines[232:236,2:3],col=y.shad[49],lwd=5)
lines(edgelines[237:239,2:3],col=y.shad[50],lwd=5)
lines(edgelines[240:242,2:3],col=y.shad[51],lwd=5)
lines(edgelines[243:246,2:3],col=y.shad[52],lwd=5)
lines(edgelines[247:251,2:3],col=y.shad[53],lwd=5)
lines(edgelines[252:254,2:3],col=y.shad[54],lwd=5)
lines(edgelines[255:267,2:3],col=y.shad[55],lwd=5)
lines(edgelines[268:273,2:3],col=y.shad[56],lwd=5)
lines(edgelines[274:276,2:3],col=y.shad[57],lwd=5)
lines(edgelines[277:283,2:3],col=y.shad[58],lwd=5)
lines(edgelines[284:291,2:3],col=y.shad[59],lwd=5)
lines(edgelines[292:333,2:3],col=y.shad[60],lwd=5)
lines(edgelines[334:347,2:3],col=y.shad[61],lwd=5)
lines(edgelines[348:360,2:3],col=y.shad[62],lwd=5)
lines(edgelines[361:365,2:3],col=y.shad[63],lwd=5)
lines(edgelines[366:368,2:3],col=y.shad[64],lwd=5)
lines(edgelines[369:378,2:3],col=y.shad[65],lwd=5)
lines(edgelines[379:382,2:3],col=y.shad[66],lwd=5)
lines(edgelines[383:385,2:3],col=y.shad[67],lwd=5)
lines(edgelines[386:402,2:3],col=y.shad[68],lwd=5)
lines(edgelines[403:407,2:3],col=y.shad[69],lwd=5)
lines(edgelines[408:410,2:3],col=y.shad[70],lwd=5)
lines(edgelines[411:413,2:3],col=y.shad[71],lwd=5)
lines(edgelines[414:416,2:3],col=y.shad[72],lwd=5)
lines(edgelines[417:420,2:3],col=y.shad[73],lwd=5)
lines(edgelines[421:423,2:3],col=y.shad[74],lwd=5)
lines(edgelines[424:428,2:3],col=y.shad[75],lwd=5)
lines(edgelines[429:431,2:3],col=y.shad[76],lwd=5)
lines(edgelines[432:434,2:3],col=y.shad[77],lwd=5)
lines(edgelines[435:439,2:3],col=y.shad[78],lwd=5)
lines(edgelines[440:444,2:3],col=y.shad[79],lwd=5)
lines(edgelines[445:447,2:3],col=y.shad[80],lwd=5)
lines(edgelines[448:450,2:3],col=y.shad[81],lwd=5)
lines(edgelines[451:453,2:3],col=y.shad[82],lwd=5)
lines(edgelines[454:456,2:3],col=y.shad[83],lwd=5)
lines(edgelines[457:459,2:3],col=y.shad[84],lwd=5)
lines(edgelines[460:470,2:3],col=y.shad[85],lwd=5)
lines(edgelines[471:483,2:3],col=y.shad[86],lwd=5)
lines(edgelines[484:489,2:3],col=y.shad[87],lwd=5)
lines(edgelines[490:496,2:3],col=y.shad[88],lwd=5)
lines(edgelines[497:500,2:3],col=y.shad[89],lwd=5)
lines(edgelines[501:507,2:3],col=y.shad[90],lwd=5)
lines(edgelines[508:516,2:3],col=y.shad[91],lwd=5)
lines(edgelines[517:519,2:3],col=y.shad[92],lwd=5)
lines(edgelines[520:522,2:3],col=y.shad[93],lwd=5)
lines(edgelines[523:526,2:3],col=y.shad[94],lwd=5)
lines(edgelines[527:529,2:3],col=y.shad[95],lwd=5)
lines(edgelines[530:534,2:3],col=y.shad[96],lwd=5)
lines(edgelines[535:537,2:3],col=y.shad[97],lwd=5)
lines(edgelines[538:540,2:3],col=y.shad[98],lwd=5)
lines(edgelines[541:543,2:3],col=y.shad[99],lwd=5)
lines(edgelines[544:548,2:3],col=y.shad[100],lwd=5)
lines(edgelines[549:551,2:3],col=y.shad[101],lwd=5)
lines(edgelines[552:554,2:3],col=y.shad[102],lwd=5)
lines(edgelines[555:557,2:3],col=y.shad[103],lwd=5)
lines(edgelines[558:560,2:3],col=y.shad[104],lwd=5)
lines(edgelines[561:563,2:3],col=y.shad[105],lwd=5)
lines(edgelines[564:567,2:3],col=y.shad[106],lwd=5)
lines(edgelines[568:570,2:3],col=y.shad[107],lwd=5)
lines(edgelines[571:573,2:3],col=y.shad[108],lwd=5)
lines(edgelines[574:576,2:3],col=y.shad[109],lwd=5)
lines(edgelines[577:579,2:3],col=y.shad[110],lwd=5)
lines(edgelines[580:582,2:3],col=y.shad[111],lwd=5)
lines(edgelines[583:585,2:3],col=y.shad[112],lwd=5)
lines(edgelines[586:591,2:3],col=y.shad[113],lwd=5)
lines(edgelines[592:594,2:3],col=y.shad[114],lwd=5)
lines(edgelines[595:597,2:3],col=y.shad[115],lwd=5)
lines(edgelines[598:602,2:3],col=y.shad[116],lwd=5)
lines(edgelines[603:605,2:3],col=y.shad[117],lwd=5)
lines(edgelines[606:610,2:3],col=y.shad[118],lwd=5)
lines(edgelines[611:619,2:3],col=y.shad[119],lwd=5)
lines(edgelines[620:634,2:3],col=y.shad[120],lwd=5)
lines(edgelines[635:637,2:3],col=y.shad[121],lwd=5)
lines(edgelines[638:640,2:3],col=y.shad[122],lwd=5)
lines(edgelines[641:643,2:3],col=y.shad[123],lwd=5)
lines(edgelines[644:646,2:3],col=y.shad[124],lwd=5)
lines(edgelines[647:649,2:3],col=y.shad[125],lwd=5)
lines(edgelines[650:654,2:3],col=y.shad[126],lwd=5)
lines(edgelines[655:663,2:3],col=y.shad[127],lwd=5)
lines(edgelines[664:666,2:3],col=y.shad[128],lwd=5)
lines(edgelines[667:669,2:3],col=y.shad[129],lwd=5)
lines(edgelines[670:674,2:3],col=y.shad[130],lwd=5)
lines(edgelines[675:679,2:3],col=y.shad[131],lwd=5)
lines(edgelines[680:684,2:3],col=y.shad[132],lwd=5)
lines(edgelines[685:687,2:3],col=y.shad[133],lwd=5)
lines(edgelines[688:690,2:3],col=y.shad[134],lwd=5)
lines(edgelines[691:694,2:3],col=y.shad[135],lwd=5)
lines(edgelines[695:699,2:3],col=y.shad[136],lwd=5)
lines(edgelines[700:702,2:3],col=y.shad[137],lwd=5)
lines(edgelines[703:707,2:3],col=y.shad[138],lwd=5)
lines(edgelines[708:710,2:3],col=y.shad[139],lwd=5)
lines(edgelines[711:717,2:3],col=y.shad[140],lwd=5)
lines(edgelines[718:722,2:3],col=y.shad[141],lwd=5)
lines(edgelines[723:728,2:3],col=y.shad[142],lwd=5)
lines(edgelines[729:733,2:3],col=y.shad[143],lwd=5)
lines(edgelines[734:750,2:3],col=y.shad[144],lwd=5)
lines(edgelines[751:758,2:3],col=y.shad[145],lwd=5)
lines(edgelines[759:761,2:3],col=y.shad[146],lwd=5)
lines(edgelines[762:764,2:3],col=y.shad[147],lwd=5)
lines(edgelines[765:771,2:3],col=y.shad[148],lwd=5)
lines(edgelines[772:774,2:3],col=y.shad[149],lwd=5)
lines(edgelines[775:777,2:3],col=y.shad[150],lwd=5)
lines(edgelines[778:780,2:3],col=y.shad[151],lwd=5)
lines(edgelines[781:783,2:3],col=y.shad[152],lwd=5)
lines(edgelines[784:786,2:3],col=y.shad[153],lwd=5)
lines(edgelines[787:789,2:3],col=y.shad[154],lwd=5)
lines(edgelines[790:794,2:3],col=y.shad[155],lwd=5)
lines(edgelines[795:797,2:3],col=y.shad[156],lwd=5)
lines(edgelines[798:800,2:3],col=y.shad[157],lwd=5)
lines(edgelines[801:803,2:3],col=y.shad[158],lwd=5)
lines(edgelines[804:806,2:3],col=y.shad[159],lwd=5)
lines(edgelines[807:809,2:3],col=y.shad[160],lwd=5)
lines(edgelines[810:812,2:3],col=y.shad[161],lwd=5)
lines(edgelines[813:816,2:3],col=y.shad[162],lwd=5)
lines(edgelines[817:819,2:3],col=y.shad[163],lwd=5)
lines(edgelines[820:825,2:3],col=y.shad[164],lwd=5)
lines(edgelines[826:828,2:3],col=y.shad[165],lwd=5)
lines(edgelines[829:831,2:3],col=y.shad[166],lwd=5)
lines(edgelines[832:836,2:3],col=y.shad[167],lwd=5)
lines(edgelines[837:839,2:3],col=y.shad[168],lwd=5)
lines(edgelines[840:844,2:3],col=y.shad[169],lwd=5)
lines(edgelines[845:847,2:3],col=y.shad[170],lwd=5)
lines(edgelines[848:850,2:3],col=y.shad[171],lwd=5)
lines(edgelines[851:855,2:3],col=y.shad[172],lwd=5)
lines(edgelines[856:859,2:3],col=y.shad[173],lwd=5)
lines(edgelines[860:862,2:3],col=y.shad[174],lwd=5)
lines(edgelines[863:865,2:3],col=y.shad[175],lwd=5)
lines(edgelines[866:868,2:3],col=y.shad[176],lwd=5)
lines(edgelines[869:871,2:3],col=y.shad[177],lwd=5)
lines(edgelines[872:874,2:3],col=y.shad[178],lwd=5)
lines(edgelines[875:877,2:3],col=y.shad[179],lwd=5)
lines(edgelines[878:880,2:3],col=y.shad[180],lwd=5)
lines(edgelines[881:885,2:3],col=y.shad[181],lwd=5)
lines(edgelines[886:888,2:3],col=y.shad[182],lwd=5)
lines(edgelines[889:893,2:3],col=y.shad[183],lwd=5)
lines(edgelines[894:896,2:3],col=y.shad[184],lwd=5)
lines(edgelines[897:901,2:3],col=y.shad[185],lwd=5)
lines(edgelines[902:909,2:3],col=y.shad[186],lwd=5)
lines(edgelines[910:913,2:3],col=y.shad[187],lwd=5)
lines(edgelines[914:916,2:3],col=y.shad[188],lwd=5)
lines(edgelines[917:926,2:3],col=y.shad[189],lwd=5)
lines(edgelines[927:931,2:3],col=y.shad[190],lwd=5)
lines(edgelines[932:936,2:3],col=y.shad[191],lwd=5)
lines(edgelines[937:939,2:3],col=y.shad[192],lwd=5)
lines(edgelines[940:944,2:3],col=y.shad[193],lwd=5)
lines(edgelines[945:958,2:3],col=y.shad[194],lwd=5)
lines(edgelines[959:962,2:3],col=y.shad[195],lwd=5)
lines(edgelines[963:968,2:3],col=y.shad[196],lwd=5)
lines(edgelines[969:974,2:3],col=y.shad[197],lwd=5)
lines(edgelines[975:977,2:3],col=y.shad[198],lwd=5)
lines(edgelines[978:980,2:3],col=y.shad[199],lwd=5)
lines(edgelines[981:983,2:3],col=y.shad[200],lwd=5)
lines(edgelines[984:992,2:3],col=y.shad[201],lwd=5)
lines(edgelines[993:1000,2:3],col=y.shad[202],lwd=5)
lines(edgelines[1001:1020,2:3],col=y.shad[203],lwd=5)
lines(edgelines[1021:1023,2:3],col=y.shad[204],lwd=5)
lines(edgelines[1024:1037,2:3],col=y.shad[205],lwd=5)
lines(edgelines[1038:1040,2:3],col=y.shad[206],lwd=5)
lines(edgelines[1041:1043,2:3],col=y.shad[207],lwd=5)
lines(edgelines[1044:1046,2:3],col=y.shad[208],lwd=5)
lines(edgelines[1047:1049,2:3],col=y.shad[209],lwd=5)
lines(edgelines[1050:1052,2:3],col=y.shad[210],lwd=5)
lines(edgelines[1053:1056,2:3],col=y.shad[211],lwd=5)
#par(old.par)
}