inst/old_files/SDEF_Biomass_Analysis_22July2017.R

In bmcnellis/SDEF.analysis: Data Analysis for 2016 San Dimas Experimental Forest Project

# 22 July 2017
# Brandon McNellis

# Load functions
source('~/docs/R_workspace/SDEF.analysis/R/CalculateBioDensity.R')
source('~/docs/R_workspace/SDEF.analysis/R/CalculateCarbonContent.R')
source('~/docs/R_workspace/SDEF.analysis/R/ConvertFrameDepth.R')

# Load data

data.root <- read.table("SDEF_2017_RootLengths.txt", sep="\t", header=TRUE)
data.root$date <- as.character(data.root$date)
data.root$date <- as.Date(data.root$date, format="%m/%d/%Y")
data.ash <- read.table("SDEF_2016_RootAshWeights.txt", sep="\t", header=TRUE)
colnames(data.root) <- ifelse(colnames(data.root) == "tile", "frame", colnames(data.root))

data.root <- ConvertFrameDepth(data=data.root)
data.ash <- CalculateCarbonContent(data=data.ash)
# Frame height 6.75mm, frame width 9.0 mm
data.root <- CalculateBioDensity(data=data.root, frame.size = c(6.75, 9.0, 0.74), scale="tube")
# end