Group SNPs; adapted from group.multi by Chris Wallace
groupmulti( SM2, snp.data, is.snpmat, min.mppi = 0.01, minsnpmppi = 0.01, r2.minmerge = 0.5 )
| SM2 | snpmod object |
|---|---|
| snp.data | SnpMatrix or SNP covariance matrix with snp ids as column names |
| is.snpmat | logical taking value TRUE when SnpMatrix is provided and FALSE when covariance matrix is given |
| min.mppi | trim snp groups with total MPPI < min.mppi in all diseases |
| minsnpmppi | only group snps with total MPPI > minsnpmppi |
| r2.minmerge | merge groups with minimum between-group r2 > r2.minmerge |
list with three components.
First is a data.frame with each row giving summary statistics for each group.
Second is a groups object, each elements ordered according to the rows of the summary
Third is the r2 matrix calculated.