Books.q

From InterSciWiki

books <-c(19076822,11325299,11000000,9919785,8065398,8061812,7000000,6978211,6578314,6326470,6000000,5899000,5563841,5490651,5473710,5416857,5390105,5363909,5089472,5000000,5000000,4988225,4916074,4873563,4835966,4828044,4810418,4712588,4637734,4466200,4400000,4375620,4186935,4171838,4007158,4000000,3913341,3910155,3858948,3821608,3800000,3768144,3662089,3646004,3642411,3633467,3593268,3581210,3577729,3549276,3501866,3501281,3500000,3499948,3475947,3445000,3430505,3426000,3389940,3316791,3283000,3262193,3220246,3216600,3208361,3190334,3171512,3170056,3157592,3156000,3142184,3130306,3100000,3087532,3079790,3075123,3007000,3000000,3000000,3000000,2974030,2950807,2940211,2925268,2911111,2833993,2830061,2825368,2820313,2816772,2816662,2816028,2814394,2786223,2778577,2773841,2763486,2760000,2750000,2739427,2719500,2716816,2702597,2700261,2680597,2672065,2667977,2655208,2635250,2634472,2632358,2625000,2621842,2611000,2609236,2600000,2595125,2585397,2560806,2530964,2527756,2527459,2524649,

+ 2521250,2505000,2499660,2472710,2471140,2469306,2461874,2437336,2435928,2433154,2414460,2409163,2403403,2403269,2400000,2400000,2370720,2358945,2341531,2340000,2336004,2330000,2318230,2316989,2312378,2311904,2310000,2299778,2286000,2282322,2266268,2262024,2253982,2253453,2251960,2246396,2244552,2237449,2230000,2229000,2219645,2214669,2213167,2192133,2184145,2183109,2160481,2153000,2152210,2149000,2147000,2146812,2136177,2134289,2133810,2121708,2111873,2106804,2105951,2102522,2100908,2089523,2087837,2087173,2080985,2075000,2073434,2068000,2063775,2054928,2053892,2051488,2050000,2033679,2029248,2024107,2011251,2007000,2007000,2006808,2005000,2000000,2000000,2000000,2000000,1983000,1968799,1968297,1956000,1941769,1933949,1932000,1930000,1929391,1928513,1926110,1920000,1919214,1918865,1900000,1894000,1888785,1886465,1874961,1863625,1863260,1862387,1853917,1843756,1842000,1836815,1834704,1830060,1829364,1827493,1826348,1820308,1819000,1818649,1816000,1805277,1805000,1801097,1794020,1793366,
+ 1783923,1778831,1777027,1758760,1755011,1755000,1754975,1751498,1750000,1748000,1747908,1740930,1738597,1736299,1730000,1728151,1725939,1724658,1720328,1715993,1710000,1706019,1703816,1702500,1700000,1696589,1692921,1691878,1690126,1686367,1682891,1676588,1671968,1661578,1660000,1653742,1652235,1651713,1642905,1640314,1640000,1631630,1626135,1625000,1620464,1617462,1616783,1612339,1611700,1611007,1609011,1601079,1600000,1594000,1591489,1589451,1588972,1586529,1586260,1584793,1584135,1583612,1583561,1582738,1581685,1578234,1578169,1575396,1575304,1573479,1571570,1568601,1563117,1560980,1558878,1558175,1550223,1550000,1550000,1548471,1547000,1546000,1544166,1543000,1540572,1539524,1536793,1536764,1534371,1531786,1529095,1525000,1525000,1520000,1516156,1515880,1515825,1507502,1507000,1500000,1500000,1500000,1500000,1496118,1493167,1488097,1480130,1480000,1479047,1473496,1469432,1469428,1468250,1467635,1464357,1462045,1457000,1455452,1455110,1450000,1436000,1436000,1434203,1433477,1432000,
+ 1430000,1425000,1422000,1421157,1421000,1415572,1415551,1412000,1409790,1406700,1405936,1400000,1398445,1396015,1396000,1390895,1388320,1388220,1383440,1381956,1378500,1377373,1375000,1373288,1369660,1366600,1363909,1362139,1357000,1356863,1354347,1354000,1351993,1351254,1350620,1349000,1338950,1336425,1333417,1327895,1325000,1322402,1320000,1319610,1318000,1317137,1316228,1314077,1312700,1311484,1310000,1308883,1308014,1303028,1300201,1300000,1300000,1296140,1291961,1290453,1290291,1286463,1285000,1282414,1280000,1279518,1279000,1279000,1278951,1275953,1272563,1269007,1266957,1266206,1265336,1265000,1262676,1261000,1260725,1259390,1253793,1250000,1250000,1250000,1247000,1246693,1244400,1243307,1241743,1241550,1239916,1237767,1237060,1235000,1235000,1229000,1226720,1225498,1221986,1220874,1217322,1216938,1212426,1211394,1210994,1210900,1209029,1203889,1200190,1200000,1200000,1200000,1200000,1200000,1200000,1200000,1200000,1194440,1194000,1193650,1193492,1193330,1193230,1191374,1190086,
+ 1188395,1183560,1183063,1177618,1177060,1175323,1173500,1170042,1168181,1165499,1163900,1153749,1152000,1151029,1150000,1148669,1148260,1146449,1143675,1142000,1139572,1138891,1137563,1136917,1132856,1131270,1131018,1128840,1126365,1123136,1121000,1120000,1119134,1118800,1117418,1115657,1111520,1110000,1108000,1107754,1107064,1104600,1100000,1100000,1100000,1100000,1100000,1100000,1100000,1099270,1099000,1098001,1097716,1095155,1095000,1094941,1090383,1090000,1089767,1088744,1087000,1083623,1081906,1081810,1079655,1079000,1077453,1076380,1070392,1068363,1068139,1068004,1068000,1063000,1059000,1058890,1057837,1056992,1056196,1055923,1055000,1054910,1054500,1053700,1050000,1049006,1048247,1047991,1045000,1040000,1039000,1034716,1033525,1032681,1031180,1028109,1028000,1027505,1026322,1021233,1018942,1018000,1017351,1016784,1011000,1010437,1010000,1010000,1010000,1009606,1008937,1007997,1007136,1005381,1005374,1005000,1005000,1002754,1001000,1000054,1000000,1000000,1000000,1000000,1000000,
+ 1000000,1000000,1000000,1000000,1000000,1000000,1000000,1000000,1000000,1000000) 

books.min <- 2400000

books.tailprob <- sum(books>=books.min)/length(books) 
books.pareto <- pareto.fit(books,books.min) 
#books.random <- rpareto(sum(books>=books.min),threshold=books.min,books.pareto$exponent)
#books.random <- rpareto(10000000,threshold=books.min,books.pareto$exponent)
books.random <- books.min*((1-ppareto(books,books.min,books.pareto$exponent,lower.tail=FALSE))^(1/(1-books.pareto$exponent)))
#books.tsal <- tsal.fit(books,xmin=1000000)
books.tsal <- tsal.fit(ifelse(books<=books.min,books,books.random),xmin=1000000)
plot.survival.loglog(books,ylab="Cumulative probability",main="Survival function of books data",sub="black=empirical, blue=Tsallis, red=Pareto") 
curve(ptsal(x,books.tsal$shape,books.tsal$scale,books.tsal$q,books.tsal$kappa,books.tsal$xmin,lower.tail=FALSE),add=TRUE,col="blue") 
curve(books.tailprob*ppareto(x,books.min,books.pareto$exponent,lower.tail=FALSE),from=books.min,col="red",add=TRUE) 
books.tsal

$type [1] "tsal"

$q [1] 1.258711

$kappa [1] 466754.9

$shape [1] 3.865321

$scale [1] 1804158

$loglike [1] -9338.826

$n [1] 633

$xmin [1] 1e+06

$method [1] "mle.equation"

>