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