Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Covariate file can be empty? #5

Open
nunofonseca opened this issue Jan 8, 2016 · 1 comment
Open

Covariate file can be empty? #5

nunofonseca opened this issue Jan 8, 2016 · 1 comment
Assignees

Comments

@nunofonseca
Copy link
Collaborator

yep.
UserWarning: loadtxt: Empty input file: "fqtl/cov.tsv"

@clody23
Copy link
Collaborator

clody23 commented Jan 12, 2016

the tsv file with covariates is used by get_covariates.py script.
For the moment the script generates an hdf5 file with two groups:
1 - array of covariates
2 - array of sample names

This hdf5 file is then used by three other scripts (sort_ids.py, runpeer.py and cis/tran_eqtl.py).
If we want to change this, should the logic be like:

get_covariates.py :
if cov_file.txt is empty
then generate a fake covariates.hdf5 file

sort_ids.py:
if covariates.hdf5 file is fake
then output a file with the same name (*sorted.hdf5) but which is still fake
else sort the file.

runpeer.py:
if (for obscure reasons since the file is empty) the user select peer_cov=yes
then check if the covariates.hdf5 file is fake and print something to say that covariates will not be considered in the model
else pass

cis/tran_eqtl.py (here is complicated, but there can be a trick):
if covariates.hdf5 is fake
then peer_cov=yes (this option disable usage of covariate in the lmm model)
else pass

???

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants