#P1(B.napus)
B.napu_sRNA <- srnapredata(sRNAseq = P1_sRNA_seq,Group = "B.napus(AACC)")
#P2(B.rapa)
B.rapa_sRNA <- srnapredata(sRNAseq = P2_sRNA_seq,Group = "B.rapa(AA)")
#F1(B.napus X B.rapa)
B.nr_sRNA <- srnapredata(sRNAseq = F1_sRNA_seq,Group = "B.napus x B.rapa(AAAACC)")
#intergrate these data for length distribution plot
sRNA_data <- rbind(B.napu_sRNA,B.rapa_sRNA,B.nr_sRNA)
#plot
lenplot(sRNAdata = sRNA_data,type = "line")
lenplot(sRNAdata = sRNA_data,type = "bar")
Run the code above in your browser using DataLab