1

我对 R 相当陌生，尽管我经历了许多教程，但在索引方面存在问题。例如，在以下情况下，我有一个带有相应重要 p 值的探针名称列表，如下所示

      unique.et.                                            pBH
1     AFFX-r2-P1-cre-5_at                                   1.646105e-14
2     AFFX-r2-P1-cre-3_at                                   1.752547e-14
3     AFFX-r2-Ec-bioD-5_at                                  2.426254e-14
127            267584_at                                    1.820809e-14
149            267562_at                                    1.646105e-14
200            267511_at                                    1.880226e-14

下面是原始矩阵，我必须从中仅提取具有显着 p 值的探针，以便获得信号强度值。

AFFX-r2-P1-cre-2_at 12.16271 12.70304 12.16271 12.70304 12.16271
AFFX-r2-P1-cre-3_at 12.70304 12.28280 12.35039 12.38397 12.36304
AFFX-r2-P1-cre-5_at 12.28280 12.35039 12.38397 12.36304 12.16271

因此，如果我考虑名为 AFFX-r2-P1-cre-5_at 且 p 值为 1.646105e-14 的探针，我将需要从其他探针列表中提取相同的探针以及强度值。我试过但做不到。

我还想问大家是否有带有清晰示例的教程，以便我理解这个概念。

4

1 回答 1

1 
#data.frame with significant cases
df1 <- read.table(text="      unique.et.                                            pBH
1     AFFX-r2-P1-cre-5_at                                   1.646105e-14
2     AFFX-r2-P1-cre-3_at                                   1.752547e-14
3     AFFX-r2-Ec-bioD-5_at                                  2.426254e-14
127            267584_at                                    1.820809e-14
149            267562_at                                    1.646105e-14
200            267511_at                                    1.880226e-14",header=T,stringsAsFactors=FALSE)

#all data
df2 <- read.table(text="AFFX-r2-P1-cre-2_at 12.16271 12.70304 12.16271 12.70304 12.16271
AFFX-r2-P1-cre-3_at 12.70304 12.28280 12.35039 12.38397 12.36304
AFFX-r2-P1-cre-5_at 12.28280 12.35039 12.38397 12.36304 12.16271",header=F,stringsAsFactors=FALSE)

#select cases of df2 that occur in df1
df2[df2$V1 %in% df1$unique.et,]
#                   V1       V2       V3       V4       V5       V6
#2 AFFX-r2-P1-cre-3_at 12.70304 12.28280 12.35039 12.38397 12.36304
#3 AFFX-r2-P1-cre-5_at 12.28280 12.35039 12.38397 12.36304 12.16271
于 2012-08-14T07:39:36.233 回答