Interlude II (simulations) [RSABE / ABEL]

posted by PharmCat  – Russia, 2020-08-19 22:26 (160 d 14:03 ago) – Posting: # 21890
Views: 2,136

Hi Helmut!

using DataFrames, CSV, ReplicateBE, LinearAlgebra

path      = dirname(@__FILE__)
dfa       = CSV.File(path*"/sim500-2x3x3.csv") |> DataFrame
dfa.logpk = log.(dfa.PK)

res = Vector{Any}(undef, 500)
pd  = Vector{Any}(undef, 500)
os  = Vector{Any}(undef, 500)
for i in 1:500
    df = filter(r -> i == r.set , dfa)
    res[i] = ReplicateBE.rbe!(df, dvar = :logpk, subject = :subject, formulation = :treatment, period = :period, sequence = :sequence, g_tol = 1e-10, singlim = 1e-12)
    pd[i]  = isposdef(Symmetric(res[i].result.H))
    os[i]  = ReplicateBE.optstat(res[i])
    println(i)
end
println("Positive definite: ", sum(pd), "(",sum(pd)/500*100,"%)")
println("Converged: ", sum(os), "(",sum(os)/500*100,"%)")


Positive definite: 332(66.4%)
Converged: 500(100.0%)

I can make table with results if necessary.

Complete thread:

Activity
 Admin contact
21,316 posts in 4,446 threads, 1,489 registered users;
online 6 (0 registered, 6 guests [including 2 identified bots]).
Forum time: Wednesday 11:30 CET (Europe/Vienna)

Nothing fails like success because you do not learn anything from it.
The only thing we ever learn from is failure.
Success only confirms our superstitions.    Kenneth E. Boulding

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5