Detailed description |
Hi Dylan, Pierre,
Im using aqp 1.7-9 (directly from R-forge), on linux.
When i use profile_compare() to calculate a dissimilarity matrix the result its 2 magnitude orders different if I use only horizon data OR if i use both horizon and site data.
The range for the first diss matrix its 0-100, and for the sencond one 0-1.
This its an issue when i try to plot the data.
CODE
data(sp1)
## init SoilProfileCollection objects from data.frame
depths(sp1) <- id ~ top + bottom
## setting site-level data
# site-level data from horizon-level data (stored in @horizons)
site(sp1) <- ~ group
d <- data.frame(id=profile_id(sp1), p=runif(n=length(sp1)), q=runif(n=length(sp1)),stringsAsFactors=FALSE)
site(sp1) <- d
sp1
d.horizon <- profile_compare(sp1, vars=c('prop', 'structure_grade'), max_d=100, k=0.01, plot.depth.matrix=TRUE)
range(d.horizon)
d.site <- profile_compare(sp1, vars=c('prop', 'structure_grade', 'p', 'q'), max_d=100, k=0.01, plot.depth.matrix=TRUE)
range(d.site)
END CODE
Regards!
|
|