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

COMPARISON BUG #4

Open
H1889 opened this issue Jul 14, 2023 · 1 comment
Open

COMPARISON BUG #4

H1889 opened this issue Jul 14, 2023 · 1 comment

Comments

@H1889
Copy link

H1889 commented Jul 14, 2023

Hi again.
I have run a comparison between two samples and when I compare the methylation and expression values that appear in the file are different from the values in the input files (plus adding 1 in expression and multiplying by 100 in methylation), for example:

CG_Gene_Body_chosen_genes.txt:
gene_ID trtmethmean trtexpmean ctrlmethmean ctrlexpmean Δmethylation Δgene expression
PHYBL_145520 0.0 17.51147464 51.250000279396765 1.0 -51.250000279396765 4.130228673018407

treatment_exp.txt
PHYBL_145520 0.0 (but the CG_Gene_Body_chosen_genes.txt it is shown 17.51147464)

control_exp.txt
PHYBL_145520 0.0 (correct, the program adds 1 to all expression values)

And same for the methylation values, the treatment values are incorrect and the control ones are ok. This occurs for all genes.

I noticed that the methylation and expression values correspond to other genes, it is as if there has been an exchange of values and it always happens in the treatment sample but not in the control sample.

@H1889
Copy link
Author

H1889 commented Jul 20, 2023

Hi.
I have been testing with the input files and I think this error occurs because in the preprocessing phase of the samples (preprocess.py) the methylation files (GeneMeant.txt, PromoterMean.txt etc) that are generated have a number of genes that differs with the number of genes that appear in the expression 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

1 participant