Entering edit mode
2.1 years ago
TM
▴
20
Hi, I got counting table from HTSeq, now I ran through the pipeline in edgeR. After I got the topTags gene, I created a table to show differential expression and I got an error like this:
Error in data.frame(..., check.names = FALSE) :
arguments imply differing number of rows: 60649, 0
This is my script:
genediiffBB$table <- cbind(genediffBB$table, FDR=p.adjust(genediffBB$tablePValue, method='fdr'))
Could you tell what is the problem and solution please? I have attached the row 60649 in the picture, Note. I have already delete the line after 60649 before taking the counting table into edgeR.
Thank you in advance!
We can't see the figure but it is likely that
FDR=p.adjust(genediffBB$tablePValue, method='fdr'))
does not work, please could you tell us the output ofstr(genediffBB)
Hi this is the output
You need
cbind(genediffBB$table, FDR=p.adjust(genediffBB$table$PValue, method='fdr'))
, pay attention to the syntax you forgot one $ after genediffBB$tableI have problem in uploading the picture but this is the line 60649 (the last row of my count table)
It is working. Thank you your much!!
code
option) to present your post better. You can use backticks for inline code (`text` becomestext
), or select a chunk of text and use the highlighted button to format it as a code block. If your code has long lines with a single command, break those lines into multiple lines with proper escape sequences so they're easier to read and still run when copy-pasted. I've done it for you this time.I am new to this platform, thank you for your guidance!!
No problem, but you've still not done #1 properly - edit your post and upload your image the right way, please.