Hi, I need to do a genome-wide in silico restriction digest for the yeast genome using R. I need to do a dual restriction digest (digest using two restriction enzymes) and then plot a histogram of size of fragments vs their frequency. I am using this R-Bloggers script as a base.

Would someone please help me modify the script so that it correctly performs a dual restriction digest (ex. EcoRI and KpnI) instead of a single restriction digest (as shown in original script)?

Thank you!