График для иллюстрации закона Зипфа (Zipf's law)

* График, иллюстрирующий действие закона Зипфа (Ципфа) - Zipf's law (см. http://www.nslij-genetics.org/wli/zipf/). 
* Требуется отобразить логарифм частоты против логарифма ранга частоты.

GET FILE='c:\\program files\\spss\\world95.sav'.

* Допустим, желаем изучить частотное распределение религий в странах мира.
AGGREGATE
  /OUTFILE=*
  /BREAK=religion
  /freq = N(region).

RANK
  VARIABLES=freq  (D) /RANK INTO rank1 /PRINT=YES
  /TIES=MEAN .

COMPUTE lnfreq=LN( freq).
COMPUTE lnrank=LN(rank1).

GRAPH
  /SCATTERPLOT(BIVAR)=lnrank WITH lnfreq
  /MISSING=LISTWISE.