非標準の密度関数からのデータをシミュレートしたいと考えています。私はすでに次のリンクを見つけました (確率関数を使用して任意の一変量確率変量をシミュレートするにはどうすればよいですか? )。ただし、これは奇妙な結果をもたらします。どういうわけか、この累積密度関数 ( cdf() ) はうまく機能しません。いくつかの値から、非常に奇妙な結果が得られます。たとえば、次のコードを見てください。
density=function(x)(25*200.7341^25/x^26*exp(-(200.7341/x)^25))
cdf<-function(x) integrate(density,1,x)[[1]]
cdf(9701)
[1] 1
cdf(9702)
[1] 6.33897e-05
私の質問は、「良い」CDF関数を作成するにはどうすればよいですか? または、より直接的に、PDF からデータをシミュレートするにはどうすればよいですか?