我有一个数据框,名为trip:
> head(trip.mutations)
Ref.y Variant.y
1 T C
2 G C
3 A C
4 T C
5 C A
6 G A
我想添加第三列 mutType,它遵循以下规则:
for (i in 1:nrow(trip)) {
if(trip$Ref.y=='G' & trip$Variant.y=='T'|trip$Ref.y=='C' & trip$Variant.y=='A') {
trip[i, 'mutType'] <- "G:C to T:A"
}
else if(trip$Ref.y=='G' & trip$Variant.y=='C'|trip$Ref.y=='C' & trip$Variant.y=='G') {
trip[i, 'mutType'] <- "G:C to C:G"
}
else if(trip$Ref.y=='G' & trip$Variant.y=='A'|trip$Ref.y=='C' & trip$Variant.y=='T') {
trip[i, 'mutType'] <- "G:C to A:T"
}
else if(trip$Ref.y=='A' & trip$Variant.y=='T'|trip$Ref.y=='T' & trip$Variant.y=='A') {
trip[i, 'mutType'] <- "A:T to T:A"
}
else if(trip$Ref.y=='A' & trip$Variant.y=='G'|trip$Ref.y=='T' & trip$Variant.y=='C') {
trip[i, 'mutType'] <- "A:T to G:C"
}
else if(trip$Ref.y=='A' & trip$Variant.y=='C'|trip$Ref.y=='T' & trip$Variant.y=='G') {
trip[i, 'mutType'] <- "A:T to C:G"
}
}
但是我遇到了错误:
Warning messages:
1: In if (trip$Ref.y == "G" & trip$Variant.y == "T" | trip$Ref.y == ... :
the condition has length > 1 and only the first element will be used
我不认为我的逻辑语句应该产生向量,但也许我漏掉了什么。trip$mutType 应该会变成这样:
mutType
A:T to G:C
G:C to C:G
A:T to C:G
A:T to G:C
G:C to T:A
G:C to A:T
有人能发现这里的问题吗?我需要使用||而不是|吗?