Skip to content

Instantly share code, notes, and snippets.

@DavidRdgz
Last active March 27, 2019 21:59
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save DavidRdgz/c84da11ed6981da013b2358145ae7b4a to your computer and use it in GitHub Desktop.
Save DavidRdgz/c84da11ed6981da013b2358145ae7b4a to your computer and use it in GitHub Desktop.
Minimal, Complete, and Verifiable Example Zipf on List[Int]
7.96 |
| · · · · ·· ·· ·
| ·· ·· · · · · · ·· · ·
| · · ·· · · · ········· ······ ··· ·· ·
| · ·· · ·· ·· ····· ······················ ·
7.62 | · · · ········································ ··· ·
| · · · ·· ·····················∘·········∘············ ·· ·· ·
|· · ·· ·················∘·∘∘∘··∘·················· · ··· ·
| · ·· ···············∘·······∘·∘∘∘·○·∘··∘∘∘················· ·
| · · · ··· ·············∘·∘∘∘∘∘··∘∘∘∘∘∘∘·∘∘··∘∘············· ·
7.28 | ·· · ·· · ··············∘∘·∘·∘∘∘∘∘∘○∘∘∘∘·∘∘∘∘∘·············· · · · ·
| ·· · ·· · ···············∘·∘·∘∘∘·∘∘○∘∘∘··∘∘∘················· ·
| · · ··························∘∘··∘∘··∘···················· ·
| · ··· ···························∘··············· ····· · ·· ·
| · · · ····· ········································ ·· · · ··
6.94 | · · ·········································· ··· ·· ·
| · · ···· ·· ············ ················ ·· ·· ··
| · ·· ··· ···· ······ · ·······
| · · · · · ·· · ··· ·
| · ·
6.60 | · · · ·
|--------|--------|--------|--------|--------|--------|--------|--------|--------
0.667 0.681 0.695 0.709 0.723 0.736 0.750 0.764 0.778
3.8 |
| ·
|
| ·
| ·
3.1 | · · ·· ·
| ·· · ·· · · · ·
| · · ··· · ········ ·· ·· ·
| · ·· ·· ························ ··
| ····· ····························· ··· ··
2.4 | · · · ·················∘···∘·∘············ ···
| · ·················∘·∘·∘∘∘∘∘··∘∘∘··∘············ ·
| · · · · ··············∘∘∘∘∘∘○∘·○∘∘∘∘∘∘∘∘∘∘··········
| ····· ·· ·······∘∘·∘∘∘∘∘∘∘○∘∘∘∘∘∘∘∘∘∘∘∘∘∘············ ·
| ···········∘··∘·∘∘∘∘∘○∘○○○∘○∘∘○∘∘∘∘∘∘············ ·
1.7 | · · ··············∘·∘∘∘∘○∘∘∘∘∘∘○∘∘∘∘∘○∘∘∘∘∘∘·········· ···
| · ·· ···············∘·∘∘∘·∘○∘∘∘∘∘∘∘∘∘··∘············ ·
| · ····················∘∘·∘·∘··∘∘········∘········ ·
|· · ·· ············································ ·· ·
| · · ·· ·· ····················· · · · · · ·
1.1 | ··· · · ·· ·· · · · ·
|--------|--------|--------|--------|--------|--------|--------|--------|--------
0.524 0.569 0.615 0.661 0.706 0.752 0.797 0.843 0.889
5.0 |
| ·
|
| ··
| · · · · · ·· · · · ·
4.1 | · · · ·· · · · · ·
| · · · · · · ······· ·· ·········· ·· ····· ·· ·· ·
| · ·· · ·································· · ·· ·
| · · ·· ···· ······································· ··· ·· ·
| · · ·········· ···············∘·∘·∘··················
3.3 | · ·· · ······················∘····∘∘∘·∘·∘·∘················· ··
| · · · ·······················∘∘·∘∘∘∘∘∘·∘∘○∘·∘·∘·∘∘················· ·
| · · ····· ···················∘∘∘∘∘○∘∘○∘∘∘○∘∘∘∘∘∘∘∘···············
| ·· ·· ··············∘·∘·∘∘∘·∘∘∘○∘○∘○∘○○∘∘∘∘∘·∘···············
| · · ········ · ··········∘∘·∘·∘∘∘∘○∘∘∘∘○○∘∘∘∘∘∘·∘∘···············
2.5 | ·· ··· ·················∘∘·∘∘∘∘∘∘∘∘·∘·∘∘∘∘∘·∘∘∘·············
| ·· · · ····························∘··∘···∘∘···∘∘················ ·· ·
| ·· · ·· ················································ · ··
| · · · ······· ··· ·································· ·· ·
|· · · · · ···· ······ ········ · · ·· ··
1.7 | · · · · · · ·
|--------|--------|--------|--------|--------|--------|--------|--------|--------
0.649 0.682 0.715 0.749 0.782 0.815 0.848 0.882 0.915
4.1 |
| ·
| · ·
| · · ·
| · · · · · ·
3.4 | ·· ·· ·· ·· ·· ·· · ······ ·· · ·
| · · ··· ······· ········· ········ · · ·
| · · · · · ·· · ······················ ··········· ·
| · · · ··· ·················∘··················· ·
| ·· ·· ·····················∘·∘······∘··········· ···
2.7 | · · ···················∘··∘··∘∘∘∘∘·∘∘·∘∘∘∘∘∘∘············ · · ·
|· ·· ·· ·················∘·∘∘∘∘∘∘○○∘∘∘∘∘○∘∘∘·∘∘······· ······
| · · ···············∘·∘····∘∘··∘∘∘○∘∘∘○∘∘∘∘·∘○∘∘············· ·
| · ·· ·· ··············∘∘·∘∘∘∘∘∘∘∘○○∘○○∘○○∘∘∘∘∘···∘············
| ··· · · ···········∘······∘··∘∘∘∘∘·∘∘○∘∘∘∘∘∘∘∘∘·∘·········· · ··
1.9 | · · · · · ··············∘··∘∘·∘·∘∘∘∘∘∘∘∘∘∘∘∘∘··················
| ···· ······ ··················∘∘·∘∘··∘∘···················· · ··
| ·· · ··· ······· ······································ ·
| · · · ···· ············ ····················· · ·
|· · · ··· · ·· · ·· ·· · ·· ·
1.2 | ·· · · · ·
|--------|--------|--------|--------|--------|--------|--------|--------|--------
0.599 0.635 0.671 0.708 0.744 0.781 0.817 0.853 0.890
9.92 |
| · ·
| ·· · ·
| · · · · · · ·
| · · · · ····· ·· ·· · · ·
9.47 | · · · ·· ············ ······· ·· · · ·· · ·
| · ··· · · ······· ······················ ··· · · · ·
| · · · ········· ··································· · ···· · ·
| · · · ····· ··· ··················∘···∘······················· ·
| · ·················∘··∘∘∘··∘∘∘∘∘∘∘··∘·∘∘∘·············· · ·· ·
9.02 | · ··········· ·······∘∘∘∘∘∘∘∘∘○○○∘∘∘∘∘∘∘∘∘·∘·∘·∘·············· ···
| · ················∘··∘∘·∘∘∘∘○∘∘∘∘∘○∘∘∘○∘∘·∘∘∘·∘················ ·
|· ·· ·· · ···············∘···∘∘·∘∘○∘○○∘∘○∘○∘∘○∘∘∘∘∘∘···∘········· · ···
| · ··· ·· ···········∘·∘·∘·∘·∘∘·∘∘∘∘○∘∘∘∘∘∘∘·∘∘∘··············· ·· · ·
| · ·················∘····∘∘∘∘∘·∘··∘∘∘∘·∘∘∘∘···∘······· ··· · ·
8.57 | · ······ ·······················∘···∘·∘························ ··
| · · ···· ··· ·········································· · · ·
| ···· ···· ·· ····· ··························· · ·
| · · ···· ····· ······ ·········· · ··
| · ·· · · · ··· · · · · ·· ·· · · ·
8.12 | · · · · · ··
|--------|--------|--------|--------|--------|--------|--------|--------|--------
0.593 0.611 0.629 0.647 0.665 0.683 0.701 0.719 0.737
val data = List(2,0,0,3,3,3,0,3,5,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,2,2,1,0,0,0,0,0,0,0,0,0,0,2,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,2,0,2,2,0,0,0,2,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,2,5,2,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,4,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,2,0,3,3,0,6,0,0,0,6,1,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,4,0,2,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,1,0,0,0,0,0,0,0,2,0,0,8,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,3,0,0,0,0,2,0,0,0,0,0,0,0,8,0,1,2,0,0,0,2,0,0,2,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,2,0,3,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,3,2,0,0,0,0,2,0,0,0,0,6,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,2,0,0,0,2,0,0,3,2,0,3,0,0,0,0,0,0,3,2,0,0,0,0,0,0,0,0,0,3,1,0,0,0,0,1,0,0,0,2,0,0,0,0,0,2,3,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,1,2,0,0,0,0,6,20,0,2,2,2,2,0,2,0,0,0,1,3,0,0,0,1,0,0,1,66,2,1,1,1,0,0,0,0,0,0,0,8,0,1,3,2,3,0,0,3,0,0,1,0,0,0,0,0,0,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,41,0,3,0,0,0,0,0,2,0,6,3,2,0,0,0,0,0,0,0,0,0,0,0,0,0,3,0,0,3,0,0,0,0,0,0,3,0,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,1,0,3,0,0,2,0,0,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,69,82,46,4,23,1,0,1,2,10,2,9,4,0,2,0,2,0,2,1,3,0,2,0,2,0,2,0,2,0,2,0,12,0,2,1,4,10,22,11,12,0,1,0,6,0,0,0,1,0,10,137,0,0,0,0,0,0,0,0,0,0,10,2,1,3,0,0,10,2,11,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,3,0,0,0,0,0,0,1,0,0,0,0,2,0,0,0,1,0,0,0,0,6,0,0,0,10,0,0,0,10,0,1,0,1,0,2,1,2,0,6,3,0,3,0,0,0,0,82,0,1,0,0,0,0,0,0,4,0,4,0,12,0,10,10,76,12,5,1,12,3,4,0,0,136,2,6,4,24,3,6,1,0,0)
val zipModel = for {
psi <- Beta(1,1).param
lambda <- Normal(5, 5).param
zip <- Poisson(lambda).zeroInflated(psi).fit(data)
} yield (psi, lambda)
plot2D(zipModel.sample())
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment