1

我想比较 Q 方法和 L 方法，我考虑了 2 个不同的对比（最后），但我不确定哪个是正确的？

有 2 种不同的方法（Q 和 L），每种方法有 2 个生物复制（L4、L6-L8 和 Q3、Q5-Q7），每个生物复制有 2 个技术复制。如下：

设计

                    biological_replicate   method
 L4_rep1                              L4       L
 L4_rep2                              L4       L
 L6_L8_rep1                        L6_L8       L
 L6_L8_rep2                        L6_L8       L
 Q3_rep1                              Q3       Q
 Q3_rep2                              Q3       Q
 Q5_Q7_rep1                        Q5_Q7       Q
 Q5_Q7_rep2                        Q5_Q7       Q                       

design$biological_replicate <- factor(design$biological_replicate, levels =  c("L4","L6_L8", "Q3", "Q5_Q7"))

design$method <- factor(design$method, levels = c("L", "Q"))

Group <- factor(paste(design$biological_replicate,design$method,sep="."))

design<- cbind(design,Group)

                    biological_replicate method   Group
L4_rep1                              L4      L    L4.L
L4_rep2                              L4      L    L4.L
L6_L8_rep1                        L6_L8      L L6_L8.L
L6_L8_rep2                        L6_L8      L L6_L8.L
Q3_rep1                              Q3      Q    Q3.Q
Q3_rep2                              Q3      Q    Q3.Q
Q5_Q7_rep1                        Q5_Q7      Q Q5_Q7.Q
Q5_Q7_rep2                        Q5_Q7      Q Q5_Q7.Q

design.matrix <- model.matrix(~0+Group,design)

colnames(design.matrix) <- levels(Group)

design.matrix

L4.L L6_L8.L Q3.Q Q5_Q7.Q
L4_rep1       1       0    0       0
L4_rep2       1       0    0       0
L6_L8_rep1    0       1    0       0
L6_L8_rep2    0       1    0       0
Q3_rep1       0       0    1       0
Q3_rep2       0       0    1       0
Q5_Q7_rep1    0       0    0       1
Q5_Q7_rep2    0       0    0       1
attr(,"assign")
[1] 1 1 1 1
attr(,"contrasts")
attr(,"contrasts")$Group
[1] "contr.treatment"

my.contrasts_1 <- makeContrasts(QvsL = (Q3.Q+Q5_Q7.Q)/2-(L4.L+L6_L8.L)/2, levels = design.matrix)

my.contrasts_2 <- makeContrasts(QvsL = (Q3.Q+Q5_Q7.Q)-(L4.L+L6_L8.L), levels = design.matrix)
4

1 回答 1

1

首先，需要对技术重复进行求和以进行差异表达分析。您只需将两个计数相加即可。

对技术复制求和后，您的数据应如下所示：

Sample <- c("L4", "L6_L8", "Q3", "Q5_Q7")
Method <- c("L", "L", "Q", "Q")
df <- data.frame(Sample, Method)

  Sample Method
1     L4      L
2  L6_L8      L
3     Q3      Q
4  Q5_Q7      Q

现在设计矩阵很简单：

design.matrix <- model.matrix(~0 + Method)
colnames(design.matrix) <- c("L","Q")
  MethodL MethodQ
1       1       0
2       1       0
3       0       1
4       0       1

对比也变得很容易生成：

my_contrast <- makeContrasts(L-Q, levels = design.matrix)

Contrasts
Levels L - Q
     L     1
     Q    -1

这不是您的问题的一部分，但是假设上表代表您的实际数据，您的生物复制数量绝对最少。不要将技术复制视为生物复制。您将增加您的类型 1 错误。

于 2018-04-11T02:17:52.210 回答