How to convert data containing species name and occurrence site into species by site matrix in r?

Question

For example, I have a data frame looks like this, the \'species\' column (total 8 species) represent species name, the \'site\' column (total 4 sites) represent which site e