Skip to content

[BUG] Error with RexProcess when first residue is not 1 #4

Description

@nlgittens

This one had been puzzling me for a while: RexProcess function will not work unless the starting residue is 1, giving the following error message:

`data("BRD4_apo")
BRD4_apo$Start <- BRD4_apo$Start + 10
BRD4_apo$End <- BRD4_apo$End + 10
BRD4_apo <- DataFrame(BRD4_apo)
BRD4_apo <- cleanHDX(res = BRD4_apo, clean = TRUE)
BRD4_apo <- data.frame(BRD4_apo) %>% filter(End < 100)
BRD4_apo <- DataFrame(BRD4_apo)
numTimepoints <- length(unique(BRD4_apo$Exposure))
Timepoints <- unique(BRD4_apo$Exposure)
numPeptides <- length(unique(BRD4_apo$Sequence))

set.seed(1)
rex_test <- rex(HdxData = BRD4_apo,
numIter = 10,
R = max(BRD4_apo$End),
density = "laplace",
numtimepoints = numTimepoints,
timepoints = Timepoints,
seed = 1L,
tCoef = c(0, rep(1, numTimepoints - 1)),
phi = 1,
BPPARAM = SerialParam())

rex_test <- RexProcess(HdxData = DataFrame(BRD4_apo),
params = rex_test,
range = 1:10,
thin = 1,
whichChains = c(1,2))`
Error in dimnames(x) <- dn :
length of 'dimnames' [2] not equal to array extent

This didn't seem to be a problem when running on a local PC, but when I compiled everything to run on hpc, I started getting this error. So wondering if it might be a version issue in one of my packages? I can provide the sessionInfo if that's useful. Manually working through RexProcess didn't seem to provide any clues because then the function worked.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions