2012-02-09 75 views
4

我知道我的问题很简单,但不适合我。这里是小数据集。长,如果其他循环和重新编码在R

mark1 <- c("AB", "BB", "AB", "BB", "BB", "AB", "--", "BB") 
mark2 <- c("AB", "AB", "AA", "BB", "BB", "AA", "--", "BB") 
mark3 <- c("BB", "AB", "AA", "BB", "BB", "AA", "--", "BB") 
mark4 <- c("AA", "AB", "AA", "BB", "BB", "AA", "--", "BB") 
mark5 <- c("AB", "AB", "AA", "BB", "BB", "AA", "--", "BB") 
mark6 <- c("--", "BB", "AA", "BB", "BB", "AA", "--", "BB") 
mark7 <- c("AB", "--", "AA", "BB", "BB", "AA", "--", "BB") 
mark8 <- c("BB", "AA", "AA", "BB", "BB", "AA", "--", "BB") 
mymark <- data.frame (mark1, mark2, mark3, mark4, mark5, mark6, mark7, mark8) 
tmymark <- data.frame (t(mymark)) 
names (tmymark) <- c("P1", "P2","I1", "I2", "I3", "I4", "KL", "MN") 

因此数据集变为:

 P1 P2 I1 I2 I3 I4 KL MN 
mark1 AB BB AB BB BB AB -- BB 
mark2 AB AB AA BB BB AA -- BB 
mark3 BB AB AA BB BB AA -- BB 
mark4 AA AB AA BB BB AA -- BB 
mark5 AB AB AA BB BB AA -- BB 
mark6 -- BB AA BB BB AA -- BB 
mark7 AB -- AA BB BB AA -- BB 
mark8 BB AA AA BB BB AA -- BB 

欲分类MARK1:基于所述P1和P2对比8,并提供一个代码,这将使得一个新的变量:

loctype <- NULL 

if (tmymark$P1 == "AB" & tmymark$P2 == "AB"){ 
     loctype = "<hkxhk>" 
     } else { 
if (tmymark$P1== "AB" & tmymark$P2 == "BB") { 
     loctype = "<lmxll>" 
     } else { 
     if (tmymark$P1 == "AA" & tmymark$P2 == "AB") { 
     loctype = "<nnxnp>" 
     } else { 
     if (tmymark$P1 == "AA" & tmymark$P2 == "BB") { 
     loctype = "MN" 
     } else { 
     if (tmymark$P1 == "BB" & tmymark$P2 == "AA"){ 
     loctype = "MN" 
     } else { 
     if (tmymark$P1 == "--" & tmymark$P2 == "AA"){ 
     loctype = "NR" 
     } else { 
if (tmymark$P1 == "AA" & tmymark$P2 == "--"){ 
      loctype = "NR" 
} else { 
    cat ("error wrong input in P1 or P2") 
    }} }}}}} 

在这里,我试图做它比较P1和P2值并生成一个新的变量。 for example,if tmymark $ P1 ==“AB”& tmymark $ P2 ==“AB”loctype应该是“”。如果不是第二个条件将是应用程序等。

这是我的错误消息。

Warning messages: 
1: In if (tmymark$P1 == "AB" & tmymark$P2 == "AB") { : 
    the condition has length > 1 and only the first element will be used 
2: In if (tmymark$P1 == "AB" & tmymark$P2 == "BB") { : 
    the condition has length > 1 and only the first element will be used 

一旦LOCTYPE向量生成我要重新编码与该变量的信息tmymark:

tmymark1 <- data.frame (loctype, tmymark)  
require(car) 
for(i in 2:length(tmymark)){ 

     if (loctype = "<hkxhk>") { 
     tmymark[[i]] <- recode (x, "AB" = "hk", "BA" = "hk", "AA" = "hh", "BB" = "kk") 
     } else { 
     if (loctype = "<lmxll>") { 
     tmymark[[i]] <- recode ((x, "AB" = "lm", "BA" = "lm", "AA" = "--", "BB" = "kk") 
     } else { 

     if (loctype = "<nnxnp>") { 
     tmymark[[i]] <- recode ((x, "AB" = "np", "BA" = "np", "AA" = "nn", "BB" = "--") 
      } else { 
     if (loctype = "MN") { 
     tmymark[[i]] <- "--" 
     } esle { 
     if (loctype = "NR") { 
     tmymark[[i]] <- "NA" 
     } else { 
     cat ("error wrong input code") 
     } } }}} 

上午我在正确的轨道?

编辑:预期输出

 loctype P1 P2 I1 I2 I3 I4 KL MN 
mark1 <lmxmm> lm mm lm mm mm lm -- mm 
mark2 <hkxhk> hk hk hh kk kk hh -- kk 
mark3 <nnxnp> nn np nn -- -- nn -- -- 
and so on 

回答

1

match肯定是要走的路我会想办法让两个数据帧作为键,这样的:

key <- data.frame(
      P1=c("AB", "AB", "AA", "AA", "BB", "--", "AA"), 
      P2=c("AB", "BB", "AB", "BB", "AA", "AA", "--"), 
      loctype=c("<hkxhk>", "<lmxll>", "<nnxnp>", "MN", "MN", "NR", "NR")) 

key2 <- cbind(
    `<hkxhk>` = c("hk","hk","hh","kk"), 
    `<lmxll>` = c("lm", "lm", "--", "kk"), 
    `<nnxnp>` = c("np", "np", "nn", "--"), 
    MN = rep("--", 4), 
    NR = rep("NA", 4)) 
rownames(key2) = c("AB","BA", "AA", "BB") 

然后用matchkey1得到loctype(贾斯汀也建议) ,以及在key2的rownames和列上以获得期望的替换,使用矩阵索引来从密钥中获得期望的值。

loctype <- key$loctype[match(with(tmymark, paste(P1, P2, sep="\b")), 
          with(key, paste(P1, P2, sep="\b")))] 
ii <- match(as.vector(as.matrix(tmymark)), rownames(key2)) 
jj <- rep(match(loctype, colnames(key2)), nrow(tmymark)) 
out <- as.data.frame(matrix(key2[cbind(ii,jj)], nrow=nrow(tmymark))) 
colnames(out) <- colnames(tmymark) 
rownames(out) <- rownames(tmymark) 
out$loctype <- loctype 

结果然后看起来像这样,其中缺少的值是因为我没有这些组合的值在我的键。

> print(out, na="") 
     P1 P2 I1 I2 I3 I4 KL MN loctype 
mark1 lm kk lm kk kk lm kk <lmxll> 
mark2 hk hk hh kk kk hh kk <hkxhk> 
mark3         
mark4 nn np nn -- -- nn -- <nnxnp> 
mark5 hk hk hh kk kk hh kk <hkxhk> 
mark6         
mark7         
mark8 -- -- -- -- -- -- --  MN 
1

发生的第一错误,因为要想(计算结果为一个或表达)的单个逻辑值。您可以使用ifelse(),而不是这是一个“矢量” if

ifelse(tmymark$P1 == "AB" & tmymark$P2 == "AB", loctype = "<hkxhk>", else clauses...) 

为了避免长if() else()结构(或ifelse()因为它是),你可以使用匹配。使P1和P2的您的预计组合的数据帧和期望LOCTYPE添加一列:

matches <- data.frame(p1p2 = c('AB AB', 'AB BB', 'AA AB', 'AA BB', 'BB AA', '-- AA', 'AA --'), 
         loctype = c('<hkxhk>', '<lmxll>', '<nnxnp>', 'MN', 'MN', 'NR', 'NR')) 
loctype <- matches$loctype[match(paste(tmymark$P1, tmymark$P2), matches$p1p2),] 

第二部分可以做多种方式,但我画一个空白上的“整洁。整洁”一个