SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Please HELP! TCGA RNAseq data and Limma Mous Bioinformatics 4 11-22-2014 02:16 AM
An stubborn error in limma wbb1813 Bioinformatics 2 08-17-2014 11:43 PM
FPKMs and Limma R package EA01 Bioinformatics 3 06-15-2014 11:55 PM
Error in duplicateCorrelation.. limma megcyan Bioinformatics 1 06-15-2014 11:51 PM
Limma removeBatchEffect question Justin AC Powell Bioinformatics 0 02-17-2014 06:13 AM

Reply
 
Thread Tools
Old 10-10-2016, 04:32 PM   #1
vsmuir
Junior Member
 
Location: Seattle, WA

Join Date: Oct 2016
Posts: 3
Default Limma selectModel output

I've been given a large RNAseq dataset on a biological process which I know little about. There is a ton of associated sample data, and I've built a ton of different linear models that include variables that correlate with my process of interest in either the literature or some basic, exploratory analyses.

I'd like to figure out which of these models is worth playing around with while I learn more about the new field. My prior RNAseq analysis experience was all in really simplistic experimental systems, so I'm wondering how to appropriately apply tools like Limma's selectModel AIC/BIC criteria.

Specifically, what is pref? How's it calculated? How should I interpret it?
The two folks that I've spoken to have opposing opinions: One thinks that pref is a broad information criterion score and that I should pick whatever model has the lowest value in this table. The other has told me that the values indicate the number of genes fit well by each model and that I should select the model with the highest value in this table.

As per the manual:
IC
matrix of information criterion scores, rows for probes and columns for models
pref
factor indicating the model with best (lowest) information criterion score
vsmuir is offline   Reply With Quote
Old 10-11-2016, 10:00 PM   #2
Dario1984
Senior Member
 
Location: Sydney, Australia

Join Date: Jun 2011
Posts: 166
Default

Quote:
The other has told me that the values indicate the number of genes fit well by each model and that I should select the model with the highest value in this table.
That is the correct explanation.
Dario1984 is offline   Reply With Quote
Reply

Tags
limma, linear model

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 10:12 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO