Nothing
## ----setup, echo = FALSE-------------------------------------------------
knitr::opts_chunk$set(echo = TRUE, message = FALSE, warning = FALSE, fig.align = 'center', fig.width = 4.65, fig.height = 4)
## ---- include=FALSE------------------------------------------------------
library("bookdown")
## ---- fig.show='hold', fig.cap = "Contaminated normal distribution Y", include=FALSE, eval=FALSE----
# #Setting $\pi = 0.3$ and $\sigma^2 = 1$ produces the following distribution for $n = 10000$:
# library("SimCorrMix")
# Nstcum <- calc_mixmoments(mix_pis = c(0.3, 0.7), mix_mus = c(0, 0),
# mix_sigmas = c(3, 1), mix_skews = c(0, 0), mix_skurts = c(0, 0),
# mix_fifths = c(0, 0), mix_sixths = c(0, 0))
# Nmix <- contmixvar1(n = 10000, "Polynomial", Nstcum[1], Nstcum[2]^2,
# mix_pis = c(0.3, 0.7), mix_mus = c(0, 0), mix_sigmas = c(3, 1),
# mix_skews = c(0, 0), mix_skurts = c(0, 0), mix_fifths = c(0, 0),
# mix_sixths = c(0, 0))
# plot_simpdf_theory(sim_y = Nmix$Y_mix[, 1], overlay = FALSE,
# legend.position = "none", title = "PDF of Y")
# qqnorm(Nmix$Y_mix[, 1], main = "Q-Q plot for Y")
# qqline(Nmix$Y_mix[, 1], col = 2)
# #The Q-Q plot shows that the assumption of homogeneous variance for the error distribution would not be valid in this situation
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.