plot(density(andmed[1001:2000]))

Tundub, et silutud ei ole piisavalt. Parema pildi saab näiteks nii.

plot(density(andmed[1001:2000]),ylim=c(0,0.03))

Silutud tihedus ei ole kuigi veenev. See on ka igati ootuspärane, sest kasutatakse rusikareeglit, mis põhineb ootusel, et jaotus on ligikaudu sümmeetriline ja kergete sabadega.

Üks võimalus on proovida veel teisigi  automaatseid silumislaiuse valikuid.

plot(density(andmed[1001:2000],bw="bcv"),ylim=c(0,0.03),from=0)

annab näiteks mõneti veenvama pildi. Lõpliku kandjaga tuum siin ilmselt olukorda parandada ei saaks (peaksime silumislaiuse määrama väga suureks)

Kindlasti oodatav -- üksikud ülejäänutest kaugel asuvad punktid jäävad paratamatult kas põhiosaga liitmata, või siis tekib üpris imelik normaaljaotuse sabaga jaotus.

Lihtne on vaadata korduvalt nt graafikut

plot(density(rgpd(1000),from=0),ylim=c(0,0.01))

Olukord on sageli kaunis sarnane.