logLik.hapassoc | R Documentation |
This function is used to return the log-likelihood at the maximum
likelihood estimates computed by hapassoc
and to return
the number of parameters fit by hapassoc
(i.e. the degrees of freedom in R
) for cohort or
cross-sectional data.
## S3 method for class 'hapassoc' logLik(object, ...)
object |
a list of class |
... |
additional arguments to the summary function (currently unused) |
See the hapassoc vignette, of the same name as the package, for details.
logLik |
log-likelihood computed at the maximum likelihood estimates
if |
df |
number of parameters in the model (i.e. regression coefficients, any dispersion parameters and haplotype frequencies). This is not the residual degrees of freedom, which is the number of subjects minus the number of parameters estimated. |
Burkett K, McNeney B, Graham J (2004). A note on inference of trait associations with SNP haplotypes and other attributes in generalized linear models. Human Heredity, 57:200-206
Burkett K, Graham J and McNeney B (2006). hapassoc: Software for Likelihood Inference of Trait Associations with SNP Haplotypes and Other Attributes. Journal of Statistical Software, 16(2):1-19
pre.hapassoc
,hapassoc
,
summary.hapassoc
.
data(hypoDatGeno) example2.pre.hapassoc<-pre.hapassoc(hypoDatGeno, numSNPs=3, allelic=FALSE) example.regr <- hapassoc(affected ~ attr + hAAA+ hACA + hACC + hCAA + pooled, example2.pre.hapassoc, family=binomial()) logLik(example.regr) # Returns: # Log Lik: -322.1558 (df=14)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.