Неверное подмножество на основе colSums - PullRequest
Новая
       3

Неверное подмножество на основе colSums

0 голосов
/ 24 августа 2018

У меня большой фрейм данных (100 x 1748), ниже уменьшенная версия 9x10: 

                Fraction Treatment Time Replicate        A        B         C        D         E         F
LL10.5AT           T        LL      10.5    A        11.11428 11.82154 10.445625 8.849699 10.373386  9.109676
LL10.5BT           T        LL      10.5    B        12.17890 11.01224 11.720548 9.405390 10.206708 10.653205
LL10.5CT           T        LL      10.5    C        10.80697 11.19782 11.175291 8.305949  9.696153  8.791403
OL10.5AT           T        OL      10.5    A        10.46481 10.81123  9.975277 7.783538  9.784773  8.640531
OL10.5BT           T        OL      10.5    B        10.75621 10.76371 10.625745 7.592059  9.820686  8.760861
OL10.5CT           T        OL      10.5    C        12.00054 11.02080 11.615536 8.903105  9.963635 10.547791
HL10.5AT           T        HL      10.5    A        10.87092 11.45102 10.780183 6.422136 10.424391  9.489396
HL10.5BT           T        HL      10.5    B        12.12334 11.29960 11.541679 9.774041  9.563639 10.532936
HL10.5CT           T        HL      10.5    C 10.21460 10.64746  9.886603 7.834040  9.828347  8.261546

Я хочу сделать его поднабором, чтобы он содержал только столбцы, для которых сумма> 100. Я использую следующий код 

dt.sub <- dt[,colSums(dt[,5:ncol(dt)]) > 100]

но я все еще получаю столбец с суммами <100. Я создал вектор для проверки: </p> 

z<-colSums(dt.sub[,5:ncol(dt.sub)])

а это его хвост: 

> tail(z)
    501         502         503         504         505
107.9368630  90.6337275   0.8724593   0.8724593   1.3497445   1.3497445

Спасибо за помощь, Кася

1 Ответ

0 голосов
/ 24 августа 2018

Вы неправильно индексируете столбцы. 

df[, c(1:4, which(colSums(df[, 5:ncol(df)]) > 100) + 4)]
#         Fraction Treatment Time Replicate        A        B
#LL10.5AT     TRUE        LL 10.5         A 11.11428 11.82154
#LL10.5BT     TRUE        LL 10.5         B 12.17890 11.01224
#LL10.5CT     TRUE        LL 10.5         C 10.80697 11.19782
#OL10.5AT     TRUE        OL 10.5         A 10.46481 10.81123
#OL10.5BT     TRUE        OL 10.5         B 10.75621 10.76371
#OL10.5CT     TRUE        OL 10.5         C 12.00054 11.02080
#HL10.5AT     TRUE        HL 10.5         A 10.87092 11.45102
#HL10.5BT     TRUE        HL 10.5         B 12.12334 11.29960
#HL10.5CT     TRUE        HL 10.5         C 10.21460 10.64746

Объяснение: which(colSums(df[, 5:ncol(df)]) > 100) возвращает индексы в пределах df[, 5:ncol(df)] (не в пределах df!), Где сумма столбцов равна >100;затем мы добавляем 4 (поскольку мы начали с индекса 5) и включаем столбцы с 1 по 4, чтобы получить индексы столбцов в df, которые мы хотим сохранить.

Образец данных 

#df <- read.table(text =
#    "                Fraction Treatment Time Replicate        A        B         C        D         E         F
#LL10.5AT           T        LL      10.5    A        11.11428 11.82154 10.445625 8.849699 10.373386  9.109676
#LL10.5BT           T        LL      10.5    B        12.17890 11.01224 11.720548 9.405390 10.206708 10.653205
#LL10.5CT           T        LL      10.5    C        10.80697 11.19782 11.175291 8.305949  9.696153  8.791403
#OL10.5AT           T        OL      10.5    A        10.46481 10.81123  9.975277 7.783538  9.784773  8.640531
#OL10.5BT           T        OL      10.5    B        10.75621 10.76371 10.625745 7.592059  9.820686  8.760861
#OL10.5CT           T        OL      10.5    C        12.00054 11.02080 11.615536 8.903105  9.963635 10.547791
#HL10.5AT           T        HL      10.5    A        10.87092 11.45102 10.780183 6.422136 10.424391  9.489396
#HL10.5BT           T        HL      10.5    B        12.12334 11.29960 11.541679 9.774041  9.563639 10.532936
#HL10.5CT           T        HL      10.5    C 10.21460 10.64746  9.886603 7.834040  9.828347  8.261546", header = T)
...