Do you really need ANOVA type III yet? No, you can’t do it with pure Julia, but you can do it with SASLMR.jl and R Project. SASLMR.jl wrap some sasLM functions (aov1, aov2, aov3).
Install:
import Pkg; Pkg.add(url="https://github.com/PharmCat/SASLMR.jl.git")
and use:
using SASLMR
using StatsModels, StatsBase, DataFrames
# data (DataFrame)
bedf = SASLMR.bedata()
# run ANOVA
beaov = SASLMR.aov(@formula(CMAX ~ PRD + TRT + SEQ + SUBJ&SEQ), bedf; beta=true, resid=true, type = "III")
# get DF
beaovdf = DataFrame(beaov.ct)
# or get coeftable
ct = coeftable(beaov)
Also functions: aov1
, aov2
, aov3
can be used in the same way:
beaov = SASLMR.aov3(@formula(CMAX ~ PRD + TRT + SEQ + SUBJ&SEQ), bedf; beta=true, resid=true)
AOVSumm
fields:
-
heading
-
ct
-
parameters
-
fitted
-
residual
For more details see: https://cran.r-project.org/web/packages/sasLM/sasLM.pdf
SASLMR.jl not cover all sasLM
functionality… so… wellcome any PR for extending.