## Not run:
# library(embryogrowth)
# m <- c(10, 14, 7, 4, 3, 0, 0)
# f <- c(0, 1, 2, 4, 15, 10, 13)
# t <- c(25, 26, 27, 28, 29, 30, 31)
# result <- tsd(males=m, females=f, temperatures=t)
# logLik(result)
# AIC(result)
# ## End(Not run)
Run the code above in your browser using DataLab