Search
Question: Error in checkFullRank(modelMatrix) : the model matrix is not full rank, ?
0
gravatar for bijaya.malla
11 weeks ago by
bijaya.malla0 wrote:

I used following command to filter out 3 groups (NC, BF, CF) out of 5 groups while remaining two groups are "AF" and "DF".

design_NCaF <- filter(design, Groups1!="NC",Groups1!="BF",Groups1!="CF")

However, at following command

NCaF <- DESeqDataSetFromMatrix(data_max_woNCaF, colData=design_NCaF, design = ~Groups1 + Patient)

It gives me following error message 

factor levels were dropped which had no samples
Error in checkFullRank(modelMatrix) : 
  the model matrix is not full rank, so the model cannot be fit as specified.
  One or more variables or interaction terms in the design formula are linear
  combinations of the others and must be removed.

Is it because of remaining two groups have different sample size (AF=6 and DF=3) ?

 

ADD COMMENTlink modified 8 weeks ago by Michael Love15k • written 11 weeks ago by bijaya.malla0
0
gravatar for Michael Love
11 weeks ago by
Michael Love15k
United States
Michael Love15k wrote:

You need to droplevels() from the factor after filtering, e.g. if you filter out samples from 'group':

dds$group <- droplevels(dds$group)
ADD COMMENTlink written 11 weeks ago by Michael Love15k
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 2.2.0
Traffic: 287 users visited in the last hour