## Not run: ------------------------------------
# dat <- sim10GDINA$simdat
# Q <- sim10GDINA$simQ
# Q[10,] <- c(0,1,0)
# mod1 <- GDINA(dat = dat, Q = Q, model = "GDINA")
# out <- Qval(mod1,eps = 0.95)
# out
# #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(out,what = "PVAF")
# #See also:
# extract(out,what = "varsigma")
# extract(out,what = "sug.Q")
#
# # Draw a mesa plot
# mesaplot(out,item=10,type="all",no.qvector=5)
## ---------------------------------------------
Run the code above in your browser using DataLab