四种方法结果的比较
以下是关于整个模拟程序的解释与程序: 加载所需要的包,生成所需要的数据:目标分布服从威布尔分布,截断分布服从指数分布。 install . packages ( "rootSolve" ) install . packages ( "Rsolnp" ) library ( Rsolnp ) library ( rootSolve ) ################################################################# ##### define a new log function log_n ### log_n < - function ( z , n ) { an = n ^ 2 ind = ( z < 1 / an ) y1 = log ( 1 / an ) - 1.5 + 2 * an * z - ( an * z ) ^ 2 / 2 y1 [ ! ind ] = 0 z [ ind ] = 1 y2 = log ( z ) y = y1 + y2 return ( y ) } ################################################################################ #### Function used to calculate log ( gamma ( N