Given a matrix and the output of noise type classification, compute the Hurst exponent of each matrix row and draw a box plot that plots for each noise type the range of Hurst exponents. Hurst exponents are computed with function hurstexp in the pracma package (simple R/S Hurst is used).

plotNoisetypesVsHurst(x, noisetypes, d = round(ncol(x)/3), header = "")

Arguments

x

a matrix

noisetypes

the noisetypes of the matrix

d

window size for Hurst exponent computation with function hurstexp

header

header string

Examples

if (FALSE) { N=50 M=500 metapop=generateAbundances(N=M, mode=5, probabs=TRUE) ts=simHubbell(N=N, M=M,I=1500,d=N, m.vector=metapop, tskip=50, tend=500) noisetypes=identifyNoisetypes(ts) noisetypesHurst(ts,noisetypes,header="Hurst exponent stratified by noise type") }