2

これは、使用後に列をより完全にマージするために書いているスクリプトの一部ですmerge()。両方のデータセットに同じ名前の列がある場合、列とがmerge()得られます。このデータをまとめて、不要な列を削除するスクリプトを作成しました (これは、場合に警告を表示するために追加した列であるとになります。また、データセット内の不要な手動操作を減らすために、 に名前を変更したいと思います。に名前を変更できませんでした。詳細については、コードを参照してください。column.xcolumn.y
column.ycolumn.x_errordat$column.x != dat$column.y)
column.xcolumncolumn.xcolumn

datを行うことで得られるdat = merge(data1,data2, by= "ID", all.x=TRUE)

#obtain a list of double columns
dubbelkol = cbind()
sorted = sort(names(dat))
for(i in as.numeric(1:length(names(dat)))) {
  if(grepl(".x",sorted[i])){
    if (grepl(".y", sorted[i+1]) && (sub(".x","",sorted[i])==sub(".y","",sorted[i+1]))){
      dubbelkol = cbind(dubbelkol,sorted[i],sorted[i+1])
    } 
  }  
}

#Check data, fill in NA in column.x from column.y if poss
temp = cbind()
for (p in as.numeric(1:(length(dubbelkol)-1))){
  if(grepl(".x",dubbelkol[p])){
    dat[dubbelkol[p]][is.na(dat[dubbelkol[p]])] = dat[dubbelkol[p+1]][is.na(dat[dubbelkol[p]])]
    temp = (dat[dubbelkol[p]] != dat[dubbelkol[p+1]])
    colnames(temp) = (paste(dubbelkol[p],"_error", sep=""))
    dat[colnames(temp)] = temp
  }
}
#If every value in "column.x_error" is TRUE or NA, delete "column.y" and "column.x_error"
#Rename "column.x" to "column"
#from here until next comment everything works
droplist= c()
for (k in as.numeric(1:length(names(dat)))) {
  if (grepl(".x_error",colnames(dat[k]))) {
    if (all(dat[k]==FALSE, na.rm = TRUE)) {
      droplist = c(droplist,colnames(dat[k]), sub(".x_error",".y",colnames(dat[k])))
#the next line doesnt work, it's supposed to turn the .x column back to "" before the .y     en .y_error columns are dropped.
      colnames(dat[sub(".x_error",".x",colnames(dat[k]))])= paste(sub(".x_error","",colnames(dat[k])))
    }
  }
}
dat = dat[,!names(dat) %in% droplist]

paste(sub(".x_error","",colnames(dat[k])))"BNR"うまくいきますがcolnames(...) = ...、の列名は変更されませんdat

何がうまくいかないのですか?

data1
+----+-------+
| ID | BNR   | 
+----+-------+
|  1 | 123   | 
|  2 | 234   |
|  3 | NA    | 
|  4 | 456   | 
|  5 | 677   |
|  6 | NA    | 
+----+-------+

data2
+----+-------+
| ID | BNR   | 
+----+-------+
|  1 | 123   | 
|  2 | 234   |
|  3 | 345   | 
|  4 | 456   | 
|  5 | 677   |
|  6 | NA    | 
+----+-------+
dat
+----+-------+-------+-----------+
| ID | BNR.x | BNR.y |BNR.x_error|
+----+-------+-------+-----------+
|  1 | 123   | NA    |FALSE      |
|  2 | 234   | 234   |FALSE      |
|  3 | NA    | 345   |FALSE      |
|  4 | 456   | 456   |FALSE      |
|  5 | 677   | 677   |FALSE      |
|  6 | NA    | NA    |NA         |
+----+-------+-------+-----------+

desired output
+----+-------+
| ID | BNR   | 
+----+-------+
|  1 | 123   |
|  2 | 234   | 
|  3 | 345   | 
|  4 | 456   | 
|  5 | 677   | 
|  6 | NA    | 
+----+-------+
4

1 に答える 1

2

交換することをお勧めします:

sub(".x_error",".x",colnames(dat[k]))]

と:

sub("\\.x_error", "\\.x", colnames(dat[k]))] 

実際の を置き換えたい場合.。でエスケープする必要が.あり\\.ます。正規表現の A.any character.

さらに良いのは、次のように置き換え.ているからです.

sub("x_error", "x", colnames(dat[k]))] 

(または)_error以外にない場合はx_error、単に:

sub("_error", "", colnames(dat[k]))] 

編集:問題は、データ形式が左右に追加の列をロードしているように見えることです。最初に必要な列を選択してから結合できます。

d1 <- read.table(textConnection("| ID | BNR   | 
|  1 | 123   | 
|  2 | 234   |
|  3 | NA    | 
|  4 | 456   | 
|  5 | 677   |
|  6 | NA    |"), sep = "|", header = TRUE, stringsAsFactors = FALSE)[,2:3]

d1$BNR <- as.numeric(d1$BNR)

d2 <- read.table(textConnection("|  1 | 123   | 
|  2 | 234   |
|  3 | 345   | 
|  4 | 456   | 
|  5 | 677   |
|  6 | NA    |"), header = FALSE, sep = "|", stringsAsFactors = FALSE)[,2:3]

names(d2) <- c("ID", "BNR")
d2$BNR <- as.numeric(d2$BNR)

# > d1
#   ID BNR
# 1  1 123
# 2  2 234
# 3  3  NA
# 4  4 456
# 5  5 677
# 6  6  NA

# > d2
#   ID BNR
# 1  1 123
# 2  2 234
# 3  3 345
# 4  4 456
# 5  5 677
# 6  6  NA

dat <- merge(d1, d2, by="ID", all=T)
> dat

#   ID BNR.x BNR.y
# 1  1   123   123
# 2  2   234   234
# 3  3    NA   345
# 4  4   456   456
# 5  5   677   677
# 6  6    NA    NA

# replace all NA values in x from y
dat$BNR.x <- ifelse(is.na(dat$BNR.x), dat$BNR.y, dat$BNR.x)

# now remove y
dat$BNR.y <- null
于 2013-02-15T12:58:14.750 に答える