rburr： Random Number Generation for Burr Distribution (R)

You can generate it by inverting the PDF function of burr.

1) randomly generate F ~ uniform
2) Solve PDF function {  F=1- (1+x**c)**(-k) }  for X.


Here is sample code in SAS.


data t1;
c=1.5; k=2;
do i=1 to 10000;
F=ranuni(123);
X= ((1-F)**(-1/K)-1) ** (1/C);
output;
end;
run;
proc univariate data=t1;
var x;
histogram x;
run;