if(interactive()){
file = file.path(path.package("BioSeqClass"), "example", "acetylation_K.pos40.pep")
seq = as.matrix(read.csv(file,header=F,sep="\t",row.names=1))[,1]
CKSAAP0 = featureCKSAAP(seq,0,elements("aminoacid"))
CKSAAP2 = featureCKSAAP(seq,2,elements("aminoacid"))
}
Run the code above in your browser using DataLab