Birds.q

From InterSciWiki

birds <- c(138705,98611,83906,77656,69928,69745,59698,57377,42162,40134,35594,32482,30790,29974,28293,27548,27344,26177,21714,21168,20120,19133,18944,17831,15604,15192,14029,13987,13965,13862,12988,12920,12253,11825,11358,11313,11163,10975,10968,10851,10127,10094,9946,9860,9826,9774,9746,9270,9236,9134,9093,8671,8623,8436,8205,8046,7914,7869,7791,7266,7202,7071,6876,6793,6733,6679,6213,6101,6043,5964,5963,5860,5734,5530,5485,5464,5341,5229,5093,5061,5041,4956,4896,4859,4814,4811,4792,4783,4665,4501,4501,4411,4384,4206,4137,3955,3911,3883,3875,3782,3744,3728,3721,3522,3520,3480,3474,3471,3462,3461,3435,3425,3398,3344,3341,3312,3287,3271,3243,3194,3139,3111,3059,3039,3037,2899,2893,2797,2727,2705,2600,2595,2558,2557,2514,2501,2467,2442,2386,2363,2356,2346,2278,2262,2229,2229,2183,2158,2109,2047,2037,2035,2000,1973,1952,1923,1908,1886,1873,1868,1853,1845,1845,1770,1760,1754,1737,1684,1670,1666,1655,1654,1619,1617,1614,1539,1530,1526,1479,1458,1438,1428,1408,1380,1372,1361,1346,1344,1336,

+ 1282,1220,1217,1198,1161,1133,1128,1116,1069,1068,1044,1028,1012,1007,1004,1000,971,945,938,937,925,924,915,898,882,879,879,877,870,869,863,857,855,854,850,849,814,813,812,782,760,752,750,750,746,728,722,715,710,696,690,685,684,677,677,656,651,632,632,628,627,
+ 622,608,603,596,587,573,561,557,552,549,548,542,539,538,538,535,535,524,519,516,507,506,500,500)
birds.min <- 6679 
birds.tailprob <- sum(birds>=birds.min)/length(birds) 
birds.pareto <- pareto.fit(birds,birds.min) 
#birds.cdf <- birds.min*((1-ppareto(birds,birds.min,birds.pareto$exponent,lower.tail=FALSE))^(1/(1-birds.pareto$exponent)))
birds.tsal <- tsal.fit(birds,xmin=500)
#birds.tsal <- tsal.fit(ifelse(birds<=birds.min,birds,birds.cdf),xmin=500)
plot.survival.loglog(birds,ylab="Cumulative probability",main="Survival function of birds data", sub="zagged=empirical, solid=Tsallis, dashed=Pareto, ") 
curve(ptsal(x,birds.tsal$shape,birds.tsal$scale,lower.tail=FALSE),add=TRUE,col="blue") 
curve(birds.tailprob*ppareto(x,birds.min,birds.pareto$exponent,lower.tail=FALSE),from=birds.min,col="red",lty="dashed",add=TRUE) 
birds.tsal

Not a good result. Had a better on earlier.