定義因素，其水平依賴於另一個變量

Question

是這樣的模擬數據：

set.seed(20120220) 
x <- c(rep("a", 4), rep("b", 4)) 
y <- c(sample(c(1, 2), 8, replace = TRUE)) 
z <- data.frame(cbind(x, y)) 

數據幀z看起來就像這樣：

x y 
1 a 1 
2 a 1 
3 a 1 
4 a 2 
5 b 2 
6 b 1 
7 b 2 
8 b 2 

我想運行到factor(z$y, levels = 1:2, labels = c("alpha", "beta"))一個類似於，但我不不想讓每個1變成alpha，並且每個2變成beta。我希望這發生只有爲x = a。如果x = b，我想1變成gamma和2變成delta。

換句話說，我希望我的數據幀，看起來像這樣：

x y 
1 a alpha 
2 a alpha 
3 a alpha 
4 a beta 
5 b delta 
6 b gamma 
7 b delta 
8 b delta 

這是我想出了這麼遠：

for (i in 1:nrow(z)) { 
    if (z$x[i] == "a") 
    z$y[i] <- factor(z$y[i], levels = 1:2, labels = c("alpha", "beta")) 
    else 
    z$y[i] <- factor(z$y[i], levels = 1:2, labels = c("gamma", "delta")) 
} 

但它給了我好幾個警告信息（爲每個i）是這樣的：

Warning messages: 
1: In `[<-.factor`(`*tmp*`, i, value = c(NA, 1L, 1L, 2L, 2L, 1L, 2L, : 
    invalid factor level, NAs generated 

然後，當我打電話z再次，數據幀亂七八糟，每y已被製作成<NA>。

我敢打賭，這是一個簡單的解決方案，但我一直在嘗試小時幾種方法無濟於事。我的頭即將爆炸！幫幫我！

Answer 1

> z$ynew <- ifelse(z$x == "a", ifelse(z$y==1, "alpha", "beta"), 
           ifelse(z$y==1, "delta", "gamma")) 
> z 
    x y ynew 
1 a 1 alpha 
2 a 1 alpha 
3 a 1 alpha 
4 a 2 beta 
5 b 2 gamma 
6 b 1 delta 
7 b 2 gamma 
8 b 2 gamma 

（我猜我換你的三角洲和伽瑪的。如果你想「ynew」是一個因素則剛：？z$ynew <- factor(z$ynew)

Answer 2

有關使用合併什麼

# define x and y to 'alpha', 'beta' etc. correspondences 
# (it's just one row for each possible factor) 
auxDf <- data.frame(x = c('a',  'a', 'b',  'b' ), 
        y = c(1,  2,  1,  2 ), 
        newy= c('alpha', 'beta', 'gamma', 'delta')) 

# merge the 2 data.frame getting a new data.frame with the factors column 
newDf <- merge(z,auxDf) 
newDf 

Answer 3

這裏有一個額外的步驟使以前的答案更快 - 您可以使用'unique'來抽出數據幀中的所有唯一組合。

auxDf=unique(z) 
auxDf$newy=c('alpha','beta','gamma','delta') 

然後，在以前的帖子

newDf <- merge(z,auxDf) 
newDf 

Answer 4

我已經成功地拿出這樣的作品，即使它是相當混亂的解決方案。

首先，對於每個x

z1 <- subset(z, x == "a") 
z2 <- subset(z, x == "b") 

創建數據幀z的子集然後，應用factor()到每個子集：

z1$y <- factor(z1$y, levels = 1:2, labels = c("alpha", "beta")) 
z2$y <- factor(z2$y, levels = 1:2, labels = c("gamma", "delta")) 

最後，團聚子集成原始對象。

z <- rbind(z1, z2)