Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- ```{r}
- ratio = function(k, alpha){
- (k / (k+1))**(1 + alpha)
- }
- x1=1:20
- x2 = 2:21
- par(mfrow=c(1, 3))
- plot(x2, zipf(x2, alpha=0.5, Cte=2.612) / zipf(x1, alpha=0.5, Cte=2.612), type='h', ylab="P(K+1) / P(K)", xlab="alpha=0.5", xlim=c(1, 20))
- lines(x1, ratio(x1, 0.5), col='red')
- plot(x2, zipf(x2, alpha=1, Cte=1.645) / zipf(x1, alpha=1, Cte=1.645), type='h', ylab="P(K+1) / P(K)", xlab="alpha=1", xlim=c(1, 20))
- lines(x1, ratio(x1, 1), col='red')
- plot(x2, zipf(x2, alpha=2, Cte=1.202) / zipf(x1, alpha=2, Cte=1.202), type='h', ylab="P(K+1) / P(K)", xlab="alpha=2", xlim=c(1, 20))
- lines(x1, ratio(x1, 2), col='red')
- ```
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement