与因子变量输出

Question

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 

Answer 1

处理因子变量A1和A2以字符变量形式给出适当的输出。

hdl.dir$new_A1 <- ifelse(hdl.dir$HDL_Effect<0, as.character(hdl.dir$A2), as.character(hdl.dir$new_A1))

Answer 2

你可以试试行/列的索引。即使没有明确地将“因素”转换为“字符”类，这也应该起作用。

df$new_a1 <- df[3:4][cbind(1:nrow(df),(df$Effect<=0)+1L)]