loafercreek(soilDB)
loafercreek()所属R语言包:soilDB
LoaferCreek Soils
LoaferCreek土壤
译者:生物统计家园网 机器人LoveR
描述----------Description----------
A SoilProfileCollection Object of Loafercreek Soils from the CA630 Soil Survey Area.
一个SoilProfileCollection对象的Loafercreek CA630土壤普查区的土壤。
用法----------Usage----------
data(loafercreek)
格式----------Format----------
实例----------Examples----------
require(reshape)
data(loafercreek)
##[#]
## aggregate major horizon types over 1cm intervals[#主要地平线总额超过1cm间隔]
##[#]
# categorize major horizon types[主要地平线类型分类]
hz.tab <- rev(sort(table(loafercreek$hzname)))
hz.tab[hz.tab > 5]
# add generalized hz name[添加广义赫兹名]
loafercreek$hz <- rep(NA, times=nrow(horizons(loafercreek)))
# generalize horizons[概括视野]
loafercreek$hz[grep('O', loafercreek$hzname)] <- 'O'
loafercreek$hz[grep('A', loafercreek$hzname)] <- 'A'
loafercreek$hz[grep('Bt', loafercreek$hzname)] <- 'Bt'
loafercreek$hz[grep('Cr', loafercreek$hzname)] <- 'Cr'
loafercreek$hz[grep('R', loafercreek$hzname)] <- 'R'
# convert generalized hz to factor[转换广义Hz至因素]
loafercreek$hz <- factor(loafercreek$hz)
loafercreek.hz.agg <- slab(loafercreek, fm= ~ hz)
# wide->long format[宽长格式]
loafercreek.hz.agg.long <- melt(loafercreek.hz.agg,
id.var=c('top','bottom', 'contributing_fraction', 'all.profiles', 'variable'),
variable_name='horizon')
# plot horizon type proportions[图地平线类型比例]
require(lattice)
xyplot(top ~ value, groups=horizon, data=loafercreek.hz.agg.long, ylim=c(150, -5),
type=c('S','g'), horizontal=TRUE, subset=value > 0, asp=1.5, ylab='', xlab='',
auto.key=list(columns=5, lines=TRUE, points=FALSE))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|