I have run coseq on my gene expression matrix from RNAseq (15000 genes* 750 samples). I run the following:
runArcsin <- coseq(counts, K=2:3, model="kmeans", transformation="clr", meanFilterCutoff=200, iter=10,GaussianModel = "Gaussian_pk_Lk_I")
I am running this with Gaussian_pk_Lk_I because using the default gausian model did not work. Also i have specified K in restricted way just for testing.
The error I get is:
Error in h(simpleError(msg, call)) : error in evaluating the argument 'object' in selecting a method for function 'DDSEextract': object 'cap' not found
could you help?