I am trying to run CDHIT , on DNAStringSet object of DNA sequences to get non-redundent set of sequences
library(Biostrings) library("BioSeqClass") ## flank is DNAStrinSet object of equal length sequences seq = as.character(flank) ## Homolog reduction of whole-length sequence by cd-hit need cd-hit program;
I have downloaded latest CD-HIT programme and named its fiolder as "cd-hit" from the given link http://www.bioinformatics.org/download/cd-hit/ and given path of the given folder to cdhit.path in above command.
> reducSeq50 = hr(seq, method="cdhit", identity=0.5, cdhit.path="cd-hit")
Error in cdhitHR(seq, identity, cdhit.path) :
cdhit.path cd-hit is not corrected
Please let me know where I am wrong.