Entering edit mode
2.9 years ago
Nai
▴
50
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]