ラプラス分布のサンプリング

誰でもスタータのラプラス分布から複数回描く方法を知っていますか？モンテカルロ分析を実行し、私のデータがラプラス分布に適合していることを知りたい。ラプラス分布のサンプリング

2017-04-07 TaxpayersMoney

'）私は、専用の機能を知らないが、私はあなたの質問を理解し、これは' generate'と（runiform 'への呼び出しの周りクォン機能をラップを使用して1行になります。 –

ありがとう、私の統計は弱いですが、これを拡大してもよろしいですか？ – TaxpayersMoney

完了。申し訳ありませんが、私は誇張しました。私はそれには2行が必要だと思います。 –

ここにサンプルスクリプトがあります。当然のことながら、尺度パラメータはそれが何であっても同じです。位置パラメータはここでは含意によってゼロです。ただでそれを追加しない場合。

clear 
version 10: set seed 2803 
set obs 10000 
scalar sigma = 1 
gen P = runiform() 
gen y = sigma * cond(P <= 0.5, log(2 * P), -log(2 * (1 - P)))

を私たちは、尾の動作が正常またはガウスとは全く異なっていることを示す参考として、通常の分位数プロットを使用することができます。

qnorm y

多くの人々は、密度推定のいくつかの種類

kdensity y, biweight bw(0.2)

を見ることを好むが、最も重要なグラフは、専用の分位数 - 分位数プロットです。これはqplotを使用します。これはStataのsearch qplotの後にStata Journalアーカイブからインストールする必要があります。 @は、ここではtypoではないことに注意してください。x軸にプロットされるもののプレースホルダーです。

qplot y, trscale(cond(@ <= 0.5, log(2 * @), -log(2 * (1 - @))))

出典

2017-04-08 08:17:14

パーフェクト、ありがとう。 – TaxpayersMoney

申し訳ありませんが、もう1つ質問：ラプラスの平均と標準デベロッパーを設定したい場合、どうすればいいですか？おそらく、シグマスカラーは関数のst devを設定しますか？私はコードをテストしたが、それはシグマデベロッパーでディストリビューションを生成するようには見えなかった。 – TaxpayersMoney

locationパラメータは平均値（および中央値とモード）です。すでに述べたように、生成された変数に追加するだけです。 SDは、任意の標準ソースに記載されているように、ルート2でシグマ乗算される。 https://en.wikipedia.org/wiki/Laplace_distribution –

ラプラス分布のサンプリング

答えて

関連する問題