pop.cdf.density: Overlay 2 PDFs as Weighted Histograms

In ejanalysis/ejanalysis: Tools for Environmental Justice (EJ) Analysis

Overlay 2 PDFs as Weighted Histograms

Description

OK but there now are better ways to compare pdf plots

Usage

pop.cdf.density(e, dcount, refcount, etxt, dtxt, brks = 10, ...)

Arguments

e	Environmental or other indicator values vector
dcount	Vector of weights for the demographic group of interest, such as population counts of Hispanics by Census tract.
refcount	Vector of weights for the reference group, such as population counts of individuals who are not Hispanic, by Census tract.
etxt	Character string to name e in graph
dtxt	Character string to name d in graph
brks	Default is 10. Passed as breaks param to plot function
...	other parameters passed to plot # used to pass to density() but this is more useful

Details

For an easy, nice, smoothed density plot, try plot(stats::density(bw = .001)); points(density(, bw = .001), col="red") to see overlay of two frequency distributions. Assumes you have weights for each and are comparing values in one group vs another.

Value

Creates a plot

See Also

pop.cdf() pop.cdf2() pop.ecdf() pop.cdf.density()

Examples

## Not run: 
  bg <- ejscreen::bg22[, c(ejscreen::names.d, 'pop', ejscreen::names.e, 'REGION')]

e <- bg$pm[!is.na(bg$pm)]
dpct <- bg$pctmin
dcount   <- bg$pop[!is.na(bg$pm)] *      dpct[!is.na(bg$pm)]
refcount <- bg$pop[!is.na(bg$pm)] * (1 - dpct[!is.na(bg$pm)])
brks <- 0:17
etxt <- 'PM2.5'
dtxt <- 'Minorities'

pop.cdf(        e, pcts = dpct, pops = bg$pop)
pop.cdf2(       e, dcount, refcount, etxt, dtxt, brks)
pop.cdf.density(e, dcount, refcount, etxt, dtxt )

  pop.cdf.density(e = e, dcount = dcount, refcount = refcount, etxt = etxt, dtxt = dtxt,
               adjust=2, brks = brks)

## End(Not run)

ejanalysis/ejanalysis documentation built on April 2, 2024, 10:12 a.m.