[R-br] Modelo Regressão Beta Inflacionado

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