Returns the minimum value of r2 in r2.vec that keeps the number of constraint matrix columns less than maxcol.

1 | ```
r2.cut.fn(block.obj, block.cood, Xa, Y, maxcol, r2.cut.min, r2.cut.max, r2.vec.length)
``` |

`block.obj` |
An object of class gwaa.data from GenABEL. |

`block.cood` |
A vector of length p+1, where p is the number of SNPs. block.cood is an indicator vector that indicates block boundaries at all p+1 SNP bounded intervals. Use find.bounds to create this vector. |

`Xa` |
If block.obj is NA then a genotype matrix must be provided. Xa is a matrix of genotype values codes as 0, 1 or 2 for homozygous major, heterozygous, or homozygous minor, respectively. |

`Y` |
If block.obj is NA then a phenotype vector Y must be provided. Y is a vector of diagnoses, where 0 is non-diseased and 1 is diseased. |

`maxcol` |
The maximum number of columns allowed in the constraint matrix. Default is 5,000. |

`r2.cut.min` |
The lower bound of the r2.cut interval |

`r2.cut.max` |
The upper bound of the r2.cut interval |

`r2.vec.length` |
The number of r2.cut values uniformly spaced in r2.cut interval |

The function r2.cut.fn returns the minimum value of r2.cut that satsisfies, ncol(A) < maxcol, where A is the constraint matrix

Samuel G. Younkin

ld_lasso

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.