Kasuta andmevektori 300 esimest vaatlust. Mäletame, et üks vaatlustest ületas teisi oluliselt. Kas aga valimi sees tehtud rekordite arv tundub vastavat IID realisatsioonidele?
recordsPlot(andmed[1:300])
Rekordite arv on küll keskmiselt oodatavast mõneti suurem, ent siiski igati oodatavas skaalas. Seega ei ole selle graafiku põhjal IID eelduses kahelda.
Mõlemad autokorrelatsionifunktsioonid mahuvad ilusasti IID korral oodatava usalduspiiridesse.
Praktiline ülesanne
Lõpetuseks väike üllatus. Kas need andmed võiksid olla pärit üldistatud Pareto jaotusest?
gpdFit(andmed[1:300],u=0)
ks.test(andmed[1:300],pgpd,xi=0.41,beta=1.65)
Võimalikke kontrollimise viise on muidugi mitmeid, ent kui juba "valesti kasutatud" Kolmogorov-Smirnovi test ütleb, et tegu ei ole sellest jaotusest pärit andmetega, siis võib selles küll üpris kindel olla.