Question: Error in running "diffExpressedVariants" function in kissDE R package
0
gravatar for c.evang
14 months ago by
c.evang10
Italy
c.evang10 wrote:

Dear all, I want to run kissDE R package on my set of SNPs, that I found using KisSplice software. I'm using this code:

snpB<-kissplice2counts("results_all_reads_mDr_all_reads_WW_k25_coherents_type_0a.fa", pairedEnd=FALSE)
conditions<-c("mDr","WW")
resB<-diffExpressedVariants(snpB, conditions, pvalue=1)

But it gives me the following error:

"Error in rowMeans(Y2[, idxA]): 'x' must be an array of at least two dimensions\n An error occured, unable to fit models on data".

Any suggestions will be appreciated.

Thanks

kissde snp kissplice • 650 views
ADD COMMENTlink modified 11 months ago by vincent.lacroix50 • written 14 months ago by c.evang10
1
gravatar for vincent.lacroix
14 months ago by
vincent.lacroix50 wrote:

Dear c.evang,

KissDE requires replicates. This is most likely the cause of your problem. I undertand you have two biological conditions (mDr, and WW). If you have replicates for each condition, then you need to run KisSplice on all replicates of all conditions. Those are 4 reads files (or possibly 8 reads file if you have paired-end reads). The command line should look like this

kissplice -r mDr1.fa -r mDr2.fa -r WW1.fa -r WW2.fa

this will produce an output file such as: results_mDr1_mDr2_WW1_WW2_type0.fa you can then run KissDE:

snpB<-kissplice2counts("results_mDr1_mDr2_WW1_WW2_type0.fa", pairedEnd=FALSE)
conditions<-c("mDr","mDr","WW","WW")
resB<-diffExpressedVariants(snpB, conditions, pvalue=1)

If you do not have replicates, then you cannot use KissDE.

ADD COMMENTlink written 14 months ago by vincent.lacroix50

Thank you very much vincent.lacroix.

So I will repeat kisSplice using replicates and then kissDE.

Thanks again Chiara

ADD REPLYlink written 14 months ago by c.evang10

Dear Vincent.lacroix,

I tried running kissplice using replicates, but another error occured while running kissDE:

[1] "Error in paste(pathToFigs, \"/models.png\", sep = \"\"): object 'pathToFigs' not found\n An error occured, unable to fit models on data."

Any suggestions?

Thanks

ADD REPLYlink modified 14 months ago • written 14 months ago by c.evang10

Hi @vincent.lacroix,

I also encountered an almost similar error:

"Error in .validate_assayDataElementReplace(obj, value): object and replacement value have different dimensions\n An error occured, unable to fit models on data."

This actually occured when I ran KissDE using my data; however, this also resulted after running the example script documented in ?diffExpressedVariants. I also successfully ran qualityControl() and got no error from it.

And I have no idea what and where to look for. Sorry. and thanks for helping out!

ADD REPLYlink modified 12 months ago • written 12 months ago by jfo0
0
gravatar for vincent.lacroix
11 months ago by
vincent.lacroix50 wrote:

Dear all,

This bug is now solved. It was related to KissDE's dependency to Biobase and DSS. These packages had been upgraded, and there were compatibiliy issues to fix. The estimation of the dispersion is now also better. Hence, the final output of KissDE 1.4.0 is expected to change compared to KissDE 1.3.0. Some new candidates can be detected. Some borderline candidates are now considered as non significant. The new version can be downloaded here (http://kissplice.prabi.fr/tools/kissDE/). We encourage users to move to version 1.4.0 of KissDE, and 2.14 of DSS. Please tell us if you experience problems.

Vincent

ADD COMMENTlink written 11 months ago by vincent.lacroix50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 962 users visited in the last hour