如何在R中利用已有的分布函数，构造新的分布函数

Python023

如何在R中利用已有的分布函数，构造新的分布函数,第1张

可以利用逆转法生成服从Pareto分布的随机数。原理：对于可逆的连续分布函数F,如果U服从[0,1]上的均匀分布，则随机变量F**(-1)(U)的分布函数为F。对于Pareto分布，生成对应随机数的代码如下： %let a=1/*the shape parameter (>0)of pareto dist*/ %let b=0.1/*the scale parameter (>0) of pareto dist*/ data pareto_distcount=1do until(count gt 100)uni_var=uniform(12345)pareto=&b/((1-uni_var)**(1/&a))if pareto ge 0 and pareto le 1 then docount+1outputendenddrop count uni_varrunproc print data=pareto_disttitle "Listing of Data Set Pareto_data"run

先把A、B的列宽设置好，用格式刷把其它列也刷一下，如果你的一页可以打印47行（是其它的数，就把公式里的47改一下）

在C1里写公式：

=INDIRECT(CHAR(65+MOD(COLUMN(B1),2))&(ROW(A1)+INT(COLUMN(B1)/2)*47))

按住C1的填充柄将公式右拖复制，再按住填充柄下拖复制。

选定打印区域，打印。

函数说明：

CHAR(65+MOD(COLUMN(B1),2)) 生成字符串 "A" 或 "B"（偶数列为"A"奇数列为"B"）。

(ROW(A1)+INT(COLUMN(B1)/2)*47 生成相应的行号。

INDIRECT() 生成字符串的引用。