学习资源
R的官方文档Writing R Extensions是目前编写R包最完整也是最好的参考文档。
本文主要参考此官方文档。
推荐中文阅读文档 R语言忍者秘笈
R包的结构
R文件夹 (其中包含R包里的code),DESCRIPTION文件
(1)DESCRIPTION
(2)Licensing
发布到CRAN的包必须要有开源许可证书。
(3)Package Dependencies
这一部分需要给出我们的包需要依赖的包。包后面可以自行选择是否添加版本信息。比如‘MASS (>= 3.1-20)’。
这一部分同样也可以给出需要依赖的R版本,比如‘R (>= 2.14.0), R (>= r56550)’
实例
写一个R包,提供以下四个方程,分别叫做 pxkcd, qxkcd, dxkcd 和 rxkcd。类似于R中其他的p, q, d, 和 r 方程,比如(qnorm,pnorm,dnorm,rnorm)
我们需要实现的分布是在xkcd comic 2118中介绍的。
假设f是正态分布的概率密度函数,我们现在有一个二元随机向量 (X,Y)均匀分布在坐标横轴和f函数之间。X的边缘分布是正态分布,给定X的Y的条件分布是(0,f(x))之间的均匀分布。Y的边缘分布就是xkcd分布,也是我们需要实现R函数的对象。
让h(y)代表从X均值到f(x)=y任意一点的距离,
Y的分布函数是:
用Mathematica算出Y的概率密度函数是
计算Y的分布是很困难的,目前看来只有数值的方法。
解答将会在后续文章中给出。
君君提示:你也可以发布优质内容,点此查看详情 >>
本文著作权归作者本人和北美省钱快报共同所有,未经许可不得转载。长文章仅代表作者看法,如有更多内容分享或是对文中观点有不同见解,省钱快报欢迎您的投稿。