[R] Biostrings

[Herve Pages]

qing wrote:
> 
> 
> Dear All，
> 
> I am a beginner in learning the package of Biostrings currently and for 
> practice doing an example.
> 
>  library(Biostrings);
>  dnaAlph<-new("BioPatternAlphabet",DNAAlphabet(),c(N="AGCT",
> + B="CGT",D="AGT",H="ACT",K="GT",M="AC",R="AG",S="CG",
> + V="ACG",W="AT",Y="CT"));
> Error in getClass(Class, where = topenv(parent.frame())) : 
>         "BioPatternAlphabet" is not a defined class


Hi Qing,


"BioPatternAlphabet" was a class defined in Biostrings 1 (Biostrings
version 1.y.z). In Biostrings 2, the class system has changed and you don't
need to create an instance of the DNA alphabet anymore:

  - To create a DNAString object, just do:

      > mydna <- DNAString("AGG-HCNTT")
      > mydna
        9-letter "DNAString" object
      Value: AGG-HCNTT

  - To get the DNA alphabet, call alphabet() on a DNAString object:

      > alphabet(mydna)
       [1] "A" "C" "G" "T" "M" "R" "S" "V" "W" "Y" "H" "K" "D" "B" "N" "-"

    or, if you don't have a DNAString instance yet:

      > DNA_ALPHABET
       [1] "A" "C" "G" "T" "M" "R" "S" "V" "W" "Y" "H" "K" "D" "B" "N" "-"

    Note that the "DNA alphabet" is the "IUPAC Extended Genetic Alphabet".

  - To get the mapping between the DNA alphabet and the set of ambiguities
    associated to each "extended" letter:

      > IUPAC_CODE_MAP
           A      C      G      T      M      R      W      S      Y      K      V
         "A"    "C"    "G"    "T"   "AC"   "AG"   "AT"   "CG"   "CT"   "GT"  "ACG"
           H      D      B      N
       "ACT"  "AGT"  "CGT" "ACGT"

  - For more details, see:

      > ?DNAString

Cheers,

H.

PS: As Ducan said, the Bioconductor mailing list is a more appropriate place
to ask help about a Bioconductor package.



> 
> I am running R 2.4.1,  Wimdows XP, Biostrings_2.2.1
> Any help or suggestions that you can provide will be greatly appreciated.
> 
> Qing
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.