[過去ﾛｸﾞ] 高校数学の質問スレ Part438 (1002ﾚｽ)
上下前次1-新
抽出解除 必死ﾁｪｯｶｰ(本家) (べ) ﾚｽ栞 あぼーん
このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ
193: 09/17(火)08:27 ID:Xm+5bLR5(1/6) AAS
>>192
PDFをsin(x)/2 [0,π]に設定して実行。
画像ﾘﾝｸ[png]:i.imgur.com
ヒストグラムが乱数
実線はy=sin(x)/2
Rのコードは兄弟スレにおいた。
196: 09/17(火)13:13 ID:Xm+5bLR5(2/6) AAS
>>195
N=2e6程度にすれば100万個以上できる。
197: 09/17(火)14:00 ID:Xm+5bLR5(3/6) AAS
乱数発生させずに100万個作るなら
sapply(seq(0,1,length=1e6),\(x) 2*asin(sqrt(x)))
で終了。
200: 09/17(火)19:20 ID:Xm+5bLR5(4/6) AAS
定義域を 0 < x < 1 確率密度関数をlog(1/x) とする。
(1) 95%信頼区間を算出せよ。
(2) この分布に従う乱数をちょうど10万個つくれ。
Rでの算出例
画像ﾘﾝｸ[png]:i.imgur.com
画像ﾘﾝｸ[png]:i.imgur.com
> hist(y,freq=FALSE,breaks = 'scott' ,col='lightblue',axes=FALSE, ann=F) ; axis(1)
> curve(log(1/x),lwd=2,add=TRUE)
> length(y)
[1] 100000
省1
201: 09/17(火)19:28 ID:Xm+5bLR5(5/6) AAS
>>198
確率密度関数 sin(x)/2 (0<x<π) で、ちょうど100万個にすると算出に時間がかかるので10万個で実行。
画像ﾘﾝｸ[png]:i.imgur.com
画像ﾘﾝｸ[png]:i.imgur.com
> y=invcdf(x)
> hist(y,freq=FALSE,breaks = 'scott' ,col='lightblue',axes=FALSE, ann=F) ; axis(1)
> curve(sin(x)/2,lwd=2,add=TRUE)
> length(y)
[1] 100000
202: 09/17(火)19:29 ID:Xm+5bLR5(6/6) AAS
尿瓶チンパンフェチのPhimoseくんによれば、確率密度関数は山の形らしいね。
東大非合格者であると推定される。
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.026s