Combinations of factors not observed, non-full rank design matrix. How to explain to investigator?
So I am helping someone with a differential expression analysis, there are only 10 samples, two variables with two levels each. Let's say Sex (M/F) and Age (Old/Young). They originally wanted to model: ~ sex + age + sex*age
However Sex = F & Age = Young does not exist in
the data (no sample with that combination observed), so model matrix is not full rank and DESEQ model can't be specified.
I warned them of this and their solution was to concat the Sex and Age variables to a new var (let's just say V3) and run the model with only ~
V3
I know this technically works... (as in the design matrix is full rank), but I also know it isn't a great idea, basically bc we are extrapolating and are now unable to make any claims about M vs For Old vs. Young.
Any tips on how to explain this to the investigators?
Answer
Answers can only be viewed under the following conditions:
- The questioner was satisfied with and accepted the answer, or
- The answer was evaluated as being 100% correct by the judge.
Kav10
1.9K
The answer is accepted.
Join Matchmaticians Affiliate Marketing
Program to earn up to a 50% commission on every question that your affiliated users ask or answer.
- answered
- 391 views
- $8.00
Related Questions
- How do we describe an intuitive arithmetic mean that gives the following? (I can't type more than 200 letters)
- Variance of Autoregressive models, AR(1)
- "Estimate of standard error of the proportion"
- Finding the probability that a roughly normal distributed will have the highest value among multiple curves
- applied probability
- Elements in a set
- How to use bootstrap techniques to criticise a linear model?
- Currently studying a grad level Statistical Inference course. I'd just like some clarification regarding how to obtain the Rao-Cramer Lower Bound for a statistic.
The offered bounty is low for the level of the question.