# NOT RUN { pdata <- data.frame(Mean = 1.0, SD = exp(1.0)) pdata <- transform(pdata, y = rposnorm(n <- 1000, m = Mean, sd = SD)) # } # NOT RUN { with(pdata, hist(y, prob = TRUE, border = "blue", main = paste("posnorm(m =", Mean[1], ", sd =", round(SD[1], 2),")"))) # } # NOT RUN { fit <- vglm(y ~ 1, posnormal, data = pdata, trace = TRUE) coef(fit, matrix = TRUE) (Cfit <- Coef(fit)) mygrid <- with(pdata, seq(min(y), max(y), len = 200)) # Add the fit to the histogram # } # NOT RUN { lines(mygrid, dposnorm(mygrid, Cfit[1], Cfit[2]), col = "orange") # }
Run the code above in your browser using DataCamp Workspace