data(lucaDat)
pen.model <- formula(d ~ I(allele.count(g, "C")) +
a + a2 + I(allele.count(g, "C")):a2)
fitDep <- luca(pen.model = pen.model, gLabel = "g",
dat = lucaDat, dep.model = formula(g ~ a))
# Summarize the results:
summary.luca(fitDep) # or just summary(fitDep)
#Returns:
#Call:
#luca(dat = lucaDat, pen.model = pen.model, gLabel = "g", dep.model =
#formula(g ~ a))
#
#Coefficients:
# Estimate Std. Error zscore Pr(>|z|)
#I(allele.count(g, "C")) 0.61738385 0.10820323 5.7057800 1.158115e-08
#a 0.11629696 0.07815014 1.4881222 1.367187e-01
#a2 -0.03087368 0.10787965 -0.2861863 7.747354e-01
#I(allele.count(g, "C")):a2 0.31879401 0.08236130 3.8706772 1.085334e-04
Run the code above in your browser using DataLab