Interlude II (simulations) [RSABE / ABEL]

posted by PharmCat  – Russia, 2020-08-20 00:26 (1785 d 15:59 ago) – Posting: # 21890
Views: 8,443

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:

UA Flag
Activity
 Admin contact
23,428 posts in 4,929 threads, 1,682 registered users;
60 visitors (0 registered, 60 guests [including 11 identified bots]).
Forum time: 16:26 CEST (Europe/Vienna)

To know that we know what we know,
and to know that we do not know what we do not know,
that is true knowledge.    Nicolaus Copernicus

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