R/print.seqhap.q

Defines functions print.seqhap

Documented in print.seqhap

#$Header: /projects/genetics/cvs/cvsroot/haplo.stats/R/print.seqhap.q,v 1.6 2011/11/22 20:29:22 sinnwell Exp $
#$Locker:  $
#$Log: print.seqhap.q,v $
#Revision 1.6  2011/11/22 20:29:22  sinnwell
#take out S exception
#
#Revision 1.5  2009/03/04 14:57:30  sinnwell
#for R use pchisq( x, p, lower=FALSE) instead of 1-pchisq
#
#Revision 1.4  2008/09/23 19:07:54  sinnwell
#use digits on signif(p, digits) for region-based p-values
#
#Revision 1.3  2007/05/25 15:38:20  sinnwell
#change inlist to scanned.loci
#
#Revision 1.2  2007/04/16 20:11:29  sinnwell
#add ...
#
#Revision 1.1  2007/04/06 19:31:34  sinnwell
#Initial revision
#
#$Author: sinnwell $
#$Date: 2011/11/22 20:29:22 $

# License: 
# 
# Copyright 2003 Mayo Foundation for Medical Education and Research. 
# 
# This program is free software; you can redistribute it and/or modify it under the terms of 
# the GNU General Public License as published by the Free Software Foundation; either 
# version 2 of the License, or (at your option) any later version.
# 
# This program is distributed in the hope that it will be useful, but WITHOUT ANY 
# WARRANTY; without even the implied warranty of MERCHANTABILITY or 
# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for 
# more details.
# 
# You should have received a copy of the GNU General Public License along with this 
# program; if not, write to the Free Software Foundation, Inc., 59 Temple Place, Suite 330, 
# Boston, MA 02111-1307 USA
# 
# For other licensing arrangements, please contact Daniel J. Schaid.
# 
# Daniel J. Schaid, Ph.D.
# Division of Biostatistics
# Harwick Building Room 775
# Mayo Clinic
# 200 First St., SW
# Rochester, MN 55905
# 
# phone: 507-284-0639
# fax:      507-284-9542
# email: schaid@mayo.edu
#

print.seqhap <- function(x, digits=max(options()$digits-2, 5), ...)

{
  
  if(x$converge==0)
    warning("EM for haplotype probabilities failed to converge")
  else
    {
      printBanner("Single-locus Chi-square Test")
      cat("Regional permuted P-value based on single-locus test is ", 
          signif(x$chi.p.region, digits),"\n")
      chi.test <- data.frame(chi.stat=round(x$chi.stat,digits), 
                     perm.point.p=x$chi.p.point, 
                     asym.point.p= round(pchisq(x$chi.stat,1,
                       lower.tail=FALSE),digits))
      row.names(chi.test) <- c(x$locus.label)
      print(chi.test)
      cat("\n\n")
      
      printBanner("Sequential Scan")     
      cat("Loci Combined in Sequential Analysis\n")
      scanned.loci <- x$scanned.loci
      for(i in 1:dim(scanned.loci)[1]){
        cat('seq-')
        cat(x$locus.label[i], scanned.loci[i,scanned.loci[i,]!=0],
            sep=" ","\n")
      }
      cat("\n\n")

      printBanner("Sequential Haplotype Test")
      cat("Regional permuted P-value based on sequential haplotype test is ", 
          signif(x$hap.p.region, digits),"\n")
      hap.test <- data.frame(hap.stat=round(x$hap.stat,digits),
                    df=x$hap.df,
                    perm.point.p=x$hap.p.point,
                    asym.point.p= round(pchisq(x$hap.stat,x$hap.df,
                      lower.tail=FALSE),digits))
      row.names(hap.test) <- paste(rep('seq-',length(x$locus.label)),
                                   x$locus.label,sep='')
      print(hap.test)
      cat("\n\n")
      
      printBanner("Sequential Sum Test")
      cat("Regional permuted P-value based on sequential sum test is ", 
            signif(x$sum.p.region, digits),"\n")
      sum.test <- data.frame(sum.stat=round(x$sum.stat,digits),
                    df=x$sum.df,
                    perm.point.p=x$sum.p.point,
                    asym.point.p = round(pchisq(x$sum.stat,x$sum.df,
                      lower.tail=FALSE),digits))
      row.names(sum.test) <- paste(rep('seq-',length(x$locus.label)),
                                   x$locus.label,sep='')
      print(sum.test)
    }

  invisible()
}

Try the haplo.stats package in your browser

Any scripts or data that you put into this service are public.

haplo.stats documentation built on May 29, 2024, 9:53 a.m.