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

The two generatestatetrees commands #18

Open
ybta2957e opened this issue Jan 3, 2022 · 2 comments
Open

The two generatestatetrees commands #18

ybta2957e opened this issue Jan 3, 2022 · 2 comments

Comments

@ybta2957e
Copy link

ybta2957e commented Jan 3, 2022

Happy New Year! Thank you for developing DeCiFer.
I found that there are two different generatestatetrees commands and I am not sure which one is appropriate.
The one of this two commands got an error and another one worked.

Here are commands of generatestatetrees

(decifer_v2.0.4) [toby@nis test_input]$ generatestatetrees -help
Usage:
  generatestatetrees [--help|-h|-help] [-S str] [-SS str] [-maxCN int]
     [-maxXY int]
Where:
  --help|-h|-help
     Print a short help message
  -S str
     Input state tree file
  -SS str
     Output state tree file
  -maxCN int
     Maximum number of copy number events (default: 2)
  -maxXY int
     Maximum number of maternal/paternal copies (default: 2)
(decifer_v2.0.4) [toby@nis test_input]$ generatestatetrees -S CCLF-020_cn_states_unique.txt -SS CCLF-020_my_state_trees_unique.txt
Line: 4. Error: invalid copy number '1,1;2,2'
(decifer_v2.0.4) [toby@nis test_input]$ generatestatetrees CCLF-020_cn_states_unique.txt > CCLF-020_my_state_trees_unique.txt
(1,1)... Done.
(1,1) (2,0)... Done.
(1,1) (2,1)... Done.
(1,1) (2,2)... Done.
(1,1) (2,2) (5,1)... Done.
(1,1) (4,1)... Done.
Processed CCLF-020_cn_states_unique.txt...
(decifer_v2.0.4) [toby@nis test_input]$ 

Thank you
Toby

@brian-arnold
Copy link
Collaborator

Thanks for catching this. For now please use generatestatetrees without the -S flag, i.e. the second way you've used it.

@melkebir
Copy link
Member

Just to add to this, -S takes as input generatestatetrees output. That way you can build a comprehensive state trees file.

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

3 participants