Plot_pdf.R
In gbutils: Utilities for Simulation, Plots, Quantile Functions and Programming

### R code from vignette source 'Plot_pdf.Rnw'

###################################################
### code chunk number 1: Plot_pdf.Rnw:61-62
###################################################
library(gbutils)


###################################################
### code chunk number 2: Plot_pdf.Rnw:77-80
###################################################
pdf1 <- function(x) dnorm(x, mean = 100, sd = 5) # pdf
qdf1 <- function(x) qnorm(x, mean = 100, sd = 5) # qf
cdf1 <- function(x) pnorm(x, mean = 100, sd = 5) # cdf


###################################################
### code chunk number 3: Plot_pdf.Rnw:85-86
###################################################
plotpdf(pdf1, qdf1)


###################################################
### code chunk number 4: Plot_pdf.Rnw:91-92
###################################################
plotpdf(pdf1, cdf = cdf1)


###################################################
### code chunk number 5: Plot_pdf.Rnw:97-98
###################################################
plotpdf(pdf1, cdf = cdf1, lq = 0.001, uq = 0.999)


###################################################
### code chunk number 6: Plot_pdf.Rnw:101-102
###################################################
plotpdf(cdf1, cdf = cdf1, lq = 0.001, uq = 0.999) # plot a cdf


###################################################
### code chunk number 7: Plot_pdf.Rnw:106-112
###################################################
pf1 <- function(x){
    0.25 * pnorm(x, mean = 3, sd = 0.2) + 0.75 * pnorm(x, mean = -1, sd = 0.5)
}
df1 <- function(x){
    0.25 * dnorm(x, mean = 3, sd = 0.2) + 0.75 * dnorm(x, mean = -1, sd = 0.5)
}


###################################################
### code chunk number 8: Plot_pdf.Rnw:116-117
###################################################
plotpdf(df1, cdf = pf1) # plot the pdf


###################################################
### code chunk number 9: Plot_pdf.Rnw:121-122
###################################################
plotpdf(pf1, cdf = pf1) # plot the cdf


###################################################
### code chunk number 10: Plot_pdf.Rnw:127-129
###################################################
plotpdf(pf1, cdf = pf1)                           # plot the cdf
plotpdf(df1, cdf = pf1, add = TRUE, col = "blue") # overlay the pdf


###################################################
### code chunk number 11: Plot_pdf.Rnw:136-137
###################################################
c(q5pc = cdf2quantile(0.05, pf1), q95pc = cdf2quantile(0.95, pf1))