按不同属性创建presence-absence表

Question

我有一个大的不同地点的物种检测数据集表（12个村庄，每个4个街区 - 每个街区9个点）。我需要为每个块（48块= 12 * 4）创建所有物种组合的存在/不存在表，但也需要为每个物种组/公会（总共20组）创建存在/不存在表。这将是总共48 * 20 = 960表。我可以通过excel Pivot表来实现这一目标，但是不是有更快的方法来实现这个吗？

我的表格列如下所示：

mydata <- structure(list(Sample = c("1", "2", "3", "7", "8", "9", "13", "14", "15"), 
                     Block = c("A", "A", "A", "B", "B","B", "A","A", "A"), 
                     Location = c("Be", "Be", "Be", "Be", "Be", "Be", "Be", "Be", "Be"),
                                  Species = c("a","d","r","a","c","r","a","c", "c"),
                     Guild1 = c("O", "H", "H","O","I","H","O","I","I")), 
                .Names = c("Sample", "Block", "Location", "Species", "Guild1"), 
                class = "data.frame", row.names = c(NA, -3L))

预期的输出数据应该是存在/不存在，行中的样本和列中的种类。但是对于按位置和公会的每个街区，应该建立一个单独的位置表。如此准备，我可以将数据读入EstimateS程序来估算物种丰富度。