Last data update: 2014.03.03

R: Matrix version of HWChisq
HWChisqMatR Documentation

Matrix version of HWChisq

Description

HWChisqMat executes the Chisquare test for HWE for each row in a matrix.

Usage

HWChisqMat(X, ...)

Arguments

X

A n times 3 matrix of genotypic counts (AA,AB,BB)

...

extra arguments that are passed on to HWChisq

Value

pvalvec

Vector with the p-values of each test

chisqvec

Vector with the chi-square statistics

Dvec

Vector with deviations from independence

Author(s)

Jan Graffelman jan.graffelman@upc.edu

See Also

HWChisq

Examples

X <- HWData(100,10)
colnames(X) <- c("MM","MN","NN")
Results <- HWChisqMat(X)
Output <- cbind(X,Results$chisqvec,Results$pvalvec)
print(Output)

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

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.

> library(HardyWeinberg)
Loading required package: mice
Loading required package: Rcpp
mice 2.25 2015-11-09
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/HardyWeinberg/HWChisqMat.Rd_%03d_medium.png", width=480, height=480)
> ### Name: HWChisqMat
> ### Title: Matrix version of HWChisq
> ### Aliases: HWChisqMat
> ### Keywords: htest
> 
> ### ** Examples
> 
> X <- HWData(100,10)
> colnames(X) <- c("MM","MN","NN")
> Results <- HWChisqMat(X)
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.610754 DF =  1 p-value =  0.4345041 D =  -0.975 f =  0.3939394 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2318594 DF =  1 p-value =  0.6301485 D =  -0.1 f =  0.04761905 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.05 DF =  1 p-value =  0.8230633 D =  0.5 f =  -0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.04444444 DF =  1 p-value =  0.8330289 D =  -0.375 f =  0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.610754 DF =  1 p-value =  0.4345041 D =  -0.975 f =  0.3939394 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2318594 DF =  1 p-value =  0.6301485 D =  -0.1 f =  0.04761905 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.161565 DF =  1 p-value =  0.2811413 D =  -1.1 f =  0.5238095 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.05 DF =  1 p-value =  0.8230633 D =  0.5 f =  -0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.6128118 DF =  1 p-value =  0.4337311 D =  0.9 f =  -0.4285714 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.02777778 DF =  1 p-value =  0.8676323 D =  -0.4 f =  0.1666667 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.161565 DF =  1 p-value =  0.2811413 D =  -1.1 f =  0.5238095 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1777778 DF =  1 p-value =  0.67329 D =  0.625 f =  -0.3333333 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2318594 DF =  1 p-value =  0.6301485 D =  -0.1 f =  0.04761905 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1111111 DF =  1 p-value =  0.7388827 D =  0.6 f =  -0.25 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1796875 DF =  1 p-value =  0.6716419 D =  -0.6 f =  0.375 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.04444444 DF =  1 p-value =  0.8330289 D =  -0.375 f =  0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  4.463218 DF =  1 p-value =  0.0346322 D =  2.025 f =  -0.8181818 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.833333 DF =  1 p-value =  0.1757343 D =  -1.4 f =  0.5833333 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3219418 DF =  1 p-value =  0.570443 D =  -0.775 f =  0.3406593 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.7125803 DF =  1 p-value =  0.3985885 D =  1.025 f =  -0.4141414 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.6128118 DF =  1 p-value =  0.4337311 D =  0.9 f =  -0.4285714 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2269156 DF =  1 p-value =  0.6338203 D =  0.025 f =  -0.01010101 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.05 DF =  1 p-value =  0.8230633 D =  0.5 f =  -0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1796875 DF =  1 p-value =  0.6716419 D =  -0.6 f =  0.375 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  9.263158 DF =  1 p-value =  0.002338095 D =  0.025 f =  -0.05263158 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.08018355 DF =  1 p-value =  0.7770488 D =  0.225 f =  -0.0989011 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.0546875 DF =  1 p-value =  0.8150986 D =  0.4 f =  -0.25 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1111111 DF =  1 p-value =  0.7388827 D =  0.6 f =  -0.25 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.04444444 DF =  1 p-value =  0.8330289 D =  -0.375 f =  0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.05 DF =  1 p-value =  0.8230633 D =  0.5 f =  -0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.08018355 DF =  1 p-value =  0.7770488 D =  0.225 f =  -0.0989011 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.08018355 DF =  1 p-value =  0.7770488 D =  0.225 f =  -0.0989011 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.6128118 DF =  1 p-value =  0.4337311 D =  0.9 f =  -0.4285714 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.161565 DF =  1 p-value =  0.2811413 D =  -1.1 f =  0.5238095 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3219418 DF =  1 p-value =  0.570443 D =  -0.775 f =  0.3406593 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1796875 DF =  1 p-value =  0.6716419 D =  -0.6 f =  0.375 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1111111 DF =  1 p-value =  0.7388827 D =  0.6 f =  -0.25 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  2.611111 DF =  1 p-value =  0.1061174 D =  1.6 f =  -0.6666667 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.02777778 DF =  1 p-value =  0.8676323 D =  -0.4 f =  0.1666667 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.08018355 DF =  1 p-value =  0.7770488 D =  0.225 f =  -0.0989011 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.0546875 DF =  1 p-value =  0.8150986 D =  0.4 f =  -0.25 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1796875 DF =  1 p-value =  0.6716419 D =  -0.6 f =  0.375 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3219418 DF =  1 p-value =  0.570443 D =  -0.775 f =  0.3406593 
warning: monomorphic marker
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  NA DF =  1 p-value =  1 D =  0 f =  0 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.05 DF =  1 p-value =  0.8230633 D =  -0.5 f =  0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.6128118 DF =  1 p-value =  0.4337311 D =  0.9 f =  -0.4285714 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2269156 DF =  1 p-value =  0.6338203 D =  0.025 f =  -0.01010101 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1796875 DF =  1 p-value =  0.6716419 D =  -0.6 f =  0.375 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1777778 DF =  1 p-value =  0.67329 D =  0.625 f =  -0.3333333 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3219418 DF =  1 p-value =  0.570443 D =  -0.775 f =  0.3406593 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.04444444 DF =  1 p-value =  0.8330289 D =  -0.375 f =  0.2 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1796875 DF =  1 p-value =  0.6716419 D =  -0.6 f =  0.375 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2318594 DF =  1 p-value =  0.6301485 D =  -0.1 f =  0.04761905 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  9.263158 DF =  1 p-value =  0.002338095 D =  0.025 f =  -0.05263158 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  9.263158 DF =  1 p-value =  0.002338095 D =  0.025 f =  -0.05263158 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3219418 DF =  1 p-value =  0.570443 D =  -0.775 f =  0.3406593 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.389204 DF =  1 p-value =  0.2385395 D =  1.225 f =  -0.5384615 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.3475586 DF =  1 p-value =  0.5554984 D =  0.225 f =  -0.1764706 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.7789312 DF =  1 p-value =  0.3774682 D =  -0.775 f =  0.6078431 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2269156 DF =  1 p-value =  0.6338203 D =  0.025 f =  -0.01010101 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.669753 DF =  1 p-value =  0.1962916 D =  0.1 f =  -0.1111111 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.08018355 DF =  1 p-value =  0.7770488 D =  0.225 f =  -0.0989011 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.02777778 DF =  1 p-value =  0.8676323 D =  -0.4 f =  0.1666667 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.1111111 DF =  1 p-value =  0.7388827 D =  0.6 f =  -0.25 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  9.263158 DF =  1 p-value =  0.002338095 D =  0.025 f =  -0.05263158 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  0.2318594 DF =  1 p-value =  0.6301485 D =  -0.1 f =  0.04761905 
Chi-square test with continuity correction for Hardy-Weinberg equilibrium (autosomal)
Chi2 =  1.389204 DF =  1 p-value =  0.2385395 D =  1.225 f =  -0.5384615 
There were 50 or more warnings (use warnings() to see the first 50)
> Output <- cbind(X,Results$chisqvec,Results$pvalvec)
> print(Output)
       MM MN NN                       
  [1,]  8  2  0 1.66975309 0.196291610
  [2,] 10  0  0         NA 1.000000000
  [3,]  4  3  3 0.61075400 0.434504055
  [4,]  5  4  1 0.23185941 0.630148518
  [5,]  2  6  2 0.05000000 0.823063274
  [6,]  6  3  1 0.04444444 0.833028894
  [7,]  3  3  4 0.61075400 0.434504055
  [8,]  5  4  1 0.23185941 0.630148518
  [9,]  2  2  6 1.16156463 0.281141287
 [10,]  8  2  0 1.66975309 0.196291610
 [11,]  0  3  7 0.34755863 0.555498394
 [12,]  8  2  0 1.66975309 0.196291610
 [13,]  7  3  0 0.34755863 0.555498394
 [14,]  0  3  7 0.34755863 0.555498394
 [15,]  2  6  2 0.05000000 0.823063274
 [16,]  8  2  0 1.66975309 0.196291610
 [17,]  0  3  7 0.34755863 0.555498394
 [18,]  0  6  4 0.61281179 0.433731081
 [19,]  2  4  4 0.02777778 0.867632335
 [20,]  6  2  2 1.16156463 0.281141287
 [21,]  0  5  5 0.17777778 0.673289980
 [22,]  5  4  1 0.23185941 0.630148518
 [23,]  0  0 10         NA 1.000000000
 [24,]  8  2  0 1.66975309 0.196291610
 [25,]  0  0 10         NA 1.000000000
 [26,]  3  6  1 0.11111111 0.738882680
 [27,]  0  3  7 0.34755863 0.555498394
 [28,]  1  2  7 0.17968750 0.671641936
 [29,]  1  3  6 0.04444444 0.833028894
 [30,]  0  9  1 4.46321804 0.034632195
 [31,] 10  0  0         NA 1.000000000
 [32,]  5  2  3 1.83333333 0.175734336
 [33,]  5  3  2 0.32194179 0.570443028
 [34,]  1  7  2 0.71258035 0.398588527
 [35,]  4  6  0 0.61281179 0.433731081
 [36,]  2  5  3 0.22691562 0.633820341
 [37,]  2  6  2 0.05000000 0.823063274
 [38,]  1  2  7 0.17968750 0.671641936
 [39,]  0  2  8 1.66975309 0.196291610
 [40,]  7  3  0 0.34755863 0.555498394
 [41,]  9  1  0 9.26315789 0.002338095
 [42,]  1  5  4 0.08018355 0.777048820
 [43,] 10  0  0         NA 1.000000000
 [44,]  0  4  6 0.05468750 0.815098647
 [45,]  3  6  1 0.11111111 0.738882680
 [46,]  0  3  7 0.34755863 0.555498394
 [47,]  1  3  6 0.04444444 0.833028894
 [48,]  2  6  2 0.05000000 0.823063274
 [49,]  4  5  1 0.08018355 0.777048820
 [50,]  8  2  0 1.66975309 0.196291610
 [51,]  1  5  4 0.08018355 0.777048820
 [52,] 10  0  0         NA 1.000000000
 [53,]  0  6  4 0.61281179 0.433731081
 [54,]  6  2  2 1.16156463 0.281141287
 [55,]  2  3  5 0.32194179 0.570443028
 [56,]  7  2  1 0.17968750 0.671641936
 [57,]  0  0 10         NA 1.000000000
 [58,]  0  2  8 1.66975309 0.196291610
 [59,]  0  0 10         NA 1.000000000
 [60,]  7  3  0 0.34755863 0.555498394
 [61,]  0  2  8 1.66975309 0.196291610
 [62,]  0  2  8 1.66975309 0.196291610
 [63,]  3  6  1 0.11111111 0.738882680
 [64,]  2  8  0 2.61111111 0.106117388
 [65,]  4  4  2 0.02777778 0.867632335
 [66,]  4  5  1 0.08018355 0.777048820
 [67,]  6  4  0 0.05468750 0.815098647
 [68,]  1  2  7 0.17968750 0.671641936
 [69,]  0  0 10         NA 1.000000000
 [70,]  2  3  5 0.32194179 0.570443028
 [71,]  0  0 10         NA 1.000000000
 [72,]  3  4  3 0.05000000 0.823063274
 [73,]  0  6  4 0.61281179 0.433731081
 [74,]  2  5  3 0.22691562 0.633820341
 [75,]  7  2  1 0.17968750 0.671641936
 [76,]  7  3  0 0.34755863 0.555498394
 [77,]  0  2  8 1.66975309 0.196291610
 [78,]  0  5  5 0.17777778 0.673289980
 [79,]  5  3  2 0.32194179 0.570443028
 [80,]  0  2  8 1.66975309 0.196291610
 [81,]  0  2  8 1.66975309 0.196291610
 [82,]  0  2  8 1.66975309 0.196291610
 [83,]  1  3  6 0.04444444 0.833028894
 [84,]  1  2  7 0.17968750 0.671641936
 [85,]  1  4  5 0.23185941 0.630148518
 [86,]  0  1  9 9.26315789 0.002338095
 [87,]  9  1  0 9.26315789 0.002338095
 [88,]  2  3  5 0.32194179 0.570443028
 [89,]  7  3  0 0.34755863 0.555498394
 [90,]  3  7  0 1.38920420 0.238539537
 [91,]  0  3  7 0.34755863 0.555498394
 [92,]  8  1  1 0.77893118 0.377468222
 [93,]  3  5  2 0.22691562 0.633820341
 [94,]  8  2  0 1.66975309 0.196291610
 [95,]  4  5  1 0.08018355 0.777048820
 [96,]  2  4  4 0.02777778 0.867632335
 [97,]  1  6  3 0.11111111 0.738882680
 [98,]  9  1  0 9.26315789 0.002338095
 [99,]  5  4  1 0.23185941 0.630148518
[100,]  0  7  3 1.38920420 0.238539537
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>