Question: deltaCq error in package NormqPCR
1
gravatar for jacorvar
3.7 years ago by
jacorvar40
European Union
jacorvar40 wrote:

Hi,

I get the following error when I type deltaCq(data, "hsa-miR-875-3p-002204"):

Error in .local(qPCRBatch, ...) : 
  given housekeeping gene, hsa.miR.875.3p.002204 not found in file. Ensure entered housekeeping genes appear in the file

However, when I do "hsa-miR-875-3p-002204" %in% featureNames(data), I get  TRUE.

Could you help me to solve it?

normqpcr qrt-pcr deltacq • 543 views
ADD COMMENTlink modified 3.7 years ago • written 3.7 years ago by jacorvar40
Answer: deltaCq error in package NormqPCR
2
gravatar for James W. MacDonald
3.7 years ago by
United States
James W. MacDonald51k wrote:

This is an infelicity in the NormqPCR package. Internally it calls make.names on the housekeeping gene name:

> make.names("hsa-miR-875-3p-002204")
[1] "hsa.miR.875.3p.002204"

And you will note that hsa.miR.875.3p.00204 is not the same as hsa-miR-875-3p-002204, hence the error. You could get around this by running make.names on all your featureNames first.

featureNames(data) <- make.names(featureNames(data))
ADD COMMENTlink written 3.7 years ago by James W. MacDonald51k

It gives me an error:

featureNames(data) <- make.names(rownames(featureData(data)))
Error in `featureNames<-`(`*tmp*`, value = character(0)) : 
  'value' length (0) must equal feature number in AssayData (384)
ADD REPLYlink modified 3.7 years ago • written 3.7 years ago by jacorvar40
Answer: deltaCq error in package NormqPCR
2
gravatar for jacorvar
3.7 years ago by
jacorvar40
European Union
jacorvar40 wrote:

Thanks for the hint. I tried it and didn't work, since it was unable to change the feature names. Then I did the following and it worked:

fData = data.frame(row.names=paste(1:length(featureNames(data)), make.names(featureNames(data)), sep="."))

featureData(data) = new("AnnotatedDataFrame", data=fData)

featureNames(data) = rownames(featureData(data))
ADD COMMENTlink written 3.7 years ago by jacorvar40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 332 users visited in the last hour