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

Integer copy number matrix output for large sample sets ? #6

Open
RodrigoGM opened this issue May 7, 2022 · 1 comment
Open

Integer copy number matrix output for large sample sets ? #6

RodrigoGM opened this issue May 7, 2022 · 1 comment

Comments

@RodrigoGM
Copy link

Hi !
I tested ACE in a data set of 3000+ samples and had a successful run, I really like the figures you have with the absolute copy number predictions. However, I noticed the runACE does not export a table with segment integer copy numbers; is there a way to export this? i.e. on a per-sample basis, have the integer copy number from the best fit exported a table ? Even if it's a text file per sample, these can be combined into a single table with a loop of some sorts.

I have not been able to find this in the documentation; for one or two samples, the manual process is OK, but for a data set of 3000+ samples, it becomes impractical to go one by one.

Thank you, any help would be most appreciated,

@jpoell
Copy link
Collaborator

jpoell commented Jan 14, 2023

I found out I did not receive notifications for issues on ACE, so I am now dealing with these very old issues. So sorry for that!
Analyzing over 3000 samples means you do not want to go check your fits manually. I think you should only use runACE if you can be quite sure that 2N fits are appropriate for the large majority of your sample. You could switch to looping through the samples and do squaremodel fits. It will take a bit longer, but shouldn't be too bad. Since you are then looping through all samples, you can run getadjustedsegments for each sample and retrieve the segment integer copies from the resulting data frame.

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