Entering edit mode
2.9 years ago
Nai
▴
50
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()
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 in file(description = filename, open = "rt", blocking = FALSE, raw = FALSE) :
cannot open the connection
In addition: Warning message:
In file(description = filename, open = "rt", blocking = FALSE, raw = FALSE) :
cannot open file 'C:/Program Files/R/R-4.0.3/library/MatrixEQTLSNP.txt': No such file or directory