morrislab / tracksigfreq Goto Github PK

Dear developers of TrackSigFreq,

I am encountering an issue when running TrackSigFreq, I sometimes have the error

Error in do.call(cbind, chunkFits) : second argument must be a list

when running the command

results <- TrackSig:::getChangepointsPELT(trinut_count_data, MU, vcaf, scoreMethod, binSize, desiredMinSegLen)

with trinut_count_data a matrix of size (96, 10) with 10 timepoints, MU the signature matrix of size (96, nSig), vcaf a description of mutations of size (1000, 3), with columns (phi, qi, bin), scoreMethod = "SigFreq", binSize = 100, desiredMinSegLen = NULL

In particular, this command runs properly with a matrix MU with nSig=47, but fails with a subset matrix with nSig=7, included in the 47.
I am happy to share the files for trinut_count_data, MU (and the extract), and vcaf so you can reproduce the error, let me know how I can share that with you.

Thank you in advance for your help
Best wishes
Judith