購入回数の確率を求めるNBDモデル(負の二項分布)の確率分布を与えるパラメータKの計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので)
参考文献によると、浸透率PenとMから
Pen = 1 - (1+M/K)^(-K)
でKが算出できるとのこと。ただし、非線形の方程式なのでここはニュートン・ラフソン法で計算しよう。
1-Pen=qとして
f(K) = (1+M/K)^(-K) - q
として、
K(n+1) = K(n) - f(K(n)) / f'(K(n))
として反復計算すればいい。ただf'を求めるのがめんどくさい。今qは定数で微分すると消えるので置いておいて、
y=(1+a/x)^(-x)を微分することを考えよう。両辺のログをとると
log(y)=-x*log(1+a/x)
両辺微分すると、
y'/y = -log(1+a/x) - x /(1+a/x) * (-a/x^2) = -log(1+a/x) + a/(x+a)
なので
y' = ((1+a/x)^(-x))*( -log(1+a/x) + a/(x+a))
となる。計算に全く自信がないのでWolframAlphaでも確かめた。
初期値はどうする?ということについてはf(K)のグラフの形を見ると
単調減少してる。ということであんまり何も考えず0に近い数(0.01とか)を入れておけばよさそう。
この計算をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUPしてみた。
リンクはこちら:
画面はこんな感じ:
参考文献:「確率思考の戦略論 USJでも実証された数学マーケティングの力」(森岡毅、今西聖貴)
« 購入回数の確率を求めるNBDモデル(負の二項分布)の計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので) | トップページ | 「名探偵の証明 蜜柑花子の栄光」(市川哲也さん)を読んだ。密室殺人ゲーム(歌野さん)+探偵がどんどん追い詰められていくようなお話で今回もラストで驚く、、、 »
「パソコン・インターネット」カテゴリの記事
- TensorFlow(Keras)でモデルをsaveで保存してload_modelで読み込むときに突然エラー(utf8で読めないとかディレクトリが存在しないとか)が出始めた。なんで?といろいろやると、単に日本語がフォルダ名に使われているときにエラーになるだけだった…(Windowsネイティブの場合)(2023.05.28)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(4) 10進31桁まで計算できることを生かして4段4次のルンゲクッタ法でローレンツ方程式を計算して図示してみる。テキストベース(アスキーアート)で!(2023.05.23)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで!今回は2次元配列の練習。(2023.05.19)
- カシオの高精度計算サイトkeisan.casio.jpが復旧してログインや自作式の作成ができるようになった。ということで作った自作式のアクセスベスト10とお勧めを紹介。(2023.05.17)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(2)話題のMojoが例題としてマンデルブロ集合を使っているそうなので、COBOLでマンデルブロ集合を描いてみよう。ただしテキストベースで! 多重ループの練習。(2023.05.10)
「学問・資格」カテゴリの記事
- TensorFlow(Keras)でモデルをsaveで保存してload_modelで読み込むときに突然エラー(utf8で読めないとかディレクトリが存在しないとか)が出始めた。なんで?といろいろやると、単に日本語がフォルダ名に使われているときにエラーになるだけだった…(Windowsネイティブの場合)(2023.05.28)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(4) 10進31桁まで計算できることを生かして4段4次のルンゲクッタ法でローレンツ方程式を計算して図示してみる。テキストベース(アスキーアート)で!(2023.05.23)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで!今回は2次元配列の練習。(2023.05.19)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース(5/18) IEEE Microwave Magazineはグラフェンのマイクロ波応用、Microwave Journalはテラヘルツ波の発生、Next G allianceの新6Gレポート、Hexa-Xの新動画、ESAが軌道上に6Gラボ、De-embeddingって何?など。(2023.05.18)
- カシオの高精度計算サイトkeisan.casio.jpが復旧してログインや自作式の作成ができるようになった。ということで作った自作式のアクセスベスト10とお勧めを紹介。(2023.05.17)
「経済・政治・国際」カテゴリの記事
- 昨日のNHKスペシャル「“中流危機”を越えて 「第1回 企業依存を抜け出せるか」を見ていて世帯所得の中央値が1994年505万円、2019年374万円という話見たので、統計しらべて1985~2020年の所得分布と中央値をGIFアニメにしてみた。(2022.09.19)
- ガンマ・ポアソン・リーセンシー・モデルに出てくるパラメータm,kの算出をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので(2021.12.14)
- ガンマ・ポアソン・リーセンシー・モデルの計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので(2021.12.11)
- デリシュレーNBDモデル(ディリクレ):カテゴリー内の自社ブランド購入確率のパラメータSの計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので(2021.12.06)
- デリシュレーNBDモデル:カテゴリー内の自社ブランド購入確率(ディリクレNBDモデル)の計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので (2021.12.03)
« 購入回数の確率を求めるNBDモデル(負の二項分布)の計算をカシオの高精度計算サイトkeisan.casio.jpにUP!(林先生の森岡毅さんインタビュー見て数学マーケティングというのに興味持って確率思考の戦略論を買って読んだので) | トップページ | 「名探偵の証明 蜜柑花子の栄光」(市川哲也さん)を読んだ。密室殺人ゲーム(歌野さん)+探偵がどんどん追い詰められていくようなお話で今回もラストで驚く、、、 »
コメント