3.13-92

alexanderrobitzsch · Mar 16, 2023 · a7cb764 · a7cb764
1 parent 70e11a5
commit a7cb764
Show file tree

Hide file tree

Showing 885 changed files with 3,284 additions and 1,206 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,8 +1,8 @@
 Package: sirt
 Type: Package
 Title: Supplementary Item Response Theory Models
-Version: 3.13-35
-Date: 2023-02-10 15:14:59
+Version: 3.13-92
+Date: 2023-03-16 12:47:17
 Author: Alexander Robitzsch [aut,cre] (<https://orcid.org/0000-0002-8226-3132>)
 Maintainer: Alexander Robitzsch <[email protected]>
 Description: 

diff --git a/NAMESPACE b/NAMESPACE
@@ -129,6 +129,7 @@ importFrom(TAM, tam_outer)
 importFrom(TAM, tam_packageinfo)
 importFrom(TAM, tam_print_call)
 importFrom(TAM, tam_rsessinfo)
+importFrom(TAM, weighted_mean)
 importFrom(TAM, weighted_sd)
 
 
@@ -250,6 +251,7 @@ export(lsem.bootstrap)
 export(lsem.estimate)
 export(lsem.MGM.stepfunctions)
 export(lsem.permutationTest)
+export(lsem.test)
 export(lsem_local_weights)
 export(marginal.truescore.reliability)
 export(mcmc.2pno)
@@ -295,6 +297,7 @@ export(plausible.value.imputation.raschtype)
 export(polychoric2)
 export(pow)
 export(predict_scale_group_means)
+export(print_digits)
 export(prior_model_parse)
 export(prmse.subscores.scales)
 export(prob.guttman)

diff --git a/R/IRT.expectedCounts.mirt.R b/R/IRT.expectedCounts.mirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.expectedCounts.mirt.R
 ## File Version: 0.07
+## File Last Change: 2019-10-29
 
 # IRT.expectedCounts.mirt
 

diff --git a/R/IRT.expectedCounts_sirt.R b/R/IRT.expectedCounts_sirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.expectedCounts_sirt.R
 ## File Version: 0.16
+## File Last Change: 2020-12-07
 
 
 

diff --git a/R/IRT.factor.scores.sirt.R b/R/IRT.factor.scores.sirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.factor.scores.sirt.R
 ## File Version: 0.16
+## File Last Change: 2019-01-02
 
 
 ########################################################################

diff --git a/R/IRT.factor.scores.xxirt.R b/R/IRT.factor.scores.xxirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.factor.scores.xxirt.R
 ## File Version: 0.07
+## File Last Change: 2018-12-30
 
 
 ###########################################################

diff --git a/R/IRT.irfprob.mirt.R b/R/IRT.irfprob.mirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.irfprob.mirt.R
 ## File Version: 0.27
+## File Last Change: 2019-10-29
 
 
 #--- irfprob singleGroup

diff --git a/R/IRT.irfprob.sirt.R b/R/IRT.irfprob.sirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.irfprob.sirt.R
 ## File Version: 0.16
+## File Last Change: 2018-12-30
 
 ########################################################
 # convert theta.k into a matrix

diff --git a/R/IRT.likelihood.mirt.R b/R/IRT.likelihood.mirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.likelihood.mirt.R
 ## File Version: 0.284
+## File Last Change: 2019-10-29
 
 
 IRT_likelihood_mirt_singleGroup <- function( object, type, ... )

diff --git a/R/IRT.likelihood_sirt.R b/R/IRT.likelihood_sirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.likelihood_sirt.R
 ## File Version: 0.18
+## File Last Change: 2022-10-13
 
 
 ################################################

diff --git a/R/IRT.mle.R b/R/IRT.mle.R
@@ -1,10 +1,12 @@
 ## File Name: IRT.mle.R
-## File Version: 0.32
+## File Version: 0.34
+## File Last Change: 2023-03-08
 
 
 #--- IRT.mle
 IRT.mle <- function(data, irffct, arg.list, theta=rep(0,nrow(data)), type="MLE",
-        mu=0, sigma=1, maxiter=20, maxincr=3, h=0.001, convP=1e-04, maxval=9, progress=TRUE)
+        mu=0, sigma=1, maxiter=20, maxincr=3, h=0.001, convP=1e-04, maxval=9,
+        progress=TRUE)
 {
     N <- length(theta)
     I <- ncol(data)

diff --git a/R/IRT.modelfit.sirt.R b/R/IRT.modelfit.sirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.modelfit.sirt.R
 ## File Version: 0.20
+## File Last Change: 2018-12-30
 
 
 ###########################################################

diff --git a/R/IRT.posterior.mirt.R b/R/IRT.posterior.mirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.posterior.mirt.R
 ## File Version: 0.21
+## File Last Change: 2019-10-29
 
 
 

diff --git a/R/IRT.posterior_sirt.R b/R/IRT.posterior_sirt.R
@@ -1,5 +1,6 @@
 ## File Name: IRT.posterior_sirt.R
 ## File Version: 0.15
+## File Last Change: 2018-12-30
 
 
 

diff --git a/R/L0_polish.R b/R/L0_polish.R
@@ -1,5 +1,6 @@
 ## File Name: L0_polish.R
 ## File Version: 0.12
+## File Last Change: 2020-04-26
 
 
 L0_polish <- function(x, tol, conv=0.01, maxiter=30, type=1, verbose=TRUE)

diff --git a/R/L0_polish_one_iteration.R b/R/L0_polish_one_iteration.R
@@ -1,5 +1,6 @@
 ## File Name: L0_polish_one_iteration.R
 ## File Version: 0.354
+## File Last Change: 2019-01-14
 
 
 L0_polish_one_iteration <- function(x, tol, eps=0.01, maxiter=30, type=1)

diff --git a/R/L1_polish.R b/R/L1_polish.R
@@ -1,5 +1,6 @@
 ## File Name: L1_polish.R
 ## File Version: 0.26
+## File Last Change: 2019-01-14
 
 
 

diff --git a/R/L2_polish.R b/R/L2_polish.R
@@ -1,5 +1,6 @@
 ## File Name: L2_polish.R
 ## File Version: 0.07
+## File Last Change: 2019-01-14
 
 L2_polish <- function(x, eps=1e-4, maxiter=20)
 {

diff --git a/R/Probtrace_sirt.R b/R/Probtrace_sirt.R
@@ -1,5 +1,6 @@
 ## File Name: Probtrace_sirt.R
 ## File Version: 0.05
+## File Last Change: 2019-04-25
 
 
 # auxiliary function for mirt package

diff --git a/R/Q3.R b/R/Q3.R
@@ -1,5 +1,6 @@
 ## File Name: Q3.R
 ## File Version: 1.21
+## File Last Change: 2018-12-30
 
 ###############################################################
 # Yen's Q3 statistic (1984)

diff --git a/R/Q3.testlet.R b/R/Q3.testlet.R
@@ -1,5 +1,6 @@
 ## File Name: Q3.testlet.R
 ## File Version: 1.21
+## File Last Change: 2018-12-30
 
 
 ################################################################################

diff --git a/R/R2conquest.R b/R/R2conquest.R
@@ -1,5 +1,6 @@
 ## File Name: R2conquest.R
 ## File Version: 1.41
+## File Last Change: 2019-05-03
 
 
 

diff --git a/R/R2noharm-utility.R b/R/R2noharm-utility.R
@@ -1,5 +1,6 @@
 ## File Name: R2noharm-utility.R
 ## File Version: 1.14
+## File Last Change: 2018-12-30
 
 
 

diff --git a/R/R2noharm.EAP.R b/R/R2noharm.EAP.R
@@ -1,5 +1,6 @@
 ## File Name: R2noharm.EAP.R
 ## File Version: 0.26
+## File Last Change: 2019-04-30
 
 R2noharm.EAP <- function( noharmobj, theta.k=seq(-6,6,len=21), print.output=TRUE )
 {

diff --git a/R/R2noharm.R b/R/R2noharm.R
@@ -1,5 +1,6 @@
 ## File Name: R2noharm.R
 ## File Version: 2.24
+## File Last Change: 2018-12-30
 
 #------------------------------------------------------------------
 # NOHARM exploratory factor analysis

diff --git a/R/R2noharm.jackknife.R b/R/R2noharm.jackknife.R
@@ -1,5 +1,6 @@
 ## File Name: R2noharm.jackknife.R
 ## File Version: 1.13
+## File Last Change: 2018-12-30
 
 
 #**** jackknife R2noharm

diff --git a/R/RcppExports.R b/R/RcppExports.R
@@ -1,5 +1,6 @@
 ## File Name: RcppExports.R
-## File Version: 3.013034
+## File Version: 3.013092
+## File Last Change: 2023-03-16
 # Generated by using Rcpp::compileAttributes() -> do not edit by hand
 # Generator token: 10BE3573-1514-4C36-9D1C-5A225CD40393
 

diff --git a/R/Rhat_sirt.R b/R/Rhat_sirt.R
@@ -1,5 +1,6 @@
 ## File Name: Rhat_sirt.R
 ## File Version: 1.12
+## File Last Change: 2022-05-16
 
 
 ####################################

diff --git a/R/amh_plot.R b/R/amh_plot.R
@@ -1,9 +1,9 @@
 ## File Name: amh_plot.R
-## File Version: 0.32
+## File Version: 0.336
+## File Last Change: 2023-03-15
 
 
-######################################################
-# plot results of objects of class amh
+#*** plot results of objects of class amh
 amh_plot <- function( x, conflevel=.95, digits=3, lag.max=.1, col.smooth="red",
     lwd.smooth=2, col.split="blue", lwd.split=2, lty.split=1,
     col.ci="orange", cex.summ=1, ask=FALSE, ... )
@@ -33,7 +33,8 @@ amh_plot <- function( x, conflevel=.95, digits=3, lag.max=.1, col.smooth="red",
         sparm.vv <- smcmcobj[ smcmcobj$parameter==parm.vv, ]
 
         #*** traceplot
-        graphics::plot( iterindex, x.vv, type="l", main=paste0( "Traceplot of ", parm.vv ),
+        graphics::plot( iterindex, x.vv, type="l",
+                            main=paste0( "Traceplot of ", parm.vv ),
                             xlab="Iterations", ylab="", ... )
         x1 <- as.numeric( x.vv )
         xmin <- min(x1)
@@ -48,7 +49,7 @@ amh_plot <- function( x, conflevel=.95, digits=3, lag.max=.1, col.smooth="red",
             i1 <- seq(g1[ii]+1, g1[ii+1] )
             m1 <- mean( x1[ i1 ] )
             graphics::lines( iterindex[ i1 ], rep( m1, length(i1) ), col=col.split,
-                        lwd=lwd.split, lty=lty.split )
+                                lwd=lwd.split, lty=lty.split )
         }
         # include moving average here!!
         l1 <- sirt_moving_average(x1, B=round( lag.max / 2 ), fill=FALSE)
@@ -72,7 +73,8 @@ amh_plot <- function( x, conflevel=.95, digits=3, lag.max=.1, col.smooth="red",
         # blue dashed line at
         bd <- .05
         ylim <- c( min( acf1, - bd), 1 )
-        graphics::plot( iter_vv, acf1, xlab="Lag", ylab="ACF", main=mtitle, type="n", ylim=ylim)
+        graphics::plot( iter_vv, acf1, xlab="Lag", ylab="ACF",
+                            main=mtitle, type="n", ylim=ylim)
         NL <- length(iter_vv)
         for (hh in 1:NL){
             graphics::lines( rep( iter_vv[hh],2), c(0, acf1[hh]) )
@@ -88,26 +90,30 @@ amh_plot <- function( x, conflevel=.95, digits=3, lag.max=.1, col.smooth="red",
         x0 <- 0
         y0 <- 0
         heights.summ=c( .05,  .15, .25,  .35, .45, .55, .65, .75)
-        graphics::text( x0 + .0015, y0 + heights.summ[8], "Posterior Mean=", cex=cex.summ, pos=4)
+        graphics::text( x0 + .0015, y0 + heights.summ[8], "Posterior Mean=",
+                                cex=cex.summ, pos=4)
         graphics::text( x0 + .5, y0 + heights.summ[8],
-                paste0( sirt_format_numb( x=mean( x1 ), digits=round.summ)  ), pos=4 )
+                paste0( sirt_format_numb( x=mean( x1 ), digits=round.summ) ), pos=4 )
         hvv <- heights.summ[7]
         graphics::text( x0 + .0015, y0 + hvv, "Posterior Mode=", cex=cex.summ, pos=4)
         graphics::text( x0 + .5, y0 + hvv,
-                paste0( sirt_format_numb( x=sparm.vv$MAP, digits=round.summ)  ), pos=4 )
+                paste0( sirt_format_numb( x=sparm.vv$MAP, digits=round.summ) ), pos=4 )
 
-        graphics::text( x0 + .0015, y0 + heights.summ[6], "Posterior SD=", cex=cex.summ, pos=4)
+        graphics::text( x0 + .0015, y0 + heights.summ[6], "Posterior SD=",
+                            cex=cex.summ, pos=4)
         graphics::text( x0 + .5, y0 + heights.summ[6],
-                paste0( sirt_format_numb( x=stats::sd( x1 ), digits=round.summ)  ), pos=4 )
+                paste0( sirt_format_numb( x=stats::sd( x1 ), digits=round.summ) ), pos=4 )
 
         hvv <- heights.summ[5]
         graphics::text( x0 + .0015, y0 + hvv,
-                            paste( round(100*conflevel ), "% Credibility Interval=",sep=""),
-                            cex=cex.summ, pos=4 )
+                        paste( round(100*conflevel ), "% Credibility Interval=",sep=""),
+                        cex=cex.summ, pos=4 )
 
         hvv <- heights.summ[4]
-        ci.lower <- sirt_format_numb( stats::quantile( x1, ( 1 - conflevel  ) / 2 ), digits=round.summ )
-        ci.upper <- sirt_format_numb( stats::quantile( x1, 1- ( 1 - conflevel  ) / 2 ), digits=round.summ )
+        ci.lower <- sirt_format_numb( stats::quantile( x1, ( 1 - conflevel  ) / 2 ),
+                                            digits=round.summ )
+        ci.upper <- sirt_format_numb( stats::quantile( x1, 1-( 1 - conflevel  ) / 2 ),
+                                            digits=round.summ )
         graphics::text( x0 + .25, y0 + hvv,
                             paste( "[", ci.lower,    ",", ci.upper, "]",  sep=""),
                             cex=cex.summ, pos=4)
@@ -117,11 +123,14 @@ amh_plot <- function( x, conflevel=.95, digits=3, lag.max=.1, col.smooth="red",
                 paste0( sirt_format_numb( x=sparm.vv$Rhat, digits=3)  ), pos=4 )
         hvv <- heights.summ[2]
         graphics::text( x0 + .0015, y0 + hvv, "SERatio=", cex=cex.summ, pos=4)
-        graphics::text( x0 + .5, y0 + hvv, paste0( sirt_format_numb( x=sparm.vv$SERatio, digits=3)  ), pos=4 )
+        graphics::text( x0 + .5, y0 + hvv,
+                        paste0( sirt_format_numb( x=sparm.vv$SERatio, digits=3)), pos=4 )
 
         hvv <- heights.summ[1]
-        graphics::text( x0 + .0015, y0 + hvv, "Effective Sample Size=", cex=cex.summ, pos=4)
-        graphics::text( x0 + .705, y0 + hvv, paste0( sirt_format_numb( x=sparm.vv$effSize, digits=1)  ), pos=4 )
+        graphics::text( x0 + .0015, y0 + hvv, "Effective Sample Size=",
+                        cex=cex.summ, pos=4)
+        graphics::text( x0 + .705, y0 + hvv,
+                        paste0( sirt_format_numb( x=sparm.vv$effSize, digits=1)), pos=4 )
         graphics::par(ask=ask)
     }
     graphics::par(mfrow=c(1,1))

diff --git a/R/anova_sirt.R b/R/anova_sirt.R
@@ -1,5 +1,6 @@
 ## File Name: anova_sirt.R
 ## File Version: 0.23
+## File Last Change: 2018-12-30
 
 ##############################################################
 # anova rasch.mml

diff --git a/R/attach.environment.sirt.R b/R/attach.environment.sirt.R
@@ -1,5 +1,6 @@
 ## File Name: attach.environment.sirt.R
 ## File Version: 0.07
+## File Last Change: 2018-12-30
 ##################################################
 # attach all elements of an object in an environment
 .attach.environment.sirt <- function( res, envir ){

diff --git a/R/automatic.recode.R b/R/automatic.recode.R
@@ -1,5 +1,6 @@
 ## File Name: automatic.recode.R
 ## File Version: 1.18
+## File Last Change: 2018-12-30
 
 #*** automatic recoding of a dataset
 automatic.recode <- function( data, exclude=NULL, pstart.min=.6,

diff --git a/R/bounds_parameters.R b/R/bounds_parameters.R
@@ -1,5 +1,6 @@
 ## File Name: bounds_parameters.R
 ## File Version: 0.08
+## File Last Change: 2018-12-30
 
 
 bounds_parameters <- function( pars, lower=NULL, upper=NULL)

diff --git a/R/brm.irf.R b/R/brm.irf.R
@@ -1,5 +1,6 @@
 ## File Name: brm.irf.R
 ## File Version: 0.11
+## File Last Change: 2019-08-04
 
 #--- item response function (discretized) beta response model
 brm.irf <- function( Theta, delta, tau, ncat, thdim=1, eps=1E-10 )

diff --git a/R/brm.sim.R b/R/brm.sim.R
@@ -1,5 +1,6 @@
 ## File Name: brm.sim.R
 ## File Version: 0.11
+## File Last Change: 2018-12-30
 
 
 ##################################################################

diff --git a/R/btm.R b/R/btm.R
@@ -1,5 +1,6 @@
 ## File Name: btm.R
 ## File Version: 1.535
+## File Last Change: 2020-04-18
 
 
 #--- Bradley-Terry model in sirt

diff --git a/R/btm_fit_combine_tables.R b/R/btm_fit_combine_tables.R
@@ -1,5 +1,6 @@
 ## File Name: btm_fit_combine_tables.R
 ## File Version: 0.07
+## File Last Change: 2019-03-09
 
 
 btm_fit_combine_tables <- function( win1, win2, ind1, ind2, TP )

diff --git a/R/btm_fit_statistics.R b/R/btm_fit_statistics.R
@@ -1,5 +1,6 @@
 ## File Name: btm_fit_statistics.R
 ## File Version: 0.209
+## File Last Change: 2022-04-17
 
 
 #**** item outfit and infit statistic

diff --git a/R/btm_sim.R b/R/btm_sim.R
@@ -1,5 +1,6 @@
 ## File Name: btm_sim.R
 ## File Version: 0.03
+## File Last Change: 2019-03-10
 
 btm_sim <- function(theta, eta=0, delta=-99, repeated=FALSE)
 {