Ошибка с подгонкой обобщенного экстремального значения (GEV) с использованием `extRemes` в R? - PullRequest
Новая
       15

Ошибка с подгонкой обобщенного экстремального значения (GEV) с использованием `extRemes` в R?

1 голос
/ 10 апреля 2019

У меня есть некоторые данные, и я хочу соответствовать распределению обобщенного экстремального значения (GEV), используя пакет extRemes в R. Однако возникает ошибка: 

library(extRemes)
Mydata = c(6,3,3,3,5,5,4,3,5,5,4,3,4,4,6,5,5,4,5,2,6,4,6,5,3,3,8,3,4,4,6,6,6,6,6,5,6,6,5,5)
fit_gev <- fevd(x=Mydata, method = "MLE", type="GEV", period.basis = "year")
summary(fit_gev)

Error in diag(cov.theta) : invalid 'nrow' value (too large or NA)
In addition: Warning message:
In diag(cov.theta) : NAs introduced by coercion

Интересно, как я могу исправить эту ошибку? Спасибо за любую помощь.

1 Ответ

1 голос
/ 30 апреля 2019

Вы можете настроить с помощью пакета EnvStats следующим образом: 

library(EnvStats)
# Data
Mydata =
c(6,3,3,3,5,5,4,3,5,5,4,3,4,4,6,5,5,4,5,2,6,4,6,5,3,3,8,3,4,4,6,6,6,6,6,5,6,6,5,5)    
# Generalized Extreme Value (EnvStats)
egevd(Mydata, method = "pwme")# (Method: probability-weighted moments)

Results of Distribution Parameter Estimation
--------------------------------------------

Assumed Distribution:            Generalized Extreme Value

Estimated Parameter(s):          location = 4.268896
                                 scale    = 1.314489
                                 shape    = 0.353434

Estimation Method:               Unbiased pwme

Data:                            Mydata

Sample Size:                     40
...