# NOT RUN {
dat <- sim10GDINA$simdat
Q <- sim10GDINA$simQ
Q[10,] <- c(0,1,0)
mod1 <- GDINA(dat = dat, Q = Q, model = "GDINA")
pvaf <- Qval(mod1,method = "PVAF",eps = 0.95)
pvaf
#If many entries are modified, you may want to check
#the PVAF plot using the function plotPVAF or
#to change eps. eps = 0.9 or 0.8 seems another two
#reasonable choices.
extract(pvaf,what = "PVAF")
#See also:
extract(pvaf,what = "varsigma")
extract(pvaf,what = "sug.Q")
# Draw a mesa plot
plot(pvaf,item=10,type="best",no.qvector=5)
stepwise <- Qval(mod1,method = "wald")
stepwise
#If many entries are modified, you may want to check
#the PVAF plot using the function plotPVAF or
#to change eps. eps = 0.9 or 0.8 seems another two
#reasonable choices.
extract(stepwise,what = "PVAF")
#See also:
extract(stepwise,what = "varsigma")
extract(stepwise,what = "sug.Q")
# Draw a mesa plot
plot(stepwise,item=10,type="best",no.qvector=5)
# }
Run the code above in your browser using DataLab