R/print.segmented.R

In segmented: Regression Models with Break-Points / Change-Points Estimation (with Possibly Random Effects)

`print.segmented` <- function(x,digits = max(3, getOption("digits") - 3),...){
#revisione 15/05/03; 24/02/04/.. 5/06/19!!! :))
        if(is.null(x$psi)) x<-x[[length(x)]]
        if(is.null(x$psi)) stop("the object does not include the 'psi' component")
        #if(!"segmented"%in%class(x)) stop("a `segmented' object is requested")
        cat( "Call: " )
        print( x$call )
        cat("\nCoefficients of the linear terms:\n")
        #print(x$coef[(1:(length(x$coef)-length(x$psi[,2])))])
        .coef <- coef(x)
        if(is.null(.coef)) .coef<- x$coef
        if(is.null(.coef)) stop(" No coefficients in the fit object?")
        iV<- -match(x$nameUV$V,names(.coef))#iV<- -grep("psi.",names(coef(x)))#indices all but V
        

        #print(x$coef[iV])
        print.default(format(.coef[iV], digits = digits), print.gap = 2, quote = FALSE)
        cat("\n")
        cat("Estimated Break-Point(s):\n")
        a<-as.vector(x$psi[,"Est."])
        names(a)<-rownames(x$psi)
        print.default(a, digits, print.gap=2)
        #print(x$psi[,"Est."], digits)
        #ng<-max(x$id.psi.group)
        #for(i in 1:ng){
        #  nomi<-x$nameUV$V[x$id.psi.group==i]
        #  cat("\t",nomi,":", format(signif(x$psi[nomi,"Est."],digits)),"\n")
        #  }
        if("glm"%in%class(x)){    cat("\nDegrees of Freedom:", x$df.null, "Total (i.e. Null); ",
                x$df.residual, "Residual\n")
            cat("Null Deviance:    ", format(signif(x$null.deviance,
                digits)), "\nResidual Deviance:", format(signif(x$deviance,
                digits)), "     AIC:", format(signif(x$aic, digits)), "\n")
            }
        if("Arima"%in%class(x)){
            cm <- x$call$method
            if (is.null(cm) || cm != "CSS")
                cat("\nsigma^2 estimated as ", format(x$sigma2, digits = digits),
                    ":  log likelihood = ", format(round(x$loglik, 2)),
                    ",  aic = ", format(round(x$aic, 2)), "\n", sep = "")
            else cat("\nsigma^2 estimated as ", format(x$sigma2, digits = digits),
                ":  part log likelihood = ", format(round(x$loglik, 2)),
                "\n", sep = "")
            }
            invisible(x)
}