mtdt: Transmission/disequilibrium test of a multiallelic marker

View source: R/mtdt.R

mtdtR Documentation

Transmission/disequilibrium test of a multiallelic marker

Description

Transmission/disequilibrium test of a multiallelic marker

Usage

mtdt(x, n.sim = 0)

Arguments

x

the data table.

n.sim

the number of simulations.

Details

This function calculates transmission-disequilibrium statistics involving multiallelic marker. Inside the function are tril and triu used to obtain lower and upper triangular matrices.

Value

It returned list contains the following components:

  • SE Spielman-Ewens Chi-square from the observed data.

  • ST Stuart or score Statistic from the observed data.

  • pSE the simulated p value.

  • sSE standard error of the simulated p value.

  • pST the simulated p value.

  • sST standard error of the simulated p value.

Author(s)

Mike Miller, Jing Hua Zhao

References

\insertRef

miller97gap

\insertRef

sham95gap

\insertRef

spielman96gap

\insertRef

zhao99gap

See Also

bt

Examples

## Not run: 
x <- matrix(c(0,0, 0, 2, 0,0, 0, 0, 0, 0, 0, 0,
              0,0, 1, 3, 0,0, 0, 2, 3, 0, 0, 0,
              2,3,26,35, 7,0, 2,10,11, 3, 4, 1,
              2,3,22,26, 6,2, 4, 4,10, 2, 2, 0,
              0,1, 7,10, 2,0, 0, 2, 2, 1, 1, 0,
              0,0, 1, 4, 0,1, 0, 1, 0, 0, 0, 0,
              0,2, 5, 4, 1,1, 0, 0, 0, 2, 0, 0,
              0,0, 2, 6, 1,0, 2, 0, 2, 0, 0, 0,
              0,3, 6,19, 6,0, 0, 2, 5, 3, 0, 0,
              0,0, 3, 1, 1,0, 0, 0, 1, 0, 0, 0,
              0,0, 0, 2, 0,0, 0, 0, 0, 0, 0, 0,
              0,0, 1, 0, 0,0, 0, 0, 0, 0, 0, 0),nrow=12)

# See note to bt for the score test obtained by SAS

mtdt(x)

## End(Not run)


gap documentation built on Aug. 26, 2023, 5:07 p.m.