This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MATRIXEQTL: Error during load SNPMATRIX.txt in MATRIXEQTL
library(MatrixEQTL)

useModel = modelLINEAR

SNP_file_name = paste("/Users/user1/Desktop/Plant_1/eqtl/209_40.txt", sep = "")

expression_file_name = paste( "/Users/user1/Desktop/Plant_1/eqtl/Gene_expression.txt.", sep = "")

covariates_file_name = character()

output_file_name = tempfile()

pvOutputThreshold = 1e-2

errorCovariance = numeric()

snps = SlicedData$new()

snps$fileDelimiter = "\t"

snps$fileOmitCharacters = "NA"

snps$fileSkipRows = 1

snps$fileSkipColumns = 1

snps$fileSliceSize = 2000

snps$LoadFile( SNP_file_name )

Error: Error in dim(data) = c(length(rowvals[[1]]), length(rowvals)) : dims [product 48375] do not match the length of object [48362]
data matrixeqtl dims snpmatrix.txt
• 618 views • link
updated by GenoMax GenoMax • written by Nai • Nai

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Help needed
    written by Nai • Nai
    var snp_id=0; var samples = header.getSampleNamesInOrder(); out.print("id"); for(var i=0;i< samples.size();++i) { out.print("\t"+samples.get(i)); } out.println(); while(iter.hasNext()) { var ctx= iter.next(); ++snp_id; out.print(new java.lang.Integer(snp_id)); for(var i=0;i< samples.size();++i) …
  • EQTL MATRIX _EQTL
    written by Nai • Nai
    ```r setwd("C:/Users/genetics/Desktop/Plant_39/eqtl/") getwd() [1] "C:/Users/genetics/Desktop/Plant_39/eqtl" library(MatrixEQTL) base.dir = find.package("MatrixEQTL") useModel = modelLINEAR SNP_file_name = paste(base.dir, "SNP.txt", sep="") expression_file_name = paste(base.dir, "GE.txt", sep="") covariates_file_name = character() …
  • Question about MatrixEQTL
    written by mms140130 mms140130
    my Question is : i have the gene expression data =genexpt.txt column1 is gene name , the rest of the columns are expression values I …