pfKruskal: Analyse composite records by a Kruskal-Wallis ANOVA

In oblarquez/paleofire: Analysis of Charcoal Records from the Global Charcoal Database

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.

Author(s)

O. Blarquez

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)

## End(Not run)

oblarquez/paleofire documentation built on Dec. 29, 2021, 11 a.m.