I did RNA-seq analysis using DESeq2. Got the results and want to draw a Kegg pathway with the enriched gene sets. So, I am trying to use gage function.

I followed the instructions in the manual and got to exp.fc as below:
head(exp.fc)
lp*0001 lp*0002 lp*0004 lp*0005 lp*0006 lp*0007
-0.3297675 0.2703983 -0.1043349 -0.4814408 1.2165759 -0.3020035

and then I extracted genesets using kegg.gsets(species ="lpl", id.type = "kegg", check.new=FALSE). this is my output

head(gsets)
$kg.sets
$kg.sets$`lpl00010 Glycolysis / Gluconeogenesis`

[1] "lp*0082" "lp*0330" "lp*0350" "lp*0440" "lp*0537" "lp*0597" "lp*0764" "lp*0789" "lp*0790"
[10] "lp*0791" "lp*0792" "lp*0826" "lp*0884" "lp*0906" "lp*1101" "lp*1245" "lp*1401" "lp*1573"
[19] "lp*1665" "lp*1731" "lp*1897" "lp*1898" "lp*1920" "lp*2151" "lp*2152" "lp*2153" "lp*2154"
[28] "lp*2349" "lp*2502" "lp*2777" "lp*2778" "lp*2873" "lp*3008" "lp*3011" "lp*3132" "lp*3170"
[37] "lp*3418" "lp*3487" "lp*3512" "lp*3520" "lp*3525" "lp*3526" "lp*3629" "lp*3662"

So my gsets have the same id as my exp.fc. However when I apply gage function, fc.kegg.p <-gage(exp.fc, gsets = gsets, ref = NULL, samp = NULL)

my fc.kegg.p is all NAN as below: head(fc.kegg.p) $greater p.geomean stat.mean p.val q.val set.size exp1 kg.sets NA NaN NA NA 3 NA sigmet.idx NA NaN NA NA 0 NA sig.idx NA NaN NA NA 0 NA met.idx NA NaN NA NA 0 NA dise.idx NA NaN NA NA 0 NA

$less p.geomean stat.mean p.val q.val set.size exp1 kg.sets NA NaN NA NA 3 NA sigmet.idx NA NaN NA NA 0 NA sig.idx NA NaN NA NA 0 NA met.idx NA NaN NA NA 0 NA dise.idx NA NaN NA NA 0 NA

$stats stat.mean exp1 kg.sets NaN NA sigmet.idx NaN NA sig.idx NaN NA met.idx NaN NA dise.idx NaN NA

I saw some suggestions as using ENTREZ ID. but the IDs are matching in both kg.sets and my exp.fc. there is one thing that I cannot figure out how gage does not match it. are the IDs really the problem? or is there something that I miss?