2015-11-23 28 views
1

我想知道是否有人可以幫助我。 我想對我的數據進行差異分析,但不斷得到一個theta.sparse錯誤,我不確定如何解決。我真的很感激你能給我的任何東西。theta.sparse錯誤與lorDIF

library(lordif) 
dat<- read.csv2("OPSO.csv",header=TRUE) 
datgender <- read.csv2("DATA.csv",header=TRUE) 
group<-datgender$Gender 
sink("outputDIFopso.txt") 
gender.difopso <- lordif(dat, group, selection = NULL, 
criterion = c("Chisqr", "R2", "Beta"), 
pseudo.R2 = c("McFadden", "Nagelkerke", "CoxSnell"), alpha = 0.01, 
beta.change = 0.1, R2.change = 0.02, maxIter = 10, minCell = 5, 
minTheta = -4, maxTheta = 4, inc = 0.1, control = list(), model = "GRM", 
anchor = NULL, MonteCarlo = FALSE, nr = 100) 
print(gender.difopso) 
summary(gender.difopso) 
sink() 
pdf("graphtestop.pdf") 
plot(gender.difopso) 
dev.off() 
dev.off() 

錯誤lordif(DAT,組,選擇= NULL,標準= C( 「Chisqr」: 對象theta.sparse'未找到

謝謝:)

回答

0

您應該檢查之前的錯誤行,輸出可能會說您沒有標記爲DIF的項目,如果是這種情況,您應該運行mirt函數並提取thetaipar對象作爲nec埃森。

作者可能compare(flags, flags.matrix)爲真時添加一些案件處理。最起碼,它似乎當有與DIF沒有項目以同樣的方式,它說

if (ndif == ni) { 
     warning("all items got flagged for DIF - stopping\n") 
    } 

並沒有辦案時(ndif == 0)雖然compare(flags, flag.matrix)評估爲TRUE省略了警告。

的影響時,項目的全或無有DIF的是,你會被合併後的樣品(無DIF)或兩個以上mirt在裝修mirt得到相同的結果(產生相同的ICC曲線,同樣的推理等)每組的模型(全部DIF)。所以這是一個正確的節省時間的程序,當所有這些發生故障時都能繞過。