Age | Commit message (Collapse) | Author |
|
The script now takes the --outdir parameter as the output directory, only using TMPDIR if nothing is passed
"output" in the output directory was also changed to "gn3" to be consistent with recent changes to GN3's code
|
|
|
|
|
|
environment
Also fixes issue where the control marker keyword was wrong
|
|
https://github.com/zsloan/genenetwork3 into feature/add_rqtl_pairscan
|
|
Now it checks for pairscan first, just in case interval ends up being
passed (which is an irrelevant parameter for pairscan)
Also added a couple more verbose prints
|
|
we need the list of markers/pseudomarkers and their positions)
|
|
step-size to 10cM for pair-scan
|
|
(pairscan) is used
- For pairscan default to using step 20 (subject to change, but some
step is required during calc.genoprob to make it run fast enough)
- Added some new verbose prints
|
|
|
|
The rqtl_wrapper script was throwing an error when only a single
categorical covariate was used. This is apparently because "covars[,name]"
throws an error in such a situation. Using just "covars" in such a
situation prevents the error. So I just added an if statement checking
the number of covariates. There might be some better way to deal with
this in R, but this is the best I could come up with.
|
|
|
|
|
|
|
|
we need the list of markers/pseudomarkers and their positions)
|
|
step-size to 10cM for pair-scan
|
|
(pairscan) is used
- For pairscan default to using step 20 (subject to change, but some
step is required during calc.genoprob to make it run fast enough)
- Added some new verbose prints
|
|
|
|
|
|
stratified permutations + covariates
|
|
|
|
Fixed issue where covariates with numerical names were being ignored
Fixed issue where the hash filename wasn't being used for the permutation output
|
|
|