高校数学の質問スレ Part410
■ このスレッドは過去ログ倉庫に格納されています
【質問者必読!!】
まず>>1-4をよく読んでね
数学@5ch掲示板用 掲示板での数学記号の書き方例と一般的な記号の使用例
http://mathmathmath.dotera.net/
・まずは教科書、参考書、web検索などで調べるようにしましょう。(特に基本的な公式など)
・問題の写し間違いには気をつけましょう。
・長い分母分子を含む分数はきちんと括弧でくくりましょう。
(× x+1/x+2 ; ○((x+1)/(x+2)) )
・丸文字、顔文字、その他は環境やブラウザによりうまく表示できない場合があります。
どうしても画像を貼る場合はPCから直接見られるところに見やすい画像を貼ってください。
ピクトはPCから見られないことがあるので避けてください。
・質問者は名前を騙られたくない場合、トリップを付けましょう。
(トリップの付け方は 名前(N)に 俺!#oretrip ←適当なトリ)
・質問者は回答者がわかるように問題を書くようにしましょう。
でないと放置されることがあります。
(変に省略するより全文書いた方がいい、また説明なく習慣的でない記号を使わないように)
・質問者は何が分からないのか、どこまで考えたのかを明記しましょう。
それがない場合、放置されることがあります。
(特に、自分でやってみたのに合わないので教えてほしい、みたいなときは必ず書くように)
・回答者も節度ある回答を心がけてください。
・970くらいになったら次スレを立ててください。
※前スレ
高校数学の質問スレ Part409
http://rio2016.5ch.net/test/read.cgi/math/1608682829/ >>706
n=242, p=0.04で二項分布B(n,p)がピークになるのは成功回数が9の時
8,9,10で確率を出してみると
> cbind(8:10,dbinom(8:10,242,0.04))
[,1] [,2]
[1,] 8 0.1208300
[2,] 9 0.1308991
[3,] 10 0.1270812
なので 正規分布N(9,σ)の確率密度が0.1308991になるようなσを求めればピークが一致するはず。
ピークの一致以外は何も考えない近似ではある。
やってみた。
> n=242
> p=0.04
> q=1-p
> nn=0:n
> (mu=nn[which.max(dbinom(nn,n,p))])
[1] 9
> (sd=uniroot(function(sd) dnorm(mu,mu,sd)-dbinom(mu,n,p),c(1,5))$root)
[1] 3.047708
検算
> dnorm(mu,mu,sd)
[1] 0.1308991
> dbinom(mu,n,p)
[1] 0.1308991
このスレにちなんで、罵倒厨の近似とでも呼ぶかなw >>710
ほら、やっぱり100が十分大きいか答えることができない。 全然、近似できてないんだな。
さいころを100回投げたとき、3の倍数の目が出る回数Xをとする。
40 <= Xとなる確率の近似値を求めよ。
https://www.ozl.jp/unit/statistics/2585.html こんな問題ができる。
二項分布B(n=242,p=0.04)をポワソン分布で近似させるとき、ピークが完全に一致するようなポワソン分布のパラメータを求めよ。 >>714
意外にも答が二個あった。
> (lambda1=uniroot(function(lambda) dpois(mu,lambda) - dbinom(mu,n,p) ,c(0,9))$root)
[1] 8.661711
> (lambda2=uniroot(function(lambda) dpois(mu,lambda) - dbinom(mu,n,p) ,c(9,10))$root)
[1] 9.34699
> dpois(mu,lambda1)
[1] 0.1308992
> dpois(mu,lambda2)
[1] 0.1308991
> dbinom(mu,n,p)
[1] 0.1308991 期待値も知らない
理由も知らない
そんなアホが今日も荒らしているのかよ >>713
Wolframに入力して二項分布の確率を加算すると
Sum[Binomial[100, k] (1/3)^k(2/3)^(100-k), {k, 40,100}]
0.0966230702となった。 >>705
期待値の定義に従ってプログラムに計算させているじゃん。
>708の1/244はβ分布の期待値の公式を使ったけど。
別の方法(JAGSでのMCMC)でそれを検算。 >633の答はまだかよ?
俺の数値解と照合したいんだが。 >>548
> 期待値の計算は
> Σ[n=0,242] n * 242Cn * p^n * (1-p)^(242-n)
>
> 手計算は大変なので
> 全部プログラム(R)が計算してくれる。
頭悪過ぎw nが十分大きいときある種の統計量の分布が正規分布に近づく、その誤差がどれくらいかの見積もりもできる
というのは統計学のイロハのイですわな 統計も期待値も分かってないのに朝から連投とはよほど恥を晒したいと見える >>722
イロハなら、これに即答できる?
n=242 p=0.04の二項分布は正規分布で近似してよいほどnは十分大きいか?
正規分布で近似したとき成功数が0の確率は(1-0.04)^242に近似していると判断してよいか? >>690
>普通に考えて、242回失敗したってだけで成功確率の推定なんかできないでしょ。
事前確率分布を想定すれば推定できる。信頼区間もだせる。
>主催者の言い分によって推定確率が変わるってのも変だし。
事前確率分布によって推定値が変わるのは別に変でもない。
主観的だといわれるが、日本人成人の平均身長を1〜2mの間に分布すると想定するのは俺には違和感はない。 正規分布って負の値も許すから、二項分布での成功回数が負の値をとる確率0でないのは本当はおかしい。 n=242 p=0.04の二項分布で成功回数が0の確率は
> 0.96^242
[1] 0.00005124345
正規分布近似で計算すると
> p=0.04
> n=242
> q=1-p
> # 1まで
> pnorm(1,n*p,sqrt(n*p*q))
[1] 0.00220399
> # 0から1まで
> pnorm(1,n*p,sqrt(n*p*q)) - pnorm(-0,n*p,sqrt(n*p*q))
[1] 0.001455908
全然、近似していない。
∴p=0.04のとき242は十分大きな数とは言えない。 >>727
【演習問題】 既述の罵倒厨の補正を用いると近似値が改善するか検討してみよ。 >>728
YESかNOで答えるだけの問題なのにw >>731
しないかどうかに難易度関係あるんwwww 医者板に書き込みないなと思ったらここで発狂してんのかプロおじ >>730
階層モデルの分散パラメータの事前分布には半コーシー分布(コーシー分布の正の方)を使う。
何故か? 分散は負の値をとらないから。
よくある試験の問題で〇〇の値は正規分布に従うという設定は−∞から+∞の値をとるか変数なのかを考えると当てはまらないものが多い。 >>736
プロおじ=トケジ=ウリュウ
これで確定だな。 >>727
pが1/2から離れると正規分布での近似は外れるので
λ=n*pを使ってポアソン分布で近似してみると
n=242 ; p=0.04
> dpois(0,n*p)
[1] 0.0000625215037748
> 0.96^242
[1] 0.0000512434540528
なので正規分布よりは近似した。
罵倒厨の補正をしてみる。
n=242 ; p=0.04
binom_pois <- function(x) sum((dbinom(0:n,n,p)-dpois(0:n,x))^2)
binom_pois=Vectorize(binom_pois)
b=optimise(binom_pois,c(0,30))$minimum
> dpois(0,b)
[1] 0.000061586237615
わずかながら近づいた。 nが十分大きければ正規分布で近似できるか?
レアアイテム排出確率が4%のガチャが242回連続して外れる確率の計算に
(1-0.04)^242を手書き計算するのは大変なので、正規分布近似で求めることにした。
二項分布の平均、標準偏差を用いてN(242*0.04,√(242*0.04*0.96))の正規分布で近似する。
求めたいのは成功0回の確率なので、この正規分布が-0.5から+0.5の確率を求めることにした。
(負の値があるのはおかしいという議論はここではしない)
n=242
p=0.04
q=1-p
pnorm(0.5,n*p,sqrt(n*p*q)) - pnorm(-0.5,n*p,sqrt(n*p*q))
> pnorm(0.5,n*p,sqrt(n*p*q)) - pnorm(-0.5,n*p,sqrt(n*p*q))
[1] 0.000880474614122
電卓で(1-0.04)^242を出すと0.00005124345405で1桁違っているのでとても近似とは言い難い。
nが十分大きな数であればいう人もいる。
n回連続して外れる確率を正規分布近似でだすにはnがどれほど大きければよいか?
有効数字1桁は合致していれば十分な近似が得られたと判断する。 >>740(実験)
nを1万にしたら十分大きいだろうか?
> p=0.04
> q=1-p
> n=10000
> q^n
[1] 5.15620761213e-178
> pnorm(0.5,n*p,sqrt(n*p*q)) - pnorm(-0.5,n*p,sqrt(n*p*q))
[1] 7.08172480364e-93
むしろ、近似が悪くなったぜぃ! >>740(脱字訂正)
nが十分大きな数であればいう人もいる。
↓
nが十分大きな数であればいいという人もいる。 暴れてる
暴れてる
期待値npを知らなかったアホが暴れてるwww 成功確率pの試行をn回繰り返し、全部失敗する確率 PB =( 1-p)^n
PN:正規分布近似で−0.5から0.5までに入る確率
確率は小さな数になるので比を対数として(同じ値に近づけば比が1になるのでその対数は0に近づく)
log(PN / PB )をp:[0,0.5] n:[10,1000]の範囲でグラフ化
https://i.imgur.com/o7rbX5X.png
nが大きくなるほど、正規分布近似で計算した全部失敗する確率が理論値( 1-p)^nから乖離するのが見て取れる。 高校数学の文字を読んで灯光やめちゃったらきっと本人の中で負けたことになっちゃうんだろうね そやね
勝ち負けが全て
今回の話も統計の教科書買ってきて読めばすぐ解決する話だけど彼はやらない
それは負けを認める事になるからな
この先もずっと我流を貫き通すつもりなんやろ 哀れだね
一生目を背けたまま妄執に囚われて生きるなんて
まあこんなところに平日の昼間から20レスもする時点で終わってるんだけどw わかすれで質問しとるwww
ベイズ統計学んだんじゃないんですかねぇ?wwwwwww >>711
N[m,σ^2] m = 9.21589, σ^2 = 9.24188
では以下のようです。
x f(x)
8 0.121142
9 0.130899
10 0.126936
9.21589 0.131229 (ピーク) >>666
AB上に, ∠AMN=15°となるように点Nをとる。MN=AN。
△MNBは∠N=30°=∠Bだから二等辺三角形。
よってMB=MN=MCになるので、三角形BNCは∠N=90°の直角三角形。特に30°定規形。
よって∠NCM=60°。よって△MCNは正三角形。よってCN=MN=AN。
よって△ANCは直角二等辺三角形。よって∠ACB=45°+60°=105°。 a,b,c,d,p,qは実数で、|ad-bc|=|pq|≠0をみたしている。
xy平面上において|ax+by|≦|p|かつ|cx+dy|≦|q|をみたす
点(x,y)全体からなる領域の面積を求めよ。
この問題なんですけど
領域は平行四辺形になるんでしょうか?
平行四辺形の頂点を求めて面積を出す方針で解けますか? >>757
Mが中点
> Calc(30,15,1/2)
[1] 105
Mが1/4内分点
> Calc(30,15,1/4)
[1] 39.89609
Mが1/3内分点
> Calc(30,15,1/3)
[1] 65.10391
交点の座標を求めて内積を使って角度を出すプログラムを書くだけ。 >>756
さっぱりイメージが沸かないので乱数発生させて1万個の点を抽出してどんな形になるかやってみた。
俺のプログラムだと台形になった。
https://i.imgur.com/JGVkMye.png
オマケ(作図のRのコード)
r=runif(5)
a=r[1];b=r[2]c=r[3];d=r[4];p=r[5]
q=(a*d-b*c)/p
c(a=a,b=b,c=c,d=d,p=p,q=q)
k=1e4
i=0
re=NULL
while(i<k){
i=i+1
x=runif(1) ; y=runif(1)
if(abs(a*x+b*y) <= abs(p) & abs(c*x+d*y) <= abs(q)){
re=rbind(re,c(x,y))
}
}
plot(re,col=2) >>760
発生させる乱数を正に限定しているというバグを発見
コードを修正したら確かに平行四辺形になった。
https://i.imgur.com/M6fQwhj.png
なので>760は撤回
オマケの修正
r=runif(5)
a=r[1];b=r[2];c=r[3];d=r[4];p=r[5]
q=(a*d-b*c)/p
c(a=a,b=b,c=c,d=d,p=p,q=q)
k=1e4
re=NULL
while(length(re)<k){
x=runif(1,-1,1) ; y=runif(1,-1,1)
if(abs(a*x+b*y) <= abs(p) & abs(c*x+d*y) <= abs(q)){
re=rbind(re,c(x,y))
}
} ド底辺の私立医でも台形じゃないってことはわかるだろうな 台形とか言っちゃうプロおじには答えて欲しくなかったんだけど… プロおじあれだけ恥を晒しておきながらまだ懲りてなかったのか?ww >>758
u = ax + by,
v = cx + dy,
を xy平面から uv平面への1次変換と考える。
この変換によって 面積は |J| 倍になる。
J = | a b | = ad - bc, (ヤコビアン)
| c d |
また uv平面では長方形
|u| ≦ |p|, |v|≦|q|
となるから、面積は 4|pq|
∴ 元のxy平面での面積は
4|pq|/|J| = 4|pq|/|ad-bc| = 4 なるほど
第二、第三の引数が範囲で省略すると[0,1]なんだな
直前に[0,1]で失敗してるのに >>768
よく意味が分からないのですが、
どこをみればUVが長方形って分かるのですか? >>771
ヤコビアンを使わないと解けないのですか?
ヤコビアンはまだ習ってないんですが モンテカルロ法で発生させる乱数の範囲が狭すぎると取りこぼしがでるし、広すぎるとヒット率が低下する。
a,b,c,d, p,q の値に応じて変化させてみたが、平行四辺形になるみたいだな。
https://i.imgur.com/vux0c9b.png >>773
> モンテカルロ法で発生させる乱数の範囲が狭すぎると取りこぼしがでるし、広すぎるとヒット率が低下する。
当たり前やん
> a,b,c,d, p,q の値に応じて変化させてみたが、平行四辺形になるみたいだな。
>
適正値評価する方が遥かに普通に解くよりめんどくさい
自己満しか得られない >>773
|ax+by|≦|p|かつ|cx+dy|≦|q|
の不等式で等号が成立するときのx,yの値を求めて、その範囲で乱数発生させることにした。
おまけのコードは長くなるのでココに置いた。
https://egg.5ch.net/test/read.cgi/hosp/1612996282/483
モンテカルロ法で面積の近似値も出せるように機能追加。
a, b, c, d, p を-1から1の範囲で乱数発生させてq=(a*d-b*c)/pとして
|ax+by|≦|p|かつ|cx+dy|≦|q|を満たす領域の例
https://i.imgur.com/Kr02E1Z.png >>774
どれか数値が間違っていたか?
間違っていたなら指摘してくれ。 >>775
不等式で等号が成立する値を範囲で乱数の境界に設定すればいいんじゃないの? >>779
どれか数値が間違っていたか?
生き恥といえば罵倒にしか喜びを見いだせない椰子のことだろ。 じゃあなんで反応した?
心当たりがあったから反応したんだろ?
アンカーもつけてないしなぁw よっぽど悔しかったみたいだねw
しばらく大人しくしてたくせにここにきて反応早すぎワロタ まぁ正規分布の話は酷すぎやからな
自分で“ベイズ統計を学んだ者”と言っておいて統計学で最も重要な定理の使い方知らんかったわけだしなぁ
他に何から始めるん?っていうくらいの1番大切な定理なのに 二項分布の正規分布近似ネタ
サイコロをn回投げて1の目の出る回数がn/10回以下である確率として
二項分布で求めた値をpb
正規分布で近似した値をpnとする。
nを増やしていいくとpb/pnは1に収束するか? 誤差を正規分布に設定するのは問題ないけど
高校生の身長を正規分布に設定するのは実はおかしい。
負の値の確率が0でないから。
これは、ある統計学の本に記載してあってなるほどと唸った。
この本:http://www.intuitivebiostatistics.com/ >>784
中心極限定理からすると>785は1に収束しそうと予測してグラフ化してみたら以外な結果だった。
pb - pn は予想とおり0に収束するグラフが得られた。
nを10から1000まで増やしたとき。
https://i.imgur.com/gNVBgWE.png
>785
二項分布での値が厳密値だからpb/pnより、分母に置いたpn/pbの方がいいな
改題
二項分布の正規分布近似ネタ
サイコロをn回投げて1の目の出る回数がn/10回以下である確率として
二項分布で求めた値をpb
正規分布で近似した値をpnとする。
nを増やしていいくとpn/pbは1に収束するか? この数学の問題を解いて欲しいです!! #知恵袋_ https://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q13240204423?fr=ios_other
なあその回答間違ってるつってんだろ
回答書き終わるまでにBAにならないようにわざわざ一旦短く回答してから返信する工夫までしたのにクソだな なんですぐわかることをプログラム使ったの?
本気でわかんなかったのかな...? >756を
|ax|+|by|≦|p|かつ|cx|+|dy|≦|q|
という条件に変更してモンテカルロで描画してみた。
菱形になるような印象。
https://i.imgur.com/PC4NkzY.png
証明は知らん。 プログラムキチガイってPC使わないと何一つ問題解けないよな
以前、補助線1本引けば解ける中学数学の問題を
わざわざPCで解いてたしなwww >>756
便宜上 p>0, q>0 とする。
↑OF = (a,b) ↑OG =(c,d) とおく。
|ax+by| ≦ p,
∴ (a,b) 方向の成分の絶対値が p/|OF| 以下。
∴ (x,y) は 幅が 2P = 2p/|OF| である帯の境界or内部にある。
|cx+dy| ≦ q,
∴ (c,d) 方向の成分の絶対値が q/|OG| 以下。
∴ (x,y) は 幅が 2Q = 2q/|OG| である帯の境界or内部にある。
2つの帯の交角θは ↑OF と ↑OG の交角だから
sinθ = 2儖FG / (|OF|・|OG|) = | ad - bc | / (|OF|・|OG|),
2つの帯の共通部分の面積は
(2P)(2Q)/sinθ = 4pq / | ad - bc | = 4. >>757
> 幾何学の王道:作図して計測
ダウト。それは図学と測量だ。
言うなれば図学版計算科学だ。
一方で幾何学は図学版純粋数学であり、
幾何学は純粋数学の範疇だ。
よって作図して計測しての解は
幾何学の王道などでは断じてないばかりか
幾何学でさえない。 白45個、赤55個の玉を無作為に1個ずつ取り出す。
どちらかの色が全て取り出されたら終了。
白が取り出されて終了した場合に取り出した玉の総数をnとする。
nの期待値と95%信頼区間を求めよ。 >>793
プログラムをする過程が楽しいんだね。
試験じゃないし。
立体図形の俯瞰図とか描けると面白い。
自分の勉強になる。
>797の期待値も定義に基づく値とシミュレーション解との近似を確認して検算。 期待値npを知らなかったバカが
また懲りずに書いているのかwww >755のような解より>757みたいに数値を変えても計算できる方が俺は好きだな。 PCに頼った数値解しか求められないアホ
期待値すら暗算出来ずにPCに頼るマヌケw >>798
お前が楽しいだけで周りは興味ないんだわ
スレタイ百回読んで他所でやれ 角度を求めよ
と言われて自分の頭では計算出来ずに
分度器使って求めるようなバカがこのスレにはいるよねwww 需要があるところでやればいいのになあ
ここでやるのは荒らしだわ >>787 のグラフは
p_n - p_b ?
小生は p_n を NORMSDIST((?-np)/√(np(1-p))) で求めた。(p=1/6)
上限を (n/10 + 1/2) にするとオーバーハングするけど
n/10 だと殆どしない… 要するに計算機に図描かせて眺めても理解できるほど数学は甘くないという事
一生理解できんやろ >>787
p_b 〜 (1.0705/√n) e^(-0.018182n)
p_n 〜 (0.7413/√n) e^(-0.0160n)
ゆえ n→∞ のとき
p_n/p_b は 1 には収束しないでしょう。
10log(10) - 9log(9) + log(p) + 9log(1-p) = - 0.018182
- (1/10 - p)^2 /{2p(1-p)} = - 0.0160
1/√(2πnp(1-p))} = 1.0705 >>807
中心極限定理からの直観だと1に収束すると思って、やってみたら意外な結果だったな。 >>799
知っていたけど、期待値の定義に従って算出しただけの話。
>797のような問題だと期待値の定義に従って計算することになると思うぞ。 >>805
p_n - p_bだったかp_b-p_nだったか忘れたが、確率の差のグラフです。 ■ このスレッドは過去ログ倉庫に格納されています