2011-05-06 57 views
2

如何使用包「廣義雙曲分佈」來使NIG分佈中的參數適合數據集?R中的廣義雙曲線分佈

+0

什麼是「NIG」分佈? – Aniko 2011-05-06 20:56:35

+1

'ghyp'軟件包的功能是'fit.NIGuv'(用於單變量數據)和'fit.NIGmv'(用於多變量)數據,並且在文檔中對這些軟件包都有非常清楚的描述。你看過或嘗試過嗎? – 2011-05-06 20:59:00

+0

嗨!謝謝您的回答。但是,當我嘗試使用此功能時,我收到無法找到此功能的消息。我必須下載哪些軟件包才能使用此功能?當我下載了一個包以使其工作時,我應該做些什麼特別的事情? – Claire 2011-05-07 10:17:29

回答

1

制定@Prasad Chalasani的評論,您需要安裝ghyp包。當我這樣做時,包裝gtools,gdata, numDerivgplots也自動安裝。但是,然後我得到了與Error: could not find function "fit.NIGuv"相同的錯誤,我通過手動安裝bitops軟件包解決了這個錯誤。

documentation給出了使用fit.NIGuv()

library(ghyp) 
data(smi.stocks) 
nig.fit <- fit.NIGuv(smi.stocks[,"SMI"], opt.pars = c(alpha.bar = FALSE), 
        alpha.bar = 1, control = list(abs.tol = 1e-8)) 
nig.fit 
summary(nig.fit) 
hist(nig.fit) 

下面的代碼的例子,其中輸出包括

Asymmetric Normal Inverse Gaussian Distribution: 

Parameters: 
    alpha.bar   mu   sigma   gamma 
1.0000000000 0.0008370731 0.0112098776 -0.0007205143 

log-likelihood: 
5495.705 

,我認爲這是你正在尋找的那種東西。