Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Error: could not find function "read.DNAStringSet" [closed]

Tags:

r

This is very awkward. I use BSgenome to write a DNAStringSet fasta file,

writeXStringSet(DNAcutout1.seq, "DNAcutout1.fa", append=FALSE,
compress=FALSE,compression_level=NA, format="fasta")

but when I run another function, it is said

library(NuPoP)
predNuPoP("DNAcutout1.fa", species=5, model=4)
[1] "The input file is not in FASTA format, or contains characters other than A/a,
C/c,G/g, T/t, N/n, stop!"

so I wanna read it back to see what goes wrong, but I can't use read.DNAString

read.DNAStringSet(file="DNAcutout1.fa", format="fasta")
Error: could not find function "read.DNAStringSet"

Please HELP! Thank you

> sessionInfo()
R version 3.1.0 (2014-04-10)
Platform: x86_64-apple-darwin13.1.0 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets 
[7] methods   base     

other attached packages:
[1] NuPoP_1.14.1                          
[2] BSgenome.Dmelanogaster.UCSC.dm3_1.3.99
[3] BSgenome_1.32.0                       
[4] Biostrings_2.32.1                     
[5] XVector_0.4.0                         
[6] GenomicRanges_1.16.3                  
[7] GenomeInfoDb_1.0.2                    
[8] IRanges_1.22.10                       
[9] BiocGenerics_0.10.0                   
[10] BiocInstaller_1.14.2                  
[11] RSQLite_0.11.4                        
[12] DBI_0.2-7                             

loaded via a namespace (and not attached):
[1] bitops_1.0-6     Rsamtools_1.16.1 stats4_3.1.0    
[4] tools_3.1.0      zlibbioc_1.10.0 
like image 411
lxcfuji Avatar asked Mar 19 '23 21:03

lxcfuji


1 Answers

As I mentioned in the comments. read.DNAStringSet has been renamed to readDNAStringSet. Use the latter function name instead

(posting here so the question doesn't appear as unanswered)

like image 137
MrFlick Avatar answered Apr 02 '23 06:04

MrFlick