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
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)
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With