Nothing
print.thetaMSAR <-
function(x,digits=4,...) {
if (!is.null(attr(x,'n_par'))) {
n_par <- attr(x,'n_par')
} else {
n_par <- 0
}
thetadim <- dim.thetaMSAR(x)
order <- attributes(x)$order
label=attributes(x)$label
M=attributes(x)$NbRegimes
d=attributes(x)$NbComp
if(label=='NULL'){label='HH'}
ncov.trans = length(x$par.trans[1,])-1
ncov.emis = length(c(x$par.emis[[1]]))
if (thetadim[2]>=2) {
x <- as.thetaMSAR(x,label=label,ncov.emis=ncov.emis,ncov.trans=ncov.trans) # this adds dimnames if they don't exist
sel <- 1:n_par
} else {
sel <- TRUE
}
class(x) <- NULL
dimname <- dimnames(x)
if(order>0){cat("AUTOREGRESSIVE COEFFICIENTS\n")
if(d==1){
print(x$A,digits)
cat("\n")
} else {
for(i in 1:M){
cat(paste("Regime",i,"\n"))
for(j in 1:order){
cat(paste("A",j,sep=""))
print(x$A[[i]][[j]],digits)
cat("\n")
}
}
}
}
if (order==0) {
cat("MEAN\n")
print(x$A0,digits)
} else {
cat("Intercept\n")
print(x$A0,digits)
# cat("MEAN\n")
# A = NULL
# for(m in 1:M){
# for (o in 1:order) {
# A[m,] = A[[m]]+x$A[m,o]
# }
# }
# print(solve(diag(d)-A[[m]])%*%x$A0,digits)
}
cat("\nVARIANCE")
if(d==1){
print(x$sigma,digits)
cat("\n")
}
else{
for(i in 1:M){
cat(paste('\nRegime',i))
print(x$sigma[[i]],digits)
}
}
cat("\nINITIAL DISTRIBUTION OF THE MARKOV CHAIN")
print(x$prior,digits)
cat("\nTRANSITION PROBABILITY MATRIX\n")
print(x$transmat,digits)
if(label=="HN" || label=="NN"){
cat('\nEMISSION PARAMETERS\n')
for(i in 1:M){
cat("Regime",i,'\n')
print(x$par.emis[[i]],digits)
cat('\n')
}
}
if(label=='NH' || label=='NN'){
cat("\nTRANSITION PARAMETERS\n")
print(x$par.trans,digits)
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.