summary.hapLong {SimHap} | R Documentation |
Summary method for objects of class snpLong
## 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"), ...)
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. |
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'. |
Pamela A. McCaskie
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.
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)