Now uses the annual inventory start years from Appendix J of the FIADB user guide to filter both raw data (in fia_tidy()) and EVALIDs (in fia_assign_strata()) to include only annual inventory years and exclude earlier periodic inventories.
Fix to internal function inter_extra_polate() which was not handling leading NAs correctly in linear interpolation (#225).
EVALIDs specific to East or West Texas won’t be matched to plots by fia_assign_strata().
Early EVALIDs in NM and WY that “don’t work” no longer assigned to plots by fia_assign_strata().
Changed behavior of fia_assign_strata() so that for each plot x year, the first EVALID that contains the YEAR (i.e. between START_INVYR and END_INVYR) is matched.
Changed package license to MIT.
Fixed a bug where for woodland species, carbon was being replaced with biomass.
As part of the workaround for woodland species (#163), any negative interpolated carbon and biomass values were set to 0.
fia_assign_strata() was added to match plots in each year of the annualized data to an EVALID, estimation unit, and stratum along with associated values necessary for population level estimation with variance calculations.
fia_tidy() now adds several columns from the POP_* tables relating to EVALIDs. Each plot in each INVYR is matched to the EVALID of each type that has the latest end year. Then, for EVALIDs ending in 01 and 02 (“EXPCURR” and “EXPVOL”, respectively), two columns are created for each of EVALID, ESTN_UNIT, STRATUMCD, P1POINTCNT, P1PNTCNT_EU, and AREA_USED suffixed with _EXPVOL or _EXPCURR. There are also two logical columns EVAL_TYPE_EXPCURR and EVAL_TYPE_EXPVOL. These get filled down and switch at the next inventory (rather than at the midpoint like everything else) (#192).
fia_tidy() no longer filters to base intensity plots and reverts to just filtering to INVYR≥2000.
Changeed the way composit IDs plot_ID and tree_ID are constructed by changing the order of STATECD and UNITCD so that now it is UNITCD_STATECD_COUNTYCD_PLOT (#189)
Fixed a bug in fia_tidy() resulting in duplicated SUBCYCLE columns (#173).
A temporary workaround was implmented to deal with missing code for calculating biomass and carbon for woodland species (#163). Eventually, we will modify fia_estimate() to be able to produce carbon and biomass estimates for woodland species and this can be reverted to avoid confusion.
fia_tidy() now keeps the CARBON_AG and DRYBIO_AG columns
fia_estimate() now overwrites CARBON_AG and DRYBIO_AG columns in the input unless the result of carbon estimation is NA (as is the case with woodland species).
The calculation of the EXPNS column no longer occurs as part of interpolate_data().
forestTIME 2.1.0
This package has been re-named from forestTIME.builder to forestTIME
fia_tidy() now only keeps base-intensity plots (INTENSITY == 1 & SUBCYCLE != 99 & SUBCYCLE != 0).
Fixed a bug that was causing TPA_UNADJ to not be populated for all trees in states that use macroplots (#160).
Fixed a bug in fia_split_composite_ids() that caused it to fail when tree_ID was NA (as it is in conditions with no observations). It now falls back on the information in plot_ID when tree_ID is in the data frame but NA (#149).
fia_estimate() now returns the additional variables DRIBIO_FOLIAGE, VOLTSGRS, and VOLTSSND in addition to DRYBIO_AG and CARBON_AG.
Code to deal with negative extrapolated values has moved to adjust_mortality(). Therefore, the results of interpolate_data() may now contain negative numbers, which are non-sensible. Use fia_annualize() whenever possible to ensure sensible results.
fia_download() arguments have changed. keep_zip now defaults to TRUE. extract options have changed from TRUE/FALSE to "forestTIME", "rFIA", "all", or "none" to extract just files needed by forestTIME, those used by rFIA (for compatibility), all files, or none.
Fixed a bug in interpolation of CONDPROP_UNAJ. Now the workflow retains all “empty” conditions (i.e. CONDIDs with no trees in them) and properly interpolates CONDPROP_UNADJ so the proportion for all conditions in a plot in a year sum to 1 (within rounding error) (#64).
fia_annualize() now adds and EXPNS column calculated as the total land area of the state in acres divided by the number of plots in the interpolated data. It should be usable in the same ways the EXPNS column in the “raw” FIA data can be used.
Renames the state_codes dataset to state_areas and adds a column for state land area in acres.