/* 方法 (1) 簡易版 PR = ...; DV = 0; IF PR >= 0.5 THEN DV = 1; 方法 (2) PR = ...; DV = 0; R = RANUNI(); IF R <= PR THEN DV = 1; */ %LET NSIM = 1000; *** 1 群 (1 用量) あたりのシミュレーション回数 ***; %LET SEED = 2601; *** 乱数の種 ***; DATA p; *** Group = 1 ~ 9, それぞれ母確率 10%, ..., 90% に対応させる ***; DO Group = 1 TO 9; DO i = 1 TO &NSIM; *** +/- 50% のばらつきを持たせる ***; pr = Group / 10 + RANUNI(&SEED) - 0.5; *** pr: 確率 probability ***; *** 0 <= pr <= 1 にする ***; IF pr < 0 THEN pr = 0; IF pr > 1 THEN pr = 1; *** 方法 (1) 簡易版 ***; DV_simple = 0; *** pr >= 0.5 だったら「有効」とする ***; IF pr >= 0.5 THEN DV_simple = 1; DV = 0; *** 確率 pr に応じて有効/無効を生成する ***; IF RANUNI(&SEED) <= pr THEN DV = 1; OUTPUT; END; END; RUN; PROC FREQ DATA=p; TITLE1 '方法 (1) 簡易版'; TABLES DV_simple * Group / NOPERCENT NOROW NOFREQ; RUN; PROC FREQ DATA=p; TITLE1 '方法 (2)'; TABLES DV * Group / NOPERCENT NOROW NOFREQ; RUN; TITLE1;
最近作成されたWikiのアクセスランキングです。見るだけでなく加筆してみよう!
atwikiでよく見られているWikiのランキングです。新しい情報を発見してみよう!
最近アクセスの多かったページランキングです。話題のページを見に行こう!
- 参加者一覧 - ストグラ まとめ @ウィキ
- 成瀬 力二 - ストグラ まとめ @ウィキ
- ギャルがアップした動画 - 検索してはいけない言葉 @ ウィキ
- マイティーストライクフリーダムガンダム - 機動戦士ガンダム EXTREME VS.2 INFINITEBOOST wiki
- 魔獣トゲイラ - バトルロイヤルR+α ファンフィクション(二次創作など)総合wiki
- 7月20日 安房鴨川駅イベントで流れた発メロをまとめたサイト - 発車メロディーwiki
- ブラック ジャックス - ストグラ まとめ @ウィキ
- 桃兎 すずめ - ストグラ まとめ @ウィキ
- ストライクフリーダムガンダム - 機動戦士ガンダム EXTREME VS.2 INFINITEBOOST wiki
- コメント/雑談・質問 - マージマンション@wiki