首先，讓我們了解copula的工作方式。
.seed（） < - < - < - mvrnorm（n，mu = rep（，m），Sigma = sigma，empirical = T）
我們使用cor()和散點圖矩陣檢查樣本相關性。
pairs.panels（Z） \[，\] \[，\] \[，\] \[，\] \[，\] - \[，\] -
建站客服QQ：88888888
圖片
pairs.panels（U）
這是包含新隨機變量的散點圖矩陣u。
點擊標題查閱往期內容
圖片
{jz:field.toptypename/}R語言多元Copula GARCH 模型時間序列預測
圖片
左右滑動查看更多
圖片
01
圖片
02
圖片
03
圖片
04
圖片
圖片
我們可以繪制矢量的3D圖表示u。
圖片
現在，作為最后一步，我們只需要選擇邊緣并應用它。我選擇了邊緣為Gamma，Beta和Student，并使用下面指定的參數。
< - qgamma（u \[，\]，shape = ，scale = ） < - qbeta（u \[，\]，，） < - qt（u \[，\]，df = ）
下面是我們模擬數據的3D圖。
圖片
< - cbind（x1，x2，x3） .panels（DF） x2 x3 . . . . . -. . -. .
米蘭體育官方網站 37， 41); text-align: start; white-space: normal; background-color: rgb(255， 255， 255);">這是隨機變量的散點圖矩陣：
圖片
使用copula讓我們使用copula復制上面的過程。
現在我們已經通過copula（普通copula）指定了相依結構并設置了邊緣，mvdc()函數生成了所需的分布。然后我們可以使用rmvdc()函數生成隨機樣本。
colnames（Z）< - （“”，“”，“”） pairs.panels（Z）
模擬數據當然非常接近之前的數據，顯示在下面的散點圖矩陣中：
圖片
簡單的應用示例現在為現實世界的例子。我們將擬合兩個股票 ，并嘗試使用copula模擬 。
讓我們在R中加載 ：
cree < - .csv（，莊閑和 = F）$ V2 yahoo < - .csv（， = F）$ V2
在直接進入copula擬合過程之前，讓我們檢查兩個股票收益之間的相關性并繪制回歸線：
我們可以看到 正相關 ：
圖片
在上面的第一個例子中，我選擇了一個正態的copula模型，但是，當將這些模型應用于實際數據時，應該仔細考慮哪些更適合數據。例如，許多copula更適合建模非對稱相關，其他強調尾部相關性等等。我對股票收益率的猜測是，t-copula應該沒問題，但是猜測肯定是不夠的。本質上， 允許我們通過函數使用BIC和AIC執行copula選擇 ：
pobs（.matrix（cbind（cree，yahoo））） selectedCopula $ PAR $ PAR2
擬合算法確實選擇了t-copula并為我們估計了參數。 讓我們嘗試擬合建議的模型，并檢查參數擬合。
.cop .seed（） < - pobs（as.matrix（cbind（cree，yahoo））） （FIT） . df . .
我們來看看我們剛估計的copula的密度
< - coef（fit）\[\] < - coef（fit）\[\]
圖片
現在我們只需要建立Copula并從中抽取3965個隨機樣本。
rCopula（，tCopula（ = ， ，df = df）） \[，\] \[，\] \[，\] \[，\]
這是包含的樣本的圖：
圖片
t-copula通常適用于在極值（分布的尾部）中存在高度相關性的現象。現在我們面臨困難：對邊緣進行建模。為簡單起見，我們將假設正態分布 。因此，我們估計邊緣的參數。
直方圖顯示如下：
圖片
圖片
現在我們在函數中應用copula，從生成的多變量分布中獲取模擬觀測值。最后，我們將模擬結果與原始數據進行比較。
這是在假設正態分布邊緣和相依結構的t-copula的情況下數據的最終散點圖：
圖片
正如您所看到的，t-copula導致結果接近實際觀察結果 。
讓我們嘗試df=1和df=8：
顯然，該參數df對于確定分布的形狀非常重要。隨著df增加，t-copula傾向于正態分布copula。
圖片
圖片
圖片
本站僅提供存儲服務，所有內容均由用戶發布，如發現有害或侵權內容，請點擊舉報。
備案號: