with(byTrophBySite[byTrophBySite$trophic == 'H', ],
plot(factor(as.character(Site), levels = c('VO', 'LA', 'KH', 'MO', 'KA')),
plot(byTrophBySite$S0, byTrophBySite$sad.z, col = factor(byTrophBySite$trophic))
Add the following code to your website.
Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.