R/plotNoisetypesVsHurst.R
plotNoisetypesVsHurst.Rd
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 = "")
x | a matrix |
---|---|
noisetypes | the noisetypes of the matrix |
d | window size for Hurst exponent computation with function hurstexp |
header | header string |
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") }