Code
── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ──
✔ ggplot2 3.3.5 ✔ purrr 0.3.4
✔ tibble 3.1.6 ✔ dplyr 1.0.7
✔ tidyr 1.1.4 ✔ stringr 1.4.0
✔ readr 2.0.2 ✔ forcats 0.5.1
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag() masks stats::lag()
Code
library (readxl)library (janitor)
Attaching package: 'janitor'
The following objects are masked from 'package:stats':
chisq.test, fisher.test
Code
library (glue)library (here)
here() starts at /Users/jgol0005/Code/report-satellites-conflict
Code
library (pins)library (themes360info)
✔ Using preferred font, ITC Franklin Gothic.
ℹ Specify a different font to use with 360info themes by calling
register_360fonts() or by setting options("themes360info.franklin.pref") to
either "itc" or "libre" (or "none" to disable automatic font loading).
Code
library (ggtext)library (nVennR)# create /data/.cache <- here ("data" , ".cache" )dir.create (data_cache, showWarnings = FALSE )# download the zip file (keeping it cached if we redo this) <- board_url (c (sat = "https://www.ucsusa.org/media/11492" ),cache = data_cache) %>% pin_download ("sat" )
Code
# import and tidy the data # (user guide: https://s3.amazonaws.com/ucs-documents/nuclear-weapons/sat-database/4-11-17-update/User+Guide+1-1-17+wAppendix.pdf) <- read_excel (sat_wb) %>% clean_names ()
New names:
* `` -> ...29
* Source -> Source...31
* Source -> Source...32
* Source -> Source...33
* Source -> Source...34
* ...
Code
%>% select (name_official = current_official_name_of_satellite,name_alt = name_of_satellite_alternate_names,country_unreg = country_org_of_un_registry,owner = operator_owner,country_owner = country_of_operator_owner,country_contractor = country_of_contractor,-> %>% group_by (country_owner, detailed_purpose) %>% summarise (num_sats = n (),total_mass_dry_kg = sum (as.numeric (dry_mass_kg), na.rm = TRUE )) ->
Warning in mask$eval_all_summarise(quo): NAs introduced by coercion
`summarise()` has grouped output by 'country_owner'. You can override using the
`.groups` argument.
Main countries and purposes
Most satellites serve multiple purposes. Since I want to look at satellites by purpose, I’m going to separate those pruposes out. We’ve gotta be careful about not double-counting satellites if we add up total satellites, though!
Code
%>% separate_rows (purpose, sep = "/" ) %>% separate_rows (country_owner, sep = "/" ) %>% mutate (country_owner = fct_lump_n (country_owner, n = 15 ,other_level = "Other countries" ),purpose = fct_infreq (purpose)) %>% group_by (country_owner, purpose) %>% summarise (num_sats = n (),total_mass_dry_kg = sum (as.numeric (dry_mass_kg), na.rm = TRUE )) ->
Warning in mask$eval_all_summarise(quo): NAs introduced by coercion
Warning in mask$eval_all_summarise(quo): NAs introduced by coercion
Warning in mask$eval_all_summarise(quo): NAs introduced by coercion
Warning in mask$eval_all_summarise(quo): NAs introduced by coercion
`summarise()` has grouped output by 'country_owner'. You can override using the
`.groups` argument.
Let’s have a look at country of ownership and purpose:
Code
%>% filter (country_owner != "USA" ) %>% ggplot (.) + aes (x = reorder (country_owner, num_sats),y = num_sats,fill = purpose) + geom_col () + # scale_fill_discrete(guide = NULL) + coord_flip () + theme_360 () + theme (plot.subtitle = element_markdown (family = "Body 360info" , face = "plain" ),panel.grid.major.y = element_blank (),panel.grid.minor.y = element_blank (),legend.position = c (0.6 , 0.2 )) + labs (y = "Country of owner organisation (No USA)" ,x = "Number of satellites" ,title = toupper ("Number of active satellites" ),subtitle = "Satellites counted more than once for multiple owners and purposes" ,fill = "Purpose" ,caption = paste ("**CHART:** James Goldie, 360info" ,"**DATA:** Union of Concerned Scientists" ,sep = "<br>" ))%>% save_360plot (here ("out" , "sat-count-separated-nousa.png" ), shape = "sdtv-portrait" ):: include_graphics (here ("out" , "sat-count-separated-nousa.png" ))
Overlap of commercial and government/military use
But what about users? Satellite users include NA. There’s potentially a lot of overlap between a few discrete groups here!
Let’s investigate that overlap:
Code
%>% mutate (rowid = 1 : n ()) %>% select (rowid, country_owner, users) %>% separate_rows (users, sep = "/" ) %>% filter (! is.na (users)) %>% # test: USA only filter (country_owner == "USA" ) %>% # get a vector of all the row ids in each users group group_by (users) %>% summarise (rows = list (c (rowid))) -> # https://cran.r-project.org/web/packages/nVennR/vignettes/nVennR.html plotVenn (us_group_overlap$ rows, sNames = us_group_overlap$ users,nCycles = 5000 )
$def
[1] "nVenn1.2" "4" "Civil" "Commercial" "Government"
[6] "Military" "0" "189" "63" "1"
[11] "2503" "35" "78" "0" "28"
[16] "6" "19" "0" "4" "0"
[21] "0" "0"
$set
[1] "_F\n0\n_L\n17.9057\n26.4125\n17.8559\n26.3751\n17.7939\n26.361\n17.7328\n26.3733\n17.6808\n26.4103\n17.6499\n26.4636\n17.6424\n26.5256\n17.6196\n26.5777\n17.5821\n26.6213\n17.5319\n26.6516\n17.4738\n26.6669\n17.4176\n26.6868\n17.3645\n26.7132\n17.3114\n26.7446\n17.2613\n26.7806\n17.2134\n26.8209\n17.1646\n26.8599\n17.1104\n26.8912\n17.0522\n26.9112\n16.9917\n26.9189\n16.931\n26.9146\n16.8714\n26.8985\n16.8161\n26.8705\n16.7678\n26.8325\n16.7286\n26.7847\n16.6996\n26.7304\n16.681\n26.672\n16.6754\n26.6101\n16.6818\n26.5482\n16.7008\n26.4893\n16.7315\n26.4352\n16.7721\n26.3877\n16.8158\n26.3436\n16.8597\n26.2996\n16.9037\n26.2565\n16.9468\n26.2125\n16.9926\n26.1693\n17.0363\n26.1263\n17.0811\n26.0833\n17.1298\n26.0425\n17.1851\n26.0213\n17.245\n26.0206\n17.303\n26.0402\n17.3541\n26.0723\n17.404\n26.1035\n17.4536\n26.1357\n17.504\n26.1691\n17.555\n26.2024\n17.6038\n26.2368\n17.6537\n26.2705\n17.7073\n26.3013\n17.7648\n26.3212\n17.8252\n26.3308\n17.8854\n26.3287\n17.9444\n26.3161\n17.9995\n26.294\n18.0558\n26.2712\n18.1105\n26.2479\n18.1681\n26.2343\n18.2285\n26.2331\n18.2867\n26.2451\n18.3408\n26.2699\n18.3892\n26.3098\n18.4426\n26.383\n18.4654\n26.4409\n18.4725\n26.502\n18.4652\n26.5643\n18.4439\n26.6226\n18.41\n26.6747\n18.3657\n26.7167\n18.3124\n26.7483\n18.2544\n26.7669\n18.1936\n26.7717\n18.1328\n26.7628\n18.0763\n26.7402\n18.0249\n26.7052\n17.984\n26.6599\n17.9544\n26.6078\n17.9375\n26.5477\n17.9351\n26.4684\n_L\n18.59\n25.4159\n18.5013\n25.584\n18.4003\n25.7418\n18.2991\n25.8988\n18.1978\n26.0571\n18.0957\n26.2147\n17.9298\n26.318\n17.735\n26.3127\n17.5428\n26.301\n17.334\n26.3243\n17.1124\n26.3336\n16.9172\n26.3128\n16.7373\n26.361\n16.5569\n26.4076\n16.3766\n26.4531\n16.1951\n26.4987\n16.011\n26.5281\n15.8239\n26.5412\n15.6364\n26.5375\n15.4511\n26.5173\n15.2679\n26.4805\n15.0894\n26.4279\n14.9171\n26.3598\n14.7516\n26.2769\n14.5939\n26.1797\n14.445\n26.0686\n14.3063\n25.9447\n14.1798\n25.8093\n14.066\n25.6638\n13.9651\n25.5084\n13.8785\n25.3451\n13.8067\n25.1746\n13.75\n24.9985\n13.7093\n24.8182\n13.6847\n24.6355\n13.6761\n24.4521\n13.6836\n24.2675\n13.7074\n24.0846\n13.7473\n23.9043\n13.8028\n23.7286\n13.8732\n23.559\n13.9584\n23.3953\n14.0581\n23.2397\n14.1713\n23.0933\n14.2968\n22.9575\n14.4335\n22.8335\n14.5804\n22.722\n14.7361\n22.6243\n14.9005\n22.5403\n15.0716\n22.4709\n15.2468\n22.4174\n15.4261\n22.3792\n15.6088\n22.3569\n15.7933\n22.3508\n15.9776\n22.3609\n16.1607\n22.3872\n16.342\n22.4297\n16.519\n22.4884\n16.6903\n22.5627\n16.854\n22.6518\n17.0081\n22.7543\n17.1525\n22.8698\n17.286\n22.9976\n17.4077\n23.1367\n17.5169\n23.2861\n17.6125\n23.4445\n17.705\n23.6034\n17.7968\n23.7624\n17.8886\n23.9227\n17.9804\n24.0836\n18.0733\n24.2439\n18.1658\n24.4046\n18.2579\n24.5654\n18.3505\n24.7259\n18.4431\n24.8869\n18.5363\n25.0486\n18.6025\n25.227\n_L\n18.3231\n26.7433\n18.2228\n26.7703\n18.1185\n26.758\n18.0213\n26.7243\n17.9244\n26.6895\n17.8279\n26.6573\n17.7164\n26.6389\n17.6515\n26.5623\n17.6493\n26.464\n17.6194\n26.3695\n17.5894\n26.275\n17.5573\n26.1796\n17.5241\n26.0837\n17.5138\n25.9822\n17.5333\n25.8819\n17.5595\n25.7859\n17.5863\n25.6878\n17.6137\n25.5897\n17.6405\n25.4926\n17.6684\n25.395\n17.6963\n25.2979\n17.7237\n25.2012\n17.7666\n25.0817\n17.8255\n24.9938\n17.9058\n24.923\n18.0006\n24.8734\n18.1037\n24.8476\n18.2106\n24.8469\n18.3141\n24.8715\n18.41\n24.9204\n18.4915\n24.9907\n18.5542\n25.0786\n18.5939\n25.1782\n18.6084\n25.283\n18.5981\n25.3888\n18.5634\n25.4889\n18.5063\n25.5774\n18.4308\n25.6502\n18.3397\n25.7038\n18.2367\n25.7346\n18.1307\n25.7403\n18.0281\n25.7213\n17.9301\n25.6782\n17.8248\n25.6676\n17.7224\n25.6906\n17.6318\n25.7446\n17.5637\n25.8231\n17.5235\n25.9168\n17.514\n26.0193\n17.5365\n26.1209\n17.5896\n26.2107\n17.668\n26.2798\n17.7631\n26.3205\n17.8667\n26.3302\n17.9683\n26.3069\n18.0569\n26.2547\n18.1263\n26.177\n18.1954\n26.1006\n18.2663\n26.024\n18.3375\n25.946\n18.4065\n25.867\n18.4772\n25.79\n18.5631\n25.7335\n18.6594\n25.7\n18.7615\n25.6917\n18.8633\n25.7084\n18.9581\n25.7506\n19.0389\n25.8143\n19.1017\n25.8961\n19.1422\n25.9912\n19.1574\n26.0934\n19.1476\n26.1961\n19.1113\n26.2945\n19.0526\n26.3794\n18.9758\n26.4463\n18.8838\n26.4929\n18.7903\n26.535\n18.6961\n26.5759\n18.6022\n26.6175\n18.5095\n26.6594\n18.4163\n26.701\n_L\n18.2553\n27.5514\n18.2206\n27.6176\n18.1779\n27.6794\n18.1281\n27.7352\n18.0717\n27.7863\n18.0096\n27.8303\n17.9438\n27.8661\n17.8753\n27.8932\n17.8036\n27.914\n17.7288\n27.9247\n17.6533\n27.9276\n17.5784\n27.9215\n17.5047\n27.9064\n17.4322\n27.8829\n17.3638\n27.8513\n17.2997\n27.8114\n17.2411\n27.7654\n17.1891\n27.7141\n17.1422\n27.6556\n17.1025\n27.5924\n17.0707\n27.5254\n17.0465\n27.4547\n17.0311\n27.3814\n17.0236\n27.3085\n17.0251\n27.2343\n17.0351\n27.1599\n17.0531\n27.0879\n17.0798\n27.0186\n17.1135\n26.9523\n17.1547\n26.8916\n17.1979\n26.8338\n17.2425\n26.7759\n17.274\n26.7095\n17.2888\n26.6373\n17.2857\n26.565\n17.2638\n26.4874\n17.273\n26.4125\n17.3097\n26.3472\n17.3616\n26.2936\n17.3905\n26.2279\n17.4197\n26.1621\n17.448\n26.0954\n17.4755\n26.0285\n17.502\n25.9614\n17.5283\n25.8952\n17.5597\n25.8293\n17.604\n25.7704\n17.6603\n25.7228\n17.7251\n25.6891\n17.7955\n25.6704\n17.8677\n25.6676\n17.9379\n25.6813\n18.0187\n25.7173\n18.0764\n25.7604\n18.1226\n25.8161\n18.1556\n25.8811\n18.1738\n25.9525\n18.1763\n26.0253\n18.1624\n26.0979\n18.1325\n26.1652\n18.0892\n26.2246\n18.0399\n26.2789\n17.9915\n26.3337\n17.9542\n26.3963\n17.9358\n26.4669\n17.9375\n26.5504\n17.961\n26.6214\n18.0041\n26.6849\n18.0539\n26.7401\n18.102\n26.797\n18.1515\n26.8529\n18.1978\n26.9121\n18.2371\n26.9766\n18.2686\n27.0443\n18.2918\n27.114\n18.3066\n27.186\n18.3145\n27.2599\n18.3122\n27.3344\n18.3016\n27.4085\n18.2824\n27.4814\n_C\n0\n4.29497e+09\n0\n16.9828\n26.6118\n0.211533\n17.2175\n26.1939\n0.079952\n0\n12\n0\n0\n18\n0\n15.7717\n24.4463\n2\n18.16\n25.2927\n0.353059\n6\n12\n0\n17.8452\n25.9989\n0.236501\n18.745\n26.1039\n0.3173\n18.2026\n26.5018\n0.174251\n18\n12\n0\n17.7922\n26.5113\n0.0548131\n17.6687\n27.2822\n0.549579\n17.4563\n26.4739\n0.0979209\n53.9925\n28.9408\n0\n"
$reg
$reg[[1]]
[1] 12 13 14 15 16 17 18 19 20 21 22 23 32 33 34
[16] 55 290 295 305 306 314 331 334 335 341 400 420 421 422 423
[31] 424 594 595 596 597 598 599 600 601 602 603 605 654 690 728
[46] 729 730 731 732 736 894 951 952 953 955 956 957 965 967 1051
[61] 1052 1053 1054 1277 1278 1279 1280 1281 1282 1313 1314 1335 1488 1489 1497
[76] 1498 1499 1524 1525 1526 1527 1535 1544 1545 1546 1547 1548 1634 1635 1636
[91] 1637 1638 1639 1640 1641 1642 1695 1697 2141 2142 2143 2144 2145 2146 2155
[106] 2156 2178 2179 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2232
[121] 2241 2242 2243 2312 2313 2314 2315 2316 2317 2318 2319 2320 2321 2322 2323
[136] 2324 2325 2326 2327 2328 2329 2336 2337 2338 2339 2340 2341 2342 2346 2372
[151] 2463 2465 2466 2467 2468 2470 2471 2472 2473 4448 4459 4460 4461 4479 4480
[166] 4488 4593 4618 4619 4620 4621 4622 4627 4628 4629 4630 4631 4632 4657 4658
[181] 4659 4660 4661 4662 4663 4664 4665 4666 4674
$reg[[2]]
[1] 24 46 299 307 310 378 379 380 381 382 383 384 385 386 387
[16] 645 727 738 815 895 896 897 898 1026 1039 1109 1186 1196 1200 1326
[31] 1327 1328 1506 1529 1530 1531 1532 1622 1623 1624 1625 1632 1672 1694 2148
[46] 2201 2236 2335 2373 2464 2485 4489 4490 4491 4492 4493 4496 4497 4498 4513
[61] 4592 4635 4636
$reg[[3]]
[1] 1543
$reg[[4]]
[1] 1 26 27 28 29 30 31 35 36 37 38 39 40 41
[15] 42 44 45 65 66 67 68 69 70 71 72 73 108 120
[29] 141 153 154 230 231 232 233 234 235 236 237 238 239 240
[43] 241 250 273 274 275 276 277 278 322 323 324 325 326 327
[57] 342 343 344 345 346 406 407 408 409 410 411 412 413 426
[71] 427 428 429 430 431 432 433 434 435 436 437 438 439 440
[85] 441 442 443 444 445 446 447 448 449 450 451 452 453 454
[99] 455 456 457 458 459 460 461 462 463 464 465 466 467 468
[113] 469 470 471 472 473 474 475 476 477 478 479 480 481 482
[127] 483 484 485 486 487 488 489 490 491 492 493 494 495 496
[141] 497 498 499 500 501 502 503 504 505 506 507 508 509 510
[155] 511 512 513 514 515 516 517 518 519 520 521 522 523 524
[169] 525 526 527 528 529 530 531 532 533 534 535 536 537 538
[183] 539 540 541 542 543 544 545 546 547 548 549 550 551 552
[197] 553 554 555 556 557 558 559 560 561 562 563 564 565 566
[211] 567 568 569 570 571 572 573 574 575 576 577 578 579 580
[225] 581 582 583 584 585 586 587 588 589 590 591 592 608 610
[239] 611 612 613 614 615 616 617 618 619 620 621 622 631 744
[253] 745 747 748 749 750 751 752 753 754 755 756 818 819 831
[267] 832 833 834 835 836 837 838 839 840 841 842 843 844 845
[281] 846 847 848 849 850 851 852 853 854 855 856 857 858 859
[295] 860 861 862 892 893 968 969 970 971 972 973 974 975 976
[309] 1017 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087
[323] 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101
[337] 1102 1103 1104 1105 1106 1107 1324 1325 1337 1338 1339 1340 1341 1342
[351] 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356
[365] 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370
[379] 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384
[393] 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398
[407] 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412
[421] 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426
[435] 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440
[449] 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454
[463] 1455 1456 1461 1462 1485 1516 1517 1549 2105 2106 2107 2108 2109 2110
[477] 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124
[491] 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135 2136 2137 2138
[505] 2139 2147 2154 2157 2159 2160 2213 2418 2419 2421 2422 2423 2425 2426
[519] 2427 2428 2429 2430 2441 2442 2443 2444 2445 2446 2447 2448 2449 2450
[533] 2451 2452 2453 2454 2455 2456 2457 2458 2459 2460 2461 2462 2487 2488
[547] 2489 2490 2491 2492 2493 2494 2495 2496 2497 2498 2499 2500 2501 2502
[561] 2503 2504 2505 2506 2507 2508 2509 2510 2511 2512 2513 2514 2515 2516
[575] 2517 2518 2519 2520 2521 2522 2523 2524 2525 2526 2527 2528 2529 2530
[589] 2531 2532 2533 2534 2535 2536 2537 2538 2539 2540 2541 2542 2543 2544
[603] 2545 2546 2547 2548 2549 2550 2551 2552 2553 2554 2555 2556 2557 2558
[617] 2559 2560 2561 2562 2563 2564 2565 2566 2567 2568 2569 2570 2571 2572
[631] 2573 2574 2575 2576 2577 2578 2579 2580 2581 2582 2583 2584 2585 2586
[645] 2587 2588 2589 2590 2591 2592 2593 2594 2595 2596 2597 2598 2599 2600
[659] 2601 2602 2603 2604 2605 2606 2607 2633 2634 2635 2636 2637 2638 2639
[673] 2640 2641 2642 2643 2644 2645 2646 2647 2648 2649 2650 2651 2652 2653
[687] 2654 2655 2656 2657 2658 2659 2660 2661 2662 2663 2664 2665 2666 2667
[701] 2668 2669 2670 2671 2672 2673 2674 2675 2676 2677 2678 2679 2680 2681
[715] 2682 2683 2684 2685 2686 2687 2688 2689 2690 2691 2692 2693 2694 2695
[729] 2696 2697 2698 2699 2700 2701 2702 2703 2704 2705 2706 2707 2708 2709
[743] 2710 2711 2712 2713 2714 2715 2716 2717 2718 2719 2720 2721 2722 2723
[757] 2724 2725 2726 2727 2728 2729 2730 2731 2732 2733 2734 2735 2736 2737
[771] 2738 2739 2740 2741 2742 2743 2744 2745 2746 2747 2748 2749 2750 2751
[785] 2752 2753 2754 2755 2756 2757 2758 2759 2760 2761 2762 2763 2764 2765
[799] 2766 2767 2768 2769 2770 2771 2772 2773 2774 2775 2776 2777 2778 2779
[813] 2780 2781 2782 2783 2784 2785 2786 2787 2788 2789 2790 2791 2792 2793
[827] 2794 2795 2796 2797 2798 2799 2800 2801 2802 2803 2804 2805 2806 2807
[841] 2808 2809 2810 2811 2812 2813 2814 2815 2816 2817 2818 2819 2820 2821
[855] 2822 2823 2824 2825 2826 2827 2828 2829 2830 2831 2832 2833 2834 2835
[869] 2836 2837 2838 2839 2840 2841 2842 2843 2844 2845 2846 2847 2848 2849
[883] 2850 2851 2852 2853 2854 2855 2856 2857 2858 2859 2860 2861 2862 2863
[897] 2864 2865 2866 2867 2868 2869 2870 2871 2872 2873 2874 2875 2876 2877
[911] 2878 2879 2880 2881 2882 2883 2884 2885 2886 2887 2888 2889 2890 2891
[925] 2892 2893 2894 2895 2896 2897 2898 2899 2900 2901 2902 2903 2904 2905
[939] 2906 2907 2908 2909 2910 2911 2912 2913 2914 2915 2916 2917 2918 2919
[953] 2920 2921 2922 2923 2924 2925 2926 2927 2928 2929 2930 2931 2932 2933
[967] 2934 2935 2936 2937 2938 2939 2940 2941 2942 2943 2944 2945 2946 2947
[981] 2948 2949 2950 2951 2952 2953 2954 2955 2956 2957 2958 2959 2960 2961
[995] 2962 2963 2964 2965 2966 2967 2968 2969 2970 2971 2972 2973 2974 2975
[1009] 2976 2977 2978 2979 2980 2981 2982 2983 2984 2985 2986 2987 2988 2989
[1023] 2990 2991 2992 2993 2994 2995 2996 2997 2998 2999 3000 3001 3002 3003
[1037] 3004 3005 3006 3007 3008 3009 3010 3011 3012 3013 3014 3015 3016 3017
[1051] 3018 3019 3020 3021 3022 3023 3024 3025 3026 3027 3028 3029 3030 3031
[1065] 3032 3033 3034 3035 3036 3037 3038 3039 3040 3041 3042 3043 3044 3045
[1079] 3046 3047 3048 3049 3050 3051 3052 3053 3054 3055 3056 3057 3058 3059
[1093] 3060 3061 3062 3063 3064 3065 3066 3067 3068 3069 3070 3071 3072 3073
[1107] 3074 3075 3076 3077 3078 3079 3080 3081 3082 3083 3084 3085 3086 3087
[1121] 3088 3089 3090 3091 3092 3093 3094 3095 3096 3097 3098 3099 3100 3101
[1135] 3102 3103 3104 3105 3106 3107 3108 3109 3110 3111 3112 3113 3114 3115
[1149] 3116 3117 3118 3119 3120 3121 3122 3123 3124 3125 3126 3127 3128 3129
[1163] 3130 3131 3132 3133 3134 3135 3136 3137 3138 3139 3140 3141 3142 3143
[1177] 3144 3145 3146 3147 3148 3149 3150 3151 3152 3153 3154 3155 3156 3157
[1191] 3158 3159 3160 3161 3162 3163 3164 3165 3166 3167 3168 3169 3170 3171
[1205] 3172 3173 3174 3175 3176 3177 3178 3179 3180 3181 3182 3183 3184 3185
[1219] 3186 3187 3188 3189 3190 3191 3192 3193 3194 3195 3196 3197 3198 3199
[1233] 3200 3201 3202 3203 3204 3205 3206 3207 3208 3209 3210 3211 3212 3213
[1247] 3214 3215 3216 3217 3218 3219 3220 3221 3222 3223 3224 3225 3226 3227
[1261] 3228 3229 3230 3231 3232 3233 3234 3235 3236 3237 3238 3239 3240 3241
[1275] 3242 3243 3244 3245 3246 3247 3248 3249 3250 3251 3252 3253 3254 3255
[1289] 3256 3257 3258 3259 3260 3261 3262 3263 3264 3265 3266 3267 3268 3269
[1303] 3270 3271 3272 3273 3274 3275 3276 3277 3278 3279 3280 3281 3282 3283
[1317] 3284 3285 3286 3287 3288 3289 3290 3291 3292 3293 3294 3295 3296 3297
[1331] 3298 3299 3300 3301 3302 3303 3304 3305 3306 3307 3308 3309 3310 3311
[1345] 3312 3313 3314 3315 3316 3317 3318 3319 3320 3321 3322 3323 3324 3325
[1359] 3326 3327 3328 3329 3330 3331 3332 3333 3334 3335 3336 3337 3338 3339
[1373] 3340 3341 3342 3343 3344 3345 3346 3347 3348 3349 3350 3351 3352 3353
[1387] 3354 3355 3356 3357 3358 3359 3360 3361 3362 3363 3364 3365 3366 3367
[1401] 3368 3369 3370 3371 3372 3373 3374 3375 3376 3377 3378 3379 3380 3381
[1415] 3382 3383 3384 3385 3386 3387 3388 3389 3390 3391 3392 3393 3394 3395
[1429] 3396 3397 3398 3399 3400 3401 3402 3403 3404 3405 3406 3407 3408 3409
[1443] 3410 3411 3412 3413 3414 3415 3416 3417 3418 3419 3420 3421 3422 3423
[1457] 3424 3425 3426 3427 3428 3429 3430 3431 3432 3433 3434 3435 3436 3437
[1471] 3438 3439 3440 3441 3442 3443 3444 3445 3446 3447 3448 3449 3450 3451
[1485] 3452 3453 3454 3455 3456 3457 3458 3459 3460 3461 3462 3463 3464 3465
[1499] 3466 3467 3468 3469 3470 3471 3472 3473 3474 3475 3476 3477 3478 3479
[1513] 3480 3481 3482 3483 3484 3485 3486 3487 3488 3489 3490 3491 3492 3493
[1527] 3494 3495 3496 3497 3498 3499 3500 3501 3502 3503 3504 3505 3506 3507
[1541] 3508 3509 3510 3511 3512 3513 3514 3515 3516 3517 3518 3519 3520 3521
[1555] 3522 3523 3524 3525 3526 3527 3528 3529 3530 3531 3532 3533 3534 3535
[1569] 3536 3537 3538 3539 3540 3541 3542 3543 3544 3545 3546 3547 3548 3549
[1583] 3550 3551 3552 3553 3554 3555 3556 3557 3558 3559 3560 3561 3562 3563
[1597] 3564 3565 3566 3567 3568 3569 3570 3571 3572 3573 3574 3575 3576 3577
[1611] 3578 3579 3580 3581 3582 3583 3584 3585 3586 3587 3588 3589 3590 3591
[1625] 3592 3593 3594 3595 3596 3597 3598 3599 3600 3601 3602 3603 3604 3605
[1639] 3606 3607 3608 3609 3610 3611 3612 3613 3614 3615 3616 3617 3618 3619
[1653] 3620 3621 3622 3623 3624 3625 3626 3627 3628 3629 3630 3631 3632 3633
[1667] 3634 3635 3636 3637 3638 3639 3640 3641 3642 3643 3644 3645 3646 3647
[1681] 3648 3649 3650 3651 3652 3653 3654 3655 3656 3657 3658 3659 3660 3661
[1695] 3662 3663 3664 3665 3666 3667 3668 3669 3670 3671 3672 3673 3674 3675
[1709] 3676 3677 3678 3679 3680 3681 3682 3683 3684 3685 3686 3687 3688 3689
[1723] 3690 3691 3692 3693 3694 3695 3696 3697 3698 3699 3700 3701 3702 3703
[1737] 3704 3705 3706 3707 3708 3709 3710 3711 3712 3713 3714 3715 3716 3717
[1751] 3718 3719 3720 3721 3722 3723 3724 3725 3726 3727 3728 3729 3730 3731
[1765] 3732 3733 3734 3735 3736 3737 3738 3739 3740 3741 3742 3743 3744 3745
[1779] 3746 3747 3748 3749 3750 3751 3752 3753 3754 3755 3756 3757 3758 3759
[1793] 3760 3761 3762 3763 3764 3765 3766 3767 3768 3769 3770 3771 3772 3773
[1807] 3774 3775 3776 3777 3778 3779 3780 3781 3782 3783 3784 3785 3786 3787
[1821] 3788 3789 3790 3791 3792 3793 3794 3795 3796 3797 3798 3799 3800 3801
[1835] 3802 3803 3804 3805 3806 3807 3808 3809 3810 3811 3812 3813 3814 3815
[1849] 3816 3817 3818 3819 3820 3821 3822 3823 3824 3825 3826 3827 3828 3829
[1863] 3830 3831 3832 3833 3834 3835 3836 3837 3838 3839 3840 3841 3842 3843
[1877] 3844 3845 3846 3847 3848 3849 3850 3851 3852 3853 3854 3855 3856 3857
[1891] 3858 3859 3860 3861 3862 3863 3864 3865 3866 3867 3868 3869 3870 3871
[1905] 3872 3873 3874 3875 3876 3877 3878 3879 3880 3881 3882 3883 3884 3885
[1919] 3886 3887 3888 3889 3890 3891 3892 3893 3894 3895 3896 3897 3898 3899
[1933] 3900 3901 3902 3903 3904 3905 3906 3907 3908 3909 3910 3911 3912 3913
[1947] 3914 3915 3916 3917 3918 3919 3920 3921 3922 3923 3924 3925 3926 3927
[1961] 3928 3929 3930 3931 3932 3933 3934 3935 3936 3937 3938 3939 3940 3941
[1975] 3942 3943 3944 3945 3946 3947 3948 3949 3950 3951 3952 3953 3954 3955
[1989] 3956 3957 3958 3959 3960 3961 3962 3963 3964 3965 3966 3967 3968 3969
[2003] 3970 3971 3972 3973 3974 3975 3976 3977 3978 3979 3980 3981 3982 3983
[2017] 3984 3985 3986 3987 3988 3989 3990 3991 3992 3993 3994 3995 3996 3997
[2031] 3998 3999 4000 4001 4002 4003 4004 4005 4006 4007 4008 4009 4010 4011
[2045] 4012 4013 4014 4015 4016 4017 4018 4019 4020 4021 4022 4023 4024 4025
[2059] 4026 4027 4028 4029 4030 4031 4032 4033 4034 4035 4036 4037 4038 4039
[2073] 4040 4041 4042 4043 4044 4045 4046 4047 4048 4049 4050 4051 4052 4053
[2087] 4054 4055 4056 4057 4058 4059 4060 4061 4062 4063 4064 4065 4066 4067
[2101] 4068 4069 4070 4071 4072 4073 4074 4075 4076 4077 4078 4079 4080 4081
[2115] 4082 4083 4084 4085 4086 4087 4088 4089 4090 4091 4092 4093 4094 4095
[2129] 4096 4097 4098 4099 4100 4101 4102 4103 4104 4105 4106 4107 4108 4109
[2143] 4110 4111 4112 4113 4114 4115 4116 4117 4118 4119 4120 4121 4122 4123
[2157] 4124 4125 4126 4127 4128 4129 4130 4131 4132 4133 4134 4135 4136 4137
[2171] 4138 4139 4140 4141 4142 4143 4144 4145 4146 4147 4148 4149 4150 4151
[2185] 4152 4153 4154 4155 4156 4157 4158 4159 4160 4161 4162 4163 4164 4165
[2199] 4166 4167 4168 4169 4170 4171 4172 4173 4174 4175 4176 4177 4178 4179
[2213] 4180 4181 4182 4183 4184 4185 4186 4187 4188 4189 4190 4191 4192 4193
[2227] 4194 4195 4196 4197 4198 4199 4200 4201 4202 4203 4204 4205 4206 4207
[2241] 4208 4209 4210 4211 4212 4213 4214 4215 4216 4217 4218 4219 4220 4221
[2255] 4222 4223 4224 4225 4226 4227 4228 4229 4230 4231 4232 4233 4234 4235
[2269] 4236 4237 4238 4239 4240 4241 4242 4243 4244 4245 4246 4247 4248 4249
[2283] 4250 4251 4252 4253 4254 4255 4256 4257 4258 4259 4260 4261 4262 4263
[2297] 4264 4265 4266 4267 4268 4269 4270 4271 4272 4273 4274 4275 4276 4277
[2311] 4278 4279 4280 4281 4282 4283 4284 4285 4286 4287 4288 4289 4290 4291
[2325] 4292 4293 4294 4295 4296 4297 4298 4299 4300 4301 4302 4303 4304 4305
[2339] 4306 4307 4308 4309 4310 4311 4312 4313 4314 4315 4316 4317 4318 4319
[2353] 4320 4321 4322 4323 4324 4325 4326 4327 4328 4329 4330 4331 4332 4333
[2367] 4334 4335 4336 4337 4338 4339 4340 4341 4342 4343 4344 4345 4346 4347
[2381] 4348 4349 4350 4351 4352 4353 4354 4355 4356 4357 4358 4359 4360 4361
[2395] 4362 4363 4364 4365 4366 4367 4368 4369 4370 4371 4372 4373 4374 4375
[2409] 4376 4377 4378 4379 4380 4381 4382 4383 4384 4385 4386 4387 4388 4389
[2423] 4390 4391 4392 4393 4394 4395 4396 4397 4398 4399 4400 4401 4402 4403
[2437] 4404 4405 4406 4407 4408 4409 4410 4411 4412 4413 4414 4415 4416 4417
[2451] 4418 4419 4420 4421 4422 4423 4424 4425 4426 4427 4428 4429 4430 4431
[2465] 4432 4433 4434 4435 4436 4437 4438 4439 4440 4441 4442 4443 4444 4445
[2479] 4446 4447 4485 4512 4518 4605 4608 4609 4610 4611 4612 4613 4614 4615
[2493] 4616 4617 4624 4642 4643 4667 4671 4672 4673 4704 4705
$reg[[5]]
[1] 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571
[16] 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586
[31] 1587 1588 1589 1590 1591
$reg[[6]]
[1] 822 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124
[16] 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139
[31] 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154
[46] 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169
[61] 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184
[76] 1185 4494 4495
$reg[[7]]
NULL
$reg[[8]]
[1] 78 80 103 104 105 109 228 247 272 329 372 376 391 392 398
[16] 740 1457 1459 1492 2140 2197 2224 2229 2230 2231 2344 4509 4638
$reg[[9]]
[1] 710 711 712 713 715 1633
$reg[[10]]
[1] 266 316 332 369 691 733 734 989 1198 1522 1673 2169 2170 2180 2484
[16] 4581 4647 4648 4649
$reg[[11]]
NULL
$reg[[12]]
[1] 330 635 636 1487
$orig
$orig$Civil
[1] 78 80 103 104 105 109 228 247 266 272 316 329 330 332 369
[16] 372 376 391 392 398 635 636 691 710 711 712 713 715 733 734
[31] 740 989 1198 1457 1459 1487 1492 1522 1633 1673 2140 2169 2170 2180 2197
[46] 2224 2229 2230 2231 2344 2484 4509 4581 4638 4647 4648 4649
$orig$Commercial
[1] 1 26 27 28 29 30 31 35 36 37 38 39 40 41
[15] 42 44 45 65 66 67 68 69 70 71 72 73 108 120
[29] 141 153 154 230 231 232 233 234 235 236 237 238 239 240
[43] 241 250 273 274 275 276 277 278 322 323 324 325 326 327
[57] 330 342 343 344 345 346 406 407 408 409 410 411 412 413
[71] 426 427 428 429 430 431 432 433 434 435 436 437 438 439
[85] 440 441 442 443 444 445 446 447 448 449 450 451 452 453
[99] 454 455 456 457 458 459 460 461 462 463 464 465 466 467
[113] 468 469 470 471 472 473 474 475 476 477 478 479 480 481
[127] 482 483 484 485 486 487 488 489 490 491 492 493 494 495
[141] 496 497 498 499 500 501 502 503 504 505 506 507 508 509
[155] 510 511 512 513 514 515 516 517 518 519 520 521 522 523
[169] 524 525 526 527 528 529 530 531 532 533 534 535 536 537
[183] 538 539 540 541 542 543 544 545 546 547 548 549 550 551
[197] 552 553 554 555 556 557 558 559 560 561 562 563 564 565
[211] 566 567 568 569 570 571 572 573 574 575 576 577 578 579
[225] 580 581 582 583 584 585 586 587 588 589 590 591 592 608
[239] 610 611 612 613 614 615 616 617 618 619 620 621 622 631
[253] 635 636 744 745 747 748 749 750 751 752 753 754 755 756
[267] 818 819 822 831 832 833 834 835 836 837 838 839 840 841
[281] 842 843 844 845 846 847 848 849 850 851 852 853 854 855
[295] 856 857 858 859 860 861 862 892 893 968 969 970 971 972
[309] 973 974 975 976 1017 1075 1076 1077 1078 1079 1080 1081 1082 1083
[323] 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097
[337] 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1111 1112 1113 1114
[351] 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128
[365] 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142
[379] 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156
[393] 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170
[407] 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184
[421] 1185 1324 1325 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347
[435] 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361
[449] 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375
[463] 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389
[477] 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403
[491] 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417
[505] 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431
[519] 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445
[533] 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1461 1462 1485
[547] 1487 1516 1517 1549 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566
[561] 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580
[575] 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 2105 2106 2107
[589] 2108 2109 2110 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121
[603] 2122 2123 2124 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135
[617] 2136 2137 2138 2139 2147 2154 2157 2159 2160 2213 2418 2419 2421 2422
[631] 2423 2425 2426 2427 2428 2429 2430 2441 2442 2443 2444 2445 2446 2447
[645] 2448 2449 2450 2451 2452 2453 2454 2455 2456 2457 2458 2459 2460 2461
[659] 2462 2487 2488 2489 2490 2491 2492 2493 2494 2495 2496 2497 2498 2499
[673] 2500 2501 2502 2503 2504 2505 2506 2507 2508 2509 2510 2511 2512 2513
[687] 2514 2515 2516 2517 2518 2519 2520 2521 2522 2523 2524 2525 2526 2527
[701] 2528 2529 2530 2531 2532 2533 2534 2535 2536 2537 2538 2539 2540 2541
[715] 2542 2543 2544 2545 2546 2547 2548 2549 2550 2551 2552 2553 2554 2555
[729] 2556 2557 2558 2559 2560 2561 2562 2563 2564 2565 2566 2567 2568 2569
[743] 2570 2571 2572 2573 2574 2575 2576 2577 2578 2579 2580 2581 2582 2583
[757] 2584 2585 2586 2587 2588 2589 2590 2591 2592 2593 2594 2595 2596 2597
[771] 2598 2599 2600 2601 2602 2603 2604 2605 2606 2607 2633 2634 2635 2636
[785] 2637 2638 2639 2640 2641 2642 2643 2644 2645 2646 2647 2648 2649 2650
[799] 2651 2652 2653 2654 2655 2656 2657 2658 2659 2660 2661 2662 2663 2664
[813] 2665 2666 2667 2668 2669 2670 2671 2672 2673 2674 2675 2676 2677 2678
[827] 2679 2680 2681 2682 2683 2684 2685 2686 2687 2688 2689 2690 2691 2692
[841] 2693 2694 2695 2696 2697 2698 2699 2700 2701 2702 2703 2704 2705 2706
[855] 2707 2708 2709 2710 2711 2712 2713 2714 2715 2716 2717 2718 2719 2720
[869] 2721 2722 2723 2724 2725 2726 2727 2728 2729 2730 2731 2732 2733 2734
[883] 2735 2736 2737 2738 2739 2740 2741 2742 2743 2744 2745 2746 2747 2748
[897] 2749 2750 2751 2752 2753 2754 2755 2756 2757 2758 2759 2760 2761 2762
[911] 2763 2764 2765 2766 2767 2768 2769 2770 2771 2772 2773 2774 2775 2776
[925] 2777 2778 2779 2780 2781 2782 2783 2784 2785 2786 2787 2788 2789 2790
[939] 2791 2792 2793 2794 2795 2796 2797 2798 2799 2800 2801 2802 2803 2804
[953] 2805 2806 2807 2808 2809 2810 2811 2812 2813 2814 2815 2816 2817 2818
[967] 2819 2820 2821 2822 2823 2824 2825 2826 2827 2828 2829 2830 2831 2832
[981] 2833 2834 2835 2836 2837 2838 2839 2840 2841 2842 2843 2844 2845 2846
[995] 2847 2848 2849 2850 2851 2852 2853 2854 2855 2856 2857 2858 2859 2860
[1009] 2861 2862 2863 2864 2865 2866 2867 2868 2869 2870 2871 2872 2873 2874
[1023] 2875 2876 2877 2878 2879 2880 2881 2882 2883 2884 2885 2886 2887 2888
[1037] 2889 2890 2891 2892 2893 2894 2895 2896 2897 2898 2899 2900 2901 2902
[1051] 2903 2904 2905 2906 2907 2908 2909 2910 2911 2912 2913 2914 2915 2916
[1065] 2917 2918 2919 2920 2921 2922 2923 2924 2925 2926 2927 2928 2929 2930
[1079] 2931 2932 2933 2934 2935 2936 2937 2938 2939 2940 2941 2942 2943 2944
[1093] 2945 2946 2947 2948 2949 2950 2951 2952 2953 2954 2955 2956 2957 2958
[1107] 2959 2960 2961 2962 2963 2964 2965 2966 2967 2968 2969 2970 2971 2972
[1121] 2973 2974 2975 2976 2977 2978 2979 2980 2981 2982 2983 2984 2985 2986
[1135] 2987 2988 2989 2990 2991 2992 2993 2994 2995 2996 2997 2998 2999 3000
[1149] 3001 3002 3003 3004 3005 3006 3007 3008 3009 3010 3011 3012 3013 3014
[1163] 3015 3016 3017 3018 3019 3020 3021 3022 3023 3024 3025 3026 3027 3028
[1177] 3029 3030 3031 3032 3033 3034 3035 3036 3037 3038 3039 3040 3041 3042
[1191] 3043 3044 3045 3046 3047 3048 3049 3050 3051 3052 3053 3054 3055 3056
[1205] 3057 3058 3059 3060 3061 3062 3063 3064 3065 3066 3067 3068 3069 3070
[1219] 3071 3072 3073 3074 3075 3076 3077 3078 3079 3080 3081 3082 3083 3084
[1233] 3085 3086 3087 3088 3089 3090 3091 3092 3093 3094 3095 3096 3097 3098
[1247] 3099 3100 3101 3102 3103 3104 3105 3106 3107 3108 3109 3110 3111 3112
[1261] 3113 3114 3115 3116 3117 3118 3119 3120 3121 3122 3123 3124 3125 3126
[1275] 3127 3128 3129 3130 3131 3132 3133 3134 3135 3136 3137 3138 3139 3140
[1289] 3141 3142 3143 3144 3145 3146 3147 3148 3149 3150 3151 3152 3153 3154
[1303] 3155 3156 3157 3158 3159 3160 3161 3162 3163 3164 3165 3166 3167 3168
[1317] 3169 3170 3171 3172 3173 3174 3175 3176 3177 3178 3179 3180 3181 3182
[1331] 3183 3184 3185 3186 3187 3188 3189 3190 3191 3192 3193 3194 3195 3196
[1345] 3197 3198 3199 3200 3201 3202 3203 3204 3205 3206 3207 3208 3209 3210
[1359] 3211 3212 3213 3214 3215 3216 3217 3218 3219 3220 3221 3222 3223 3224
[1373] 3225 3226 3227 3228 3229 3230 3231 3232 3233 3234 3235 3236 3237 3238
[1387] 3239 3240 3241 3242 3243 3244 3245 3246 3247 3248 3249 3250 3251 3252
[1401] 3253 3254 3255 3256 3257 3258 3259 3260 3261 3262 3263 3264 3265 3266
[1415] 3267 3268 3269 3270 3271 3272 3273 3274 3275 3276 3277 3278 3279 3280
[1429] 3281 3282 3283 3284 3285 3286 3287 3288 3289 3290 3291 3292 3293 3294
[1443] 3295 3296 3297 3298 3299 3300 3301 3302 3303 3304 3305 3306 3307 3308
[1457] 3309 3310 3311 3312 3313 3314 3315 3316 3317 3318 3319 3320 3321 3322
[1471] 3323 3324 3325 3326 3327 3328 3329 3330 3331 3332 3333 3334 3335 3336
[1485] 3337 3338 3339 3340 3341 3342 3343 3344 3345 3346 3347 3348 3349 3350
[1499] 3351 3352 3353 3354 3355 3356 3357 3358 3359 3360 3361 3362 3363 3364
[1513] 3365 3366 3367 3368 3369 3370 3371 3372 3373 3374 3375 3376 3377 3378
[1527] 3379 3380 3381 3382 3383 3384 3385 3386 3387 3388 3389 3390 3391 3392
[1541] 3393 3394 3395 3396 3397 3398 3399 3400 3401 3402 3403 3404 3405 3406
[1555] 3407 3408 3409 3410 3411 3412 3413 3414 3415 3416 3417 3418 3419 3420
[1569] 3421 3422 3423 3424 3425 3426 3427 3428 3429 3430 3431 3432 3433 3434
[1583] 3435 3436 3437 3438 3439 3440 3441 3442 3443 3444 3445 3446 3447 3448
[1597] 3449 3450 3451 3452 3453 3454 3455 3456 3457 3458 3459 3460 3461 3462
[1611] 3463 3464 3465 3466 3467 3468 3469 3470 3471 3472 3473 3474 3475 3476
[1625] 3477 3478 3479 3480 3481 3482 3483 3484 3485 3486 3487 3488 3489 3490
[1639] 3491 3492 3493 3494 3495 3496 3497 3498 3499 3500 3501 3502 3503 3504
[1653] 3505 3506 3507 3508 3509 3510 3511 3512 3513 3514 3515 3516 3517 3518
[1667] 3519 3520 3521 3522 3523 3524 3525 3526 3527 3528 3529 3530 3531 3532
[1681] 3533 3534 3535 3536 3537 3538 3539 3540 3541 3542 3543 3544 3545 3546
[1695] 3547 3548 3549 3550 3551 3552 3553 3554 3555 3556 3557 3558 3559 3560
[1709] 3561 3562 3563 3564 3565 3566 3567 3568 3569 3570 3571 3572 3573 3574
[1723] 3575 3576 3577 3578 3579 3580 3581 3582 3583 3584 3585 3586 3587 3588
[1737] 3589 3590 3591 3592 3593 3594 3595 3596 3597 3598 3599 3600 3601 3602
[1751] 3603 3604 3605 3606 3607 3608 3609 3610 3611 3612 3613 3614 3615 3616
[1765] 3617 3618 3619 3620 3621 3622 3623 3624 3625 3626 3627 3628 3629 3630
[1779] 3631 3632 3633 3634 3635 3636 3637 3638 3639 3640 3641 3642 3643 3644
[1793] 3645 3646 3647 3648 3649 3650 3651 3652 3653 3654 3655 3656 3657 3658
[1807] 3659 3660 3661 3662 3663 3664 3665 3666 3667 3668 3669 3670 3671 3672
[1821] 3673 3674 3675 3676 3677 3678 3679 3680 3681 3682 3683 3684 3685 3686
[1835] 3687 3688 3689 3690 3691 3692 3693 3694 3695 3696 3697 3698 3699 3700
[1849] 3701 3702 3703 3704 3705 3706 3707 3708 3709 3710 3711 3712 3713 3714
[1863] 3715 3716 3717 3718 3719 3720 3721 3722 3723 3724 3725 3726 3727 3728
[1877] 3729 3730 3731 3732 3733 3734 3735 3736 3737 3738 3739 3740 3741 3742
[1891] 3743 3744 3745 3746 3747 3748 3749 3750 3751 3752 3753 3754 3755 3756
[1905] 3757 3758 3759 3760 3761 3762 3763 3764 3765 3766 3767 3768 3769 3770
[1919] 3771 3772 3773 3774 3775 3776 3777 3778 3779 3780 3781 3782 3783 3784
[1933] 3785 3786 3787 3788 3789 3790 3791 3792 3793 3794 3795 3796 3797 3798
[1947] 3799 3800 3801 3802 3803 3804 3805 3806 3807 3808 3809 3810 3811 3812
[1961] 3813 3814 3815 3816 3817 3818 3819 3820 3821 3822 3823 3824 3825 3826
[1975] 3827 3828 3829 3830 3831 3832 3833 3834 3835 3836 3837 3838 3839 3840
[1989] 3841 3842 3843 3844 3845 3846 3847 3848 3849 3850 3851 3852 3853 3854
[2003] 3855 3856 3857 3858 3859 3860 3861 3862 3863 3864 3865 3866 3867 3868
[2017] 3869 3870 3871 3872 3873 3874 3875 3876 3877 3878 3879 3880 3881 3882
[2031] 3883 3884 3885 3886 3887 3888 3889 3890 3891 3892 3893 3894 3895 3896
[2045] 3897 3898 3899 3900 3901 3902 3903 3904 3905 3906 3907 3908 3909 3910
[2059] 3911 3912 3913 3914 3915 3916 3917 3918 3919 3920 3921 3922 3923 3924
[2073] 3925 3926 3927 3928 3929 3930 3931 3932 3933 3934 3935 3936 3937 3938
[2087] 3939 3940 3941 3942 3943 3944 3945 3946 3947 3948 3949 3950 3951 3952
[2101] 3953 3954 3955 3956 3957 3958 3959 3960 3961 3962 3963 3964 3965 3966
[2115] 3967 3968 3969 3970 3971 3972 3973 3974 3975 3976 3977 3978 3979 3980
[2129] 3981 3982 3983 3984 3985 3986 3987 3988 3989 3990 3991 3992 3993 3994
[2143] 3995 3996 3997 3998 3999 4000 4001 4002 4003 4004 4005 4006 4007 4008
[2157] 4009 4010 4011 4012 4013 4014 4015 4016 4017 4018 4019 4020 4021 4022
[2171] 4023 4024 4025 4026 4027 4028 4029 4030 4031 4032 4033 4034 4035 4036
[2185] 4037 4038 4039 4040 4041 4042 4043 4044 4045 4046 4047 4048 4049 4050
[2199] 4051 4052 4053 4054 4055 4056 4057 4058 4059 4060 4061 4062 4063 4064
[2213] 4065 4066 4067 4068 4069 4070 4071 4072 4073 4074 4075 4076 4077 4078
[2227] 4079 4080 4081 4082 4083 4084 4085 4086 4087 4088 4089 4090 4091 4092
[2241] 4093 4094 4095 4096 4097 4098 4099 4100 4101 4102 4103 4104 4105 4106
[2255] 4107 4108 4109 4110 4111 4112 4113 4114 4115 4116 4117 4118 4119 4120
[2269] 4121 4122 4123 4124 4125 4126 4127 4128 4129 4130 4131 4132 4133 4134
[2283] 4135 4136 4137 4138 4139 4140 4141 4142 4143 4144 4145 4146 4147 4148
[2297] 4149 4150 4151 4152 4153 4154 4155 4156 4157 4158 4159 4160 4161 4162
[2311] 4163 4164 4165 4166 4167 4168 4169 4170 4171 4172 4173 4174 4175 4176
[2325] 4177 4178 4179 4180 4181 4182 4183 4184 4185 4186 4187 4188 4189 4190
[2339] 4191 4192 4193 4194 4195 4196 4197 4198 4199 4200 4201 4202 4203 4204
[2353] 4205 4206 4207 4208 4209 4210 4211 4212 4213 4214 4215 4216 4217 4218
[2367] 4219 4220 4221 4222 4223 4224 4225 4226 4227 4228 4229 4230 4231 4232
[2381] 4233 4234 4235 4236 4237 4238 4239 4240 4241 4242 4243 4244 4245 4246
[2395] 4247 4248 4249 4250 4251 4252 4253 4254 4255 4256 4257 4258 4259 4260
[2409] 4261 4262 4263 4264 4265 4266 4267 4268 4269 4270 4271 4272 4273 4274
[2423] 4275 4276 4277 4278 4279 4280 4281 4282 4283 4284 4285 4286 4287 4288
[2437] 4289 4290 4291 4292 4293 4294 4295 4296 4297 4298 4299 4300 4301 4302
[2451] 4303 4304 4305 4306 4307 4308 4309 4310 4311 4312 4313 4314 4315 4316
[2465] 4317 4318 4319 4320 4321 4322 4323 4324 4325 4326 4327 4328 4329 4330
[2479] 4331 4332 4333 4334 4335 4336 4337 4338 4339 4340 4341 4342 4343 4344
[2493] 4345 4346 4347 4348 4349 4350 4351 4352 4353 4354 4355 4356 4357 4358
[2507] 4359 4360 4361 4362 4363 4364 4365 4366 4367 4368 4369 4370 4371 4372
[2521] 4373 4374 4375 4376 4377 4378 4379 4380 4381 4382 4383 4384 4385 4386
[2535] 4387 4388 4389 4390 4391 4392 4393 4394 4395 4396 4397 4398 4399 4400
[2549] 4401 4402 4403 4404 4405 4406 4407 4408 4409 4410 4411 4412 4413 4414
[2563] 4415 4416 4417 4418 4419 4420 4421 4422 4423 4424 4425 4426 4427 4428
[2577] 4429 4430 4431 4432 4433 4434 4435 4436 4437 4438 4439 4440 4441 4442
[2591] 4443 4444 4445 4446 4447 4485 4494 4495 4512 4518 4605 4608 4609 4610
[2605] 4611 4612 4613 4614 4615 4616 4617 4624 4642 4643 4667 4671 4672 4673
[2619] 4704 4705
$orig$Government
[1] 24 46 266 299 307 310 316 332 369 378 379 380 381 382 383
[16] 384 385 386 387 645 691 727 733 734 738 815 822 895 896 897
[31] 898 989 1026 1039 1109 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120
[46] 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135
[61] 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150
[76] 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165
[91] 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180
[106] 1181 1182 1183 1184 1185 1186 1196 1198 1200 1326 1327 1328 1506 1522 1529
[121] 1530 1531 1532 1543 1622 1623 1624 1625 1632 1672 1673 1694 2148 2169 2170
[136] 2180 2201 2236 2335 2373 2464 2484 2485 4489 4490 4491 4492 4493 4494 4495
[151] 4496 4497 4498 4513 4581 4592 4635 4636 4647 4648 4649
$orig$Military
[1] 12 13 14 15 16 17 18 19 20 21 22 23 32 33 34
[16] 55 290 295 305 306 314 331 334 335 341 400 420 421 422 423
[31] 424 594 595 596 597 598 599 600 601 602 603 605 654 690 710
[46] 711 712 713 715 728 729 730 731 732 736 894 951 952 953 955
[61] 956 957 965 967 1051 1052 1053 1054 1277 1278 1279 1280 1281 1282 1313
[76] 1314 1335 1488 1489 1497 1498 1499 1524 1525 1526 1527 1535 1543 1544 1545
[91] 1546 1547 1548 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568
[106] 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583
[121] 1584 1585 1586 1587 1588 1589 1590 1591 1633 1634 1635 1636 1637 1638 1639
[136] 1640 1641 1642 1695 1697 2141 2142 2143 2144 2145 2146 2155 2156 2178 2179
[151] 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2232 2241 2242 2243
[166] 2312 2313 2314 2315 2316 2317 2318 2319 2320 2321 2322 2323 2324 2325 2326
[181] 2327 2328 2329 2336 2337 2338 2339 2340 2341 2342 2346 2372 2463 2465 2466
[196] 2467 2468 2470 2471 2472 2473 4448 4459 4460 4461 4479 4480 4488 4593 4618
[211] 4619 4620 4621 4622 4627 4628 4629 4630 4631 4632 4657 4658 4659 4660 4661
[226] 4662 4663 4664 4665 4666 4674
$svg
[1] "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"700\" height=\"500\">\n<defs>\n<style type=\"text/css\"><![CDATA[\n .borderLine {\n\t stroke: none;\n\t fill-opacity: 0.4 ;\n }\n .outLine {\n\t stroke-width: 1 ;\n\t fill: none;\n }\n .circle {\n\t stroke: #888888;\n\t stroke-width: 0.5;\n\t fill: none;\n pointer-events: all;\n }\n .nLabel {\n\t font-family: Arial;\n pointer-events: none;\n\t font-size: 10px;\n\t text-anchor: middle;\n\t alignment-baseline: central;\n }\n .belong {\n\t font-family: Arial;\n pointer-events: none;\n\t font-size: 5px;\n\t text-anchor: middle;\n\t alignment-baseline: central;\n }\n .legend {\n font-family: Arial;\n font-size: 15px;\n }\n #legendBox {\n display: inline ;\n }\n .p0{\n stroke: none;\n fill: #e6194b;\n }\n .q0{\n fill: none;\n stroke: #e6194b;\n }\n .p1{\n stroke: none;\n fill: #3cb44b;\n }\n .q1{\n fill: none;\n stroke: #3cb44b;\n }\n .p2{\n stroke: none;\n fill: #ffe119;\n }\n .q2{\n fill: none;\n stroke: #ffe119;\n }\n .p3{\n stroke: none;\n fill: #0082c8;\n }\n .q3{\n fill: none;\n stroke: #0082c8;\n }\n]]>\n</style>\n<symbol id=\"bl0\">\n<path d=\"M 380.15 364.48 L 379.53 362.90 L 378.73 361.41 L 377.72 360.05 L 376.53 358.84 L 375.18 357.81 L 373.68 356.99 L 372.09 356.38 L 370.41 356.00 L 368.71 355.86 L 367.00 355.96 L 365.32 356.30 L 363.70 356.86 L 362.18 357.66 L 360.79 358.65 L 359.54 359.82 L 358.47 361.16 L 357.57 362.63 L 356.83 364.23 L 356.18 365.92 L 355.56 367.72 L 354.86 369.59 L 354.00 371.48 L 352.95 373.32 L 351.72 375.03 L 350.36 376.54 L 348.91 377.81 L 347.41 378.84 L 345.87 379.69 L 344.30 380.38 L 342.69 380.99 L 341.07 381.56 L 339.46 382.14 L 337.86 382.74 L 336.27 383.38 L 334.71 384.08 L 333.16 384.83 L 331.62 385.64 L 330.10 386.49 L 328.59 387.38 L 327.09 388.33 L 325.62 389.32 L 324.18 390.35 L 322.75 391.41 L 321.34 392.51 L 319.93 393.62 L 318.52 394.73 L 317.10 395.83 L 315.66 396.89 L 314.19 397.91 L 312.68 398.86 L 311.13 399.72 L 309.53 400.50 L 307.91 401.18 L 306.25 401.74 L 304.56 402.19 L 302.85 402.53 L 301.13 402.75 L 299.40 402.86 L 297.67 402.85 L 295.94 402.72 L 294.21 402.49 L 292.50 402.14 L 290.81 401.68 L 289.15 401.11 L 287.54 400.42 L 285.98 399.63 L 284.47 398.74 L 283.03 397.74 L 281.65 396.65 L 280.36 395.47 L 279.15 394.20 L 278.04 392.86 L 277.01 391.44 L 276.08 389.95 L 275.25 388.41 L 274.52 386.82 L 273.90 385.19 L 273.38 383.51 L 272.98 381.81 L 272.70 380.07 L 272.53 378.32 L 272.48 376.57 L 272.55 374.81 L 272.74 373.05 L 273.04 371.31 L 273.47 369.60 L 274.00 367.92 L 274.65 366.28 L 275.42 364.69 L 276.28 363.15 L 277.24 361.65 L 278.28 360.21 L 279.39 358.82 L 280.56 357.47 L 281.76 356.15 L 282.99 354.85 L 284.24 353.57 L 285.49 352.29 L 286.75 351.03 L 288.02 349.77 L 289.28 348.52 L 290.54 347.27 L 291.80 346.02 L 293.06 344.76 L 294.32 343.50 L 295.59 342.24 L 296.87 340.99 L 298.16 339.74 L 299.45 338.50 L 300.73 337.26 L 302.01 336.02 L 303.29 334.79 L 304.59 333.57 L 305.90 332.36 L 307.25 331.18 L 308.63 330.06 L 310.05 329.02 L 311.52 328.10 L 313.05 327.34 L 314.62 326.75 L 316.24 326.34 L 317.89 326.13 L 319.56 326.10 L 321.23 326.26 L 322.89 326.59 L 324.52 327.09 L 326.11 327.72 L 327.66 328.46 L 329.17 329.27 L 330.65 330.14 L 332.11 331.03 L 333.56 331.93 L 335.00 332.84 L 336.44 333.77 L 337.88 334.69 L 339.32 335.63 L 340.77 336.57 L 342.21 337.52 L 343.66 338.48 L 345.11 339.44 L 346.56 340.40 L 348.01 341.37 L 349.44 342.35 L 350.87 343.33 L 352.31 344.30 L 353.75 345.26 L 355.21 346.21 L 356.69 347.13 L 358.20 348.01 L 359.75 348.83 L 361.33 349.58 L 362.94 350.25 L 364.59 350.82 L 366.26 351.29 L 367.96 351.66 L 369.67 351.92 L 371.40 352.08 L 373.12 352.14 L 374.84 352.09 L 376.56 351.94 L 378.25 351.68 L 379.94 351.34 L 381.60 350.90 L 383.23 350.39 L 384.85 349.83 L 386.46 349.22 L 388.07 348.59 L 389.68 347.95 L 391.29 347.31 L 392.89 346.68 L 394.49 346.09 L 396.11 345.56 L 397.75 345.10 L 399.41 344.74 L 401.08 344.49 L 402.77 344.34 L 404.47 344.31 L 406.17 344.41 L 407.85 344.63 L 409.50 344.98 L 411.12 345.45 L 412.71 346.05 L 414.25 346.78 L 415.75 347.64 L 417.18 348.64 L 418.55 349.79 L 419.86 351.05 L 421.09 352.44 L 422.23 353.91 L 423.29 355.46 L 424.22 357.05 L 425.02 358.69 L 425.68 360.35 L 426.19 362.04 L 426.56 363.76 L 426.77 365.50 L 426.84 367.25 L 426.77 369.01 L 426.57 370.76 L 426.23 372.49 L 425.76 374.19 L 425.16 375.85 L 424.43 377.45 L 423.59 378.98 L 422.63 380.45 L 421.57 381.84 L 420.41 383.13 L 419.15 384.34 L 417.80 385.44 L 416.37 386.43 L 414.87 387.33 L 413.31 388.10 L 411.70 388.75 L 410.06 389.28 L 408.37 389.68 L 406.66 389.94 L 404.94 390.08 L 403.22 390.09 L 401.49 389.97 L 399.78 389.71 L 398.09 389.33 L 396.43 388.82 L 394.82 388.18 L 393.25 387.42 L 391.75 386.54 L 390.31 385.55 L 388.95 384.46 L 387.69 383.28 L 386.53 382.00 L 385.46 380.64 L 384.51 379.21 L 383.67 377.72 L 382.95 376.15 L 382.35 374.53 L 381.86 372.86 L 381.48 371.18 L 381.18 369.48 L 380.91 367.79 L 380.58 366.11 Z\" />\n</symbol>\n<symbol id=\"bl1\">\n<path d=\"M 436.22 257.39 L 436.55 262.59 L 436.22 267.79 L 435.24 272.92 L 433.68 277.94 L 431.64 282.85 L 429.24 287.64 L 426.62 292.33 L 423.85 296.95 L 421.01 301.52 L 418.12 306.08 L 415.22 310.61 L 412.30 315.15 L 409.37 319.68 L 406.41 324.18 L 403.38 328.63 L 400.20 332.94 L 396.79 337.02 L 393.05 340.73 L 388.93 343.94 L 384.40 346.52 L 379.53 348.41 L 374.37 349.62 L 369.02 350.26 L 363.54 350.51 L 357.96 350.56 L 352.31 350.59 L 346.55 350.69 L 340.69 350.91 L 334.73 351.22 L 328.67 351.55 L 322.57 351.83 L 316.50 352.00 L 310.51 352.08 L 304.66 352.11 L 298.97 352.23 L 293.43 352.56 L 288.03 353.18 L 282.73 354.08 L 277.48 355.20 L 272.27 356.44 L 267.07 357.74 L 261.87 359.06 L 256.68 360.37 L 251.48 361.67 L 246.28 362.96 L 241.06 364.21 L 235.83 365.40 L 230.58 366.49 L 225.31 367.46 L 220.02 368.29 L 214.70 368.98 L 209.36 369.51 L 204.01 369.88 L 198.64 370.09 L 193.27 370.14 L 187.91 370.04 L 182.55 369.77 L 177.21 369.35 L 171.89 368.77 L 166.60 368.03 L 161.33 367.13 L 156.08 366.08 L 150.88 364.87 L 145.72 363.51 L 140.61 362.00 L 135.56 360.35 L 130.56 358.54 L 125.62 356.58 L 120.74 354.49 L 115.93 352.25 L 111.19 349.87 L 106.53 347.35 L 101.93 344.70 L 97.42 341.91 L 92.99 338.99 L 88.64 335.93 L 84.38 332.75 L 80.22 329.44 L 76.16 326.01 L 72.20 322.46 L 68.34 318.80 L 64.61 315.02 L 60.99 311.15 L 57.49 307.17 L 54.11 303.08 L 50.85 298.91 L 47.72 294.64 L 44.71 290.28 L 41.83 285.84 L 39.08 281.30 L 36.47 276.69 L 34.00 272.01 L 31.66 267.26 L 29.47 262.44 L 27.42 257.56 L 25.51 252.62 L 23.75 247.62 L 22.14 242.58 L 20.68 237.50 L 19.37 232.37 L 18.21 227.21 L 17.20 222.02 L 16.35 216.81 L 15.65 211.58 L 15.10 206.34 L 14.71 201.08 L 14.47 195.81 L 14.38 190.54 L 14.45 185.27 L 14.67 179.98 L 15.05 174.71 L 15.58 169.46 L 16.27 164.22 L 17.11 159.00 L 18.11 153.81 L 19.25 148.65 L 20.55 143.53 L 21.99 138.46 L 23.58 133.43 L 25.32 128.45 L 27.20 123.53 L 29.22 118.66 L 31.39 113.85 L 33.69 109.11 L 36.14 104.43 L 38.72 99.81 L 41.44 95.28 L 44.30 90.82 L 47.29 86.45 L 50.40 82.17 L 53.64 77.98 L 57.01 73.88 L 60.48 69.89 L 64.08 66.00 L 67.78 62.22 L 71.59 58.55 L 75.51 55.00 L 79.52 51.56 L 83.64 48.25 L 87.84 45.05 L 92.14 41.99 L 96.52 39.05 L 100.99 36.24 L 105.53 33.57 L 110.16 31.02 L 114.86 28.62 L 119.63 26.35 L 124.46 24.22 L 129.35 22.23 L 134.29 20.39 L 139.28 18.70 L 144.30 17.16 L 149.37 15.76 L 154.47 14.51 L 159.61 13.42 L 164.78 12.47 L 169.99 11.67 L 175.22 11.03 L 180.47 10.54 L 185.73 10.21 L 191.01 10.03 L 196.30 10.00 L 201.58 10.13 L 206.86 10.41 L 212.14 10.85 L 217.40 11.44 L 222.65 12.18 L 227.88 13.09 L 233.09 14.15 L 238.28 15.36 L 243.43 16.73 L 248.54 18.25 L 253.61 19.92 L 258.63 21.74 L 263.59 23.71 L 268.49 25.83 L 273.33 28.09 L 278.08 30.49 L 282.77 33.03 L 287.37 35.71 L 291.88 38.51 L 296.30 41.44 L 300.64 44.50 L 304.87 47.68 L 309.01 50.98 L 313.05 54.41 L 316.98 57.95 L 320.81 61.60 L 324.53 65.36 L 328.13 69.23 L 331.62 73.21 L 335.00 77.29 L 338.25 81.47 L 341.39 85.74 L 344.41 90.09 L 347.33 94.52 L 350.16 99.00 L 352.92 103.51 L 355.64 108.05 L 358.33 112.60 L 360.99 117.16 L 363.64 121.73 L 366.29 126.30 L 368.93 130.88 L 371.57 135.47 L 374.22 140.06 L 376.86 144.67 L 379.50 149.28 L 382.15 153.90 L 384.80 158.51 L 387.47 163.12 L 390.13 167.73 L 392.79 172.34 L 395.46 176.96 L 398.12 181.57 L 400.78 186.19 L 403.44 190.81 L 406.09 195.43 L 408.75 200.05 L 411.41 204.68 L 414.08 209.30 L 416.74 213.92 L 419.40 218.55 L 422.05 223.19 L 424.67 227.86 L 427.22 232.56 L 429.65 237.33 L 431.87 242.19 L 433.79 247.16 L 435.27 252.23 Z\" />\n</symbol>\n<symbol id=\"bl2\">\n<path d=\"M 422.42 384.38 L 419.71 385.52 L 416.97 386.60 L 414.20 387.57 L 411.39 388.38 L 408.53 388.97 L 405.63 389.30 L 402.71 389.36 L 399.79 389.14 L 396.89 388.68 L 394.02 388.03 L 391.18 387.23 L 388.36 386.34 L 385.55 385.41 L 382.75 384.46 L 379.94 383.52 L 377.11 382.60 L 374.25 381.71 L 371.38 380.84 L 368.54 379.90 L 365.80 378.81 L 363.30 377.42 L 361.11 375.62 L 359.28 373.32 L 357.82 370.52 L 356.62 367.26 L 355.57 363.66 L 354.55 359.87 L 353.49 356.05 L 352.40 352.33 L 351.32 348.80 L 350.27 345.51 L 349.27 342.42 L 348.33 339.48 L 347.46 336.63 L 346.68 333.80 L 346.05 330.97 L 345.60 328.11 L 345.36 325.24 L 345.34 322.36 L 345.53 319.48 L 345.92 316.62 L 346.44 313.77 L 347.07 310.95 L 347.76 308.14 L 348.49 305.34 L 349.24 302.53 L 350.01 299.72 L 350.78 296.91 L 351.56 294.09 L 352.34 291.28 L 353.12 288.47 L 353.90 285.66 L 354.68 282.86 L 355.47 280.05 L 356.26 277.25 L 357.05 274.44 L 357.85 271.63 L 358.66 268.83 L 359.48 266.02 L 360.30 263.19 L 361.15 260.33 L 362.07 257.41 L 363.07 254.41 L 364.19 251.36 L 365.44 248.32 L 366.85 245.37 L 368.41 242.58 L 370.14 239.98 L 372.05 237.58 L 374.14 235.38 L 376.40 233.37 L 378.81 231.55 L 381.36 229.94 L 384.03 228.53 L 386.80 227.35 L 389.66 226.39 L 392.58 225.66 L 395.55 225.17 L 398.55 224.91 L 401.57 224.89 L 404.57 225.11 L 407.55 225.58 L 410.49 226.28 L 413.36 227.21 L 416.16 228.37 L 418.86 229.75 L 421.44 231.34 L 423.89 233.13 L 426.19 235.12 L 428.33 237.29 L 430.29 239.61 L 432.05 242.09 L 433.62 244.69 L 434.97 247.41 L 436.09 250.22 L 436.98 253.10 L 437.64 256.03 L 438.06 259.00 L 438.24 261.99 L 438.19 264.99 L 437.90 267.98 L 437.37 270.93 L 436.62 273.83 L 435.64 276.65 L 434.45 279.39 L 433.04 282.01 L 431.43 284.52 L 429.64 286.90 L 427.66 289.12 L 425.52 291.18 L 423.21 293.07 L 420.76 294.77 L 418.18 296.28 L 415.47 297.59 L 412.67 298.68 L 409.78 299.54 L 406.83 300.18 L 403.85 300.58 L 400.86 300.75 L 397.87 300.68 L 394.89 300.39 L 391.95 299.88 L 389.02 299.21 L 386.12 298.43 L 383.21 297.63 L 380.30 296.93 L 377.36 296.42 L 374.40 296.16 L 371.45 296.19 L 368.53 296.52 L 365.67 297.16 L 362.89 298.10 L 360.23 299.32 L 357.72 300.82 L 355.38 302.57 L 353.25 304.55 L 351.34 306.72 L 349.67 309.08 L 348.27 311.59 L 347.14 314.22 L 346.29 316.95 L 345.73 319.76 L 345.47 322.63 L 345.51 325.51 L 345.84 328.38 L 346.47 331.20 L 347.40 333.95 L 348.61 336.59 L 350.09 339.10 L 351.82 341.43 L 353.80 343.57 L 355.98 345.49 L 358.34 347.17 L 360.87 348.59 L 363.52 349.73 L 366.29 350.59 L 369.13 351.15 L 372.02 351.42 L 374.92 351.38 L 377.80 351.04 L 380.63 350.40 L 383.37 349.48 L 386.00 348.28 L 388.51 346.83 L 390.88 345.14 L 393.13 343.26 L 395.27 341.23 L 397.34 339.12 L 399.38 336.95 L 401.40 334.75 L 403.41 332.55 L 405.43 330.35 L 407.46 328.13 L 409.50 325.91 L 411.53 323.68 L 413.57 321.44 L 415.60 319.19 L 417.62 316.94 L 419.64 314.71 L 421.68 312.49 L 423.76 310.34 L 425.89 308.27 L 428.12 306.33 L 430.45 304.55 L 432.88 302.96 L 435.43 301.57 L 438.08 300.40 L 440.81 299.46 L 443.60 298.75 L 446.45 298.27 L 449.33 298.03 L 452.23 298.03 L 455.12 298.28 L 457.99 298.76 L 460.80 299.48 L 463.55 300.43 L 466.22 301.61 L 468.78 303.01 L 471.20 304.61 L 473.49 306.41 L 475.61 308.40 L 477.57 310.55 L 479.34 312.87 L 480.91 315.32 L 482.26 317.89 L 483.40 320.57 L 484.31 323.34 L 484.98 326.17 L 485.42 329.05 L 485.62 331.96 L 485.57 334.88 L 485.28 337.78 L 484.75 340.66 L 483.97 343.48 L 482.96 346.24 L 481.73 348.90 L 480.28 351.43 L 478.63 353.83 L 476.79 356.08 L 474.78 358.17 L 472.61 360.08 L 470.29 361.83 L 467.85 363.43 L 465.31 364.89 L 462.70 366.26 L 460.06 367.55 L 457.38 368.80 L 454.70 370.02 L 452.00 371.23 L 449.30 372.42 L 446.60 373.62 L 443.90 374.81 L 441.21 376.01 L 438.52 377.21 L 435.84 378.41 L 433.16 379.61 L 430.48 380.81 L 427.80 382.01 L 425.11 383.20 Z\" />\n</symbol>\n<symbol id=\"bl3\">\n<path d=\"M 410.61 451.69 L 409.91 453.72 L 409.13 455.73 L 408.28 457.70 L 407.35 459.63 L 406.35 461.53 L 405.27 463.39 L 404.12 465.20 L 402.90 466.97 L 401.61 468.69 L 400.25 470.35 L 398.82 471.96 L 397.32 473.52 L 395.77 475.02 L 394.16 476.48 L 392.48 477.86 L 390.76 479.19 L 388.99 480.45 L 387.18 481.63 L 385.33 482.73 L 383.44 483.76 L 381.53 484.71 L 379.59 485.57 L 377.62 486.37 L 375.61 487.09 L 373.58 487.74 L 371.51 488.31 L 369.43 488.80 L 367.31 489.21 L 365.18 489.53 L 363.03 489.77 L 360.87 489.92 L 358.71 490.00 L 356.55 490.00 L 354.40 489.91 L 352.25 489.73 L 350.11 489.47 L 347.98 489.12 L 345.86 488.69 L 343.76 488.18 L 341.67 487.59 L 339.61 486.92 L 337.58 486.17 L 335.58 485.34 L 333.61 484.43 L 331.69 483.44 L 329.81 482.38 L 327.97 481.25 L 326.18 480.06 L 324.45 478.80 L 322.78 477.48 L 321.16 476.11 L 319.60 474.68 L 318.09 473.19 L 316.64 471.65 L 315.25 470.04 L 313.91 468.38 L 312.64 466.66 L 311.44 464.90 L 310.30 463.09 L 309.23 461.24 L 308.24 459.35 L 307.33 457.43 L 306.49 455.47 L 305.72 453.47 L 305.03 451.45 L 304.42 449.40 L 303.89 447.33 L 303.45 445.25 L 303.07 443.15 L 302.79 441.04 L 302.57 438.94 L 302.44 436.83 L 302.40 434.71 L 302.44 432.59 L 302.56 430.46 L 302.77 428.33 L 303.06 426.21 L 303.42 424.11 L 303.86 422.02 L 304.38 419.95 L 304.98 417.91 L 305.66 415.89 L 306.41 413.90 L 307.24 411.94 L 308.14 410.02 L 309.11 408.13 L 310.14 406.28 L 311.23 404.47 L 312.38 402.70 L 313.56 400.97 L 314.76 399.26 L 315.99 397.57 L 317.22 395.87 L 318.43 394.17 L 319.60 392.43 L 320.69 390.66 L 321.68 388.83 L 322.55 386.93 L 323.28 384.98 L 323.85 382.99 L 324.27 380.95 L 324.52 378.88 L 324.60 376.79 L 324.54 374.69 L 324.36 372.55 L 324.10 370.39 L 323.86 368.22 L 323.72 366.05 L 323.76 363.92 L 324.02 361.82 L 324.52 359.78 L 325.25 357.83 L 326.18 355.95 L 327.27 354.15 L 328.45 352.41 L 329.65 350.69 L 330.81 348.96 L 331.89 347.19 L 332.88 345.37 L 333.79 343.52 L 334.66 341.64 L 335.51 339.75 L 336.35 337.85 L 337.17 335.94 L 337.99 334.03 L 338.80 332.11 L 339.60 330.18 L 340.39 328.25 L 341.17 326.33 L 341.94 324.40 L 342.72 322.48 L 343.49 320.56 L 344.28 318.65 L 345.08 316.75 L 345.92 314.85 L 346.82 312.98 L 347.78 311.14 L 348.84 309.34 L 350.00 307.61 L 351.27 305.95 L 352.66 304.39 L 354.15 302.93 L 355.74 301.58 L 357.42 300.36 L 359.19 299.26 L 361.03 298.31 L 362.92 297.49 L 364.87 296.81 L 366.87 296.28 L 368.89 295.90 L 370.93 295.67 L 372.98 295.60 L 375.05 295.69 L 377.12 295.94 L 379.22 296.38 L 381.33 296.98 L 383.44 297.75 L 385.52 298.68 L 387.52 299.74 L 389.39 300.91 L 391.12 302.18 L 392.70 303.55 L 394.15 305.02 L 395.47 306.60 L 396.67 308.28 L 397.74 310.04 L 398.67 311.89 L 399.48 313.79 L 400.13 315.76 L 400.65 317.78 L 401.02 319.82 L 401.23 321.89 L 401.30 323.96 L 401.22 326.04 L 400.97 328.11 L 400.58 330.16 L 400.03 332.17 L 399.34 334.14 L 398.50 336.05 L 397.54 337.91 L 396.45 339.70 L 395.27 341.44 L 394.01 343.12 L 392.69 344.76 L 391.34 346.38 L 389.98 347.98 L 388.63 349.59 L 387.33 351.23 L 386.09 352.91 L 384.96 354.66 L 383.96 356.47 L 383.10 358.36 L 382.41 360.33 L 381.90 362.38 L 381.58 364.51 L 381.44 366.70 L 381.50 368.92 L 381.76 371.11 L 382.22 373.25 L 382.87 375.32 L 383.70 377.30 L 384.70 379.20 L 385.83 381.03 L 387.07 382.80 L 388.38 384.50 L 389.74 386.17 L 391.13 387.81 L 392.52 389.44 L 393.92 391.07 L 395.32 392.70 L 396.71 394.33 L 398.10 395.97 L 399.47 397.63 L 400.81 399.32 L 402.11 401.03 L 403.36 402.79 L 404.55 404.60 L 405.67 406.44 L 406.72 408.32 L 407.70 410.22 L 408.60 412.16 L 409.42 414.12 L 410.17 416.10 L 410.83 418.11 L 411.42 420.14 L 411.93 422.19 L 412.36 424.26 L 412.73 426.34 L 413.01 428.44 L 413.21 430.56 L 413.33 432.68 L 413.36 434.82 L 413.31 436.95 L 413.17 439.08 L 412.94 441.21 L 412.64 443.33 L 412.25 445.44 L 411.79 447.55 L 411.24 449.63 Z\" />\n</symbol>\n</defs>\n<!-- isDone: 0 -->\n<rect width=\"700\" height=\"500\" style=\"fill:#fff;stroke-width:0\" />\n<use class=\"p0 borderLine\" xlink:href=\"#bl0\"/>\n<use class=\"p1 borderLine\" xlink:href=\"#bl1\"/>\n<use class=\"p2 borderLine\" xlink:href=\"#bl2\"/>\n<use class=\"p3 borderLine\" xlink:href=\"#bl3\"/>\n<use class=\"q0 outLine\" xlink:href=\"#bl0\"/>\n<use class=\"q1 outLine\" xlink:href=\"#bl1\"/>\n<use class=\"q2 outLine\" xlink:href=\"#bl2\"/>\n<use class=\"q3 outLine\" xlink:href=\"#bl3\"/>\n<circle onclick=\"fromCircle(8)\" class=\"circle\" cx=\"298.5911\" cy=\"376.7458\" r=\"18.2849\" />\n<text class=\"nLabel\" x=\"298.59\" y=\"371.75\">28</text>\n<text class=\"belong\" x=\"298.59\" y=\"381.75\">(1)</text>\n<circle onclick=\"fromCircle(12)\" class=\"circle\" cx=\"318.8787\" cy=\"340.6224\" r=\"6.9111\" />\n<text class=\"nLabel\" x=\"318.88\" y=\"335.62\">4</text>\n<text class=\"belong\" x=\"318.88\" y=\"345.62\">(1, 2)</text>\n<circle onclick=\"fromCircle(4)\" class=\"circle\" cx=\"193.9036\" cy=\"189.5599\" r=\"172.8801\" />\n<text class=\"nLabel\" x=\"193.90\" y=\"184.56\">2503</text>\n<text class=\"belong\" x=\"193.90\" y=\"194.56\">(2)</text>\n<circle onclick=\"fromCircle(6)\" class=\"circle\" cx=\"400.3484\" cy=\"262.7227\" r=\"30.5184\" />\n<text class=\"nLabel\" x=\"400.35\" y=\"257.72\">78</text>\n<text class=\"belong\" x=\"400.35\" y=\"267.72\">(2, 3)</text>\n<circle onclick=\"fromCircle(5)\" class=\"circle\" cx=\"373.1371\" cy=\"323.7666\" r=\"20.4432\" />\n<text class=\"nLabel\" x=\"373.14\" y=\"318.77\">35</text>\n<text class=\"belong\" x=\"373.14\" y=\"328.77\">(2, 4)</text>\n<circle onclick=\"fromCircle(2)\" class=\"circle\" cx=\"450.9160\" cy=\"332.8430\" r=\"27.4274\" />\n<text class=\"nLabel\" x=\"450.92\" y=\"327.84\">63</text>\n<text class=\"belong\" x=\"450.92\" y=\"337.84\">(3)</text>\n<circle onclick=\"fromCircle(10)\" class=\"circle\" cx=\"404.0308\" cy=\"367.2375\" r=\"15.0623\" />\n<text class=\"nLabel\" x=\"404.03\" y=\"362.24\">19</text>\n<text class=\"belong\" x=\"404.03\" y=\"372.24\">(1, 3)</text>\n<circle onclick=\"fromCircle(3)\" class=\"circle\" cx=\"368.5558\" cy=\"368.0585\" r=\"4.7124\" />\n<text class=\"nLabel\" x=\"368.56\" y=\"363.06\">1</text>\n<text class=\"belong\" x=\"368.56\" y=\"373.06\">(3, 4)</text>\n<circle onclick=\"fromCircle(1)\" class=\"circle\" cx=\"357.8804\" cy=\"434.6952\" r=\"47.5057\" />\n<text class=\"nLabel\" x=\"357.88\" y=\"429.70\">189</text>\n<text class=\"belong\" x=\"357.88\" y=\"439.70\">(4)</text>\n<circle onclick=\"fromCircle(9)\" class=\"circle\" cx=\"339.5206\" cy=\"364.8257\" r=\"8.4643\" />\n<text class=\"nLabel\" x=\"339.52\" y=\"359.83\">6</text>\n<text class=\"belong\" x=\"339.52\" y=\"369.83\">(1, 4)</text>\n<g id=\"legendBox\">\n<rect class=\"p0 borderLine\" x=\"500.00\" y=\"50.00\" width=\"30.00\" height=\"15.00\" />\n<rect class=\"q0\" x=\"500.00\" y=\"50.00\" width=\"30.00\" height=\"15.00\" />\n<text class=\"legend\" x=\"540.00\" y=\"65.00\">Civil</text>\n<rect class=\"p1 borderLine\" x=\"500.00\" y=\"90.00\" width=\"30.00\" height=\"15.00\" />\n<rect class=\"q1\" x=\"500.00\" y=\"90.00\" width=\"30.00\" height=\"15.00\" />\n<text class=\"legend\" x=\"540.00\" y=\"105.00\">Commercial</text>\n<rect class=\"p2 borderLine\" x=\"500.00\" y=\"130.00\" width=\"30.00\" height=\"15.00\" />\n<rect class=\"q2\" x=\"500.00\" y=\"130.00\" width=\"30.00\" height=\"15.00\" />\n<text class=\"legend\" x=\"540.00\" y=\"145.00\">Government</text>\n<rect class=\"p3 borderLine\" x=\"500.00\" y=\"170.00\" width=\"30.00\" height=\"15.00\" />\n<rect class=\"q3\" x=\"500.00\" y=\"170.00\" width=\"30.00\" height=\"15.00\" />\n<text class=\"legend\" x=\"540.00\" y=\"185.00\">Military</text>\n</g>\n</svg>\n"
attr(,"class")
[1] "list" "nVennObj"
Code
# https://observablehq.com/@simulmedia/basic-venn-diagram
This starts to give us an idea of the kinds of overlaps we can explore! In the case of the US, for example, there’s a significant overlap between commercial use and government/military uses.
Overlaps
If we naturally group our data by users when some are specified jointly, we’ll end with single-use totals that are exclusive of joint use. So for example, it may say that 3 sats are specified as Commercial use and 3 more are specified as Government + Commercial use.
In order to visualise this data with venn.js, we need single-use totals to be inclusive : that 8 are Commercial use and that 3 of those 8 are Government + Commercial.
Code
# do the basic aggregation %>% select (country_owner, users) %>% # make sure we also catch any cases where there's both B/A and A/B listed mutate (users = map_chr (users,~ paste (sort (unlist (str_split (.x, "/" ))), collapse = "/" ))) %>% # sum up group_by (country_owner, users) %>% summarise (size = n ()) %>% ungroup () %>% rename (sets = users) ->
`summarise()` has grouped output by 'country_owner'. You can override using the
`.groups` argument.
Code
# in order to use this data with venn.js, data for single uses needs to be # _inclusive_ of overlaps, not exclusive. so for example, if 3 sats are civil # only and 3 more are gov+civil, cviil should be specified as 8, not 3! %>% separate_rows (sets, sep = "/" ) %>% group_by (country_owner, sets) %>% summarise (size = sum (size, na.rm = TRUE )) %>% ->
`summarise()` has grouped output by 'country_owner'. You can override using the
`.groups` argument.
Code
# now combine exclusive and inclusive use together again %>% filter (str_detect (sets, "/" )) -> sat_joint_use%>% filter (! str_detect (sets, "/" )) -> sat_exclusive_use# rows (country + user) that appear in `sat_inclusive_use` but *not* # `sat_exclusive_use` are totally overlapped # let's mark them so that we can remove their labels later %>% anti_join (sat_exclusive_use, by = c ("country_owner" , "sets" )) %>% mutate (overlapped = TRUE ) -> %>% semi_join (sat_exclusive_use, by = c ("country_owner" , "sets" )) %>% mutate (overlapped = FALSE ) -> %>% mutate (overlapped = FALSE ) %>% bind_rows (sat_inclusive_overlapped) %>% bind_rows (sat_inclusive_notoverlapped) -> write_csv (sat_use_all, here ("data" , "sat-user-overlap-data.csv" ))write_csv (tibble (country = unique (sat_group_data$ country_owner)),here ("data" , "sat-user-overlap-countries.csv" ))
Code
import { countrySelect, overlapChart} from "embed-user-overlaps.qmd" ; ;
Code
Aside: what kinds of purposes do different users put satellites to?
Code
%>% mutate (rowid = 1 : n ()) %>% select (rowid, users, purpose) %>% separate_rows (users, sep = "/" ) %>% separate_rows (purpose, sep = "/" ) %>% filter (! is.na (users), ! is.na (purpose), users != "Earth Observation" ) %>% mutate (purpose = fct_infreq (purpose)) %>% group_by (users, purpose) %>% summarise (n = n ()) %>% ggplot (.) + aes (x = users, y = n, fill = purpose) + geom_col (position = "fill" ) + scale_fill_brewer (type = "qual" ) + scale_y_continuous (labels = scales:: label_percent ()) + coord_flip () + theme_minimal ()
`summarise()` has grouped output by 'users'. You can override using the
`.groups` argument.
Warning in RColorBrewer::brewer.pal(n, pal): n too large, allowed maximum for palette Accent is 8
Returning the palette you asked for with that many colors
