[R] fit a distribution

[Simone Gabbriellini]

hello list,

I have this array of distribution of frequencies:

freq<-c(45,18,10,8,13,5,9,3,4,1,2,2,2,3,4,4,1,1,2,1,2,1,1,1,1,1,2,3,1,1,1,1,1,1,1,1,2,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)

for this array of values:

values<-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,16,18,19,20,21,24,25,26,27,28,29,32,34,35,37,39,40,42,43,44,45,46,54,61,66,70,72,75,78,80,87,90,95,101,110,169,175,193,204,209,218,269,284,290,299,310,345,403,506,520,523,842,981,1034)

I would like to fit this distribution with a theoretical model, and i need a help to select the best one... is this an example of Zipf's law? 

thanks in advance for any help.

best regards,
Simone