Displaying 2 results from an estimated 2 matches for "thres1".
Did you mean:
thresh
2008 Jan 28
0
Grid search: avoid "for" loops thanks to apply&co?
...ues one has to compute the
sum of squares of all models, which takes n x n/2 operations for a model
with two regimes (for my case: 300x150=45'000).
Transcription of the Matlab code into R is inefficient because of the
loops (very slow), a schema of the actual code gives:
for(in in 1:length(thres1))
firstThresh<-tresh1[i]
for(j in 1:length(thres1))
if(j>i)
secondThresh<-tresh1[j]
regime1<-ifelse(x<thresh1)*x
regime2<-ifelse(x>thresh2)*x
...(some matrix algebra in order to obtai...
2012 Nov 08
0
mirt vs. eRm vs. ltm vs. winsteps
...d(mirt[order(mirt$loc,decreasing=TRUE),],2)
mirt #I get the following order of item parameters: x6,x5,x1,x4,x7,x3,x2
library(eRm)
res<-PCM(pcmdat)
pres<-person.parameter(res)
thres<-thresholds(res)
erm<-thres$threshtable
erm<-as.data.frame(erm)
names(erm)<-c("loc","thres1","thres2","thres3")
rownames(erm)<-c("x1","x2","x3","x4","x5","x6","x7")
erm<-round(erm[order(erm$loc,decreasing=TRUE),],2)
erm<-erm[c(2:4,1)]
erm #I get the following order of item parameters: x...