計算方差-協(xié)方差矩陣

《計算方差-協(xié)方差矩陣》由會員上傳分享，免費在線閱讀，更多相關內容在教育資源-天天文庫。

1、lecture10FINANCIALMODELING金融建模第10章計算方差-協(xié)方差矩陣要計算有效投資組合，我們就必須計算股票收益數(shù)據(jù)的方差-協(xié)方差矩陣。本章中，我們將討論在Excel中怎樣實現(xiàn)這個計算。其中最顯而易見的計算為樣本方差-協(xié)方差矩陣：這是直接由歷史收益計算而得的矩陣。我們介紹幾種計算方差-協(xié)方差的方法，包括在電子表中用超額收益矩陣直接計算、VBA實現(xiàn)該方法計算。即使樣本方差-協(xié)方差矩陣看起來像一個很明顯的選擇，但我們將用大量的文字說明它也許不是方差與協(xié)方差最好的估計。樣本方差-協(xié)方差矩陣有兩個不盡人意的缺陷：一是它常使用不現(xiàn)實的參數(shù)，二是它難以用于預測。這些將主要在10.5和10

2、.6節(jié)中討論。作為樣本矩陣的替換，10.9和10.10節(jié)將討論用于優(yōu)化方差-協(xié)方差矩陣估計的“壓縮”方法。在開始本章之前，你應先閱讀第34章數(shù)組函數(shù)的內容。里面有一些Excel函數(shù)，其參數(shù)是向量和矩陣；它們的實施與標準Excel函數(shù)略有不同。本章重點討論這些數(shù)組函數(shù)Transpose()和MMult()，還有“自制”的數(shù)組函數(shù)的使用。10.1引言我們用我們的數(shù)字例子來說明計算方差-協(xié)方差矩陣的矩陣方法。我們通過減去資產(chǎn)各自的平均收益，得到超額收益矩陣（接下來的電子表中的42-52行）。在55-61行中我們計算樣本方差-協(xié)方差矩陣。10.2.1一個稍微更有效率的替代方法正如你所期望那樣，的確存在

3、其他計算方差-協(xié)方差矩陣可選方法。這里所介紹的方法跳過了超額收益的計算，并且直接使用單元格B71：G76中的公式進行計算。它通過使用數(shù)組函數(shù)=MMULT(TRANSPOSE(B23：G33-B35：G35)，B23：G33-B35：G35)/10。通過寫入B23：G33-B35我們直接將每項收益減去均值得到超額收益向量：10.3我們應該除以M還是M-1？Excel與統(tǒng)計量？在前面的計算中我們除以M-1而非M，以此得到無偏的方差和協(xié)方差的估計。不過這個選擇看起來幾乎沒有多大影響。我們引用主流的教科書：“對于為什么要用M-1取代M這兒有一段很長的歷史。如果你從來沒有聽說過，你可以參考任何一本好的統(tǒng)

4、計教材。這里我們主要想提醒你，如果你在計算一個分布的方差時，這個分布存在已知的先驗的均值，而不需要從歷史數(shù)據(jù)估計的時候，那么M-1應該變回M。（我們同樣想說關于在分母上用M-1替代M上，我們認為對你是已知的，但這卻是對你不負責任的——例如，試圖用圖例說明去證明一個充滿疑問的假設）”Excel本身某程度上在除以M還是M-1這個問題上也有些混亂。在下面的電子表中我們給出幾種計算均值，方差，標準差和協(xié)方差的方法。Excel區(qū)分總體方差（Varp，除以M）、樣本方差（Var，除以M-1），以及總體和樣本標準差（分別為Stdevp和Stdev）。但是Excel并沒有在協(xié)方差函數(shù)Covar中作此區(qū)分。你可

5、以看到B30中Covar除以M，和Varp一樣。如果你想得到一個相應的除以M-1的協(xié)方差函數(shù)，那么你得像單元格B33那樣用Covar乘以，或者你需要使用像單元格B32那樣的數(shù)組函數(shù)=MMULT(TRANSPOSE(B3：B13-B16)，C3：C13-C16)/10。如果Excel是完全合理的，它應該有兩個函數(shù)：Covarp，它除以M（對應Varp或Stdevp），以及Covar，它除以M-1（對應Var或Stdev）。困惑了嗎？沒關系！正如該部分開始的教科書引用指出的那樣，它不是一個至關重要的問題。10.4計算方差-協(xié)方差矩陣的其他方法在這一節(jié)中，我們介紹兩種替代計算方差-協(xié)方差矩陣的方法。

6、第一種是使用一個VBA數(shù)組函數(shù)，它可以直接計算出樣本方差-協(xié)方差矩陣。第二種是使用Excel的Offset函數(shù)。10.4.1一個計算方差-協(xié)方差矩陣的VBA函數(shù)我們的第一種替代方法是用一個VBA函數(shù)：我們的第一種替代方法是用一個VBA函數(shù)：FunctionVarCovar(rng.Range)AsVariantDimiAsIntegerDimjAsIntegerDimnumColsAsIntegernumCols=rng.Columns.CountDimmatrix()AsDoubleReDimmatrix(numCols-1,numCols-1)Fori=jTonumColsForj=1To

7、numColsmatrix(i-1,j-1)=Application.WorksheetFunction.Covar(rng.Columns(i),rng.Columns(j))NextjVarCovar=matrixEndFunction這個函數(shù)是一個數(shù)組函數(shù)。（意味著它必須使用[Ctrl]-[Shift]-[Enter]）。下面有一個例子：應注意以上6只股票的GMVP包含兩個空頭（BA和MSF