Re: Log of non positive number
Posted: Mon Jul 06, 2015 12:15 pm
Thank you for your question. Here is the model:
:eq_rpcetotto_jpn
:eq_ratepubtot_jpn
:eq_rateopprv_jpn
@IDENTITY rpcepubto_jpn = ratepubtot_jpn * rpcetotto_jpn
@IDENTITY rpceprvto_jpn = (1 - ratepubtot_jpn) * rpcetotto_jpn
@IDENTITY rateopprv_jpn_2 = @recode(rateopprv_jpn + rateinprv_jpn >.99 , (.99 / (rateopprv_jpn + rateinprv_jpn)) * rateopprv_jpn , rateopprv_jpn)
@IDENTITY rateinprv_jpn_2 = @recode(rateopprv_jpn + rateinprv_jpn >.99 , (.99 / (rateopprv_jpn + rateinprv_jpn)) * rateinprv_jpn , rateinprv_jpn)
@IDENTITY rpceprvop_jpn = rateopprv_jpn_2 * rpceprvto_jpn
@IDENTITY rpceprvin_jpn = rateinprv_jpn_2 * rpceprvto_jpn
@IDENTITY rpceprvng_jpn = rpceprvto_jpn - (rpceprvop_jpn + rpceprvin_jpn)
:eq_rateinprv_jpn
:eq_rpcetotinv_jpn
@IDENTITY rpcetotcur_jpn = rpcetotto_jpn - rpcetotinv_jpn
:eq_rpcetotper_jpn
@IDENTITY rpcetotcol_jpn = rpcetotcur_jpn - rpcetotper_jpn
:eq_rpcetotperda_jpn
:eq_rpcetotperho_jpn
:eq_rpcetotperinp_jpn
:eq_rpcetotperout_jpn
:eq_rpcetotpermg_jpn
@IDENTITY rpcetotperan_jpn = rpcetotper_jpn - (rpcetotperda_jpn + rpcetotperho_jpn + rpcetotperinp_jpn + rpcetotpermg_jpn)
:eq_rpcetotpermgdu_jpn
@IDENTITY rpcetotpermgnd_jpn = rpcetotpermg_jpn - rpcetotpermgdu_jpn
:eq_rpcetotto_jpn
:eq_ratepubtot_jpn
:eq_rateopprv_jpn
@IDENTITY rpcepubto_jpn = ratepubtot_jpn * rpcetotto_jpn
@IDENTITY rpceprvto_jpn = (1 - ratepubtot_jpn) * rpcetotto_jpn
@IDENTITY rateopprv_jpn_2 = @recode(rateopprv_jpn + rateinprv_jpn >.99 , (.99 / (rateopprv_jpn + rateinprv_jpn)) * rateopprv_jpn , rateopprv_jpn)
@IDENTITY rateinprv_jpn_2 = @recode(rateopprv_jpn + rateinprv_jpn >.99 , (.99 / (rateopprv_jpn + rateinprv_jpn)) * rateinprv_jpn , rateinprv_jpn)
@IDENTITY rpceprvop_jpn = rateopprv_jpn_2 * rpceprvto_jpn
@IDENTITY rpceprvin_jpn = rateinprv_jpn_2 * rpceprvto_jpn
@IDENTITY rpceprvng_jpn = rpceprvto_jpn - (rpceprvop_jpn + rpceprvin_jpn)
:eq_rateinprv_jpn
:eq_rpcetotinv_jpn
@IDENTITY rpcetotcur_jpn = rpcetotto_jpn - rpcetotinv_jpn
:eq_rpcetotper_jpn
@IDENTITY rpcetotcol_jpn = rpcetotcur_jpn - rpcetotper_jpn
:eq_rpcetotperda_jpn
:eq_rpcetotperho_jpn
:eq_rpcetotperinp_jpn
:eq_rpcetotperout_jpn
:eq_rpcetotpermg_jpn
@IDENTITY rpcetotperan_jpn = rpcetotper_jpn - (rpcetotperda_jpn + rpcetotperho_jpn + rpcetotperinp_jpn + rpcetotpermg_jpn)
:eq_rpcetotpermgdu_jpn
@IDENTITY rpcetotpermgnd_jpn = rpcetotpermg_jpn - rpcetotpermgdu_jpn