search

TNC-NMFO / NWLAND

carbon accounting model
0 stars 0 forks source link

asRan | caland | missing values are not allowed in subscripted assignments of data frames #118

Open sbassett opened 2 years ago

sbassett commented 2 years ago 
Start CALAND at Fri Jan 14 12:10:12 2022 

Starting year  2021 ...
 Error in `[<-.data.frame`(`*tmp*`, man_dead_df$Management == "Afforestation" |  : 
  missing values are not allowed in subscripted assignments of data frames 
4.
stop("missing values are not allowed in subscripted assignments of data frames") 
3.
`[<-.data.frame`(`*tmp*`, man_dead_df$Management == "Afforestation" | 
    man_dead_df$Management == "Reforestation" | man_dead_df$Management == 
    "Restoration" | man_dead_df$Management == "Prescribed_burn" | 
    man_dead_df$Management == "Prescribed_burn_med_slash_util" |  ... at CALAND.r#1754
2.
`[<-`(`*tmp*`, man_dead_df$Management == "Afforestation" | man_dead_df$Management == 
    "Reforestation" | man_dead_df$Management == "Restoration" | 
    man_dead_df$Management == "Prescribed_burn" | man_dead_df$Management == 
    "Prescribed_burn_med_slash_util" | man_dead_df$Management ==  ... at CALAND.r#1754
1.
CALAND(scen_file_arg = "CONM_BAU_v01_testA_rcp85_climNA.xls", 
    c_file_arg = "2022_01_13_input_conm_testA_rcp85_noOcean.xls", 
    indir = "2022_01_13_input_conm_testA_rcp85_noOcean", outdir = "2022_01_13_output_conm_testA_noOcean", 
    start_year = 2021, end_year = 2051, value_col_dens = 7, ADD_dens = TRUE,  ...
sbassett commented 2 years ago 
> head(man_dead_df$deadcfracXarea)
  2021_frac 2021_frac.1
1        NA          NA
2        NA          NA
3        NA          NA
4        NA          NA
5        NA          NA
6        NA          NA
> head(man_dead_df)
  Land_Cat_ID Region Land_Type   Ownership Management man_area man_area_sum tot_area man_area_agg_pre excess_area_pre man_area_agg
1      100002 C08001     Water         DOD       <NA>       NA           NA       NA               NA              NA           NA
2      100003 C08001     Water         FWS       <NA>       NA           NA       NA               NA              NA           NA
3      100004 C08001     Water       Local       <NA>       NA           NA       NA               NA              NA           NA
4      100006 C08001     Water     Private       <NA>       NA           NA       NA               NA              NA           NA
5      100007 C08001     Water State_Other       <NA>       NA           NA       NA               NA              NA           NA
6      100008 C08001     Water State_Trust       <NA>       NA           NA       NA               NA              NA           NA
  man_area_sum_agg_extra excess_sum_area man_area_sum_agg SoilCaccum_frac Above_harvested_frac StandDead_harvested_frac
1                     NA              NA               NA              NA                   NA                       NA
2                     NA              NA               NA              NA                   NA                       NA
3                     NA              NA               NA              NA                   NA                       NA
4                     NA              NA               NA              NA                   NA                       NA
5                     NA              NA               NA              NA                   NA                       NA
6                     NA              NA               NA              NA                   NA                       NA
  Harvested2Wood_frac Harvested2Energy_frac Harvested2SawmillDecay_frac Harvested2Slash_frac Under2Slash_frac DownDead2Slash_frac
1                  NA                    NA                          NA                   NA               NA                  NA
2                  NA                    NA                          NA                   NA               NA                  NA
3                  NA                    NA                          NA                   NA               NA                  NA
4                  NA                    NA                          NA                   NA               NA                  NA
5                  NA                    NA                          NA                   NA               NA                  NA
6                  NA                    NA                          NA                   NA               NA                  NA
  Litter2Slash_frac Slash2Energy_frac Slash2Wood_frac Slash2Burn_frac Slash2Decay_frac Under2DownDead_frac Soil2Atmos_frac
1                NA                NA              NA              NA               NA                  NA              NA
2                NA                NA              NA              NA               NA                  NA              NA
3                NA                NA              NA              NA               NA                  NA              NA
4                NA                NA              NA              NA               NA                  NA              NA
5                NA                NA              NA              NA               NA                  NA              NA
6                NA                NA              NA              NA               NA                  NA              NA
  Above2StandDead_frac Below2Atmos_frac Below2Soil_frac VegCuptake_frac DeadCaccum_frac high_sev_frac med_sev_frac low_sev_frac
1                   NA               NA              NA              NA              NA            NA           NA           NA
2                   NA               NA              NA              NA              NA            NA           NA           NA
3                   NA               NA              NA              NA              NA            NA           NA           NA
4                   NA               NA              NA              NA              NA            NA           NA           NA
5                   NA               NA              NA              NA              NA            NA           NA           NA
6                   NA               NA              NA              NA              NA            NA           NA           NA
  start_urban_forest_fraction current_urban_forest_fraction deadc_frac_in.2021_frac deadc_frac_in.2021_frac.1 adjDeadCfrac.2021_frac
1                          NA                            NA                       0                         0                     NA
2                          NA                            NA                       0                         0                     NA
3                          NA                            NA                       0                         0                     NA
4                          NA                            NA                       0                         0                     NA
5                          NA                            NA                       0                         0                     NA
6                          NA                            NA                       0                         0                     NA
  adjDeadCfrac.2021_frac.1 deadcfracXarea.2021_frac deadcfracXarea.2021_frac.1
1                       NA                       NA                         NA
2                       NA                       NA                         NA
3                       NA                       NA                         NA
4                       NA                       NA                         NA
5                       NA                       NA                         NA
6                       NA                       NA                         NA
sbassett commented 2 years ago

selecting undefined columns ?

> man_dead_df$deadcfracXarea[man_dead_df$Management == "Afforestation" | man_dead_df$Management == "Reforestation" | man_dead_df$Management == "Restoration" | 
+                                man_dead_df$Management == "Prescribed_burn" | man_dead_df$Management == "Prescribed_burn_med_slash_util" | man_dead_df$Management == "Prescribed_burn_hi_slash_util"]
Error in `[.data.frame`(man_dead_df$deadcfracXarea, man_dead_df$Management ==  : 
  undefined columns selected
sbassett commented 2 years ago

two values for 2021 frac becasue of increased mortality? mortality table has two values for 2021 (baseline year, and first year of increased mortality period)

will try getting rid of the increased mortality period

this worked!

sbassett commented 2 years ago

to avoid this, have made the increased mortality exclude the start year, and made it a very small increase
forest_mort_fact = 1.01, #tiny increase to test hypothesis that this is mucking up mortality code later forest_mort_adj_first = 2023, forest_mort_adj_last = 2024,