adding modMAratio + function documentation #23
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
|
note: any example using # have to do it in this order to ensure we use my dev version of FIESTAutils
devtools::install_github("joshyam-k/FIESTA")
devtools::install_github("joshyam-k/FIESTAutils")
MApopdat <- FIESTA::modMApop(popTabs = list(tree = FIESTA::WYtree,
cond = FIESTA::WYcond),
pltassgn = FIESTA::WYpltassgn,
pltassgnid = "CN",
unitarea = FIESTA::WYunitarea,
unitvar = "ESTN_UNIT",
unitzonal = FIESTA::WYunitzonal,
prednames = c("dem", "tcc", "tpi", "tnt"),
predfac = "tnt")
#> INTENSITY not in dataset.. assuming single intensity plots
#> DESIGNCD not in dataset.. assuming DESIGNCD = 1 (annual inventory)
#> For FIA estimation, adjustment factors are calculated to account for plots with partial nonresponse.
#> ...there are 14 nonsampled forest conditions in the dataset.
#> filter removed 14 records: COND_STATUS_CD != 5
out <- FIESTA::modMAratio(MApopdat,
estvarn = "VOLCFNET",
estvarn.filter = "STATUSCD == 1",
rowvar = "STDSZCD",
prednames = c("dem", "tcc", "tpi"))
#> filter removed 2620 records: COND_STATUS_CD == 1
#> multiplying VOLCFNET by TPA
#> using the following predictors...dem, tcc, tpi
#> getting output...
str(out)
#> List of 4
#> $ est :'data.frame': 5 obs. of 3 variables:
#> ..$ Stand-size class : chr [1:5] "1" "2" "3" "5" ...
#> ..$ Estimate : chr [1:5] "17665.1" "14011.4" "1907.7" "129.7" ...
#> ..$ Percent Sampling Error: chr [1:5] "4.34" "10.4" "8.62" "8.32" ...
#> $ raw :List of 15
#> ..$ unit_totest :'data.frame': 23 obs. of 16 variables:
#> .. ..$ ESTN_UNIT : num [1:23] 1 3 5 7 9 11 13 15 17 19 ...
#> .. ..$ rhat : num [1:23] 346.6 336.6 69.1 1060.3 259.7 ...
#> .. ..$ rhat.var : num [1:23] 2559 92075 2546 8214 6048 ...
#> .. ..$ NBRPLT : int [1:23] 24 12 14 36 18 26 53 NA 8 22 ...
#> .. ..$ NBRPLT.gt0: int [1:23] 23 9 12 34 14 26 46 NA 5 21 ...
#> .. ..$ ACRES : num [1:23] 2757613 2021729 3072988 5096959 2729653 ...
#> .. ..$ AREAUSED : num [1:23] 2757613 2021729 3072988 5096959 2729653 ...
#> .. ..$ rhat.se : num [1:23] 50.6 303.4 50.5 90.6 77.8 ...
#> .. ..$ rhat.cv : num [1:23] 0.146 0.9015 0.7307 0.0855 0.2994 ...
#> .. ..$ pse : num [1:23] 14.6 90.15 73.07 8.55 29.94 ...
#> .. ..$ CI99left : num [1:23] 216.3 0 0 826.9 59.4 ...
#> .. ..$ CI99right : num [1:23] 477 1118 199 1294 460 ...
#> .. ..$ CI95left : num [1:23] 247 0 0 883 107 ...
#> .. ..$ CI95right : num [1:23] 446 931 168 1238 412 ...
#> .. ..$ CI68left : num [1:23] 296.3 34.8 18.9 970.2 182.4 ...
#> .. ..$ CI68right : num [1:23] 397 638 119 1150 337 ...
#> ..$ totest :'data.frame': 1 obs. of 14 variables:
#> .. ..$ TOTAL : num 1
#> .. ..$ rhat : num 12985
#> .. ..$ rhat.var : num 247783
#> .. ..$ NBRPLT.gt0: int 459
#> .. ..$ AREAUSED : num 62600430
#> .. ..$ rhat.se : num 498
#> .. ..$ rhat.cv : num 0.0383
#> .. ..$ pse : num 3.83
#> .. ..$ CI99left : num 11703
#> .. ..$ CI99right : num 14267
#> .. ..$ CI95left : num 12009
#> .. ..$ CI95right : num 13961
#> .. ..$ CI68left : num 12490
#> .. ..$ CI68right : num 13480
#> ..$ unit_rowest :'data.frame': 73 obs. of 17 variables:
#> .. ..$ ESTN_UNIT : num [1:73] 1 1 1 1 3 3 3 3 5 5 ...
#> .. ..$ Stand-size class: Factor w/ 4 levels "1","2","3","5": 1 2 3 4 1 2 3 4 1 3 ...
#> .. ..$ rhat : num [1:73] 649 2070 181 NA 1060 ...
#> .. ..$ rhat.var : num [1:73] 13950 1692805 322 NA 308078 ...
#> .. ..$ NBRPLT : int [1:73] 24 24 24 NA 12 NA NA NA 14 14 ...
#> .. ..$ NBRPLT.gt0 : int [1:73] 9 6 6 NA 8 NA NA NA 8 2 ...
#> .. ..$ ACRES : num [1:73] 2757613 2757613 2757613 2757613 2021729 ...
#> .. ..$ AREAUSED : num [1:73] 2757613 2757613 2757613 2757613 2021729 ...
#> .. ..$ rhat.se : num [1:73] 118.1 1301.1 17.9 NA 555 ...
#> .. ..$ rhat.cv : num [1:73] 0.1819 0.6286 0.0993 NA 0.5235 ...
#> .. ..$ pse : num [1:73] 18.19 62.86 9.93 NA 52.35 ...
#> .. ..$ CI99left : num [1:73] 345 0 134 0 0 ...
#> .. ..$ CI99right : num [1:73] 954 5421 227 NA 2490 ...
#> .. ..$ CI95left : num [1:73] 418 0 146 0 0 ...
#> .. ..$ CI95right : num [1:73] 881 4620 216 NA 2148 ...
#> .. ..$ CI68left : num [1:73] 532 776 163 0 508 ...
#> .. ..$ CI68right : num [1:73] 767 3364 199 NA 1612 ...
#> ..$ rowest :'data.frame': 4 obs. of 14 variables:
#> .. ..$ Stand-size class: Factor w/ 4 levels "1","2","3","5": 1 2 3 4
#> .. ..$ rhat : num [1:4] 17665 14011 1908 130
#> .. ..$ rhat.var : num [1:4] 588535 2122472 27032 116
#> .. ..$ NBRPLT.gt0 : int [1:4] 244 92 88 4
#> .. ..$ AREAUSED : num [1:4] 62600430 44542176 45447111 50085619
#> .. ..$ rhat.se : num [1:4] 767.2 1456.9 164.4 10.8
#> .. ..$ rhat.cv : num [1:4] 0.0434 0.104 0.0862 0.0832
#> .. ..$ pse : num [1:4] 4.34 10.4 8.62 8.32
#> .. ..$ CI99left : num [1:4] 15689 10259 1484 102
#> .. ..$ CI99right : num [1:4] 19641 17764 2331 158
#> .. ..$ CI95left : num [1:4] 16162 11156 1585 109
#> .. ..$ CI95right : num [1:4] 19169 16867 2230 151
#> .. ..$ CI68left : num [1:4] 16902 12563 1744 119
#> .. ..$ CI68right : num [1:4] 18428 15460 2071 140
#> ..$ domdat :'data.frame': 590 obs. of 19 variables:
#> .. ..$ PLT_CN : chr [1:590] "40404729010690" "40404730010690" "40404737010690" "40404738010690" ...
#> .. ..$ CONDID : num [1:590] 1 1 2 1 1 1 2 1 1 1 ...
#> .. ..$ CONDPROP_UNADJ : num [1:590] 1 1 0.25 1 1 ...
#> .. ..$ SUBPPROP_UNADJ : num [1:590] 1 1 0.25 1 1 ...
#> .. ..$ MICRPROP_UNADJ : num [1:590] 1 1 0.25 1 1 1 0.5 1 1 1 ...
#> .. ..$ MACRPROP_UNADJ : logi [1:590] NA NA NA NA NA NA ...
#> .. ..$ ESTN_UNIT : int [1:590] 1 1 1 1 1 1 1 1 1 1 ...
#> .. ..$ dem : int [1:590] 2382 2274 2576 2434 2202 2177 2011 1861 2161 2124 ...
#> .. ..$ tcc : int [1:590] 5 45 16 27 1 16 5 19 0 3 ...
#> .. ..$ tpi : num [1:590] 0.75 -12.12 2.5 -0.75 -8.62 ...
#> .. ..$ tnt2 : num [1:590] 1 0 0 0 1 1 1 0 1 1 ...
#> .. ..$ cadjfac : num [1:590] 1 1 1 1 1 ...
#> .. ..$ ADJ_FACTOR_SUBP : num [1:590] 1 1 1 1 1 ...
#> .. ..$ ADJ_FACTOR_MACR : num [1:590] 0 0 0 0 0 0 0 0 0 0 ...
#> .. ..$ ADJ_FACTOR_MICR : num [1:590] 1 1 1 1 1 ...
#> .. ..$ CONDPROP_ADJ : num [1:590] 1 1 0.25 1 1 ...
#> .. ..$ STDSZCD : int [1:590] 1 1 5 3 3 1 1 1 1 5 ...
#> .. ..$ TOTAL : num [1:590] 1 1 1 1 1 1 1 1 1 1 ...
#> .. ..$ VOLCFNET_TPA_ADJ: num [1:590] 41.8 1563 14.9 137.5 448.7 ...
#> ..$ estvarn : chr "VOLCFNET"
#> ..$ estvarn.filter: chr "STATUSCD == 1"
#> ..$ estvard : chr "CONDPROP_ADJ"
#> ..$ module : chr "MA"
#> ..$ esttype : chr "RATIO"
#> ..$ MAmethod : chr "gregRatio"
#> ..$ predselect : chr [1:3] "dem" "tcc" "tpi"
#> ..$ rowvar : chr "STDSZCD"
#> ..$ colvar : chr "NONE"
#> ..$ areaunits : chr "acres"
#> $ statecd: int 56
#> $ invyr : int [1:3] 2011 2012 2013
out$est
#> Stand-size class Estimate Percent Sampling Error
#> 1 1 17665.1 4.34
#> 2 2 14011.4 10.4
#> 3 3 1907.7 8.62
#> 4 5 129.7 8.32
#> 5 Total 12984.9 3.83 |
joshyam-k
pushed a commit
that referenced
this pull request
Mar 20, 2024
Update modGBarea to check for Forest Type
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
These changes rely on FIESTAutils#15