summary.hapLong {SimHap}R Documentation

Summarizing single SNP analysis models for longitudinal data

Description

Summary method for objects of class snpLong

Usage

## S3 method for class 'hapLong':
summary(object, ...)
## S3 method for class 'summary.hapLong':
print(x, digits = max(3, getOption("digits") - 3), 
        signif.stars = getOption("show.signif.stars"), ...)

Arguments

object an object of class hapLong, a result of a call to haplo.long.
x an object of class summary.hapLong, the result of a call to
summary.hapLong.
digits the number of significant digits to use when printing.
signif.stars logical. If TRUE, ``significance stars" are printed for each coefficient.
... further arguments passed to or from other methods.

Value

summary.hapLong returns an object of class summary.hapLong. A list with components

fixed_formula fixed effects formula used in haplo.long.
random_formula random effects formula used in haplo.long.
coefficients summarized results from fitted model, including coefficients, standard errors and p-values.
empiricalResults a list containing the coefficients, standard errors and p-values calculated at each simulation of haplo.long.
AIC Akaike Information Criterion for the model fitted in haplo.long.
corStruct the correlation structure specified in haplo.long.
effect the haplotypic effect modelled: `ADDITIVE', `DOMINANT' or `RECESSIVE'.

Author(s)

Pamela A. McCaskie

References

McCaskie, P.A., Carter, K.W. Hazelton, M., Palmer, L.J. (2007) SimHap: A comprehensive modeling framework for epidemiological outcomes and a multiple-imputation approach to haplotypic analysis of population-based data, [online] www.genepi.org.au/simhap.

See Also

haplo.long

Examples


data(SNPlong.dat)

# convert SNP.dat to format required by infer.haplos
longHaplo.dat <- SNP2Haplo(SNPlong.dat)

data(longPheno.dat)

myinfer<-infer.haplos(longHaplo.dat) 
myinfer$hap.freq 

myhaplo<-make.haplo.rare(myinfer,min.freq=0.05)
mymodel <- haplo.long(fixed=fev1f~h.ACV2, random=~1|id, 
        pheno=longPheno.dat, haplo=myhaplo, form=~year|id, sim=10)
summary(mymodel)


[Package SimHap version 1.0.0 Index]