# Determine intersection of compounds in features and responses | |
commonMols <- intersect(responses$CID, | |
molFeats$CID) | |
# Subset features and responses accordingly | |
responses %<>% filter(CID %in% commonMols) | |
molFeats %<>% filter(CID %in% commonMols) | |
# Compute median pleasantness across the population | |
medianPlsnt <- responses %>% | |
group_by(CID) %>% | |
dplyr::summarise(pleasantness = median(`VALENCE/PLEASANTNESS`, na.rm = T)) | |
all(medianPlsnt$CID == molFeats$CID) # TRUE - rownames match | |
# Concatenate predictors (molFeats) and population pleasantness | |
X <- mutate(molFeats, Y = medianPlsnt$pleasantness) %>% | |
select(-CID) |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment