zhqingit / giremi

GIREMI is a method that can identify RNA editing sites using one RNA-seq data set without requiring genome sequence data.
42 stars 15 forks source link

GIREMI stopped working, Rout Error attached #49

Closed LastGenomeBender closed 6 years ago

LastGenomeBender commented 6 years ago

Dr. Qing,

First of all thank you for a great tool. While running GIREMI we encountered an error. Below is the terminal error:

[mpileup] 1 samples in 1 input files Calculating the MI values... meanMI:0.245166 sdMI:0.330423 Estimating the allelic ratios... R CMD BATCH --no-save --no-restore '--args finput="RNAE.lst" fout="RNAE.lst.res"' giremi.r RNAE.lst.Rout Analysis DONE!

The Rout file error is the following:

dat <- read.table(file=finput,header=F,row.names=1,comment.char="")

dat <- read.table(file=finput,header=T,comment.char="") Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : line 98 did not have 24 elements Calls: read.table -> scan Execution halted

Our input .txt file did not include multiple SNPs for a single site. We would be happy if you could help.

Regards, Farid Ahadli