@@ -1656,6 +1656,7 @@ mirt2traditional <- function(x, vcov, nfact){
1656
1656
opar <- par <- x @ par
1657
1657
which.a <- which(x @ par [1L : nfact ] != 0 )
1658
1658
if (length(which.a ) != 1L ) return (x )
1659
+ a.nms <- if (nfact == 1L ) ' a' else paste0(' a' , 1L : nfact )
1659
1660
if (cls != ' GroupPars' )
1660
1661
ncat <- x @ ncat
1661
1662
if (cls == ' dich' ){
@@ -1683,7 +1684,7 @@ mirt2traditional <- function(x, vcov, nfact){
1683
1684
par [nfact + 1L ] <- - par [nfact + 1L ]/ par [which.a ]
1684
1685
par [nfact + 2L ] <- plogis(par [nfact + 2L ])
1685
1686
par [nfact + 3L ] <- plogis(par [nfact + 3L ])
1686
- names(par ) <- c(paste0( ' a ' , 1L : nfact ) , ' b' , ' g' , ' u' )
1687
+ names(par ) <- c(a.nms , ' b' , ' g' , ' u' )
1687
1688
} else if (cls == ' graded' ){
1688
1689
fns <- vector(' list' , ncat + nfact - 1L )
1689
1690
for (i in 2L : ncat - 1L ){
@@ -1702,7 +1703,7 @@ mirt2traditional <- function(x, vcov, nfact){
1702
1703
par [i + nfact ] <- - par [i + nfact ]/ par [which.a ]
1703
1704
delta_index [[i + nfact ]] <- c(which.a , i + nfact )
1704
1705
}
1705
- names(par ) <- c(paste0( ' a ' , 1L : nfact ) , paste0(' b' , 2 : ncat - 1L ))
1706
+ names(par ) <- c(a.nms , paste0(' b' , 2 : ncat - 1L ))
1706
1707
} else if (cls == ' gpcm' ){
1707
1708
fns <- vector(' list' , ncat + nfact )
1708
1709
for (i in 2L : ncat - 1L ){
@@ -1739,7 +1740,7 @@ mirt2traditional <- function(x, vcov, nfact){
1739
1740
}
1740
1741
delta_index [[nfact + 1L ]] <- c(which.a , ncat + nfact + 2L )
1741
1742
par <- c(x @ par [1 : nfact ], newd )
1742
- names(par ) <- c(paste0( ' a ' , 1 : nfact ) , paste0(' b' , 1 : length(newd )))
1743
+ names(par ) <- c(a.nms , paste0(' b' , 1 : length(newd )))
1743
1744
x @ est <- x @ est [c(1 : nfact , (ncat + nfact + 2L ): length(x @ est ))]
1744
1745
} else if (cls == ' nominal' ){
1745
1746
if (nfact > 1L ) return (x )
@@ -1766,7 +1767,7 @@ mirt2traditional <- function(x, vcov, nfact){
1766
1767
ds <- par [(ncat + 2 ): length(par )]
1767
1768
ds <- ds - mean(ds )
1768
1769
par <- c(as , ds )
1769
- names(par ) <- c(paste0( ' a ' , 1 : ncat ) , paste0(' c' , 1 : ncat ))
1770
+ names(par ) <- c(a.nms , paste0(' c' , 1 : ncat ))
1770
1771
x @ est <- rep(TRUE , ncat * 2 )
1771
1772
x @ SEpar <- rep(as.numeric(NA ), ncat * 2 )
1772
1773
} else if (cls == ' nestlogit' ){
@@ -1812,7 +1813,7 @@ mirt2traditional <- function(x, vcov, nfact){
1812
1813
par1 [2 ] <- - par1 [2 ]/ par1 [1 ]
1813
1814
par1 [3 ] <- plogis(par1 [3 ])
1814
1815
par1 [4 ] <- plogis(par1 [4 ])
1815
- names(par1 ) <- c(' a ' , ' b' , ' g' , ' u' )
1816
+ names(par1 ) <- c(a.nms , ' b' , ' g' , ' u' )
1816
1817
par2 <- par [5 : length(par )]
1817
1818
as <- par2 [1 : (ncat - 1 )]
1818
1819
as <- as - mean(as )
0 commit comments