【自编脚本】R语言告诉你：离三黄蛋还有多远？

动机
我也不知道我的脑子大半夜在抽什么疯，总之就是突然想想看看自己离一帖生三魂还有多远，就简单写了一小段R。
发完帖回来编辑：这帖一个都没生（悲）

使用方法
修改m和n参数，其中m代表含【2%发帖+1灵魂】词条勋章的数量，同理n则是【1%发帖+1灵魂】。
改完参数之后，运行代码即可输出（个人感觉四黄蛋概率太过渺茫、没有统计学意义，所以只到了三黄蛋，如果有需要，可以依葫芦画瓢自行添加一行）。

一帖至少3魂的概率为: 0.000752377 平均主题数: 1329.121
一帖至少2魂的概率为: 0.01471589 平均主题数: 67.95375
一帖至少1魂的概率为: 0.1745055 平均主题数: 5.730479

寄，1329帖感觉要发到下辈子（这是可以说的吗

代码

1. # 至少成功k次
   
2. KSouls <- function(m, n, k) {
   
3.   # 成功率为2%的实验次数和1%的实验次数
   
4.   prob_m <- 0.02
   
5.   prob_n <- 0.01
   
6.   # 计算至少成功k次的概率
   
7.   total_prob <- 0
   
8.   for (i in 0:m) {
   
9.     for (j in 0:n) {
   
10.       if ((i + j) >= k) {
    
11.         # 计算m次成功i次的概率和n次成功j次的概率的乘积
    
12.         prob_i <- dbinom(i, m, prob_m)
    
13.         prob_j <- dbinom(j, n, prob_n)
    
14.         total_prob <- total_prob + (prob_i * prob_j)
    
15.       }
    
16.     }
    
17.   }
    
18.   return(total_prob)
    
19. }
    
20. 
    
21. # 参数m和n分别是2%魂和1%魂勋章的数量 执行前在此处修改
    
22. m <- 8
    
23. n <- 3
    
24. cat("一帖至少3魂的概率为:",KSouls(m, n, 3),"平均主题数:",1/KSouls(m, n, 3),"\n")
    
25. cat("一帖至少2魂的概率为:",KSouls(m, n, 2),"平均主题数:",1/KSouls(m, n, 2),"\n")
    
26. cat("一帖至少1魂的概率为:",KSouls(m, n, 1),"平均主题数:",1/KSouls(m, n, 1),"\n")

复制代码

R语言虽然比较冷门，但解决概率论问题的时候还是蛮舒服的。而且本身相对轻量，用什么包就装什么包。
我个人是用部署在本地的R和R Studio，但在线工具也很好找。随便在Google搜索，比如关键词“在线运行R”，就能搜出一大片。
我这里也是在搜的几个中，随便挑了一个界面看着好看点儿的点击直达，效果如下：





来自群组: 星象占卜

一贴至少3魂也太难了惹，这想都不敢想

感觉自己没这运气，第一点灵魂只能是血液献祭了

大膽一點!來算一下一帖10魂?

对于灵魂目前已经没有更多需求了，现在都是随缘

概率是概率实际是实际靴靴

感觉实际上可能会有更大的概率发生帖子数多于预期的情况

概率99%，实际超绝命中1%

这个从旅程触发的就可以窥见一二
叠0.1的旅程，实际20贴都不一定加1……

谢谢，到现在都靠着风雪之家赌命，一胎都没生过

我没有在收集灵魂向勋章，发帖能有+1灵魂就已经谢天谢地了

瘫倒ing反正俺都40帖了一个都还没触发咧，吐血ing

0-0 一下三胞胎嘛。。感觉这个发帖数只能说努力定会成功惹。。

虽然但是，概率也只是概率，运气差可能发完数量也不行，运气好可能提前有

我连一个灵魂都没触发过，三个想都不敢想

1贴3魂目前还是想都不敢想

概率这东西。欧皇是欧，非到极致也是欧

是我做梦都不敢想的事情

不是太喜欢发贴，还是回帖容易一些，而且我基本上没发帖徽章

用阳寿回帖感觉也很难触发啊