From InterSciWiki
cities <- c(8008654,3694742,2896047,1953633,1517550,1321190,1223429,1188589,1151305,951270,895193,781864,776733,735617,711265,656562,651154,650100,596974,589141,572059,563657,563376,557834,553693,545535,541099,529184,506129,487341,484674,479639,477472,461522,448948,441545,428873,425257,416441,407018,399484,397776,393120,391019,382747,371657,362437,360988,351150,348189,337977,334563,332969,331285,328071,313782,303463,292648,286840,282956,277496,275923,272537,260512,260283,256207,255175,248408,243771,243082,242790,240055,234403,228520,226419,225638,224047,222008,220483,219773,218831,217070,215794,208903,203413,202596,201607,200172,199572,199191,199184,198768,197800,197790,196019,195629,195182,194973,193556,191615,190914,189627,189453,188864,187316,186291,185975,185787,185237,184256,183135,181767,180900,180150,177322,176643,175406,174501,173627,173618,173553,172648,170359,166196,165196,161029,158625,158518,158011,155681,154198,152082,151722,151404,151131,150349,149473,149431,149222,147854,
+ 146866,146437,146435,144137,143593,143148,143072,142685,142379,141674,140031,139715,139530,139529,138319,138247,137946,137415,136991,133936,133747,131844,131784,128937,128821,128409,128283,127743,127427,126003,125251,125056,124943,124523,124506,124471,124121,123626,122882,121872,121582,121271,120568,119371,118718,118716,117549,117083,117005,116829,116760,116273,115965,115930,115488,114117,113880,113866,113288,113019,112580,112495,111532,111473,111373,111365,109949,109578,108902,108896,108756,108685,107919,107323,107271,107006,106656,106632,106431,105441,105167,105080,104323,104197,103717,103625,103478,102767,102746,102743,102673,102361,102286,102238,102129,101355,100998,100920,100565,100545,100316,100266,99216,98359,97775,97255,96650,96375,96178,95694,94996,94911,94892,94666,94624,94580,94492,94304,94301,93768,93493,93102,92998,92779,92472,92335,91957,91938,91776,91628,90943,90607,90532,90405,90205,90032,89730,89692,89621,89050,88967,88774,88642,88443,88207,88116,88025,87933,87097,86911,
+ 86605,86602,86319,85979,85932,85808,85794,85779,85655,85616,85403,85172,85071,85025,84662,84362,84329,84112,84110,84084,83829,83629,83259,83048,82951,82934,82549,82376,82111,82026,81991,81855,81207,81111,80959,80874,80734,80537,80431,80268,80142,80114,80056,79928,79904,79554,79452,79426,79345,79269,79062,78721,78306,78296,77962,77753,77737,77685,77478,77423,77248,77145,77021,76947,76760,76415,76042,76021,75837,75720,75702,75402,75402,74848,74764,74239,73990,73539,73402,73344,72958,72878,72739,72684,72670,72664,72259,72182,72069,72043,71538,71393,71329,71283,71127,71118,71021,70705,70702,70517,70296,70126,69845,69824,69567,69473,69368,68991,68904,68825,68747,68652,68556,68393,68381,68315,68263,68015,67861,67624,67561,67504,67388,67249,67173,67088,67083,66869,66787,66715,66455,65936,65920,65894,65868,65660,65273,65269,64903,64583,64296,64249,64196,64112,64029,63900,63870,63828,63677,63603,63557,63428,63325,63317,63296,63285,63261,63096,63011,62749,62698,62590,62406,62226,62150,61968,61842,
+ 61821,61792,61714,61607,61348,61254,61182,60603,60578,60552,60525,60523,60448,60432,60389,60308,60220,60213,60062,60062,60051,60032,59880,59740,59702,59686,59684,59577,59226,58974,58969,58912,58898,58812,58757,58596,58590,58271,58266,58264,58244,58062,57955,57923,57870,57755,57746,57697,57695,57585,57247,57221,57104,56933,56929,56769,56694,56646,56521,56456,56372,56352,56340,56287,56265,56259,56255,56063,56052,55981,55976,55919,55901,55825,55765,55635,55593,55588,55566,55513,55502,55419,55266,55256,55245,55147,54978,54977,54901,54749,54653,54593,54552,54514,54370,54260,54239,54088,54059,54016,53948,53909,53505,53421,53364,53296,53205,53127,53109,53077,53071,53054,52975,52946,52913,52793,52717,52715,52648,52613,52524,52466,52457,52360,52029,52018,51988,51884,51765,51659,51599,51548,51508,51488,51460,51443,51104,51102,50868,50866,50790,50769,50767,50719,50673,50644,50624,50594,50365,50284,50211,49936,49816,49739,49693,49672,49553,49523,49435,49415,49399,49374,49366,49325,49286,49286,49170,
+ 49155,49121,48971,48950,48806,48702,48688,48501,48472,48411,48187,48151,48129,48129,48089,48048,48028,47996,47829,47821,47782,47662,47656,47637,47425,47388,47386,47380,47303,47283,47271,47214,47152,47045,46837,46822,46783,46707,46699,46614,46606,46549,46509,46506,46488,46463,46449,46332,46277,46258,46005,45946,45866,45793,45768,45697,45563,45527,45504,45444,45430,45409,45293,45287,45255,45175,45054,45027,45014,44925,44897,44870,44795,44752,44712,44681,44536,44503,44382,44292,44193,44102,44054,44003,43957,43858,43789,43783,43768,43638,43576,43567,43407,43370,43364,43250,43224,43128,43123,43081,43062,42928,42805,42776,42687,42677,42647,42584,42514,42471,42412,42358,42249,42236,42223,42097,42068,42059,41938,41872,41862,41773,41721,41715,41633,41633,41581,41550,41531,41464,41375,41303,41303,41207,41155,41145,41138,41103,41063,40866,40862,40853,40786,40687,40687,40580,40517,40513,40507,40453,40424,40407,40377,40340,40238,40236,40165,40105,40099,40073,40072,40018,40008,39968,39862,39838,39824,
+ 39804,39803,39695,39695,39573,39515,39503,39394,39375,39199,39196,39139,39116,39102,39074,39065,39065,39040,39018,38978,38915,38829,38816,38723,38649,38586,38577,38533,38446,38420,38392,38357,38327,38278,38224,38217,38216,38176,38141,38138,38101,38064,38041,38037,38025,37984,37882,37835,37784,37775,37688,37674,37648,37647,37644,37414,37336,37326,37259,37257,37251,37208,36944,36929,36854,36817,36722,36688,36664,36582,36474,36458,36430,36395,36368,36346,36330,36316,36282,36264,36255,36234,36189,36145,36117,36081,36028,35958,35918,35901,35866,35863,35838,35814,35811,35794,35690,35597,35521,35476,35462,35450,35444,35404,35361,35335,35318,35278,35258,35202,35189,35135,35110,35102,35094,35086,35083,35080,35048,34980,34950,34789,34748,34727,34701,34665,34649,34527,34496,34424,34310,34282,34211,34129,34113,34104,34077,34014,33998,33904,33892,33867,33862,33858,33857,33852,33826,33784,33769,33740,33714,33712,33690,33686,33588,33556,33496,33464,33436,33377,33262,33204,33147,33097,33058,33049,33007,
+ 32986,32914,32913,32900,32884,32861,32837,32797,32793,32787,32732,32667,32652,32647,32645,32512,32414,32393,32338,32264,32199,32190,32139,32102,32066,32052,32046,32010,31976,31909,31881,31877,31836,31814,31776,31730,31727,31719,31711,31637,31626,31615,31580,31575,31542,31492,31473,31436,31432,31411,31340,31283,31282,31258,31229,31226,31101,31077,30999,30989,30966,30965,30940,30906,30904,30785,30734,30711,30708,30588,30582,30560,30504,30423,30399,30388,30308,30273,30262,30256,30224,30211,30136,30124,30115,30093,30068,30047,30027,30007,30004,30000,29933,29926,29923,29884,29871,29855,29793,29786,29780,29771,29745,29712,29700,29693,29644,29592,29560,29549,29506,29494,29437,29405,29381,29376,29344,29339,29271,29172,29172,29072,29037,28981,28981,28978,28930,28842,28803,28782,28775,28735,28657,28648,28641,28633,28605,28600,28590,28574,28537,28520,28463,28462,28461,28456,28334,28332,28311,28307,28259,28224,28144,28135,28128,28083,28073,28031,28006,27939,27933,27931,27906,27902,27899,27867,27844,
+ 27819,27819,27812,27802,27772,27765,27735,27718,27696,27674,27656,27651,27619,27599,27588,27508,27449,27433,27433,27407,27400,27387,27368,27362,27345,27324,27217,27215,27195,27152,27144,27136,27134,27109,27069,27044,26997,26992,26987,26984,26955,26940,26884,26865,26847,26839,26809,26782,26705,26653,26622,26614,26588,26588,26549,26530,26518,26502,26500,26475,26463,26451,26411,26409,26408,26407,26386,26365,26359,26325,26320,26309,26309,26263,26247,26237,26232,26225,26206,26203,26200,26186,26185,26156,26128,26118,26107,26052,26038,26009,26001,25974,25946,25945,25927,25924,25919,25912,25897,25894,25869,25866,25849,25756,25754,25737,25720,25709,25645,25642,25639,25619,25614,25605,25600,25573,25496,25478,25470,25464,25421,25410,25405,25363,25325,25267,25233,25216,25205,25178,25171,25144,25104,25099,25018,25011,24998,24948,24936,24927,24918,24881,24863,24848,24747,24731,24668,24662,24604,24598,24568,24555,24516,24508,24498,24477,24426,24331,24325,24312,24303,24302,24244,24243,24208,24183,24157,
+ 24157,24146,24127,24100,24071,24069,24040,24031,24015,24010,23971,23970,23936,23929,23908,23908,23853,23779,23749,23748,23727,23720,23694,23669,23635,23609,23585,23574,23564,23562,23562,23546,23545,23537,23532,23526,23516,23510,23503,23470,23463,23462,23453,23376,23370,23347,23329,23327,23324,23283,23243,23242,23232,23222,23208,23203,23171,23171,23154,23135,23133,23122,23113,23083,23060,23033,23029,23007,22979,22976,22917,22904,22856,22825,22801,22795,22782,22771,22763,22725,22722,22701,22697,22666,22643,22641,22628,22621,22564,22563,22551,22542,22527,22505,22501,22464,22451,22442,22441,22439,22439,22406,22387,22376,22338,22336,22302,22279,22277,22265,22261,22242,22228,22206,22193,22152,22147,22141,22128,22107,22105,22097,22091,22076,22072,22060,22052,22036,22035,22027,22021,21987,21967,21961,21946,21945,21908,21869,21859,21845,21842,21810,21802,21786,21779,21751,21735,21675,21672,21662,21659,21631,21602,21576,21556,21536,21530,21522,21519,21511,21511,21498,21488,21479,21447,21445,21403,
+ 21400,21400,21353,21345,21329,21274,21263,21256,21253,21205,21182,21182,21179,21170,21155,21144,21131,21125,21119,21115,21103,21088,21077,21054,21043,21043,20990,20986,20976,20966,20942,20932,20920,20909,20905,20905,20878,20870,20843,20829,20821,20790,20784,20784,20770,20761,20756,20746,20742,20735,20726,20718,20714,20709,20692,20678,20664,20645,20642,20568,20550,20537,20535,20527,20512,20510,20504,20501,20490,20436,20421,20411,20410,20403,20403,20382,20353,20339,20323,20318,20308,20307,20302,20300,20281,20275,20255,20239,20238,20235,20229,20211,20208,20205,20196,20191,20173,20173,20170,20127,20112,20099,20084,20065,20057,20046,20035,20033,20026,20025,19999,19998,19973,19968,19953,19935,19911,19900,19869,19837,19805,19725,19712,19688,19684,19666,19652,19628,19615,19584,19577,19543,19522,19520,19515,19489,19488,19478,19472,19469,19463,19456,19434,19405,19386,19377,19363,19355,19310,19310,19307,19306,19297,19291,19279,19254,19234,19201,19196,19161,19130,19124,19085,19068,19064,19055,19051,
+ 19039,19012,19010,19001,18990,18988,18985,18963,18944,18931,18928,18925,18925,18919,18918,18886,18872,18822,18816,18792,18786,18740,18714,18708,18687,18681,18678,18675,18667,18635,18596,18566,18559,18555,18555,18554,18544,18520,18516,18515,18510,18505,18494,18491,18464,18464,18460,18429,18425,18420,18397,18393,18376,18370,18355,18353,18300,18291,18273,18260,18218,18216,18188,18165,18165,18162,18153,18147,18145,18130,18127,18112,18110,18110,18101,18096,18080,18076,18019,18002,17976,17970,17962,17899,17887,17875,17862,17861,17823,17806,17783,17757,17720,17720,17712,17707,17707,17665,17642,17634,17617,17599,17598,17582,17581,17572,17534,17503,17499,17477,17452,17450,17447,17444,17423,17422,17414,17411,17367,17367,17361,17357,17350,17342,17334,17320,17320,17312,17308,17304,17299,17298,17278,17268,17254,17243,17241,17225,17189,17181,17179,17172,17147,17127,17093,17081,17080,17073,17071,17048,17026,17024,17004,16992,16985,16955,16948,16930,16925,16897,16865,16849,16838,16824,16822,16812,16810,
+ 16799,16799,16791,16759,16754,16746,16743,16737,16718,16706,16687,16676,16654,16651,16629,16608,16592,16582,16538,16495,16461,16456,16444,16436,16426,16420,16413,16411,16400,16391,16390,16381,16380,16367,16354,16348,16347,16340,16337,16332,16330,16322,16318,16290,16286,16276,16263,16261,16256,16247,16212,16169,16157,16142,16136,16136,16123,16117,16116,16111,16103,16101,16087,16073,16071,16054,16051,16033,16013,16000,15994,15967,15967,15942,15935,15933,15931,15928,15922,15919,15889,15888,15886,15881,15879,15877,15876,15859,15851,15851,15833,15804,15798,15770,15761,15741,15733,15713,15696,15688,15681,15680,15643,15642,15640,15608,15605,15605,15600,15581,15559,15550,15549,15544,15542,15531,15527,15522,15521,15507,15503,15503,15497,15485,15477,15469,15460,15445,15443,15441,15416,15413,15411,15408,15404,15378,15371,15369,15362,15358,15355,15347,15347,15345,15341,15335,15322,15322,15315,15312,15297,15288,15279,15268,15266,15265,15264,15259,15244,15241,15225,15225,15200,15181,15177,15166,15157,
+ 15128,15123,15121,15109,15074,15055,15054,15053,15049,15034,15026,14966,14961,14955,14951,14949,14940,14940,14913,14911,14895,14880,14879,14879,14863,14852,14849,14835,14831,14830,14828,14827,14816,14789,14788,14781,14779,14774,14774,14734,14733,14701,14692,14681,14674,14665,14643,14619,14617,14610,14601,14600,14599,14597,14583,14582,14582,14579,14573,14572,14560,14556,14535,14532,14523,14518,14502,14458,14453,14444,14440,14431,14424,14418,14416,14405,14389,14377,14361,14355,14354,14354,14351,14344,14330,14326,14324,14315,14302,14263,14254,14248,14241,14233,14226,14223,14214,14212,14211,14209,14200,14194,14187,14173,14166,14163,14151,14146,14127,14123,14122,14114,14108,14099,14054,14018,14008,14004,14001,13999,13977,13948,13934,13926,13908,13907,13897,13895,13894,13876,13868,13864,13859,13855,13848,13844,13841,13835,13822,13815,13806,13804,13799,13783,13782,13780,13777,13772,13763,13753,13747,13728,13717,13712,13710,13710,13682,13682,13668,13664,13643,13620,13617,13613,13607,13601,13600,
+ 13595,13594,13589,13575,13568,13556,13546,13532,13532,13531,13528,13521,13518,13485,13484,13472,13469,13453,13452,13442,13422,13417,13416,13382,13372,13368,13367,13365,13347,13340,13337,13336,13329,13309,13295,13292,13281,13276,13272,13264,13263,13257,13255,13250,13249,13225,13222,13214,13213,13203,13184,13176,13162,13158,13156,13152,13151,13150,13147,13140,13134,13131,13129,13122,13098,13085,13082,13069,13067,13038,13034,13025,13023,13022,13019,12989,12988,12987,12979,12978,12968,12956,12949,12946,12945,12945,12938,12937,12935,12934,12918,12911,12910,12881,12873,12872,12871,12868,12866,12860,12841,12840,12811,12811,12808,12808,12790,12781,12780,12738,12732,12731,12727,12721,12717,12716,12710,12708,12703,12697,12680,12666,12659,12659,12648,12620,12613,12608,12604,12603,12600,12600,12594,12589,12585,12576,12575,12568,12568,12559,12555,12545,12534,12533,12532,12530,12527,12513,12512,12507,12501,12493,12485,12481,12473,12471,12466,12462,12461,12449,12443,12437,12433,12429,12423,12422,12419,
+ 12411,12401,12396,12396,12391,12378,12377,12377,12374,12364,12359,12331,12327,12311,12288,12268,12268,12264,12259,12255,12253,12235,12206,12204,12193,12189,12186,12179,12178,12172,12170,12156,12153,12145,12142,12128,12120,12095,12084,12083,12080,12067,12066,12065,12064,12049,12035,12032,12023,12020,12014,12009,11998,11977,11968,11967,11949,11947,11941,11939,11931,11929,11928,11923,11919,11912,11907,11901,11895,11893,11893,11887,11887,11878,11871,11864,11857,11855,11853,11850,11844,11842,11835,11820,11808,11807,11798,11797,11788,11779,11775,11773,11769,11762,11756,11746,11743,11742,11730,11698,11695,11692,11686,11686,11682,11681,11668,11662,11660,11659,11647,11635,11626,11624,11621,11620,11619,11617,11615,11614,11606,11586,11585,11583,11567,11566,11558,11556,11546,11536,11528,11526,11524,11517,11514,11507,11504,11500,11499,11491,11484,11477,11476,11458,11448,11446,11440,11438,11434,11427,11427,11423,11418,11415,11396,11392,11391,11385,11381,11380,11375,11373,11369,11364,11327,11320,11318,
+ 11313,11312,11311,11306,11306,11303,11301,11296,11293,11292,11289,11283,11275,11273,11262,11260,11260,11253,11251,11248,11237,11229,11212,11209,11204,11201,11188,11182,11168,11155,11148,11142,11139,11138,11135,11124,11119,11108,11102,11096,11090,11087,11060,11056,11047,11044,11036,11035,11034,11008,11006,10999,10987,10987,10986,10979,10978,10976,10970,10955,10952,10946,10945,10944,10943,10941,10934,10930,10927,10927,10922,10919,10918,10903,10902,10895,10889,10880,10876,10874,10871,10867,10866,10861,10858,10854,10853,10851,10849,10848,10842,10839,10835,10830,10825,10807,10800,10796,10783,10774,10767,10764,10763,10762,10748,10735,10723,10718,10714,10706,10705,10697,10668,10654,10640,10637,10637,10635,10622,10622,10610,10597,10592,10586,10585,10582,10582,10581,10578,10565,10563,10558,10550,10545,10540,10539,10535,10534,10529,10523,10515,10510,10507,10500,10498,10498,10486,10477,10473,10472,10467,10466,10458,10458,10445,10442,10437,10436,10434,10427,10424,10423,10422,10419,10411,10401,10389,
+ 10387,10380,10377,10375,10375,10374,10366,10360,10355,10350,10344,10343,10340,10331,10311,10307,10302,10300,10299,10297,10296,10295,10290,10278,10266,10260,10259,10255,10255,10255,10254,10242,10238,10232,10228,10207,10192,10182,10175,10174,10167,10163,10157,10155,10146,10122,10117,10116,10105,10103,10098,10087,10085,10082,10076,10074,10070,10063,10056,10050,10049,10042,10036,10034,10029,10011,10010,10001,9995,9992,9991,9991,9989,9986,9983,9975,9974,9952,9952,9929,9925,9923,9920,9917,9916,9896,9892,9888,9887,9885,9880,9873,9862,9861,9860,9860,9859,9850,9850,9848,9847,9844,9843,9839,9835,9826,9821,9813,9812,9810,9808,9800,9798,9798,9796,9792,9788,9774,9771,9764,9763,9761,9745,9743,9737,9736,9735,9733,9730,9728,9727,9722,9710,9708,9705,9703,9686,9684,9682,9679,9676,9670,9667,9664,9657,9653,9652,9652,9652,9640,9637,9628,9628,9622,9621,9621,9618,9614,9609,9602,9602,9600,9596,9592,9591,9577,9558,9557,9552,9548,9544,9541,9538,9532,9523,9523,9520,9519,9513,9511,9501,9501,9488,9488,9481,9476,9471,
+ 9470,9460,9451,9449,9445,9443,9441,9435,9426,9424,9424,9420,9417,9415,9413,9409,9405,9400,9400,9389,9388,9377,9371,9365,9363,9352,9338,9333,9331,9323,9321,9320,9311,9311,9305,9293,9291,9282,9280,9279,9279,9279,9278,9275,9275,9268,9266,9264,9257,9254,9247,9245,9242,9239,9239,9237,9237,9231,9228,9224,9222,9215,9214,9214,9212,9207,9206,9198,9189,9189,9188,9182,9175,9168,9157,9156,9149,9149,9146,9146,9144,9139,9131,9105,9100,9099,9099,9098,9097,9094,9090,9089,9086,9083,9081,9081,9056,9052,9042,9038,9029,9027,9022,9019,9017,9015,9012,9008,9008,9006,9006,9000,8998,8994,8993,8989,8987,8984,8979,8973,8969,8968,8960,8959,8958,8934,8931,8923,8914,8912,8910,8901,8895,8883,8881,8877,8876,8870,8870,8869,8860,8842,8838,8836,8835,8833,8829,8828,8811,8811,8806,8803,8803,8798,8795,8793,8789,8771,8770,8762,8758,8757,8751,8749,8748,8747,8737,8737,8733,8733,8731,8721,8716,8716,8714,8711,8708,8705,8704,8702,8695,8694,8693,8683,8674,8673,8670,8669,8667,8664,8662,8662,8658,8656,8654,8647,8647,8629,8626,8625,
+ 8620,8618,8610,8610,8609,8604,8600,8596,8595,8592,8591,8586,8586,8583,8578,8576,8571,8570,8562,8558,8556,8551,8551,8542,8540,8533,8531,8524,8516,8506,8498,8496,8495,8491,8489,8489,8488,8488,8487,8485,8463,8460,8442,8442,8438,8437,8432,8428,8426,8419,8413,8410,8409,8408,8405,8402,8399,8397,8389,8385,8385,8383,8382,8382,8378,8377,8376,8368,8364,8363,8357,8356,8354,8348,8346,8342,8342,8339,8334,8329,8314,8313,8312,8312,8303,8297,8297,8288,8286,8276,8276,8267,8264,8261,8258,8252,8250,8249,8244,8243,8242,8237,8227,8209,8206,8195,8193,8193,8192,8192,8188,8186,8183,8180,8178,8177,8175,8174,8170,8166,8165,8155,8154,8154,8153,8152,8144,8139,8137,8133,8132,8130,8129,8127,8123,8123,8104,8103,8098,8098,8089,8082,8080,8053,8052,8047,8047,8044,8042,8040,8038,8037,8035,8035,8029,8024,8023,8016,8014,8012,8009,8006,8000,7996,7996,7991,7989,7988,7985,7978,7975,7967,7962,7957,7957,7951,7948,7943,7943,7939,7936,7923,7920,7916,7916,7915,7913,7912,7901,7897,7896,7894,7888,7880,7870,7862,7861,7859,7857,7856,
+ 7853,7846,7845,7845,7844,7843,7837,7835,7829,7829,7825,7824,7823,7820,7817,7816,7812,7807,7804,7804,7799,7796,7793,7783,7782,7780,7780,7776,7774,7774,7772,7769,7765,7765,7761,7761,7759,7746,7744,7741,7741,7734,7734,7731,7728,7725,7723,7723,7723,7722,7715,7714,7712,7709,7698,7696,7696,7695,7694,7693,7678,7676,7676,7676,7675,7668,7666,7661,7661,7657,7655,7644,7643,7640,7638,7636,7633,7618,7614,7610,7610,7609,7609,7608,7606,7606,7605,7605,7603,7601,7597,7596,7595,7591,7589,7589,7587,7584,7579,7578,7576,7572,7570,7568,7563,7561,7554,7553,7547,7544,7538,7527,7524,7523,7516,7513,7511,7508,7506,7505,7502,7498,7495,7494,7490,7487,7476,7467,7465,7459,7459,7456,7452,7450,7450,7442,7438,7437,7436,7435,7435,7433,7430,7424,7422,7420,7407,7405,7403,7401,7400,7396,7393,7389,7388,7381,7375,7372,7371,7370,7356,7349,7347,7345,7341,7340,7338,7337,7336,7334,7332,7330,7325,7320,7319,7319,7315,7307,7305,7290,7289,7289,7288,7286,7281,7280,7277,7274,7271,7271,7270,7270,7266,7263,7261,7259,7257,7247,7246,7245,
+ 7244,7239,7231,7228,7228,7226,7224,7212,7210,7204,7197,7196,7194,7190,7174,7174,7174,7173,7172,7172,7170,7160,7159,7149,7149,7148,7145,7145,7141,7139,7139,7137,7135,7135,7132,7129,7127,7126,7125,7115,7111,7109,7106,7096,7095,7093,7093,7091,7091,7090,7090,7084,7080,7068,7064,7063,7057,7047,7042,7039,7036,7034,7034,7032,7027,7025,7022,7021,7015,7012,7009,7006,7000,6997,6994,6991,6989,6985,6982,6980,6978,6975,6973,6970,6964,6959,6956,6951,6950,6947,6944,6942,6928,6927,6924,6920,6918,6913,6911,6905,6903,6901,6895,6894,6890,6887,6884,6882,6880,6880,6871,6868,6867,6866,6863,6862,6862,6862,6861,6849,6847,6847,6846,6845,6842,6837,6835,6834,6834,6834,6831,6831,6830,6830,6828,6826,6825,6823,6822,6821,6820,6819,6818,6818,6815,6814,6812,6809,6808,6807,6806,6802,6802,6796,6796,6792,6789,6789,6789,6781,6781,6779,6773,6772,6765,6764,6762,6760,6757,6753,6749,6744,6744,6741,6737,6733,6733,6730,6729,6720,6716,6713,6712,6711,6710,6710,6709,6706,6704,6702,6699,6698,6696,6695,6692,6688,6686,6682,6681,6681,
+ 6676,6674,6673,6670,6668,6664,6662,6660,6658,6656,6653,6650,6650,6649,6645,6644,6644,6641,6635,6635,6632,6631,6631,6628,6626,6626,6622,6619,6619,6607,6602,6595,6595,6594,6592,6589,6588,6586,6585,6579,6576,6571,6571,6571,6569,6568,6563,6562,6561,6558,6554,6551,6543,6540,6537,6537,6537,6536,6535,6534,6530,6522,6518,6518,6508,6507,6502,6492,6486,6485,6485,6482,6480,6480,6480,6478,6476,6475,6473,6472,6471,6470,6465,6463,6460,6459,6457,6456,6452,6451,6448,6446,6440,6438,6437,6433,6433,6433,6433,6432,6426,6426,6419,6415,6414,6412,6410,6408,6407,6400,6396,6396,6395,6395,6392,6391,6390,6388,6386,6385,6383,6381,6380,6380,6372,6369,6364,6363,6361,6355,6355,6350,6349,6348,6347,6345,6338,6334,6331,6324,6323,6319,6318,6314,6314,6313,6311,6310,6310,6310,6310,6306,6303,6302,6299,6299,6298,6293,6288,6286,6284,6281,6272,6270,6266,6260,6260,6256,6256,6252,6249,6249,6239,6235,6231,6229,6224,6223,6220,6220,6217,6216,6216,6216,6212,6211,6208,6208,6205,6204,6198,6198,6189,6185,6183,6180,6179,6179,6170,6170,
+ 6168,6164,6164,6163,6160,6151,6149,6149,6148,6148,6145,6143,6138,6137,6131,6125,6120,6116,6116,6115,6115,6110,6106,6099,6098,6098,6098,6097,6097,6093,6080,6078,6078,6077,6076,6076,6075,6070,6065,6064,6062,6056,6049,6048,6045,6040,6038,6035,6034,6033,6033,6025,6023,6022,6020,6018,6017,6016,6015,6014,6013,6005,6004,6002,6002,5997,5996,5995,5993,5992,5991,5990,5989,5989,5987,5986,5985,5984,5982,5981,5980,5969,5966,5965,5963,5958,5958,5958,5955,5952,5952,5951,5950,5947,5942,5939,5937,5936,5929,5925,5924,5917,5917,5916,5915,5910,5909,5907,5906,5905,5904,5900,5900,5899,5897,5887,5877,5873,5869,5868,5866,5864,5863,5863,5860,5858,5858,5857,5856,5853,5852,5852,5851,5850,5847,5844,5843,5836,5836,5822,5817,5813,5813,5812,5809,5807,5806,5805,5802,5800,5795,5788,5786,5784,5780,5779,5778,5776,5774,5767,5766,5765,5764,5759,5751,5745,5741,5740,5740,5739,5735,5730,5728,5727,5723,5720,5717,5717,5716,5714,5714,5712,5711,5700,5699,5696,5693,5691,5690,5688,5687,5686,5682,5681,5676,5676,5676,5676,5675,5674,
+ 5670,5669,5667,5666,5665,5665,5659,5659,5655,5654,5653,5649,5645,5641,5638,5636,5634,5634,5634,5633,5630,5628,5626,5626,5624,5624,5623,5621,5620,5614,5613,5612,5607,5606,5605,5603,5600,5599,5596,5590,5589,5589,5588,5586,5586,5585,5584,5583,5582,5580,5576,5576,5575,5575,5574,5573,5572,5567,5562,5561,5558,5556,5549,5548,5546,5539,5538,5537,5537,5537,5536,5533,5533,5528,5527,5523,5523,5514,5502,5502,5500,5500,5498,5498,5496,5490,5489,5488,5488,5486,5485,5483,5482,5482,5478,5471,5470,5470,5470,5469,5468,5468,5466,5464,5463,5462,5460,5459,5458,5458,5456,5456,5455,5450,5448,5445,5441,5439,5438,5436,5436,5436,5434,5433,5430,5429,5424,5422,5419,5414,5413,5409,5409,5408,5408,5406,5395,5391,5389,5384,5383,5379,5376,5376,5376,5376,5373,5371,5369,5368,5359,5356,5352,5344,5344,5341,5336,5333,5332,5332,5332,5331,5330,5322,5319,5318,5314,5314,5312,5308,5307,5307,5305,5303,5302,5301,5298,5297,5295,5295,5294,5294,5292,5291,5289,5288,5287,5286,5285,5282,5280,5278,5276,5274,5273,5264,5262,5257,5255,5250,
+ 5248,5247,5245,5236,5236,5234,5233,5231,5228,5227,5221,5219,5218,5218,5216,5216,5213,5213,5204,5203,5203,5203,5197,5193,5192,5192,5190,5189,5189,5188,5187,5184,5182,5181,5180,5178,5177,5174,5173,5166,5160,5157,5156,5156,5156,5152,5151,5148,5146,5144,5143,5138,5137,5136,5136,5135,5135,5133,5132,5132,5131,5130,5128,5127,5127,5126,5126,5124,5121,5119,5119,5117,5117,5114,5111,5109,5107,5107,5105,5105,5105,5104,5102,5102,5099,5099,5098,5097,5092,5092,5088,5088,5088,5087,5083,5083,5080,5080,5078,5076,5073,5072,5070,5067,5066,5065,5062,5061,5059,5058,5057,5057,5056,5056,5054,5054,5049,5043,5039,5035,5034,5032,5032,5032,5029,5028,5021,5019,5019,5017,5005,5004)
cities.min <- 52640
cities.tailprob <- sum(cities>=cities.min)/length(cities)
cities.pareto <- pareto.fit(cities,cities.min)
#cities.random <- rpareto(sum(cities>=cities.min),threshold=cities.min,cities.pareto$exponent)
#cities.random <- rpareto(1000000,threshold=cities.min,cities.pareto$exponent)
cities.random <- cities.min*((1-ppareto(cities,cities.min,cities.pareto$exponent,lower.tail=FALSE))^(1/(1-cities.pareto$exponent)))
#cities.tsal <- tsal.fit(cities,xmin=5000)
cities.tsal <- tsal.fit(ifelse(cities<=cities.min,cities,cities.random),xmin=5000)
plot.survival.loglog(cities,ylab="Cumulative probability",main="Survival function of cities data",sub="zagged=empirical, solid=Tsallis, dotted=Pareto, xq-min=5000, xp-min=52640")
curve(ptsal(x,cities.tsal$shape,cities.tsal$scale,cities.tsal$q,cities.tsal$kappa,cities.tsal$xmin,lower.tail=FALSE),add=TRUE,col="blue")
curve(cities.tailprob*ppareto(x,cities.min,cities.pareto$exponent,lower.tail=FALSE),from=cities.min,col="red",lty="dotted",add=TRUE)
cities.tsal
$type
[1] "tsal"
$q
[1] 1.792614
$kappa
[1] 5022.802
$shape
[1] 1.261648
$scale
[1] 6337.007
$loglike
[1] -47365.01
$n
[1] 4347
$xmin
[1] 5000
$method
[1] "mle.equation"
