Amazonでアプリを個人開発しています

Amazonでアプリを個人開発しています

ホームピグアメブロ

芸能人ブログ人気ブログ

前ページ
次ページ

10Apr
- ExcelのVBAで指数分布の和がガンマ分布（アーラン分布）に従う事を乱数で示す
  本日は、ExcelのVBAで指数分布の和がガンマ分布（アーラン分布）に従う事を乱数で示したいと思います。ガンマ分布の形状母数が整数のとき、指数分布の和と一致し、アーラン分布と呼びます。ExcelのVBAプログラムは----------Sub Macro1()'' Macro1 Macro'' Keyboard Shortcut: Ctrl+r'simula = 10000 '乱数の個数alpha = 5 '形状母数beta = 2 '尺度母数For j = 1 To simula Gamma = 0 For i = 1 To alpha u = Rnd Gamma = Gamma - beta * Log(1 - u) '指数分布の和の分布を計算 Next i Cells(j, 1) = Gamma 'アーラン分布に従う乱数をＡ列へNext jEnd Sub------A列に10000個、ガンマ分布の特殊なケース（指数分布の和）に従う乱数が出力されます。コメントにも書いてありますが、simulaが乱数の個数、alphaが形状母数（指数分布の和の個数）でbetaは尺度母数になります。A列をExce1上で検証したいセルG2からG22まで選択しFREQUENCY(A:A,F2:F22)と入力しCtrl+Shift+Enterを押します。セルH2=GAMMA.DIST(E2,$D$1,$D$2,FALSE)*10000セルH22までコピーします。１から2の乱数は3個存在し理論値は0.634個という意味です。グラフにするととなります。一致してますね。参考にしてください。Excelでマクロを実行するには１．開発タブが無い『ExcelでVBAを使うには（開発タブを表示する）』本日は、ExcelでVBAを使うために開発タブを表示させる方法を解説したい。まず、開発タブが表示されているか確認します。表示されていなければ、「ファイル」…ameblo.jp２．マクロを有効にする『Excelでマクロを有効にするには』本日は、Excelでマクロを有効にする設定変更について解説します。「ファイル」を押します。「オプション」を押します。ポップアップ画面で「トラストセンター…ameblo.jp３．上記ソースコードの使い方『Excelでveq02020のVBAソースコードを使うには』本日は、Excelでveq02020のVBAソースコードを使い方を解説したい。「開発」タブに移動します。「マクロの記録」をクリックします。ポップアップ画…ameblo.jp西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
09Apr
- Excelで標準コーシー分布の算術平均は標準コーシー分布であることを示す
  本日は、Excelで標準コーシー分布の算術平均は標準コーシー分布であることを示したい。Excelのシートで下記のように式を組みましたセルA1=TAN(PI()*(RAND()-0.5))セルA1:E10000までコピーします。セルG1=AVERAGE(A1:E1)セルG10000までコピーします。セルL2からL22まで選択しFREQUENCY(G:G,K2:K22)と入力しCtrl+shift+Enterを押す。セルM3=1/PI()/(1+J3^2)*10000セルM22までコピーします。セルL2とM2の意味は、-10から-9まで算術平均が36個ありその理論値が34.88個という意味です。G列が標準コーシー分布に従っていることがわかると思います。。これをグラフにしました。標準コーシー分布の算術平均は標準コーシー分布に従っていることがわかると思います。参考にしてください。西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
08Apr
- ExcelのVBAで複素数を扱うプログラムを作成しました（セルから値を取得）
  本日は、ExcelのVBAで複素数の扱い方を解説したい。セルA1とセルB1から２つの複素数を取得し和を計算しC1に結果を出力するプログラムは下記のようになります。---------Sub Macro1()'' Macro1 Macro'' Keyboard Shortcut: Ctrl+r' ' IMSUM(A1, B1) の結果を C1 に出力 Cells(1, 3) = WorksheetFunction.ImSum(Cells(1, 1), Cells(1, 2))End Sub-----セルが複素数であってもCells関数で値が取得できるみたいです。参考にしてください。Excelでマクロを実行するには１．開発タブが無い『ExcelでVBAを使うには（開発タブを表示する）』本日は、ExcelでVBAを使うために開発タブを表示させる方法を解説したい。まず、開発タブが表示されているか確認します。表示されていなければ、「ファイル」…ameblo.jp２．マクロを有効にする『Excelでマクロを有効にするには』本日は、Excelでマクロを有効にする設定変更について解説します。「ファイル」を押します。「オプション」を押します。ポップアップ画面で「トラストセンター…ameblo.jp３．上記ソースコードの使い方『Excelでveq02020のVBAソースコードを使うには』本日は、Excelでveq02020のVBAソースコードを使い方を解説したい。「開発」タブに移動します。「マクロの記録」をクリックします。ポップアップ画…ameblo.jp西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
07Apr
- ExcelのVBAで複素数を扱うプログラムを作成しました（内部のみでの処理）
  本日は、ExcelのVBAで複素数の扱い方を解説したい。内部で複素数を定義し２つの値を足してA1に結果を出力するプログラムは下記のようになります。--------Sub Macro2()'' Macro2 Macro'' Keyboard Shortcut: Ctrl+r' ' 複素数を作成 x = WorksheetFunction.Complex(1, 1) y = WorksheetFunction.Complex(2, 2) ' IMSUM(x, y) を計算 result = WorksheetFunction.ImSum(x, y) ' 結果を A1 に出力 Cells(1, 1) = resultEnd Sub-----WorksheetFunctionを使えば簡単に扱えるみたいです。参考にしてください。Excelでマクロを実行するには１．開発タブが無い『ExcelでVBAを使うには（開発タブを表示する）』本日は、ExcelでVBAを使うために開発タブを表示させる方法を解説したい。まず、開発タブが表示されているか確認します。表示されていなければ、「ファイル」…ameblo.jp２．マクロを有効にする『Excelでマクロを有効にするには』本日は、Excelでマクロを有効にする設定変更について解説します。「ファイル」を押します。「オプション」を押します。ポップアップ画面で「トラストセンター…ameblo.jp３．上記ソースコードの使い方『Excelでveq02020のVBAソースコードを使うには』本日は、Excelでveq02020のVBAソースコードを使い方を解説したい。「開発」タブに移動します。「マクロの記録」をクリックします。ポップアップ画…ameblo.jp西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
03Apr
- ExcelのVBAで超幾何分布に従う乱数を生成するプログラムを作成しました
  本日は、超幾何分布に従う乱数を生成するプログラムをExcelのVBAで作成しました。超幾何分布は成功状態をもつ母集団から非復元抽出したときに成功状態がいくつあるかという確率を与える離散分布である。ExcelのVBAプログラムは------------Sub Macro1()'' Macro1 Macro'' Keyboard Shortcut: Ctrl+r'simula = 10000 '乱数の個数N = 100 '母集団の標本数K = 30 '母集団の成功数a = 50 '抜取の標本数For j = 1 To simula x = 0 b = K '母集団の成功数のカウンター c = N '母集団の標本数のカウンター For i = 1 To a If x < 30 Then If Rnd <= b / c Then x = x + 1 b = b - 1 End If End If c = c - 1 Next i Cells(j, 1) = x '乱数Next jEnd Sub--------------------------------------となります。コメントにも書いてありますがA列に幾何分布に従う乱数がSimula個生成されます。Nは母集団の標本数、K は母集団の成功数a は抜取の標本数、b は母集団の成功数のカウンターとc は母集団の標本数のカウンターになります。A列をExce1上で検証したいセルG2=COUNTIF(A:A,E2)セルG22までコピーします。セルH2=HYPGEOM.DIST(F2,$D$3,$D$2,$D$1,FALSE)*10000セルH22までコピーします。成功数が１0の乱数は147個存在し理論値は164.82個という意味です。グラフにするととなります。一致してますね。参考にしてください。Excelでマクロを実行するには１．開発タブが無い『ExcelでVBAを使うには（開発タブを表示する）』本日は、ExcelでVBAを使うために開発タブを表示させる方法を解説したい。まず、開発タブが表示されているか確認します。表示されていなければ、「ファイル」…ameblo.jp２．マクロを有効にする『Excelでマクロを有効にするには』本日は、Excelでマクロを有効にする設定変更について解説します。「ファイル」を押します。「オプション」を押します。ポップアップ画面で「トラストセンター…ameblo.jp３．上記ソースコードの使い方『Excelでveq02020のVBAソースコードを使うには』本日は、Excelでveq02020のVBAソースコードを使い方を解説したい。「開発」タブに移動します。「マクロの記録」をクリックします。ポップアップ画…ameblo.jp西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
02Apr
- 標準コーシー分布が自由度が１のt分布に従うことを乱数で示す
  本日は、標準コーシー分布と自由度が１のt分布の関係について解説したい。標準コーシー分布は期待値と分散を持たず自由度が１のt分布と一致します。Excelのシートに式を組みました。セルA1=TAN(PI()*(RAND()-0.5))セルA10000までコピーします。セルF2:F22まで選択しFREQUENCY(A1:A10000,E2:E22)と入力しCtrl+shift+Enterを押す。セルG2=1/PI()/(1+(D2^2))*10000セルG22までコピーします。セルH2=T.DIST(D2,1,FALSE)*10000セルH22までコピーします。A列が標準コーシー分布に従う乱数になります。これをグラフにしました。乱数、標準コーシー分布と自由度が１のt分布が一致していることがわかると思います。参考にしてください。西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
01Apr
- ExcelのVBAでリーマン・ゼータ関数の値を求めるカスタム関数
  ExcelのVBAでリーマン・ゼータ関数の値を求めるカスタム関数の作り方を解説したい。リーマンゼータ関数は、数学の分野のひとつである解析的整数論において素数分布の研究をはじめとした重要な研究対象である。ExcelのVBAプログラムは---------------Option ExplicitFunction myzeta(n As Integer) myzeta = cumfnc(n)End FunctionFunction cumfnc(n As Integer)Dim euler, tmp, tmp2 As DoubleDim i, m As Longm = 1000000tmp = 0For i = 1 To m tmp2 = (i ^ n) tmp2 = 1 / tmp2 tmp = tmp + tmp2Next icumfnc = tmpEnd Function--------------------------となります。古い仕様が残っているみたいでカスタム関数を使う時は変数の宣言をしないと動きません。Excelのセルに例えばmyzeta(2)と入力すれば実数2のゼータ関数の値が求められる。Rの関数とmyzeta関数を比較した。セルC2=myzeta(A2)セルC20までコピーします。実数2の時のゼータ関数の誤差は少し大きいが他の値の時は十分なレベルであると考えられる。参考にしてください。Excelでマクロを実行するには１．開発タブが無い『ExcelでVBAを使うには（開発タブを表示する）』本日は、ExcelでVBAを使うために開発タブを表示させる方法を解説したい。まず、開発タブが表示されているか確認します。表示されていなければ、「ファイル」…ameblo.jp２．マクロを有効にする『Excelでマクロを有効にするには』本日は、Excelでマクロを有効にする設定変更について解説します。「ファイル」を押します。「オプション」を押します。ポップアップ画面で「トラストセンター…ameblo.jp３．上記ソースコードの使い方『Excelでveq02020のVBAソースコードを使うには』本日は、Excelでveq02020のVBAソースコードを使い方を解説したい。「開発」タブに移動します。「マクロの記録」をクリックします。ポップアップ画…ameblo.jp西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
31Mar
- Excelで複素数の指数を計算する方法
  本日は、Excelで複素数の指数を計算する方法を解説したい。Excelのシートに式を組みました。セルA2に任意の複素数の実部を入力します。セルB2に任意の複素数の虚部を入力します。セルB4=COMPLEX(A2,B2)とし複素数を構成します。実数の関数のみで複素数の指数を計算するための関数を構成しようと考えました。セルB7=COMPLEX(EXP(A2)*COS(B2),EXP(A2)*SIN(B2))となります。Excelでは複素数の指数を計算する関数が実装されていてセルB8=IMEXP(B4)となります。参考にしてください。西洋の敗北　日本と世界に何が起きるのか (文春e-book)Amazon（アマゾン）ウクライナの問題に絡めて欧米の衰退を描く。アサヒ生ビールアサヒビール350ml24本【飲食店で愛され続けたまろやかな味わいマルエフ】Amazon（アマゾン）苦みが余りなく若い人向けのビールかなアサヒ生ビール【芳ばしい香りとまろやかなうまみ】アサヒ黒生 350ml×24本Amazon（アマゾン）マルエフの黒版アサヒ　スーパードライドライクリスタルアサヒビール生ビール 350ml×24本Amazon（アマゾン）これも若い人向けに作ってるかなよろしくお願いいたします。以上です！
27Mar
- ExcelのVBAで幾何分布に従う乱数を生成するプログラムを作成しました
  本日は、幾何分布に従う乱数を生成するプログラムをExcelのVBAで作成しました。幾何分布は負の二項分布の特殊なケースで1回成功するために何回失敗するかの分布です。ExcelのVBAプログラムは----------Sub Macro1()'' Macro1 Macro'' Keyboard Shortcut: Ctrl+r'Simula = 10000 '乱数の個数N = 1 '成功回数p = 0.3 '成功率For i = 1 To Simula k = 0 j = 0 Do While k < N If Rnd < p Then k = k + 1 Else j = j + 1 End If Loop Cells(i, 1) = j 'A列に幾何分布に従う乱数を出力Next iEnd Sub----------となります。コメントにも書いてありますがA列に幾何分布に従う乱数がSimula個生成されます。Nは成功回数、pは成功率です。A列をExce1上で検証したいセルF2=COUNTIF(A:A,E2)セルF22までコピーします。セルG2=NEGBINOM.DIST(E2,$D$3,$D$2,FALSE)*10000セルG22までコピーします。失敗数が１の乱数は2097個存在し理論値は2100個という意味です。グラフにするととなります。一致してますね。参考にしてください。Excelでマクロを実行するには１．開発タブが無い『ExcelでVBAを使うには（開発タブを表示する）』本日は、ExcelでVBAを使うために開発タブを表示させる方法を解説したい。まず、開発タブが表示されているか確認します。表示されていなければ、「ファイル」…ameblo.jp２．マクロを有効にする『Excelでマクロを有効にするには』本日は、Excelでマクロを有効にする設定変更について解説します。「ファイル」を押します。「オプション」を押します。ポップアップ画面で「トラストセンター…ameblo.jp３．上記ソースコードの使い方『Excelでveq02020のVBAソースコードを使うには』本日は、Excelでveq02020のVBAソースコードを使い方を解説したい。「開発」タブに移動します。「マクロの記録」をクリックします。ポップアップ画…ameblo.jp［図解］大学4年間の経済学が10時間でざっと学べるAmazon（アマゾン）大学4年間の経済学が10時間でざっと学べる (角川文庫)Amazon（アマゾン）veq02020はレベルが上がった。バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ・エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017年正規品辛口赤ワイン 750ml 【バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ】エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017 ★希少なバックヴィンテージAmazon（アマゾン）ムートンのチリでコスパ最高。シャンパーニュバロンドロスチャイルドコンコルディアブリュット箱なし 750ml シャンパーニュフランス白泡コク辛口 ((VABRB6Z0)) ヴェリタスAmazon（アマゾン）高くなったなあ。よろしくお願いいたします。以上です！
26Mar
- Excelで複素数の自然対数を計算する方法
  本日は、Excelで複素数の自然対数を計算する方法を解説したい。Excelのシートに式を組みました。セルA2に任意の複素数の実部を入力します。セルB2に任意の複素数の虚部を入力します。セルB4=COMPLEX(A2,B2)とし複素数を構成します。実数の関数のみで複素数の自然対数を計算するための関数を構成しようと考えました。セルB7=COMPLEX(LN(SQRT(A2^2+B2^2)),IF(SIGN(A2)=1, ASIN(B2/SQRT(A2^2+B2^2)),IF(SIGN(B2)=1, ACOS(A2/SQRT(A2^2+B2^2)),-ACOS(A2/SQRT(A2^2+B2^2)))))ASINとACOSの計算方法が使いづらいので複雑になりました。Excelに実装された関数で計算するにはセルB8=COMPLEX(LN(IMABS(B4)),IMARGUMENT(B4))となりシンプルになりました。Excelでは複素数の自然対数を計算する関数が実装されていてセルB9=IMLN(B4)となります。参考にしてください。［図解］大学4年間の経済学が10時間でざっと学べるAmazon（アマゾン）大学4年間の経済学が10時間でざっと学べる (角川文庫)Amazon（アマゾン）veq02020はレベルが上がった。バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ・エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017年正規品辛口赤ワイン 750ml 【バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ】エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017 ★希少なバックヴィンテージAmazon（アマゾン）ムートンのチリでコスパ最高。シャンパーニュバロンドロスチャイルドコンコルディアブリュット箱なし 750ml シャンパーニュフランス白泡コク辛口 ((VABRB6Z0)) ヴェリタスAmazon（アマゾン）高くなったなあ。よろしくお願いいたします。以上です！
25Mar
- Excelで複素数の3乗根を計算する方法
  本日は、Excelで複素数の3乗根を計算する方法を解説したいと思います。Excelのシートに下記のように式を組みました。埋め込まれた式はそれぞれセルB2=IMREAL(B1)セルC2=IMAGINARY(B1)セルD2=COMPLEX(B2,C2)としD17までコピーします。セルE2=IMPOWER(D2,3)としE17までコピーします。セルF2=IMSUB(E2,$B$1)としF17までコピーします。セルG2=COMPLEX(3*(B2^2),6*B2*C2)としG17までコピーします。セルH2=IMDIV(F2,G2)としH17までコピーします。セルI2=IMSUB(D2,H2)としI17までコピーします。セルB3=IMREAL(I2)としB17までコピーします。セルC3=IMAGINARY(I2)としC17までコピーします。セルB19=COMPLEX(B9,C9)としこれが3乗根になります。他の2つの根はちょっと複雑ですがセルB20=IMPRODUCT(COMPLEX(IMABS(B19),0),COMPLEX(COS(IMARGUMENT(B19)+2*PI()/3),SIN(IMARGUMENT(B19)+2*PI()/3)))セルB21=IMPRODUCT(COMPLEX(IMABS(B19),0),COMPLEX(COS(IMARGUMENT(B19)+2*PI()*2/3),SIN(IMARGUMENT(B19)+2*PI()*2/3)))です。セルC19=IMPOWER(B11,3)セルC20=IMPOWER(B11,3)セルC21=IMPOWER(B11,3)とし3乗すると元の2+2iになっている事がわかると思います。参考にしてください。［図解］大学4年間の経済学が10時間でざっと学べるAmazon（アマゾン）大学4年間の経済学が10時間でざっと学べる (角川文庫)Amazon（アマゾン）veq02020はレベルが上がった。バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ・エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017年正規品辛口赤ワイン 750ml 【バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ】エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017 ★希少なバックヴィンテージAmazon（アマゾン）ムートンのチリでコスパ最高。シャンパーニュバロンドロスチャイルドコンコルディアブリュット箱なし 750ml シャンパーニュフランス白泡コク辛口 ((VABRB6Z0)) ヴェリタスAmazon（アマゾン）高くなったなあ。よろしくお願いいたします。以上です！
24Mar
- ブラック・ショールズ式で最高気温と最低気温を予測する
  余り上手くいかなかったがブラック・ショールズ式で最高気温と最低気温を予測した。Rでプログラムを組んだ---------------------library(forecast)library(derivmkts)## 1. データ読み込み（Book1.csv は 1 列・ヘッダなしを想定）dat <- read.csv("Book1.csv", header = TRUE)price <- as.numeric(dat[[1]]) # ベクトル化ts_price <- ts(price) # 時系列オブジェクト（必要なら frequency を指定）## 2. ARIMA モデルで 20 ステップ先を予測fit <- auto.arima(ts_price)fc <- forecast(fit, h = 20)S_pred <- as.numeric(fc$mean) # 予測された将来の「株価」## 3. ブラック–ショールズ用パラメータ設定## ここは用途に応じて書き換えてくださいK <- 11.3 # 行使価格（前日の気温）r_annual <- 0.00 # 年率無リスク金利（0）temp <- ts(price) # 気温データret <- diff(temp) # 日次変化sigma_annual <- sd(ret) * sqrt(365)#sigma_annual <- 0.00 # 年率ボラティリティ（0%）q_div <- 0.00 # 配当利回り（なければ 0）## 4. 予測ステップごとの「残存期間 T」を設定## 例：1〜20営業日先を想定し、1年=252営業日として T = h/252horizon <- 1:20T_vec <- horizon / 365## 5. 予測株価 × 残存期間ごとに bscall / bsput で価格計算call_prices <- mapply( function(S, T) bscall(s = S, k = K, v = sigma_annual, r = r_annual, tt = T, d = q_div), S = S_pred, T = T_vec)put_prices <- mapply( function(S, T) bsput(s = S, k = K, v = sigma_annual, r = r_annual, tt = T, d = q_div), S = S_pred, T = T_vec)## 6. 結果をまとめるresult <- data.frame( horizon = horizon, # 何ステップ先か（1〜20） T_years = T_vec, # 残存期間（年） S_pred = S_pred, # 予測株価 Call_BS = call_prices, # BS コール価格 Put_BS = put_prices # BS プット価格)print(result)## 必要なら CSV に出力write.csv(result, "BS_forecast_20steps_260324.csv", row.names = FALSE)---------------------------------------ブラック・ショールズのパラメータ設定を気温に置き換えて予測した。最高気温と最低気温の計算の方法はセルF2=11.3-E2セルG2=11.3+D2です。考え方としては最高気温が行使価格が11.3で、3/23日のコール価格が8.37。11.3+8.37でないと利益が出ない最低気温が行使価格が11.3で、3/23日のプット価格が7.98。11.3-7.98でないと利益が出ないこれを気温の問題に置き換えて予測しました。3/23の最低気温が9.8℃、最高気温が17.5℃で最低気温が上手く予測できない。予測出来ても２、３日か。参考にしてください。［図解］大学4年間の経済学が10時間でざっと学べるAmazon（アマゾン）大学4年間の経済学が10時間でざっと学べる (角川文庫)Amazon（アマゾン）veq02020はレベルが上がった。バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ・エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017年正規品辛口赤ワイン 750ml 【バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ】エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017 ★希少なバックヴィンテージAmazon（アマゾン）ムートンのチリでコスパ最高。シャンパーニュバロンドロスチャイルドコンコルディアブリュット箱なし 750ml シャンパーニュフランス白泡コク辛口 ((VABRB6Z0)) ヴェリタスAmazon（アマゾン）高くなったなあ。よろしくお願いいたします。以上です！
19Mar
- ExcelのVBAで負の二項分布に従う乱数を生成するプログラムを作成しました
  本日は、負の二項分布に従う乱数を生成するプログラムをExcelのVBAで作成しました。負の二項分布はk回成功するために何回失敗するかの分布です。ExcelのVBAプログラムは----------Sub Macro1()'' Macro1 Macro'' Keyboard Shortcut: Ctrl+r'Simula = 10000 '乱数の個数N = 5 '成功回数p = 0.3 '成功率For i = 1 To Simula k = 0 j = 0 Do While k < N If Rnd < p Then k = k + 1 Else j = j + 1 End If Loop Cells(i, 1) = j 'A列に負の二項分布に従う乱数を出力Next iEnd Sub----------となります。コメントにも書いてありますがA列に負の二項分布に従う乱数がSimula個生成されます。Nは成功回数、pは成功率です。A列をExce1上で検証したいセルF2=COUNTIF(A:A,E2)セルF22までコピーします。セルG2=COMBIN($D$3+E2-1,E2)*((1-$D$2)^E2)*(($D$2)^$D$3)*10000セルG22までコピーします。失敗数が１の乱数は84個存在し理論値は85.05個という意味です。グラフにするととなります。一致してますね。参考にしてください。Excelでマクロを実行するには１．開発タブが無い『ExcelでVBAを使うには（開発タブを表示する）』本日は、ExcelでVBAを使うために開発タブを表示させる方法を解説したい。まず、開発タブが表示されているか確認します。表示されていなければ、「ファイル」…ameblo.jp２．マクロを有効にする『Excelでマクロを有効にするには』本日は、Excelでマクロを有効にする設定変更について解説します。「ファイル」を押します。「オプション」を押します。ポップアップ画面で「トラストセンター…ameblo.jp３．上記ソースコードの使い方『Excelでveq02020のVBAソースコードを使うには』本日は、Excelでveq02020のVBAソースコードを使い方を解説したい。「開発」タブに移動します。「マクロの記録」をクリックします。ポップアップ画…ameblo.jp［図解］大学4年間の経済学が10時間でざっと学べるAmazon（アマゾン）大学4年間の経済学が10時間でざっと学べる (角川文庫)Amazon（アマゾン）veq02020はレベルが上がった。バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ・エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017年正規品辛口赤ワイン 750ml 【バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ】エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017 ★希少なバックヴィンテージAmazon（アマゾン）ムートンのチリでコスパ最高。シャンパーニュバロンドロスチャイルドコンコルディアブリュット箱なし 750ml シャンパーニュフランス白泡コク辛口 ((VABRB6Z0)) ヴェリタスAmazon（アマゾン）高くなったなあ。よろしくお願いいたします。以上です！
18Mar
- Excelで複素数の平方根を計算する方法
  本日は、Excelで複素数の平方根を計算する方法を解説したいと思います。Excelのシートに下記のように式を組みました。埋め込まれた式はそれぞれセルB2=IMREAL(B1)セルC2=IMAGINARY(B1)セルD2=COMPLEX(B2,C2)としD8までコピーします。セルE2=IMPOWER(D2,2)としE8までコピーします。セルF2=IMSUB(E2,$B$1)としF8までコピーします。セルG2=COMPLEX(2*B2,2*C2)としG8までコピーします。セルH2=IMDIV(F2,G2)としH8までコピーします。セルI2=IMSUB(D2,H2)としI8までコピーします。セルB3=IMREAL(I2)としB9までコピーします。セルC3=IMAGINARY(I2)としC9までコピーします。セルB11=COMPLEX(B9,C9)としこれが平方根になります。セルB12=IMPOWER(B11,2)とし元の2+2iになっている事がわかると思います。参考にしてください。［図解］大学4年間の経済学が10時間でざっと学べるAmazon（アマゾン）大学4年間の経済学が10時間でざっと学べる (角川文庫)Amazon（アマゾン）veq02020はレベルが上がった。バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ・エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017年正規品辛口赤ワイン 750ml 【バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ】エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017 ★希少なバックヴィンテージAmazon（アマゾン）ムートンのチリでコスパ最高。シャンパーニュバロンドロスチャイルドコンコルディアブリュット箱なし 750ml シャンパーニュフランス白泡コク辛口 ((VABRB6Z0)) ヴェリタスAmazon（アマゾン）高くなったなあ。
17Mar
- ExcelでLN(-1)を計算する方法
  本日は、ExcelでLN(-1)を計算する方法を解説したい。Excelで普通にLN(-1)と計算すると#NUM!とエラーを返します。数を複素数まで拡張するとLN(-1)は計算できます。IMLN(COMPLEX(-1,0))=3.14159265358979iとなります。LN(COMPLEX(-1,0))としないように注意。参考にしてください。［図解］大学4年間の経済学が10時間でざっと学べるAmazon（アマゾン）大学4年間の経済学が10時間でざっと学べる (角川文庫)Amazon（アマゾン）veq02020はレベルが上がった。バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ・エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017年正規品辛口赤ワイン 750ml 【バロン・フィリップ・ド・ロスチャイルド・マイポ・チリ】エスクード・ロホ・ヴィンテージ・コレクション [ボックス付] 2017 ★希少なバックヴィンテージAmazon（アマゾン）ムートンのチリでコスパ最高。シャンパーニュバロンドロスチャイルドコンコルディアブリュット箱なし 750ml シャンパーニュフランス白泡コク辛口 ((VABRB6Z0)) ヴェリタスAmazon（アマゾン）高くなったなあ。
13Mar
- Excelにおける混合二項分布の確率質量関数と乱数の生成方法
  本日は、Excelにおける混合二項分布の確率質量関数と乱数の生成方法について解説したい。混合率をa、成功率をp1、p2、成功回数をk、試行回数をnとすると混合二項分布の確率質量関数はとなります。(　)は二項係数です。次に混合二項分布に従う乱数を生成し上記確率質量関数をチェックしたい。二項分布に従う乱数の生成方法については『ExcelのVBAで二項分布に従う乱数を生成するプログラムを作成しました』本日は、二項分布に従う乱数を生成するプログラムをExcelのVBAで作成しました。ExcelのVBAプログラムは----------Sub Macro1()…ameblo.jpを参考にしてください。VBAプログラムのNとpを変更すれば、A列に二項分布に従う乱数が生成されます。二項分布に従う乱数を混合二項分布に従う乱数に変換します。上記VBAプログラムを利用してp=0.3,N=20として実行しA列にコピーしp=0.7,N=20として実行しB列にコピーします。セルC2=RAND()とし、C10001までコピーします。混合率を0.3とし、セルD1に入力します。セルD2=IF(C2>1-$D$1,A2,B2)としD10001までコピーします。D2からD10001までに、混合二項分布に従う乱数が生成されます。乱数の個数と理論値を計算したい。セルG4=COUNTIF($D$2:$D$10001,F4)としてG24までコピーします。セルH4=($D$1*BINOM.DIST(F4,$F$3,$F$2,FALSE)+(1-$D$1)*BINOM.DIST(F4,$G$3,$G$2,FALSE))*10000としてH24までコピーします。これをグラフにまとめると下記のようになります。上記乱数の生成法が理論値とほぼ一致した。注）p=0.3,N=20として実行しA列にコピーしp=0.7,N=30として実行しB列にコピーすると面白いことが起きる。高校生のための経済学入門［新版］ (ちくま新書１７７９)Amazon（アマゾン）マクロ経済学やミクロ経済学の違いがわからないレベルでも読めたブレンデッドスコッチウイスキーバランタインファイネスト40度 [ ウイスキー日本 1750ml ]Amazon（アマゾン）21年物は一万越えジムビーム 700ml [JimBeam バーボンウイスキーアメリカ合衆国][SUNTORY'S WHISKY]Amazon（アマゾン）バーボンウイスキーの定番響(HIBIKI) JAPANESE HARMONY 43°700mlAmazon（アマゾン）ジャパニーズウィスキーは入手困難ヨサソーダ炭酸水無糖 190ml 20本缶炭酸飲料強炭酸Amazon（アマゾン）これで割ればアルコール度数がだいたい５%よろしくお願いいたします。以上です！
12Mar
- ExcelでNORM.DISTやNORM.S.DISTをBETA.DISTで近似する方法
  本日は、ExcelでNORM.DIST関数やNORM.S.DIST関数をBETA.DIST関数で近似する方法を解説したい。平均0、標準偏差1における正規分布をベータ分布で近似するにはz=1.96（z値）とすると標準正規分布の累積分布関数はNORM.S.DIST(z,TRUE)=NORM.S.DIST(1.96,TRUE)=0.97500210485178と計算できる。平均0、標準偏差1における正規分布の累積分布関数はNORM.DIST(z,0,1,TRUE)=NORM.DIST(1.96,0,1,TRUE)=0.97500210485178と計算できる。正規分布をベータ分布で近似するにはBETA.DIST(z*SQRT(0.5*0.5/20000000000)+0.5,10000000000,10000000000,TRUE)=BETA.DIST(1.96*SQRT(0.5*0.5/20000000000)+0.5,10000000000,10000000000,TRUE)=0.97500210485653と計算できる。完全に一致はしないが、小数点以下12桁程度まで合っている。平均と標準偏差が存在する正規分布を近似するには統計量=z=1.96平均=μ=1.1標準偏差=σ=2.1とすると標準正規分布の累積分布関数はNORM.S.DIST((z-μ)/σ,TRUE)=NORM.S.DIST((1.96-1.1)/2.1,TRUE)=0.658922351248734平均1.1、標準偏差2.1に対する正規分布の累積分布関数はNORM.DIST(z,μ,σ,TRUE)=NORM.DIST(1.96,1.1,2.1,TRUE)=0.658922351248734正規分布をベータ分布で近似するにはBETA.DIST((z-μ)/σ*SQRT(0.5*0.5/20000000000)+0.5,10000000000,10000000000,TRUE)=BETA.DIST((1.96-1.1)/2.1*SQRT(0.5*0.5/20000000000)+0.5,10000000000,10000000000,TRUE)=0.65892235125072完全に一致はしないが、小数点以下10桁程度まで合っている。参考にしてください。高校生のための経済学入門［新版］ (ちくま新書１７７９)Amazon（アマゾン）マクロ経済学やミクロ経済学の違いがわからないレベルでも読めたブレンデッドスコッチウイスキーバランタインファイネスト40度 [ ウイスキー日本 1750ml ]Amazon（アマゾン）21年物は一万越えジムビーム 700ml [JimBeam バーボンウイスキーアメリカ合衆国][SUNTORY'S WHISKY]Amazon（アマゾン）バーボンウイスキーの定番響(HIBIKI) JAPANESE HARMONY 43°700mlAmazon（アマゾン）ジャパニーズウィスキーは入手困難ヨサソーダ炭酸水無糖 190ml 20本缶炭酸飲料強炭酸Amazon（アマゾン）これで割ればアルコール度数がだいたい５%よろしくお願いいたします。以上です！
11Mar
- Excelで複素数の計算をする方法
  本日は、Excelで複素数の計算をする方法を解説したい。Excelでは複素数を扱った計算が可能です。セルB1=COMPLEX(2,2) 　　　複素数AセルB2=COMPLEX(2,-2)　　　複素数BセルB4=IMSUM(B1,B2)　　　A＋BセルB5=IMSUB(B1,B2)　　　A-BセルB6=IMPRODUCT(B1,B2)　A×BセルB7=IMDIV(B1,B2)　　　　A／BセルB8=IMPOWER(B1,3)　　　Aの３乗セルB9=IMABS(B1)　　　　　Aの絶対値セルB10=IMLN(B1)　　　　　Aの対数セルB11=IMEXP(B1)　　　　Aの指数セルB12=IMSQRT(B1)　　　　Aの平方根セルB13=IMREAL(B2)　　　　Bの実部セルB14=IMAGINARY(B2)　　Bの虚部主な関数は上記です。計算をするときアセンブラでプログラムを組むみたいになりますが非常に便利です。参考にしてください。高校生のための経済学入門［新版］ (ちくま新書１７７９)Amazon（アマゾン）マクロ経済学やミクロ経済学の違いがわからないレベルでも読めたブレンデッドスコッチウイスキーバランタインファイネスト40度 [ ウイスキー日本 1750ml ]Amazon（アマゾン）21年物は一万越えジムビーム 700ml [JimBeam バーボンウイスキーアメリカ合衆国][SUNTORY'S WHISKY]Amazon（アマゾン）バーボンウイスキーの定番響(HIBIKI) JAPANESE HARMONY 43°700mlAmazon（アマゾン）ジャパニーズウィスキーは入手困難ヨサソーダ炭酸水無糖 190ml 20本缶炭酸飲料強炭酸Amazon（アマゾン）これで割ればアルコール度数がだいたい５%よろしくお願いいたします。以上です！
10Mar
- ExcelでNORM.DISTやNORM.S.DISTをGAMMA.DISTで近似する方法
  本日は、ExcelでNORM.DIST関数やNORM.S.DIST関数をGAMMA.DIST関数で近似する方法を解説したい。平均0、標準偏差1における正規分布をガンマ分布で近似するにはz=1.96（z値）とすると標準正規分布の累積分布関数はNORM.S.DIST(z,TRUE)=NORM.S.DIST(1.96,TRUE)=0.97500210485178と計算できる。平均0、標準偏差1における正規分布の累積分布関数はNORM.DIST(z,0,1,TRUE)=NORM.DIST(1.96,0,1,TRUE)=0.97500210485178と計算できる。正規分布をガンマ分布で近似するにはGAMMA.DIST(z*100000+10000000000,10000000000,1,TRUE)=GAMMA.DIST(1.96*100000+10000000000,10000000000,1,TRUE)=0.975001551302255と計算できる。完全に一致はしないが、小数点以下5桁程度まで合っている。平均と標準偏差が存在する正規分布を近似するには統計量=z=1.96平均=μ=1.1標準偏差=σ=2.1とすると標準正規分布の累積分布関数はNORM.S.DIST((z-μ)/σ,TRUE)=NORM.S.DIST((1.96-1.1)/2.1,TRUE)=0.658922351248734平均1.1、標準偏差2.1に対する正規分布の累積分布関数はNORM.DIST(z,μ,σ,TRUE)=NORM.DIST(1.96,1.1,2.1,TRUE)=0.658922351248734正規分布をガンマ分布で近似するにはGAMMA.DIST((z-μ)/σ*100000+10000000000,10000000000,1,TRUE)=GAMMA.DIST((1.96-1.1)/2.1*100000+10000000000,10000000000,1,TRUE)=0.658923369011137完全に一致はしないが、小数点以下5桁程度まで合っている。高校生のための経済学入門［新版］ (ちくま新書１７７９)Amazon（アマゾン）マクロ経済学やミクロ経済学の違いがわからないレベルでも読めたブレンデッドスコッチウイスキーバランタインファイネスト40度 [ ウイスキー日本 1750ml ]Amazon（アマゾン）21年物は一万越えジムビーム 700ml [JimBeam バーボンウイスキーアメリカ合衆国][SUNTORY'S WHISKY]Amazon（アマゾン）バーボンウイスキーの定番響(HIBIKI) JAPANESE HARMONY 43°700mlAmazon（アマゾン）ジャパニーズウィスキーは入手困難ヨサソーダ炭酸水無糖 190ml 20本缶炭酸飲料強炭酸Amazon（アマゾン）これで割ればアルコール度数がだいたい５%よろしくお願いいたします。以上です！
06Mar
- Excelにおける混合ポアソン分布の確率質量関数と乱数の生成方法
  本日は、Excelにおける混合ポアソン分布の確率質量関数と乱数の生成方法について解説したい。混合率をp、平均をλ1、λ2とすると混合ポアソン分布の確率質量関数はとなります。次に混合ポアソン分布に従う乱数を生成し上記確率質量関数をチェックしたい。生成方法については『ExcelのVBAでポアソン過程とポアソン分布の関係を示す』本日は、ポアソン過程とポアソン分布の関係をExcelのVBAで示したい。ポアソン過程はある確率でイベント数が１だけ起きるような確率過程です。同時にイベントが…ameblo.jpを参考にしてください。VBAプログラムのlamdaとNumberを変更すれば、A列にポアソン分布に従う乱数が生成されます。ポアソン分布に従う乱数を混合ポアソン分布に従う乱数に変換します。上記VBAプログラムを利用してlamda=5,Number=10000として実行しA列にコピーしlamda=10,Number=10000として実行しB列にコピーします。セルC2=RAND()とし、C10001までコピーします。混合率を0.3とし、セルI1に入力します。セルD2=IF(C2>1-$I$1,A2,B2)としD10001までコピーします。D2からD10001までに、混合ポアソン分布に従う乱数が生成されます。x=kにおける乱数の個数と理論値を計算したい。セルG4=COUNTIF($D$2:$D$10001,F4)としてG20までコピーします。セルH4=($I$1*POISSON.DIST(F4,$F$2,FALSE)+(1-$I$1)*POISSON.DIST(F4,$G$2,FALSE))*10000としてH20までコピーします。これをグラフにまとめると下記のようになります。上記乱数の生成法が理論値とほぼ一致した。高校生のための経済学入門［新版］ (ちくま新書１７７９)Amazon（アマゾン）マクロ経済学やミクロ経済学の違いがわからないレベルでも読めたブレンデッドスコッチウイスキーバランタインファイネスト40度 [ ウイスキー日本 1750ml ]Amazon（アマゾン）21年物は一万越えジムビーム 700ml [JimBeam バーボンウイスキーアメリカ合衆国][SUNTORY'S WHISKY]Amazon（アマゾン）バーボンウイスキーの定番響(HIBIKI) JAPANESE HARMONY 43°700mlAmazon（アマゾン）ジャパニーズウィスキーは入手困難ヨサソーダ炭酸水無糖 190ml 20本缶炭酸飲料強炭酸Amazon（アマゾン）これで割ればアルコール度数がだいたい５%よろしくお願いいたします。以上です！

前ページ
次ページ