### Example 1: simple atrificial data
ss <- data.frame(
cbind(
c("sample1","sample1","sample2","sample2","sample3","sample4"),
c("species1","species1","species1","species2","species3","zero.count"),
c("male","female","male","female","male","male")
),
c(1, 2, 10, 3, 4, 1)
)
colnames(ss) <- c("sample.id", "species.id", "gender", "catch")
xc1 <- xcount(sscount(ss, zc="zero"))
xc1
xc2 <- xcount(sscount(ss, zc="zero.count"))
xc2
xc3 <- xcount(sscount(ss, zc="zero.count"), 0)
xc3
xc4 <- xcount(sscount(ss, zc="zero.count"), 1)
xc4
### Example 2: field data of the dolina
data(dol.count, dol.sample, landsnail)
dssc <- sscount(dol.count, zc="zero.count", fill=TRUE)
dssc
(xc <- xcount(dssc))
plot(xc, type="hist")
plot(xc, type="rank", logscale = TRUE)
plot(xc, type="biplot")
xcount(dssc, "broken")
xcount(dssc, c("adult", "juvenile"))
xcount(dssc, c(1,4), "mixed")
xcount(dssc, c(1:3))
Run the code above in your browser using DataLab