Description Usage Arguments Details Value Author(s) Examples

Performs multiple association tests.

1 |

`y` |
numeric vector with phenotype status: 0=controls, 1=cases. No missing data allowed |

`X` |
numeric matrix or data frame with genotype data coded as 0, 1, 2. Missing data is allowed |

`tests` |
character vector with names of the tests to be applied |

`maf` |
numeric value indicating the minor allele frequency threshold for rare variants ( |

`perm` |
positive integer indicating the number of permutations (100 by default) |

`weights` |
optional vector of weights for the variants ( |

`c.param` |
Optional value to specify the |

The available tests are: `"WSS", "ORWSS", "RWAS", "CMC", "CMAT", "CALPHA", "RBT", "SCORE", "SUM", "SSU", "SSUW", "UMINP", "BST", "WST", "RVT1", "RVT2", "VT"`

There is no imputation for the missing data. Missing values are simply ignored in the computations.

A data frame with test statistics and permutated p-values

Gaston Sanchez

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 | ```
## Not run:
# number of cases
cases = 250
# number of controls
controls = 250
# total (cases + controls)
total = cases + controls
# phenotype vector
phenotype = c(rep(1, cases), rep(0, controls))
# genotype matrix with 10 variants (random data)
set.seed(1234)
genotype = matrix(rbinom(total*10, 2, 0.051), nrow=total, ncol=10)
# apply MULTI with "BST", "CMC", "RWAS" and 100 permutations
mymulti1 = MULTI(phenotype, genotype, c("BST", "CMC", "RWAS"), perm=100)
# this is what we get
mymulti1
# create list with the following tests
test_list = c("BST", "CMC", "CMAT", "CALPHA", "ORWSS", "RWAS",
"RBT", "SCORE", "SUM", "SSU", "SSUW", "UMINP", "WSS", "WST")
# apply MULTI with 100 permutations
mymulti2 = MULTI(phenotype, genotype, test_list, perm=100)
# this is what we get
mymulti2
## End(Not run)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.