hapassoc (version 0.2)

summary.EM: Summarize results of the EM function

Description

Summary function for reporting the results of the EM function in a similar style to the lm and glm summaries.

Usage

summary.EM(object, ...)

Arguments

object

a list of class EM

...

additional arguments to the summary function (currently unused)

Value

• coefficients
  Table of estimated coefficients, standard errors and Wald tests for each variable
• frequencies
  Table of estimated haplotype frequencies and standard errors
• dispersion
  Estimate of dispersion parameter (Moment estimator for gamma model)

See Also

PreEM,EM.

Examples

data(hypoDat)
example.preEM<-PreEM(hypoDat, 3)
example.regr <- EM(affected ~ attr + h000+ h010 + h011 + h100 + pooled,
                     example.preEM, family=binomial())

# Summarize the results:
summary.EM(example.regr) # or just summary(example.regr)

# Results:
#$coefficients
#               Estimate Std. Error      zscore   Pr(>|z|)
#(Intercept) -1.24114270  0.7820977 -1.58694079 0.11252606
#attr         0.74036920  0.2918205  2.53707057 0.01117844
#h000         1.14968352  0.5942542  1.93466627 0.05303126
#h010        -0.59318434  0.6569672 -0.90291311 0.36657201
#h011        -0.03615243  0.9161959 -0.03945928 0.96852422
#h100        -0.85329292  1.0203105 -0.83630709 0.40298217
#pooled       0.38516864  0.8784283  0.43847478 0.66104215
#
#$frequencies
#         Estimate Std. Error
#f.h000 0.26716394 0.03933158
#f.h001 0.25191674 0.03866739
#f.h010 0.21997138 0.03881578
#f.h011 0.10094795 0.02949617
#f.h100 0.09507014 0.02371878
#f.h101 0.02584918 0.01411881
#f.h110 0.01779455 0.01386080
#f.h111 0.02128613 0.01247265
#
#$dispersion
#[1] 1