Question: No differential expression with limmaCtData possible
gravatar for h.roest
12 months ago by
h.roest0 wrote:

Dear All,

Tried rerunning an R file that I made about a year ago to analyse TLDA array cards. However, when arriving at determining differentially expressed miRNAs I get an error where this returned a perfect set of differentially expressed miRNAs a year ago.

design <- model.matrix(~0 + sampleinfo$Function)
colnames(design) <- c("EAD", "IF")
contrasts <- makeContrasts(EAD - IF,levels = design)
colnames(contrasts) <- c("EAD - IF")
contrasts<[order(featureNames(] <- limmaCtData(,design=design, contrasts=contrasts,  stringent = FALSE, sort=TRUE)

Error in data.frame(..., check.names = FALSE) : 
  arguments imply differing number of rows: 192, 0

I checked the object but 

 @ featureData

      @ data indicates a 'data.frame' of 192 obs. of 1 variable:

Being a beginner in R it took me a lot of time to get this package to work.

> sessionInfo()
R version 3.4.0 (2017-04-21)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

Matrix products: default

[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United States.1252    

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] HTqPCR_1.30.0       limma_3.32.10       RColorBrewer_1.1-2  Biobase_2.36.2      BiocGenerics_0.22.1

loaded via a namespace (and not attached):
 [1] gtools_3.5.0          bitops_1.0-6          affy_1.54.0           stats4_3.4.0          KernSmooth_2.23-15    BiocInstaller_1.26.1 
 [7] gplots_3.0.1          zlibbioc_1.22.0       gdata_2.18.0          affyio_1.46.0         preprocessCore_1.38.1 tools_3.4.0          
[13] compiler_3.4.0        caTools_1.17.1 


ADD COMMENTlink modified 11 months ago • written 12 months ago by h.roest0
gravatar for Axel Klenk
12 months ago by
Axel Klenk920
Axel Klenk920 wrote:

Dear h.roest,

looking at your code, I see that being used without ever being read from a data file

or otherwise created.  Is it possible that you're running your R script line by line and accidentally

skipped one or more lines?  Just a guess.


 - axel


ADD COMMENTlink written 12 months ago by Axel Klenk920
gravatar for h.roest
11 months ago by
h.roest0 wrote:

Dear Axel,

Thanks for your comment. Unfortunately, the "" object was created by filtering out the targets that were more then 4 out of 8 Undetermined. It looks that the problem actually starts here. limma of the unfiltered HTqPCR object can nicely pass the differential expression using limmaCtData. Remarkably, mutest and ttest work fine on the object. Maybe you have some more suggestions on how to solve this problem?



ADD COMMENTlink written 11 months ago by h.roest0
Please log in to add an answer.


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