search

PalEON-Project / PalEON-FIA

Code to take the FIA composition, density, and biomass data and aggregate it to the PalEON grid
0 stars 2 forks source link

Evergreen PFT fails to produce model - no spcd's in table #6

Open SimonGoring opened 8 years ago

SimonGoring commented 8 years ago

@mdietze here's code that replicates the failure for Evergreen:

AllomAve(list(Evergreen = list(data.frame(acronym = c("ILOP", "MAVI2"), 
                                          spcd = c(591, 653))),
         components = 2,
         outdir     = "data/output/PEcAn_allom/",
         ngibbs     = 10000,
         dmin       = 17,
         dmax       = 151)

And the output I get is:

[1] "writing output to"        "data/output/PEcAn_allom/"
[[1]]
  acronym spcd
1    ILOP  591
2   MAVI2  653

[[2]]
[1] 2

[[1]]
[1] "READ.ALLOM.DATA: ** Warning no match of PFT to tally data **"

[[2]]
  acronym spcd
1    ILOP  591
2   MAVI2  653

[1] "READ.ALLOM.DATA: ** No allometric data available **"
$Evergreen
list()

Warning message:
In nu(allom$parm$Component.ID) : NAs introduced by coercion

This affects 50 ILOP & 4 MAVI2. Looking at the source data (in data/Table3_GTR-NE-319.v2.csv) it might be that neither of those taxa have allometry data entered, so that the model can't be generated. If that's the case, that we can't build an allometric model for these trees using the package, what's our option? Should we see if we can add in other taxa to at least come close? This is a very small proportion of all trees.

@paciorek also tagged.