opencobra / cobratoolbox

The COnstraint-Based Reconstruction and Analysis Toolbox. Documentation:
https://opencobra.github.io/cobratoolbox
Other
252 stars 317 forks source link

mapExpressionToReactions returning NaN array #1942

Open anilk991 opened 2 years ago

anilk991 commented 2 years ago

The [expressionRxns, parsedGPR, gene_used] = mapExpressionToReactions(model, expData); is giving expressionRxns variable as entirely NaN values. But parsedGPR and gene_used return cell array. The genes in expData are in entrez format without version, on the other hand genes in model.genes are in entrez format with version.

I hereby confirm that I have:

(Note: You may replace [ ] with [X] to check the box)

rmtfleming commented 2 years ago

expressionRxns(j) = NaN when there is no expression data for the corresponding gene.

On Wed, 23 Feb 2022 at 11:33, anilk991 @.***> wrote:

The [expressionRxns, parsedGPR, gene_used] = mapExpressionToReactions(model, expData); is giving expressionRxns variable as entirely NaN values. But parsedGPR and gene_used return cell array. The genes in expData are in entrez format without version, on the other hand genes in model.genes are in entrez format with version.

I hereby confirm that I have:

  • [ X] Tried to solve the issue on my own
  • [ X] Retried to run my code with the latest version of The COBRA Toolbox
  • [ X] Checked that a similar issue has not already been opened

(Note: You may replace [ ] with [X] to check the box)

— Reply to this email directly, view it on GitHub https://github.com/opencobra/cobratoolbox/issues/1942, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAQMEOVCAGBVNRMG2ZUO3UDU4TAXTANCNFSM5PECQ55Q . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you are subscribed to this thread.Message ID: @.***>

--

Mr. Ronan MT Fleming B.V.M.S. Dip. Math. Ph.D.

Associate Professor, School of Medicine, National University of Ireland, Galway. & Assistant Professor, Division of Systems Biomedicine and Pharmacology, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University. https://www.universiteitleiden.nl/en/staffmembers/ronan-fleming & H2020 Project Coordinator, Systems Medicine of Mitochondrial Parkinson’s Disease, http://sysmedpd.eu

Peer-reviewed publications: https://goo.gl/FZPG23 Mobile: +353 852 109 806 Skype: ronan.fleming

(This message is confidential and may contain privileged information. It is intended for the named recipient only. If you receive it in error please notify me and permanently delete the original message and any copies.)

anilk991 commented 2 years ago

The issue is that all the reactions are coming out as NaN.

On Wed, Feb 23, 2022, 18:49 Ronan M.T. Fleming @.***> wrote:

expressionRxns(j) = NaN when there is no expression data for the corresponding gene.

On Wed, 23 Feb 2022 at 11:33, anilk991 @.***> wrote:

The [expressionRxns, parsedGPR, gene_used] = mapExpressionToReactions(model, expData); is giving expressionRxns variable as entirely NaN values. But parsedGPR and gene_used return cell array. The genes in expData are in entrez format without version, on the other hand genes in model.genes are in entrez format with version.

I hereby confirm that I have:

  • [ X] Tried to solve the issue on my own
  • [ X] Retried to run my code with the latest version of The COBRA Toolbox
  • [ X] Checked that a similar issue has not already been opened

(Note: You may replace [ ] with [X] to check the box)

— Reply to this email directly, view it on GitHub https://github.com/opencobra/cobratoolbox/issues/1942, or unsubscribe < https://github.com/notifications/unsubscribe-auth/AAQMEOVCAGBVNRMG2ZUO3UDU4TAXTANCNFSM5PECQ55Q

. Triage notifications on the go with GitHub Mobile for iOS < https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675

or Android < https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub .

You are receiving this because you are subscribed to this thread.Message ID: @.***>

--

Mr. Ronan MT Fleming B.V.M.S. Dip. Math. Ph.D.


Associate Professor, School of Medicine, National University of Ireland, Galway. & Assistant Professor, Division of Systems Biomedicine and Pharmacology, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University. https://www.universiteitleiden.nl/en/staffmembers/ronan-fleming & H2020 Project Coordinator, Systems Medicine of Mitochondrial Parkinson’s Disease, http://sysmedpd.eu


Peer-reviewed publications: https://goo.gl/FZPG23 Mobile: +353 852 109 806 Skype: ronan.fleming


(This message is confidential and may contain privileged information. It is intended for the named recipient only. If you receive it in error please notify me and permanently delete the original message and any copies.)

— Reply to this email directly, view it on GitHub https://github.com/opencobra/cobratoolbox/issues/1942#issuecomment-1048773417, or unsubscribe https://github.com/notifications/unsubscribe-auth/AI2IQMKIVWQQL2RY6KQZFRTU4TNFPANCNFSM5PECQ55Q . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you authored the thread.Message ID: @.***>

rmtfleming commented 2 years ago

what is the input? are you sure those genes correspond to the reactions?

On Wed, 23 Feb 2022 at 13:34, anilk991 @.***> wrote:

The issue is that all the reactions are coming out as NaN.

On Wed, Feb 23, 2022, 18:49 Ronan M.T. Fleming @.***> wrote:

expressionRxns(j) = NaN when there is no expression data for the corresponding gene.

On Wed, 23 Feb 2022 at 11:33, anilk991 @.***> wrote:

The [expressionRxns, parsedGPR, gene_used] = mapExpressionToReactions(model, expData); is giving expressionRxns variable as entirely NaN values. But parsedGPR and gene_used return cell array. The genes in expData are in entrez format without version, on the other hand genes in model.genes are in entrez format with version.

I hereby confirm that I have:

  • [ X] Tried to solve the issue on my own
  • [ X] Retried to run my code with the latest version of The COBRA Toolbox
  • [ X] Checked that a similar issue has not already been opened

(Note: You may replace [ ] with [X] to check the box)

— Reply to this email directly, view it on GitHub https://github.com/opencobra/cobratoolbox/issues/1942, or unsubscribe <

https://github.com/notifications/unsubscribe-auth/AAQMEOVCAGBVNRMG2ZUO3UDU4TAXTANCNFSM5PECQ55Q

. Triage notifications on the go with GitHub Mobile for iOS <

https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675

or Android <

https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub

.

You are receiving this because you are subscribed to this thread.Message ID: @.***>

--

Mr. Ronan MT Fleming B.V.M.S. Dip. Math. Ph.D.


Associate Professor, School of Medicine, National University of Ireland, Galway. & Assistant Professor, Division of Systems Biomedicine and Pharmacology, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University. https://www.universiteitleiden.nl/en/staffmembers/ronan-fleming & H2020 Project Coordinator, Systems Medicine of Mitochondrial Parkinson’s Disease, http://sysmedpd.eu


Peer-reviewed publications: https://goo.gl/FZPG23 Mobile: +353 852 109 806 Skype: ronan.fleming


(This message is confidential and may contain privileged information. It is intended for the named recipient only. If you receive it in error please notify me and permanently delete the original message and any copies.)

— Reply to this email directly, view it on GitHub < https://github.com/opencobra/cobratoolbox/issues/1942#issuecomment-1048773417 , or unsubscribe < https://github.com/notifications/unsubscribe-auth/AI2IQMKIVWQQL2RY6KQZFRTU4TNFPANCNFSM5PECQ55Q

. Triage notifications on the go with GitHub Mobile for iOS < https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675

or Android < https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub .

You are receiving this because you authored the thread.Message ID: @.***>

— Reply to this email directly, view it on GitHub https://github.com/opencobra/cobratoolbox/issues/1942#issuecomment-1048787545, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAQMEOWKS7EPUAY3ZKIZHQTU4TO5JANCNFSM5PECQ55Q . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you commented.Message ID: @.***>

--

Mr. Ronan MT Fleming B.V.M.S. Dip. Math. Ph.D.

Associate Professor, School of Medicine, National University of Ireland, Galway. & Assistant Professor, Division of Systems Biomedicine and Pharmacology, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University. https://www.universiteitleiden.nl/en/staffmembers/ronan-fleming & H2020 Project Coordinator, Systems Medicine of Mitochondrial Parkinson’s Disease, http://sysmedpd.eu

Peer-reviewed publications: https://goo.gl/FZPG23 Mobile: +353 852 109 806 Skype: ronan.fleming

(This message is confidential and may contain privileged information. It is intended for the named recipient only. If you receive it in error please notify me and permanently delete the original message and any copies.)

anilk991 commented 2 years ago

expData is in form of a struct with expData.gene as double array with entrez gene ids and expData.value as logical array.