mAPCA: Functions for mAPCA

In kralljr/share: Estimates Regional Sources of Air Pollution

Description

mAPCA Apply mAPCA to pollution data

Usage

mAPCA(data, tots = NULL, nsources = NULL, adjust = NULL,
  mdl = NULL, cut = 1, lim = 50)

Arguments

data	list of data frames of daily constituent concentrations with date as first column
tots	list of totals (PM2.5) corresponding to each monitor
adjust	method of handles adjustment
mdl	list of MDLs for each monitor if applicable
cut	cutoff for eigenvalues (see nsource), default is 1.
lim	Number of days to include data
nsource	number of sources. If null, uses number of eigenvalues of the correlation matrix greater than cut.

Examples

library(handles)
library(share)
data(consConc)
data <- lapply(consConc, function(x) x[, -which(colnames(x) == "PM25_SPEC")])
pm <- lapply(consConc, function(x) x$PM25_SPEC)
mapca1 <- mAPCA(data, pm)
# Look at source loadings
mapca1$apca$vmax$loadings 

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