tests/niche_test.R

In quicklizard99/Basingstoke: Food-web structural models

TestNicheModelLinks <- function() {
    AssertRaises(NicheModelLinks(S=0))
    AssertRaises(NicheModelLinks(S=1, n=0))
    AssertRaises(NicheModelLinks(S=1, n=1, C=0))
    AssertRaises(NicheModelLinks(S=1, n=1, C=1))

    # 100 sets of links with directed connectance of 0.15 taken from from a 
    # pool of 10 species.
    res <- NicheModelLinks(S=10, n=500, C=0.15)
    AssertEqual(500, length(res))
    AssertEqual(15, mean(sapply(res, nrow)), tolerance=1, scale=1)
    AssertEqual(FALSE, any(sapply(res, function(r) any(duplicated(r)))))

    # 100 sets of links with directed connectance of 0.5 taken from from a 
    # pool of 10 species.
    res <- NicheModelLinks(S=10, n=500, C=0.45)
    AssertEqual(500, length(res))
    AssertEqual(45, mean(sapply(res, nrow)), tolerance=6, scale=1)
    AssertEqual(FALSE, any(sapply(res, function(r) any(duplicated(r)))))
}

TestProbabilisticNicheModelLinks <- function() {
    # 100 sets of links with directed connectance of 0.15 taken from from a 
    # pool of 10 species.
    res <- NicheModelLinks(S=10, n=500, C=0.15, probabilistic=TRUE)
    AssertEqual(500, length(res))
    AssertEqual(15, mean(sapply(res, nrow)), tolerance=3, scale=1)
    AssertEqual(FALSE, any(sapply(res, function(r) any(duplicated(r)))))

    # 100 sets of links with directed connectance of 0.5 taken from from a 
    # pool of 10 species.
    res <- NicheModelLinks(S=10, n=500, C=0.45, probabilistic=TRUE)
    AssertEqual(500, length(res))
    AssertEqual(45, mean(sapply(res, nrow)), tolerance=7, scale=1)
    AssertEqual(FALSE, any(sapply(res, function(r) any(duplicated(r)))))
}