Skip to contents

Summary method for lawbl objects

Source: R/summary.lawbl.R
summary.lawbl.Rd

Provide summaries of posterior information for a lawbl object, .

Usage

# S3 method for lawbl
summary(
  object,
  what = "basic",
  med = FALSE,
  SL = 0.05,
  detail = FALSE,
  digits = 4,
  istart = 1,
  iend = -1,
  ...
)

Arguments

object

A lawbl object

what

A list of options for what to summarize.

  • basic: Basic information about the model and posteriors.

  • lambda: Loading estimates.

  • qlambda: Loading estimates in pattern/Q-matrix format.

  • eigen: Factorial eigen value.

  • dpsx: Diagonal elements in the residual covariance matrix PSX.

  • offpsx: Off-diagonal elements in PSX; local dependence terms.

  • phi: Factorial correlations.

  • thd: Threshold estimates.

  • int: Intercept estimates (for pcirm only).

  • shrink: (Ave) shrinkage for each factor's loadings and LD (if LD in pcfa = T).

  • factor: Are the factors true of not (for EFA).

  • all: All above information.

med

logical; if the posterior median (TRUE) or mean (FALSE) is used as the estimate.

SL

Significance level for interval estimate. The default is .05.

detail

logical; if only significant (FALSE) or all (TRUE) estimates are presented.

digits

Number of significant digits to print when printing numeric values.

istart

Starting point of the Markov chain for summary.

iend

Ending point of the Markov chain for summary; -1 for the actual final point.

...

additional arguments

Value

A list or matrix containing the summarized information based on the option what.

Examples

# \donttest{
dat <- sim18cfa0$dat
J <- ncol(dat)
K <- 3
Q<-matrix(-1,J,K);
Q[1:2,1]<-Q[7:8,2]<-Q[13:14,3]<-1

m0 <- pcfa(dat = dat, Q = Q, LD = FALSE,burn = 1000, iter = 1000)
summary(m0) # summarize basic information
#> $NJK
#> [1] 1000   18    3
#> 
#> $`Miss%`
#> [1] 0
#> 
#> $`LD Allowed`
#> [1] FALSE
#> 
#> $`Burn in`
#> [1] 1000
#> 
#> $Iteration
#> [1] 1000
#> 
#> $`No. of sig lambda`
#> [1] 24
#> 
#> $Selected
#> [1] TRUE TRUE TRUE
#> 
#> $`Auto, NCONV, MCONV`
#> [1]  0  0 10
#> 
#> $EPSR
#>      Point est. Upper C.I.
#> [1,]     1.3778     2.5388
#> [2,]     1.0232     1.0579
#> [3,]     1.3948     2.2712
#> 
#> $`DIC, BIC, AIC`
#> [1] 4859.694 2843.553 2475.471
#> 
#> $Time
#>    user  system elapsed 
#>    8.03    0.00    8.05 
#> 
summary(m0, what = 'lambda') #summarize significant loadings
#>       Item F     est     sd   lower   upper sig
#>  [1,]    1 1  0.6814 0.0488  0.5736  0.7669   1
#>  [2,]    2 1  0.6971 0.0468  0.6117  0.7996   1
#>  [3,]    3 1  0.7011 0.0477  0.6117  0.7983   1
#>  [4,]    4 1  0.7240 0.0469  0.6305  0.8133   1
#>  [5,]    5 1  0.6980 0.0553  0.5878  0.7918   1
#>  [6,]    6 1  0.6711 0.0551  0.5548  0.7588   1
#>  [7,]   16 1 -0.1978 0.0947 -0.3649 -0.0192   1
#>  [8,]   17 1  0.1647 0.0922  0.0042  0.3265   1
#>  [9,]    5 2  0.2922 0.0693  0.1639  0.4138   1
#> [10,]    6 2  0.3096 0.0673  0.1870  0.4250   1
#> [11,]    7 2  0.7324 0.0409  0.6534  0.8110   1
#> [12,]    8 2  0.7277 0.0428  0.6527  0.8154   1
#> [13,]    9 2  0.7393 0.0417  0.6618  0.8217   1
#> [14,]   10 2  0.7208 0.0413  0.6499  0.8113   1
#> [15,]   11 2  0.7230 0.0462  0.6425  0.8247   1
#> [16,]   12 2  0.7131 0.0429  0.6304  0.7991   1
#> [17,]   11 3  0.2850 0.0647  0.1487  0.4047   1
#> [18,]   12 3  0.2689 0.0618  0.1488  0.3875   1
#> [19,]   13 3  0.7351 0.0491  0.6341  0.8276   1
#> [20,]   14 3  0.7115 0.0503  0.6214  0.8205   1
#> [21,]   15 3  0.7481 0.0518  0.6532  0.8524   1
#> [22,]   16 3  0.7429 0.0540  0.6358  0.8493   1
#> [23,]   17 3  0.7402 0.0512  0.6355  0.8329   1
#> [24,]   18 3  0.7443 0.0511  0.6369  0.8335   1
summary(m0, what = 'qlambda') #summarize significant loadings in pattern/Q-matrix format
#>           1      2      3
#> I1   0.6814 0.0000 0.0000
#> I2   0.6971 0.0000 0.0000
#> I3   0.7011 0.0000 0.0000
#> I4   0.7240 0.0000 0.0000
#> I5   0.6980 0.2922 0.0000
#> I6   0.6711 0.3096 0.0000
#> I7   0.0000 0.7324 0.0000
#> I8   0.0000 0.7277 0.0000
#> I9   0.0000 0.7393 0.0000
#> I10  0.0000 0.7208 0.0000
#> I11  0.0000 0.7230 0.2850
#> I12  0.0000 0.7131 0.2689
#> I13  0.0000 0.0000 0.7351
#> I14  0.0000 0.0000 0.7115
#> I15  0.0000 0.0000 0.7481
#> I16 -0.1978 0.0000 0.7429
#> I17  0.1647 0.0000 0.7402
#> I18  0.0000 0.0000 0.7443
summary(m0, what = 'offpsx') #summarize significant LD terms
#> NULL
# }