# What steps should I follow for DNA analysis, for a classification problem?

I have a bed file which contains DNA sequences information as follow:

**

track name="194" description="194 methylation (sites)" color=0,60,120 useScore=1
chr1    15864   15866   FALSE   894 +
chr1    534241  534243  FALSE   921 -
chr1    710096  710098  FALSE   729 +
chr1    714176  714178  FALSE   12  -
chr1    720864  720866  FALSE   988 -


**

I loaded the bed file in R and named the matrix DataSet. I used the follow code to get the sequences:

mydataSet_Test1<-dataSet[,1:3]

library(BSgenome.Hsapiens.UCSC.hg19)
genome <- BSgenome.Hsapiens.UCSC.hg19
chr<-as.matrix(as.character(mydataSet_Test1[,1]))

#50
start<-as.matrix(as.integer(as.character(mydataSet_Test1[,2]))-50)
end<-as.matrix(as.integer(as.character(mydataSet_Test1[,3]))+50)
Seqs50_Test1<-getSeq(genome,chr,start=start,end=end)


I'm having DNA sequences and for each sequence, I have its Methylation Number (i.e. 800). I have split the data to hypomethylated (methylation < 200, class 0) and hypermethylated (methylation > 800, class 1) and I want to build a model that can predict if a sequence (or a group of sequences) is hypo or hyper - methylated.

I now have a matrix of 358.367 data and 2 columns. The first column is the DNA sequences and the second column is the class for each sequence (0 or 1).

I want to build a classification model in R, using SVM algorithm and 83 features (dinucleotides, trinucleotides, etc.).

Before calculating the features and building the model, are there any steps I should follow for my DNA sequence data (each row of the matrix is a sequence from the human genome), like filtering CpG islands or sequences, homolog reduction or overlapping sequences etc?

As you understand I have a lack in DNA analysis.

• Hi Giannis, could you please describe the classification problem you are trying to solve? Jun 3 '20 at 14:11
• I'm having DNA sequences and for each sequence, I have its Methylation Number (i.e. 800). I have split the data to hypomethylated (methylation < 200, class 0) and hypermethylated (methylation > 800, class 1) and I want to build a model that can predict if a sequence (or a group of sequences) is hypo or hyper - methylated Jun 4 '20 at 5:43
• @GiannisLazaridis Please edit the question to include the last comment, and include the code or guides you have tried
– llrs
Jun 4 '20 at 7:50
• I edited the question Jun 4 '20 at 11:23