data(trlremeasu)
head(trlremeasu)
df<-trlremeasu
df$fe<-10000/1600
df$garb.ha<- (pi/40000)*df$dbh^2*df$fe
gha.t<-tapply(df$garb.ha, df$year, sum)
nha.t<-tapply(df$fe, df$year, sum);
time<-as.numeric(rownames(gha.t))
plot(nha.t~time, type="b",las=1)
plot(gha.t~time, type="b",las=1)
Run the code above in your browser using DataLab