paleofire (version 1.2.4)

pfKruskal: Analyse composite records by a Kruskal-Wallis ANOVA

Description

The function applies a Kruskal-Wallis ANOVA on binned data issued from a "pfComposite" object (of directly on "pfTransform" objects), in order to test the difference in biomass burning activity between different time periods.

Usage

pfKruskal(data, p.adj = "none", alpha = 0.05, bins = NULL, verbose = TRUE)

Arguments

data

An object returned by pfComposite or pfTransform.

p.adj

Method for adjusting p values (see p.adjust). Includes: "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr" and "none" (default).

alpha

Numeric, confidence level.

bins

Numeric, bins to use if a "pfTransform" object is provided.

verbose

Logical, verbose or not...

Value

A "pfKruskal" object containing multiple comparison results.

See Also

plot.pfKruskal,kruskal

Examples

# NOT RUN {
## Composite charcoal record for Western Boreal North America:
ID=pfSiteSel(continent=="North America", long<(-100) & l12==1)
plot(ID)
## Transform data
res3=pfTransform(ID,method=c("MinMax","Box-Cox","Z-Score"),BasePeriod=c(200,4000))

## Composite
comp=pfComposite(res3,bins=seq(from=-500,to=12500,by=1000))
plot(comp)

## Kruskal Wallis Anova
comparison=pfKruskal(comp)
# }