Last chance! 50% off unlimited learning
Sale ends in
alignSequences
, getMutations
and mapMutations
in order to execute all these methods at once.
setup(object, repos = NULL, clustalo_filename=NULL
, mail=NULL , perlCommand="perl", use_hmm=FALSE, datum=FALSE)
alignSequences
, getMutations
and mapMutations
methods.
alignSequences
getMutations
mapMutations
#Create an object of class LowMACA for RAS domain family
lm <- newLowMACA(pfam="PF00071" , genes=c("KRAS" , "NRAS" , "HRAS"))
#Select a few tumor types
lmParams(lm)$tumor_type <- c("skcm" , "brca" , "coadread")
#Align sequences, get mutation data and map them on consensus
lm <- setup(lm)
#Same as above, but using web service
lm <- setup(lm , mail="lowmaca@gmail.com")
#Use HMM to align
lm <- setup(lm , use_hmm=TRUE)
#Use "datum"
lm <- setup(lm , datum=TRUE)
Run the code above in your browser using DataLab