делает что-то вроде этой справки:
library(tidyverse)
x <- data.frame(stringsAsFactors=FALSE,
Levels = c("1_at", "10_at", "100_at", "1000_at", "10000_at", "100009613_at"),
Tarca_001_P1A01 = c(6.062215, 3.796484, 5.849338, 3.567779, 6.166815,
4.443027),
Tarca_003_P1A03 = c(6.125023, 3.805305, 6.191562, 3.452524, 5.678373,
4.773199),
Tarca_004_P1A04 = c(5.875502, 3.450245, 6.550525, 3.316134, 6.185059,
4.393488),
Tarca_005_P1A05 = c(6.126131, 3.628411, 6.421877, 3.432451, 5.633757,
4.623783)
)
y <- data.frame(stringsAsFactors=FALSE,
gene = c("Tarca_001_P1A01", "Tarca_013_P1B01", "Tarca_025_P1C01",
"Tarca_037_P1D01", "Tarca_049_P1E01", "Tarca_061_P1F01"),
SampleID = c("Tarca_001_P1A01", "Tarca_013_P1B01", "Tarca_025_P1C01",
"Tarca_037_P1D01", "Tarca_049_P1E01", "Tarca_061_P1F01"),
GA = c(11, 15.3, 21.7, 26.7, 31.3, 32.1),
Batch = c(1, 1, 1, 1, 1, 1),
Set = c("PRB_HTA", "PRB_HTA", "PRB_HTA", "PRB_HTA", "PRB_HTA", "PRB_HTA")
)
x %>% gather(SampleID, value, -Levels) %>%
left_join(., y, by = "SampleID") %>%
group_by(SampleID) %>%
filter(value == max(value)) %>%
spread(SampleID, value)