PharmCat
PharmCat
> If you want them removed from the table, then you should remove the columns from the data before fitting. Hi! I can't remove data from a table. Also if...
> I think we're talking past each other: > > * Your F-test MWE works for me and does not error. Sorry, I really don't understand: For this models I...
> versioninfo() ``` julia> versioninfo() Julia Version 1.7.3 Commit 742b9abb4d (2022-05-06 12:58 UTC) Platform Info: OS: Windows (x86_64-w64-mingw32) CPU: AMD Ryzen 9 5950X 16-Core Processor WORD_SIZE: 64 LIBM: libopenlibm LLVM:...
> Note that nesting and rank deficiency is one of those problems that is well-defined in theory, but not always in practice due to the vagaries of floating point. Yes,...
> Note that nesting and rank deficiency is one of those problems that is well-defined in theory, but not always in practice due to the vagaries of floating point. May...
> dropcollinear=false In this case results is definitly wrong: ``` olscmax_noseq = lm(@formula(log(Cmax) ~ Period+Formulation+Subject), df; dropcollinear=false) olscmax = lm(@formula(log(Cmax) ~ Sequence+Period+Formulation+Subject), df; dropcollinear=false) ftest(olscmax.model, olscmax_noseq.model) ``` ``` F-test: 2...
> dropcollinear=false ``` using CategoricalArrays, DataFrames, CSV, GLM df = download("https://raw.githubusercontent.com/PharmCat/edu/main/csv/be.csv") |> CSV.File |> DataFrame transform!(df, :Subject => categorical, renamecols=false) transform!(df, :Period => categorical, renamecols=false) transform!(df, :Sequence => categorical, renamecols=false)...
> Now when you consider the a likelihood-ratio or F-test, how does this impact the (denominator) degrees of freedom for the comparison? I think for current case DOF will be...
> Thank you so much! Also I try to represent vector of units as the values. Is this approach is correct? ``` uvec = ones(typeof(1u"mg"), 10) T = typeof(1u"mg") vec...
> It will work if `uvec` is constructed like this, but for a general vector it might not. Basically, for this to work, the `eltype` of `uvec` must be a...