Prezados,

Ospina sugere a construção de bandas de confiança (envelopes) simuladas com o seguinte comando:

#==== normal probability plots ====

envelope.BIc(fit, type="quantil", main="quantal residual")

mas não estou conseguindo executá-lo no R.

Para executar o modelo utilizei as seguintes bibliotecas:

library(gamlss)

library(gamlss.dist)

library(gamlss.nl)

library(gamlss.mx)

library(gamlss.cens)

library(gamlss.tr)

library(foreign)

library(sas7bdat)

Alguém pode me ajudar?

Att., Arlene