share: share.

In kralljr/share: Estimates Regional Sources of Air Pollution

Description

share estimates major sources of pollution using the share method

Usage

share(data, cut = 1, nmsources = NULL, thres = pi/4, ndays = 50)

Arguments

data	data frame of daily constituent concentrations with date as first column
cut	cutoff for eigenvalues (see nmsource), default is 1.
nmsources	number of major sources. If null, uses number of eigenvalues of the correlation matrix greater than cut.
thres	cutoff for angle between local and major sources
ndays	number of days of data necessary to apply PCA for each monitor

Details

This function estimates major sources of pollution across multiple ambient monitors. Other outputs of share include a matrix to guide pooling short-term health effects of pollution to conduct regional and national studies. Works on list where each element of the list is a dataframe corresponding to one ambient monitor. The first column of each dataframe is date and all subsequent columns are concentrations of chemical constituents.

Examples

library(handles)
library(share)
data(consConc)
# remove PM2.5
data <- lapply(consConc, function(x) x[, -which(colnames(x) == "PM25_SPEC")])
share1 <- share(data)
# Look at major sources
round(share1$major.sig, 2)

kralljr/share documentation built on May 20, 2019, 1:25 p.m.