QidiPeng / eQTLMAPT

2 stars 2 forks source link

gmap.gpd and gmap.ac.gpd do not work #2

Open juanb001 opened 1 year ago

juanb001 commented 1 year ago

Hello,

I've had success running the gmap() function without issues.

However, I have been unable to get the above functions to work properly.

When I run gmap.gpd() using the same genotype, expression, etc. inputs as the gmap() function, I get this error:

[1] "call Pgpd" Error in if (abs(x) < REALMIN) return(2^(-1074)) : the condition has length > 1

Similarly, when I try to run the gmap.ac.gpd function, I get this error:

Error in svd(x, nu = 0, nv = k) : infinite or missing values in 'x'

I've double checked, and there are no NAs in my expression matrix.

For reference, I am running R 4.2.1

stormlovetao commented 1 year ago

Sorry for the late reply. We have fixed the bug. Please retry!

juanb001 commented 1 year ago

Hello!

I tried re-installing eqtlmapt today and tried the gmap.gpd function. However, the same issue remains:

[1] "call Pgpd" Error in if (abs(x) < REALMIN) { : the condition has length > 1

Any suggestions?