search for: davood

...rforms analysis, and computes a statistic 1000 times. Therefore, for each condition, I will have 1000 statistic values. My question is what's the best way to store the 1000 statistic for each condition. Any suggestion on how to manage such simulation studies is greatly appreciated. Thanks, -- Davood Tofighi Department of Psychology Arizona State University [[alternative HTML version deleted]]

Davood, I developed an MC simulation model for wind hazard analysis last year. I found three important issues to increase efficiency: 1) Reuse most variables in each loop 2) Write results (1000 stats) to external files, perhaps one file for each condition. (a good ID for each file can be implemented usin...

Hello, II used the logLik() function to get the log-likelihood estimate of an object. The function also prints the degrees of freedom. How can I extract the degrees of freedom and assign it to a variable. Below is the output: > logLik(fit2pl) 'log Lik.' -4842.912 (df=36) Thanks, Davood Tofighi [[alternative HTML version deleted]]

...en I ran the same code, I got the following error: > dbCreate("simdb") #create simdb database [1] TRUE > db<-dbInit("simdb") #initiate an object of database Error in sprintf(gettext(fmt, domain = domain), ...) : object "datafile" not found > Thanks, Davood [[alternative HTML version deleted]]

Hello All, I was wondering how I can get the overall Pearson chi^2 test of model fit with its df and p value in the LTM package for the 2PL models. Thanks, -- Davood Tofighi Department of Psychology Arizona State University [[alternative HTML version deleted]]

...s of the other variables? require(adapt) fn<-function(z,y,x) {exp(x^3+y^4+x^2*y+x*z^2+y/z)} x<-runif(200);y<-runif(200);z<-runif(200); c1<-.5;c2<-5;M<-100; #M to represent infinity i1<-adapt(3,lo=c(.0001,0,0),up=c(M,c1*z,c2/y),functn=fn)$value print(i1); Thanks, -- Davood Tofighi [[alternative HTML version deleted]]