Error in interpp.old missing values and Infs not allowed

[qvhaelen]

Hello,

I try to use your method FUNNEL with a time series data (8 time points) and pathway database. I tried with the main function result <- FUNNEL.GSEA(X, tt, genesets), this example works well and works also when replacing genesets by my own pathway list "pathway_gene_list": FUNNEL.GSEA(X, tt, pathway_gene_list). however, when using my own transcription data, I got an error:

Error in interpp.old(x, y, z, xo, yo, ncp = 0, extrap = FALSE, duplicate = duplicate, :
missing values and Infs not allowed

the vector tt is adapted to contain the right number of time points. I have checked if the format of my matrix is the same than the one provided as an example and it is same data type:

str(X)
num [1:11189, 1:16] 2.26 9.38 1.96 8.19 7.53 ...

• attr(, "dimnames")=List of 2
  ..$ : Named chr [1:11189] "NAT2" "ADA" "CDH2" "AKT3" ...
  .. ..- attr(, "names")= chr [1:11189] "10" "100" "1000" "10000" ...
  ..$ : chr [1:16] "t1" "t2" "t3" "t4" ...

str(x3array)
num [1:20340, 1:8] 0.366 0 1.077 0.169 0.139 ...

• attr(, "dimnames")=List of 2
  ..$ : Named chr [1:20340] "A1BG" "A1BG-AS1" "A1CF" "A2LD1" ...
  .. ..- attr(, "names")= chr [1:20340] "10" "20" "30" "40" ...
  ..$ : chr [1:8] "t1" "t2" "t3" "t4" ...

there is no NA in the file. I do not know what to do or where to look to fix this error. Any advise or suggestion will be helpful?

thank you

Quentin

[yunzhang813]

Hi Quentin, Thank you for bringing this issue to us. This error message is because interpp.old() is deprecated in a package called ‘akima’ that FUNNEL depends on. Please see the link below. https://www.rdocumentation.org/packages/akima/versions/0.6-2/topics/interpp.old For a quick fix, you may look for an old version of akima, or locate the interpp.old() function and replace it with interpp() as noted in the above page. Your message has brought me series attention on this. However, I am not able to work on this in coming weekends. If you can wait, I plan to submit an updated version in a month or so. Please let us know if you have any other feedbacks for the package. Thanks a lot, Yun From: qvhaelen <[email protected]> Reply-To: yunzhang813/FUNNEL-GSEA-R-Package <[email protected]> Date: Friday, September 21, 2018 at 5:25 AM To: yunzhang813/FUNNEL-GSEA-R-Package <[email protected]> Cc: Subscribed <[email protected]> Subject: [yunzhang813/FUNNEL-GSEA-R-Package] Error in interpp.old missing values and Infs not allowed (#2) Hello, I try to use your method FUNNEL with a time series data (8 time points) and pathway database. I tried with the main function result <- FUNNEL.GSEA(X, tt, genesets), this example works well and works also when replacing genesets by my own pathway list "pathway_gene_list": FUNNEL.GSEA(X, tt, pathway_gene_list). however, when using my own transcription data, I got an error: Error in interpp.old(x, y, z, xo, yo, ncp = 0, extrap = FALSE, duplicate = duplicate, : missing values and Infs not allowed the vector tt is adapted to contain the right number of time points. I have checked if the format of my matrix is the same than the one provided as an example and it is same data type: str(X) num [1:11189, 1:16] 2.26 9.38 1.96 8.19 7.53 ... * attr(, "dimnames")=List of 2 ..$ : Named chr [1:11189] "NAT2" "ADA" "CDH2" "AKT3" ... .. ..- attr(, "names")= chr [1:11189] "10" "100" "1000" "10000" ... ..$ : chr [1:16] "t1" "t2" "t3" "t4" ... str(x3array) num [1:20340, 1:8] 0.366 0 1.077 0.169 0.139 ... * attr(, "dimnames")=List of 2 ..$ : Named chr [1:20340] "A1BG" "A1BG-AS1" "A1CF" "A2LD1" ... .. ..- attr(, "names")= chr [1:20340] "10" "20" "30" "40" ... ..$ : chr [1:8] "t1" "t2" "t3" "t4" ... there is no NA in the file. I do not know what to do or where to look to fix this error. Any advise or suggestion will be helpful? thank you Quentin — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub<#2>, or mute the thread<https://github.com/notifications/unsubscribe-auth/ANi8MmUXrEJ2It2oZuM0SA2OdGTzY11zks5udNqigaJpZM4W0DKS>.

[qvhaelen]

hello

Thank you for the answer and information provided. nevertheless, I have to say I am still confused because of the following. When running FUNNEL using the example provided in the package, result <- FUNNEL.GSEA(X, tt, genesets), the program does indicate a warning saying that
Warning messages:
1: In interpp.old(x, y, z, xo, yo, ncp = 0, extrap = FALSE, duplicate = duplicate, :
interpp.old() is deprecated, future versions will only provide interpp()

But it is working anyway. The error message which stops the calculation (Error in interpp.old(x, y, z, xo, yo, ncp = 0, extrap = FALSE, duplicate = duplicate, :
missing values and Infs not allowed) appears only when running the same function using my own data X rather than the one provided in the FUNNEL package. Modifying the lsit of pathways for example does not generate error message.

thank you

quentin