r - Trinity failed with differential expression analysis test -


i tried running de test comes trinity after installing that's mentioned on website, ran error.

warning: edger comparison failed...  cmd: r --vanilla -q < trinity_trans.counts.matrix.diauxic_shift_vs_plateau.edger.rscript > library(edger) > > data = read.table("/usr/local/trinityrnaseq-r2013.02.25/sample_data/test_full_edger_pipeline/trinity_trans.counts.matrix", header=t, row.names=1, com='') > col_ordering = c(4,1) > rnaseqmatrix = data[,col_ordering] > rnaseqmatrix = round(rnaseqmatrix) > rnaseqmatrix = rnaseqmatrix[rowsums(rnaseqmatrix)>=10,] > conditions = factor(c(rep("plateau", 1), rep("diauxic_shift", 1))) > > exp_study = dgelist(counts=rnaseqmatrix, group=conditions) calculating library sizes column totals. > exp_study = calcnormfactors(exp_study) > et = exacttest(exp_study, dispersion=0.1) comparison of groups:  plateau - diauxic_shift > ttags = toptags(et,n=null) > write.table(ttags[ttags$table$pvalue <= 0.05,], file='trinity_trans.counts.matrix.diauxic_shift_vs_plateau.edger.de_results', sep='   ', quote=f, row.names=t) > source("/usr/local/trinityrnaseq-r2013.02.25/analysis/differentialexpression/r/rnaseq_plot_funcs.r") > pdf("trinity_trans.counts.matrix.diauxic_shift_vs_plateau.edger.de_results.ma_n_volcano.pdf") > result_table = ttags$table > plot_ma_and_volcano(result_table$logcpm, result_table$logfc, result_table$fdr) error in xy.coords(x, y, xlabel, ylabel, log) :   'x' , 'y' lengths differ calls: plot_ma_and_volcano -> plot_ma -> plot -> plot.default -> xy.coords execution halted error, cmd: r --vanilla -q < trinity_trans.counts.matrix.diauxic_shift_vs_plateau.edger.rscript died ret (256)  @ /usr/local/trinityrnaseq-r2013.02.25/util/..//analysis/differentialexpression/run_de_analysis.pl line 416.


Comments

Post a Comment

Popular posts from this blog

linux - Does gcc have any options to add version info in ELF binary file? -

i mean whether gcc can insert source code version infor elf binary section or similar. not want change source file, add info gcc option in makefile. you can emit version info text file, turn text file object file statically link executable. the first step simple have write code: script or write version info in format plain text file. write makefile rule produce version.o version.txt, using objcopy . you'll have object file 2 useful symbols defined in it: beginning , end of textual version info. add generated object executable, , you'll able access version 2 ways: running strings on binary, or writing code in application print version string (you'll need declare start , end symbols variables in header file).
Read more

android - send complex objects as post php java -

is possible send contact numbers via post php class inside eclipse ? or post limited primitive types , have find work arounds send them string ? i propose: send them json in post request, befor start, have @ gson project , keep in mind set content application -> json. luck
Read more

charts - What graph/dashboard product is facebook using in Dashboard: PUE & WUE -

i love facebook dashboard: pue & wue. beautifull can point out if there open source or commercial product offer similar style? here link https://www.facebook.com/forestcitydatacenter?sk=app_288655784601722&app_data so read through 29,888 lines of javascript ya.. biggest js file i've ever seen. looks make graph using variety of tools, including tween.js . it'd fascinating chat whoever made it. code wheel graph object starts @ line 27,266 text var graphcomponent = class.extend({ . can find source @ https://www.fbpuewue.com/assets/application-3a3af20f7fa3591229c79442ae3d1c26.js . sorry if dissapoints it's not awesome plugin can use :/ share pain. let me know if have questions :)
Read more