Ballgown object construction error
3
0
Entering edit mode
spaque57 • 0
@spaque57-13166
Last seen 6.8 years ago

Hello,

 

I have an error that I can't resolve when I'm creating a ballgown object:

> pheno_data = read.csv("HighLow_phenodata_2.csv", sep=";", header=T)
> bg = ballgown(dataDir = "ballgown", samplePattern = "SNBN", pData=pheno_data)

Thu Jun  1 17:38:29 2017
Thu Jun  1 17:38:29 2017: Reading linking tables
Thu Jun  1 17:38:29 2017: Reading intron data files
Thu Jun  1 17:38:30 2017: Merging intron data
Thu Jun  1 17:38:30 2017: Reading exon data files
Thu Jun  1 17:38:31 2017: Merging exon data
Thu Jun  1 17:38:31 2017: Reading transcript data files
Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  : 
  scan() attendait 'a real' et a reçu '.'

 

> traceback()
4: scan(file = file, what = what, sep = sep, quote = quote, dec = dec, 
       nmax = nrows, skip = 0, na.strings = na.strings, quiet = TRUE, 
       fill = fill, strip.white = strip.white, blank.lines.skip = blank.lines.skip, 
       multi.line = FALSE, comment.char = comment.char, allowEscapes = allowEscapes, 
       flush = flush, encoding = encoding, skipNul = skipNul)
3: read.table(file, header = TRUE, sep = "\t", colClasses = cc, 
       quote = "")
2: .readTranscript(f, meas)
1: ballgown(dataDir = "ballgown", samplePattern = "SNBN", pData = pheno_data)

 

 

If anyone has an idea to solve it or know a package to do gene expression analysis like ballgown :/

Thank you !

 

ballgown scan error package • 1.8k views
ADD COMMENT
1
Entering edit mode
spaque57 • 0
@spaque57-13166
Last seen 6.8 years ago

I fund what was the problem. For those who would have this problem:
- I used a gtf file that I created alone
- it was the source of the problem because of the last column

For a GTF file, the values in the last column must be spaced because the columns are separated by tabs.
When generating a GTF file, I had everything tabbed, that's why the scan problem!

ADD COMMENT
0
Entering edit mode
Haiying.Kong ▴ 110
@haiyingkong-9254
Last seen 5.0 years ago
Germany

If it is csv file, why is sep semicolon?

"HighLow_phenodata_2.csv", sep=";"

ADD COMMENT
0
Entering edit mode
spaque57 • 0
@spaque57-13166
Last seen 6.8 years ago

My excel software put ";" but i try with "," and i also get the error :/

ADD COMMENT
0
Entering edit mode

I never used this function, but i would check transcript data files in the ballgown folder.

probably you can find out how ballgown reads in transcript data, and see if you can read the data in same way.

ADD REPLY
0
Entering edit mode

Thanks for this advice but a can't find any tips on the ballgown script which could occur this error :/

ADD REPLY
0
Entering edit mode

i am not 100% sure.

but, most likely, scan error means read error.

something in your transcript data is not compatible with the function with which ballgown reads in the data.

ADD REPLY
0
Entering edit mode

Thank you for your help !

I try with different data and phenotype file, and the code works, so i had a problem with my data.

Now i know what to do ;)

ADD REPLY

Login before adding your answer.

Traffic: 747 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6