与因子变量输出
问题描述:
ifelse语句我有一个数据帧df
与因子变量输出
SNP Effect A1 A2
8 rs12145743 -0.0203 T G
31 rs4650994 -0.0210 A G
26 rs2602836 0.0192 A G
34 rs4976033 0.0215 A G
10 rs12525163 -0.0215 T C
1 rs10282707 -0.0250 T C
33 rs4917014 -0.0222 T G
23 rs2293889 -0.0312 T G
32 rs4871137 -0.0209 T G
43 rs970548 -0.0258 A C
27 rs2923084 0.0256 A G
24 rs2303975 0.0279 A G
36 rs499974 -0.0263 A C
35 rs4983559 -0.0197 A G
13 rs16942887 0.0831 A G
16 rs17695224 -0.0290 A G
18 rs181362 -0.0379 T C
我创建一个新的变量new_A1
将取决于Effect
的价值是不同的。
如果Effect
> 0,new_A1
应该包含值A1
。
如果Effect
< 0,new_A1
应该包含值A2
。
我试图做到这一点使用下面的ifelse语句:
df$new_A1 <- ifelse(df$Effect>0, df$A1, df$A2)
但是,并不如预期的输出。 new_A1
包含1S和2S,而不是维持原有的变量因素值A1
和A2
SNP Effect A1 A2 new_a1
8 rs12145743 -0.0203 T G 2
31 rs4650994 -0.0210 A G 2
26 rs2602836 0.0192 A G 1
34 rs4976033 0.0215 A G 1
10 rs12525163 -0.0215 T C 1
1 rs10282707 -0.0250 T C 1
33 rs4917014 -0.0222 T G 2
23 rs2293889 -0.0312 T G 2
32 rs4871137 -0.0209 T G 2
43 rs970548 -0.0258 A C 1
27 rs2923084 0.0256 A G 1
24 rs2303975 0.0279 A G 1
36 rs499974 -0.0263 A C 1
35 rs4983559 -0.0197 A G 2
13 rs16942887 0.0831 A G 1
16 rs17695224 -0.0290 A G 2
18 rs181362 -0.0379 T C 1
答
处理因子变量A1
和A2
以字符变量形式给出适当的输出。
hdl.dir$new_A1 <- ifelse(hdl.dir$HDL_Effect<0, as.character(hdl.dir$A2), as.character(hdl.dir$new_A1))
答
你可以试试行/列的索引。即使没有明确地将“因素”转换为“字符”类,这也应该起作用。
df$new_a1 <- df[3:4][cbind(1:nrow(df),(df$Effect<=0)+1L)]
将列(A1,A2)转换为字符类。现在,这将是因子类。然后'ifelse'可以工作。 – akrun