# #——初始化 ------------------------------------------------------------------ 图书馆(profileScoreDist) # #——印度卢比 ------------------------------------------------------------------- 数据(INR) INR # #——规范 ------------------------------------------------------------ 印度卢比。reg <- regularizeMatrix(INR) INR。reg # #——dist参数 ----------------------------------------------------------- 粒度gcgran < < - 0.05 - 0.01 gcmin gcmax < < - 0.01 - 0.99 # # gc分数考虑gcs < - seq (gcgran * (gcmin / gcgran), gcgran * (gcmax / gcgran), gcgran) # #——dist ------------------------------------------------------------------ ## 计算概率分布distlist < -拉普(gcs、功能(x) computeScoreDist(印度卢比。reg, x, granularity)) ## ----plot, fig.cap=“Rahmann et al.文章中图1的再现。”——distlist [[50]] plotDist (distlist [[50 ]]) ## ---- 短裤 --------------------------------------------------------------- ab5 < - scoreDistCutoffs (distlist[[50]], 500年,1,c = 1, 0.05) # # 5%罗斯福ab5 cutoffa # #美元5% FNR ab5 $ cutoffb # #罗斯福= FNR ab5 cutoffopt美元