R Under development (unstable) (2025-09-04 r88794 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> ## Purpose: Test script for pedgene package
> ## Testing a simple 10-variant gene on both the X chrom and autosome
> ## Authors: Dan Schaid and Jason Sinnwell
> ## Created: 19-AUG-2013
> ## Updated: 11/7/2013
> 
> require(CompQuadForm)
Loading required package: CompQuadForm
> 
> 
> davies(144, c(20,10, 5, 1, .4, .1), acc=1e-4)
$trace
[1]   5.08345178 128.00000000   1.00000000   0.01675413   2.13515773
[6]   1.03615893  30.00000000

$ifault
[1] 0

$Qq
[1] 0.01383433

> 
> davies(144, c(20,10, 5, 1, .4, .1), acc=1e-6)
$trace
[1] 5.365851e+00 1.291000e+03 1.000000e+00 1.103389e-02 1.423438e+01
[6] 1.036159e-01 3.300000e+01

$ifault
[1] 0

$Qq
[1] 0.01382884

> 
> davies(144, c(15,8, 5, 1, .4, .1), acc=1e-6)
$trace
[1]   6.24993456 937.00000000   1.00000000   0.01651383  15.44929251
[6]   0.10524124  33.00000000

$ifault
[1] 0

$Qq
[1] 0.004081417

> 
> 
> 
> proc.time()
   user  system elapsed 
   0.15    0.01    0.15