adding arguments for icens and cens.method (to be passed through to autonpde()) and something like blq_col to identify BLQ flag column
include BLQ records in simulations: this would mean augmenting the analysis dataset if they're ignored in the model. These would be excluded by nm_join() with .superset = FALSE
I believe this would entail:
icensandcens.method(to be passed through toautonpde()) and something likeblq_colto identify BLQ flag columnnm_join()with.superset = FALSE