Nothing
setMethod("summary",
signature(object = "bild"),
function(object,cov=FALSE,cor=FALSE)
{
cat("\nCall:\n")
print(object@call)
cat("\nNumber of profiles in the dataset: ", object@ni.cases)
cat("\nNumber of profiles used in the fit: ", object@n.cases)
cat("\nDependence structure:", object@call$dependence)
cat("\nLog likelihood: ", round(object@log.likelihood, 4))
cat("\nAIC: ", round(object@aic, 4),"\n")
coef <- object@coefficients
nas <- is.na(coef[, 1])
cnames <- names(coef[, 1][!nas])
coef <- matrix(rep(coef[, 1][!nas], 4), ncol = 4)
coef.aux<-matrix(rep(coef[, 1][!nas], 4), ncol = 4)
# coef[, 1] <- 1:dim(coef)[[1]]
coef[, 2] <- object@se[, 1][!nas]
coef[, 3] <- round(coef[, 1]/coef[, 2], 3)
coef.aux[,1]<-pnorm(coef[, 1]/coef[, 2])
coef.aux[,2]<-1-pnorm(coef[, 1]/coef[, 2])
for (i in 1:dim(coef)[[1]])
{coef.aux[i,3]<-min(coef.aux[i, 1],coef.aux[i, 2])}
coef[, 4] <- round(2*coef.aux[,3],6)
dimnames(coef) <- list(cnames, c("Estimate", "Std. Error", "z value", "p-value"))
# for dependence="ind"
if(all(is.na(match(cnames, "log.psi1"))) && all(is.na(match(cnames, "omega"))))
{
cat("\nFixed effects:\t\n")
print(coef[ , ])
}
# for dependence="MC1"
else if(!all(is.na(match(cnames, "log.psi1"))) && all(is.na(match(cnames, "log.psi2"))) &&
all(is.na(match(cnames, "omega"))))
{
cat("\nFixed effects:\t\n")
print(coef[1:(dim(coef)[[1]]-1), ])
cat("\nDependence parameter: \t\n")
table.r0<-matrix(coef[dim(coef)[[1]], ], nrow=1)
dimnames(table.r0)<-list("log.psi1", c("Estimate", "Std. Error", "z value", "p-value"))
print(table.r0)
}
# for dependence="MC2"
else if(!all(is.na(match(cnames, "log.psi1"))) && !all(is.na(match(cnames, "log.psi2"))) &&
all(is.na(match(cnames, "omega"))))
{
cat("\nFixed effects:\t\n")
print(coef[1:(dim(coef)[[1]]-2) , ])
cat("\nDependence parameter: \t\n")
table.r2<-data.frame(coef[(dim(coef)[[1]]-1):dim(coef)[[1]], ])
dimnames(table.r2)<-list(c("log.psi1","log.psi2"), c("Estimate", "Std. Error", "z value", "p-value"))
print(table.r2)
}
# for dependence="indR"
else if (all(is.na(match(cnames, "log.psi1"))) && !all(is.na(match(cnames, "omega"))))
{
cat("\nFixed effects:\t\n")
print(coef[-dim(coef)[[1]], ])
coef[dim(coef)[[1]], 1] <- exp(coef[dim(coef)[[1]], 1])
cat("\nRandom effects:\t\n")
table.r1<-data.frame(coef[dim(coef)[[1]], 1])
dimnames(table.r1)<-list("(Intercept)", "Variance")
print(table.r1)
}
# for dependence="MC1R"
else if (!all(is.na(match(cnames, "log.psi1"))) && all(is.na(match(cnames, "log.psi2"))) &&
!all(is.na(match(cnames, "omega"))))
{
cat("\nFixed effects:\t\n")
print(coef[1:(dim(coef)[[1]]-2), ])
cat("\nDependence parameter: \t\n")
table.r0<-matrix(coef[dim(coef)[[1]]-1, ], nrow=1)
dimnames(table.r0)<-list("log.psi1", c("Estimate", "Std. Error", "z value", "p-value"))
print(table.r0)
### for random coef
coef[dim(coef)[[1]], 1] <- exp(coef[dim(coef)[[1]], 1])
cat("\nRandom effect:\t\n")
table.r1<-data.frame(coef[dim(coef)[[1]], 1])
dimnames(table.r1)<-list("(Intercept)", "Variance")
print(table.r1)
}
# for dependence="MC2R"
else if (!all(is.na(match(cnames, "log.psi1"))) && !all(is.na(match(cnames, "log.psi2"))) &&
!all(is.na(match(cnames, "omega"))))
{
cat("\nFixed effects:\t\n")
print(coef[1:(dim(coef)[[1]]-3), ])
cat("\nDependence parameter: \t\n")
table.r2<-data.frame(coef[(dim(coef)[[1]]-2):(dim(coef)[[1]]-1), ])
dimnames(table.r2)<-list(c("log.psi1","log.psi2"),
c("Estimate", "Std. Error", "z value", "p-value"))
print(table.r2)
### for random coef
coef[dim(coef)[[1]], 1] <- exp(coef[dim(coef)[[1]], 1])
cat("\nRandom effect:\t\n")
table.r1<-data.frame(coef[dim(coef)[[1]], 1])
dimnames(table.r1)<-list("(Intercept)", "Variance")
print(table.r1)
}
if (cov)
{
cat("\nCovariance of Coefficients: \n")
print(object@covariance, digits = 2)}
if (cor){
cat("\nCorrelation of Coefficients: \n")
print(object@correlation, digits = 2)}
cat("\nMessage: ", object@message,"\n")
})
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.