update PKPDeleveld

src/PKPDeleveld.jl

@@ -127,10 +127,10 @@ ts = 0.0:Ts:3600 # s
 ## PD model
 ke0 = ϕ[2]/60 # [1/s]
 Ce50 = ϕ[1]*f_ageing(ϕ[7])*1000 # [μg/l]
-BIS_baseline = ϕ[3] # [BIS index]
+BIS_baseline = 100 # [BIS index]
 γ = ϕ[4]
 #BIS(Ce) = BIS_baseline*(1 - (Ce^γ)/(Ce50^γ + Ce^γ)) + ϕ[5]  # This residual error gives 93+8 = 101 as start?
-BIS(Ce) = 100*(1 - (Ce^γ)/(Ce50^γ + Ce^γ))
+BIS(Ce) = BIS_baseline*(1 - (Ce^γ)/(Ce50^γ + Ce^γ))
 
 ## PKPD model
 PD_lc = ss(tf(ke0,[1, ke0]))