S1

R Markdown

this is sample r markdown documentation for DNA sequence will be published in online platforms

generate_dna_sequence <- function(length) 
  {
  nucleotides <- c("A", "T", "C", "G")
  sequence <- sample(nucleotides, length, replace = TRUE)
  return(paste(sequence, collapse = ""))
}

# Generate a random DNA sequence
dna_sequence <- generate_dna_sequence(100)

summary(dna_sequence)

##    Length     Class      Mode 
##         1 character character

# Count nucleotide frequencies
nucleotide_counts <- table(strsplit(dna_sequence, "")[[1]])
nucleotide_counts

## 
##  A  C  G  T 
## 20 30 28 22

# Create a bar plot to visualize nucleotide frequencies
barplot(nucleotide_counts, 
        main = "DNA Nucleotide Frequencies", 
        xlab = "Nucleotide", 
        ylab = "Frequency", 
        col = "orange")

plot(nucleotide_counts)

plot(nucleotide_counts,
     main = "Scatter Plot DNA Nucleotide Frequencies",
     xlab = "Nucleotide",
     ylab = "Frequency"
     )