- Notifications
You must be signed in to change notification settings - Fork14
Open
Description
Consider profile-specific charge density for method = "E", or adaptive selection of "E" vs. "S".
Some options for profile-specific settings:
profileInformationIndex()- "number of horizons per 100cm of depth"
library(aqp)library(soilDB)s<- c('inks' ,'pardee','clarksville','palau','hao','inks','eheuiki','puaulu','zook','cecil')x<- fetchOSD(s)par(mar= c(0,0,0,3)).args<-list(width=0.3,name.style='center-center',hz.depths=TRUE,cex.names=1)options(.aqp.plotSPC.args=.args)plotSPC(x,fixOverlapArgs=list(method='E',q=1.25),max.depth=151)plotSPC(x,fixOverlapArgs=list(method='E',q=1),max.depth=151)plotSPC(x,fixOverlapArgs=list(method='E',q=0.75),max.depth=151)plotSPC(x,fixOverlapArgs=list(method='E',q=0.5),max.depth=151)plotSPC(x,fixOverlapArgs=list(method='E',q=0.25),max.depth=151)plotSPC(x,fixOverlapArgs=list(method='S'),max.depth=151)# profile information indexprofileInformationIndex(x,vars='top',baseline=TRUE,method='median')# number of horizons / 100cm depthhzdensity<- profileApply(x,function(i) { (nrow(i)/ max(i))*100})q<-aqp:::.rescaleRange(hzdensity,0.25,1.25)# ...
The latter is simpler and probably more relevant to label placement.
Metadata
Metadata
Assignees
Labels
No labels