你好:
MS給了我們一個線索。從 R 提示符:
>圖書館(VGAM)
>rgengamma
和
>?rgengamma
揭示了 rgengamma 返回一個修改後的基本 R rgamma 函數
y = rgamma(n, k)
規模*y^(1/d)
在 JMP 中,我們還可以使用基本“隨機 Gamma”函數的修改,例如
當前數據表()<<新列(“JMP GenGamma”,
公式(1*隨機伽瑪(1)^(1/2)));
當前數據表() << 運行公式;
(當然,必須擺弄參數才能獲得匹配)
匹配隨機數生成器,任何人——也許是通過固定種子?
此致,
-馬特
最好的
我認為您無法使用內置生成器獲得三參數廣義伽馬分佈的隨機數。 如果可能的話,我也想知道該怎麼做。
但是,JMP 9 可以從 R 中獲取隨機數。請參閱下面的 JSL 示例。
// 返回隨機數的函數
// 生成廣義伽馬分佈
//通過R包VGAM
rgengamma =功能( { n = 100 , b = 1 , d = 1 , ķ = 1 } ,
R初始化() ;
R 提交( "
庫(VGAM)
rggam <- rgengamma(100,1,2,3)
" ) ;
格甘= R 獲取(格甘) ;
R 術語() ;
格甘
) ;
// 如何使用函數填充 a 的示例
// n=100個隨機數的數據表
dt =新表( “廣義伽瑪” ,新列( “隨機的” ) ) ;
柱子( dt , “隨機的” ) <<設定值( rgengamma ( 100 , 1 , 2 , 3 ) ) ;