Question about MatrixEQTL
0
0
Entering edit mode
7.1 years ago
mms140130 ▴ 60

my Question is :

i have the gene expression data =genexpt.txt
column1 is gene name , the rest of the columns are expression values

I followed the steps in the MatrixEQTL website
the following is specifying the data in my working directory

expression_file_name = paste(currwd, "/genexpr.txt",sep="")

then I loaded the expression_file_name as follows:

> gene = SlicedData$new() 
  gene$fileDelimiter = ""      # the TAB
> character gene$fileOmitCharacters = "NA" # denote missing values;
> gene$fileSkipRows = 1          # one row of column labels
> gene$fileSkipColumns = 2       # one column of row labels
> gene$fileSliceSize = 5000      # read file in slices of 2,000 rows
> gene$LoadFile(expression_file_name)

    it gave me the following error
 Error: nzchar(output_file_name) is not TRUE

Can anyone help with that please
gene R snp • 1.9k views
ADD COMMENT
0
Entering edit mode

Your delimiter setting seems to be weird. Would you try changing:

gene$fileDelimiter = ""      # the TAB

into

gene$fileDelimiter = "\t"      # the TAB

and try again? If you used another delimiter, please change "\t" accordingly.

Also,

 gene$fileSkipColumns = 2       # one column of row labels

means that you have TWO columns of row labels. Trying changing it to 1 also?

ADD REPLY

Login before adding your answer.

Traffic: 1595 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