vba functionプロシージャの扱いと計算

Question

前提・実現したいこと

乱数をn個生成し、その和Fを求め
Fe＝(√(12/n))*(F-(n/2))
の結果を関数値で返すfunctionプロシージャを作る
Feを複数回呼び出しFeの値を複数発生させ
さらにそれらFeの合計と平均を計算するサブプロシージャを作成する

平均はゼロに近づいているが分散は8.何々など１に近づいている様子がない

該当のソースコード

Sub heikin()

Dim n As Long, mmm As Double, ma As Double, aver As Double

Dim f As Long

Randomize
n = InputBox("n")
v = InputBox("v")

For y = 1 To v
mmm = sh(n) + mmm
Next y

ma = ma + mmm
  aver = ma / v

mau = mau + mmm ^ 2
  aver2 = mau / v

variance = aver2 - (aver) ^ 2

MsgBox "平均; " & aver & vbCrLf _
       & "分散; " & variance
End Sub



Function sh(n As Long) As Double
Dim f As Double, fe As Double

For o = 1 To n

    f = f + Rnd

Next o

fe = Sqr(12 / n) * (f - (n / 2))
sh = fe

End Function

###試したこと

functionプロシージャ内のfor文をサブプロシージャに変えたり
式を変えたりしてみた
漠然としていてすみません

Comments

[m.ts10806]

コードはマークダウンのcode機能をご利用ください。
また、「考えられること」は今のところ質問者さんしか知りません。
書かれていないことを読み取るのは不可能ですので、可能な範囲で記述いただけませんか？

[m.ts10806]

「うまくいきません」では起きている問題が分かるわけではないので、起きている現象、どういう結果となったか、エラーが出ていないか　等々　詳細記述してください。

[SORAshining]

ごめんなさい。初めての利用でよくわかっていませんでした。

[hatena19]

コードの分散の計算方が間違っています。
分散は平均と各値の差の2乗の平均なので、ループで回して計算する必要があります。

Answer 1

コード中の分散の計算法が間違っています。
分散は平均と各値の差の2乗の平均なので、ループで回して計算する必要があります。

まず、sh(n) の計算結果をループさせて配列にいれます。
その配列からループで平均を求めます。
さらに配列をループさせて平均と各値の差の2乗の平均を求めます。これが分散です。

下記を参考に。

Excel VBAの「配列にn個のデータを読み込み、その平均と分散... - Yahoo!知恵袋

---

コード例

Sub heikin()

    Randomize
    
    Dim n As Long, v As Long
    n = InputBox("n")
    v = InputBox("v")
    
    Dim ary() As Double  '結果格納用配列
    ReDim ary(1 To v) '配列のサイズを決定
    
    Dim y As Long
    For y = 1 To v
        ary(y) = sh(n) 'shの結果を配列に格納。
    Next y

    Dim ma As Double '総計
    For y = 1 To v
        ma = ma + ary(v)
    Next y
    Dim aver As Double '平均
    aver = ma / v

    Dim variance As Double '分散
    For y = 1 To v
        variance = variance + (ary(v) - aver) ^ 2
    Next y
    variance = variance / v

    MsgBox "平均; " & aver & vbCrLf _
           & "分散; " & variance
End Sub

Answer 2

Doubleを使用しているから、精度誤差によるものじゃないでしょうか。

Doubleの代わりにCurrencyを使ってみてはどうでしょう。
それでも精度に問題があるなら、整数部と小数部を分けるしか無いかと。

Comments

[SORAshining]

試してみます。ありがとうございます。

Answer 3

たくさん問題点がありますが、一番の問題は乱数を使っていない事です。

For o = 1 To n
　f = f + Rnd
Next o

で乱数の和を計算したつもりでしょうが、乱数を発生する関数の呼び出しは Rnd()であって、Rndではありません。
f=f+Rndは、変数fに変数fと変数Rndの値の和を代入していますが、変数Rndは定義されていないので値は0です。　f=f+0を何度繰り返しても、fは0のままです。

このようなエラーは、コードの最初に"Option Explicit"と指定する事で避けられます。
このオプションは、定義されていない変数が使われた時に警告を出します。

Comments

[SORAshining]

初歩的なミスですみません。ありがとうございます

[hatena19]

VBAの場合、引数のない場合は、() は省略できるので、
　f = f + Rnd
と
　f = f + Rnd()
は同じです。

[SORAshining]

あ、だから乱数取得は普通にできてたんですね。ありがとうございます。