set_morris.r

# parameter settings for Morris method
testpars <- c("V_plasma",
            "V_inter",
            "V_muscle",
            "kgut",
            "Km",
            "Vmax",
            "m_K_ALDO",
            "P_ECF",
            "FF",
            "GFR",
            "dtKsec_eq",
            "A_dtKsec",
            "B_dtKsec",
            "cdKsec_eq",
            "A_cdKsec",
            "B_cdKsec",
            "A_cdKreab",
            "A_insulin",
            "B_insulin",
            "KMuscleBase",
            "Kecf_base",
            "ALD_eq",
            "etapsKreab"
            )

parsbinf <- c(0.50 * p$V_plasma, # V_plasma
            0.50 * p$V_inter, # V_inter
            0.50 * p$V_muscle, # V_muscle
            0.50 * p$kgut, # kgut
            0.8, # Km, Cheng gave
            0.50 * p$Vmax, # Vmax
            0.50 * p$m_K_ALDO, # m_K_ALDO
            0.50 * p$P_ECF, # P_ECF
            0.50 * p$FF, # FF
            100 / 1000, # GFR, normal healthy range
            0.75 * p$dtKsec_eq, # dtKsec_eq
            0.50 * p$A_dtKsec, # A_dtKsec
            0.50 * p$B_dtKsec, # B_dtKsec
            0.75 * p$cdKsec_eq, # cdKsec_eq
            0.50 * p$A_cdKsec, # A_cdKsec
            0.50 * p$B_cdKsec, # B_cdKsec
            0.75 * p$A_cdKreab, # A_cdKreab
            0.50 * p$A_insulin, # A_insulin
            0.50 * p$B_insulin, # B_insulin
            120, # KMuscleBase, low end of regular K
            3.5, # Kecf_base, low normal range
            31, #ALD_eq, low normal ALD range
            0.85 # etapsKreab
            )

parsbsup <- c(1.50 * p$V_plasma, # V_plasma
            1.50 * p$V_inter, # V_inter
            1.50 * p$V_muscle, # V_muscle
            1.50 * p$kgut, # kgut
            1.5, # Km, Cheng gave
            1.50 * p$Vmax, # Vmax
            1.50 * p$m_K_ALDO, # m_K_ALDO
            1.50 * p$P_ECF, # P_ECF
            1.50 * p$FF, # FF
            130 / 1000, # GFR, normal healthy range
            1.25 * p$dtKsec_eq, # dtKsec_eq
            1.50 * p$A_dtKsec, # A_dtKsec
            1.50 * p$B_dtKsec, # B_dtKsec
            1.25 * p$cdKsec_eq, # cdKsec_eq
            1.50 * p$A_cdKsec, # A_cdKsec
            1.50 * p$B_cdKsec, # B_cdKsec
            1.25 * p$A_cdKreab, # A_cdKreab
            1.50 * p$A_insulin, # A_insulin
            1.50 * p$B_insulin, # B_insulin
            140, # KMuscleBase, high normal intracellular
            5.0, # Kecf_base, high normal plasma K
            350, # ALD_eq
            0.95 # etapsKreab
            )