У меня есть фрейм данных со значениями Fst для разных окон SNP.
Я хотел бы суммировать среднее значение и медиану значений Fst для каждого окна гена:
gene chr start Nsites Fst01 Fst02 Fst12
1 AC002310.7 16 30591000 2 1 8. 5
2 ZNF785 16 30591000 2 1 30 1
3 AC002310.7 16 30593000 2 2 4. 10
4 ZNF785 16 30593000 2 3 20. 2
5 AC002310.7 16 30595000 2 3 2. 20
6 ZNF785 16 30595000 2 6. 10. 3
Вывод:
gene n Fst01 Fst02 Fst12
1 AC002310.7 3 2. 4. 10
2 ZNF785 3 3. 20 2
Я написал что-то вроде этого:
df %>% summarise_at(group_by(df, gene), vars(Fst01, Fst02,Fst12), funs(n(),mean))
, но это дает мне следующую ошибку:
Error: `.vars` must be a character/numeric vector or a `vars()` object, not list
Call `rlang::last_error()` to see a backtrace
Любая помощь, высоко ценится, ура