pfKruskal: Analyse composite records by a Kruskal-Wallis ANOVA

In 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)

Example output

Loading required package: GCD\nThis is paleofire v1.2.4 | GCD v4.0.7\nWarning messages:\n1: In get(method, envir = home) : input string \'[\x87]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\n2: In get(method, envir = home) : input string \'[\x85]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\n3: In get(method, envir = home) : input string \'[\x82]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\n4: In get(method, envir = home) : input string \'[\x8a]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\n5: In get(method, envir = home) : input string \'[\x87]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\n6: In get(method, envir = home) : input string \'[\x85]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\n7: In get(method, envir = home) : input string \'[\x82]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\n8: In get(method, envir = home) : input string \'[\x8a]\' cannot be translated to UTF-8, is it valid in \'ANSI_X3.4-1968\'?\nLoading and preparing data...\nTransforming...\nPercentage done: \r  |                                                                            \r  |                                                                      |   0%\r  |                                                                            \r  |===================================                                   |  50%\r  |                                                                            \r  |======================================================================| 100%\r  |                                                                            \r  |                                                                      |   0%\r  |                                                                            \r  |===================================                                   |  50%\r  |                                                                            \r  |======================================================================| 100%\r  |                                                                            \r  |                                                                      |   0%\r  |                                                                            \r  |===================================                                   |  50%\r  |                                                                            \r  |======================================================================| 100%\n$statistics\n     Chisq Df    p.chisq\n  22.48359 12 0.03244408\n\n$parameters\n            test p.ajusted name.t ntr alpha\n  Kruskal-Wallis      none data$x  13  0.05\n\n$means\n            data.y      rank       std  r        Min        Max         Q25\n0     -0.008673841 104.25806 0.7501623 31 -1.1939011  1.5804262 -0.56776866\n1000  -0.080565767 104.67857 0.2959661 28 -0.6351988  0.6564250 -0.26396701\n10000 -0.647049478  63.22222 0.9946368  9 -2.0351710  1.2126036 -1.04413317\n11000 -0.990827915  41.25000 0.8940749  4 -2.1264410 -0.1845970 -1.48839521\n12000 -1.682388842  27.50000 1.7090567  2 -2.8908744 -0.4739033 -2.28663163\n2000   0.141128665 135.30769 0.3177323 26 -0.6790635  0.6749508  0.04529697\n3000  -0.066994258 105.39130 0.3392909 23 -0.7804535  0.4326644 -0.22830267\n4000   0.054090211 119.47368 0.4629597 19 -0.9271397  0.7923088 -0.17119768\n5000   0.155506625 128.05263 0.7047172 19 -1.4033288  1.4329732 -0.27947213\n6000  -0.081812700 113.64706 1.0492837 17 -2.1255491  1.2861596 -0.85012222\n7000  -0.122085023 112.12500 1.1155425 16 -2.2584230  1.1269163 -0.73679914\n8000  -0.209352784 108.15385 1.4582505 13 -3.4761013  1.8868693 -0.73086407\n9000  -0.366749472  78.55556 0.9222576  9 -1.8184255  1.4018422 -0.72606816\n              Q50         Q75\n0     -0.08865951  0.08325314\n1000  -0.07337137  0.09864429\n10000 -0.68361637 -0.16786714\n11000 -0.82613683 -0.32856954\n12000 -1.68238884 -1.07814605\n2000   0.11998518  0.35787440\n3000  -0.14066682  0.30081137\n4000   0.06069577  0.43846081\n5000   0.14963591  0.58773242\n6000   0.20830457  0.85237441\n7000  -0.01823844  0.81110168\n8000  -0.18500951  0.80163551\n9000  -0.55342708  0.05268799\n\n$comparison\nNULL\n\n$groups\n     trt     means    M\n1   2000 135.30769    a\n2   5000 128.05263    a\n3   4000 119.47368   ab\n4   6000 113.64706   ab\n5   7000 112.12500  abc\n6   8000 108.15385 abcd\n7   3000 105.39130 abcd\n8   1000 104.67857 abcd\n9      0 104.25806 abcd\n10  9000  78.55556  bcd\n11 10000  63.22222   cd\n12 11000  41.25000    d\n13 12000  27.50000    d\n\nattr(,"class")\n[1] "group"\n

paleofire documentation built on Jan. 11, 2020, 9:44 a.m.