R/summary.traj.r

In adehabitat: Analysis of Habitat Selection by Animals

"summary.traj" <- function(object, id=levels(object$id), date=NULL, ...)
{
    ## Verifications
    x<-object
    if (!inherits(x, "traj"))
        stop("x should be an object of class traj")

    ## Selection of dates
    if (!is.null(date))
        x<-x[(x$date>=date[1])&(x$date<date[2]),]

    ## Selection of animals
    li<-split(x, factor(x$id))
    x<-do.call("rbind", li[id])
    x$id<-factor(x$id)
    x$burst<-factor(x$burst)

    ## Output
    ll<-list()
    for (i in id) {
        if (!is.na(match(i, levels(x$id)))) {
            cat("Animal ", i, ": ",
                nlevels(factor(li[[i]]$burst)),
                "circuits. \nNumber of relocations per circuit:")
            o<-li[[i]]
            print(table(factor(o$burst)))
            ll[[i]]<-table(factor(o$burst))
            cat("\n")
        }
    }
    invisible(ll)
}