Эффективно преобразовать стебель и лист в вектор в R - PullRequest
0 голосов
/ 05 июня 2018

Мне нужно проверить сводную статистику (среднее значение, стандартное отклонение и т. Д.) Для многих графиков стеблей и листьев, поэтому я написал несколько функций, чтобы попытаться преобразовать график стеблей и листьев в вектор, так как статистику по векторам легко получить.в R.

График ствола и листа может быть введен в виде матрицы или фрейма данных, где каждая строка является строкой."|"Символ представляет разделитель для десятичного разряда.Например, график ствола и листа ниже

100 | 9
102 | 601
104 | 0678
106 | 5
108 | 649
110 | 3857
112 | 56
114 | 29

может быть введен как

> example.stem = rbind("100|9", "102|601", "104|0678", "106|5", "108|649", "110|3857", "112|56", "114|29")

Две мои функции, которые выполняют преобразование этого графика ствола и листа:

## Convert a single row into a vector
> convert.row = function(current){

  temp.split = as.vector(strsplit(current, split="|", fixed=TRUE)[[1]])

  int = temp.split[1]

  dec = temp.split[2]
  dec = (strsplit(dec, ""))[[1]]

  temp.string = NULL

  for(i in 1:length(dec)){
    temp.string[i] = paste(int, dec[i], sep=".")
  }

  result = as.numeric(temp.string)
  return(result)
  }


## Convert matrix or dataframe with a stem and leaf plot into a vector
> stem.to.vec = function(df){  

  df = data.frame(df, stringsAsFactors = F)  

  result.vec = NULL

  for(i in 1:nrow(df)){
    current = df[i, ]
    result.vec = c(result.vec, convert.row(current))
  }

  return(result.vec)
  }    

Мы можем убедиться, что это работает, потому что мы знаем решение:

> solution = c(100.9, 102.6,102.0,102.1,104.0,104.6,104.7,104.8,106.5,108.6,108.4,108.9,110.3,110.8,110.5,110.7,112.5,112.6, 114.2, 114.9)
> stem.to.vec(example.stem) == solution

Хотя это решение работает, оно не элегантно и не эффективно.Мы преобразуем фрейм матрицы / данных со строками в числовое значение, затем обратно в строку и затем снова в числовое значение.Поэтому работа с очень большими участками стеблей и листьев может быть медленной.

Кто-нибудь может предложить лучшее и более эффективное решение с меньшим количеством конверсий?

Ответы [ 3 ]

0 голосов
/ 05 июня 2018

Метод разделения и вставки:

Разделение на элементы списка, а затем разделение второго элемента элементов списка.Наконец, вставьте два вектора в элемент списка.

x <- sapply(strsplit(example.stem, "[|]"), 
            function(x) { paste(x[1], unlist(strsplit(x[2], "")), sep= ".") })

as.numeric(unlist(x))                           

# [1] 100.9 102.6 102.0 102.1 104.0 104.6 104.7 104.8 106.5 108.6 108.4 108.9 
# [13] 110.3 110.8 110.5 110.7 112.5 112.6 114.2 114.9
0 голосов
/ 05 июня 2018

Вот вложенная lapply версия, которая выполняет все объединения в виде строк, а затем преобразует выходное значение в числовое значение:

out <- unlist(lapply(strsplit(example.stem, "[|]"), function(x){
  lapply(x[2], function(y){
    as.numeric(paste(x[1], unlist(strsplit(y, "")), sep = "."))
  })
}))

> identical(solution, out)
[1] TRUE
0 голосов
/ 05 июня 2018

Это далеко не красиво, но я думаю, что вам все равно придется выполнить какое-то обратное преобразование.

Используйте read.table для ввода данных, а затем разделите правую часть на10 и добавьте к каждому значению с левой стороны.

out <- read.table(text=example.stem, sep="|", colClasses=c("numeric","character"))
res <- unlist(Map(`+`, out$V1, lapply(strsplit(out$V2,""), function(x) as.numeric(x)/10)))
res
# [1] 100.9 102.6 102.0 102.1 104.0 104.6 104.7 104.8 106.5 108.6 108.4 108.9
# [13] 110.3 110.8 110.5 110.7 112.5 112.6 114.2 114.9

identical(solution,res)
#[1] TRUE
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...