reported bug by Leonardo Saravia ‏@larysar

test_siber.hull.metrics_AJ.r

require(plyr)
require(siar)

t1 <- data.frame(SizeClass=c("Meiofauna","Smaller Macrofauna","Bigger Macrofauna"),
           num=rpois(3,10),
           md13C=rnorm(3,-10,10),
           sd13C=rlnorm(3,1,1),
           md15N=rnorm(3,5,5) ,
           sd15N=rlnorm(3,1,1))
t1

hh <- function(x) {
  with(x,data.frame(SizeClass=SizeClass,
             #             d13C=rnorm(num,md13C,sd13C),d15N=rlnorm(num,log(md15N),log(sd15N)))
             d13C=rnorm(num,md13C,sd13C),d15N=rnorm(num,md15N,sd15N)))
}
t2 <- ddply(t1, .(SizeClass), hh)

t3 <- with(t2,siber.hull.metrics(d13C,d15N,SizeClass,R=10^4))

# ---------------------------------------------------------
# Andrew Jackson - 16/6/2015

# try long-handed
t4 <- siber.hull.metrics(t2$d13C, t2$d15N, t2$SizeClass, R = 10^4)


summary(t3)
summary(t4)

# ---------------------------------------------------------


#
# Error: not a matrix
#