[R] Simulating dataset using Parallel Latent CTT model?

William Revelle lists at revelle.net
Sat Dec 27 18:00:57 CET 2008


As I said in my previous response, if you want 20 items with equal 
factor loadings you can just specify the loading: e.g., for all 
loadings of .6
GenData <- congeneric.sim(N=500, loads = rep(.6,20), short = FALSE)

The error variances will all be equal to 1 - loading^2  where loading 
is the factor loading you specified.

This is R.  If you want to know how a function works,  list the source code.


At 9:47 AM -0500 12/27/08, Nidhi Kohli wrote:
>Hi Bill,
>Thank you very much for your response. You are right, I want to 
>simulate data set for 500 examinees across 20 items using Parallel 
>latent CTT model. As you know in Parallel Latent CTT model all the 
>error variances and factor loadings are equal across all the items. 
>Could you please let me know how can I incorporate common error 
>variance for 20 items in this R-program?
>Nidhi Kohli
>Nidhi Kohli, M.Ed.
>Doctoral Student
>Department of Measurement, Statistics
>  and Evaluation
>University of Maryland
>1230 Benjamin Building
>College Park, MD 20742-1115
>e-mail: nidhik at umd.edu
>---- Original message ----
>  >Date: Fri, 26 Dec 2008 11:31:10 -0600
>  >From: William Revelle <lists at revelle.net> 
>>Subject: Re: [R] Simulating dataset using Parallel Latent CTT model? 
>>To: Nidhi Kohli <nidhik at umd.edu>, r-help at r-project.org
>>Presumably, you are trying to simulate 20 items all sharing one
>>general factor but having some error.
>>The model as you specified it has no error.  Thus all the
>>correlations will be 1 and the factors will not make any sense.
>>Most items have loadings on a general factor of  the  order of about
>>.4 to .6.  You might try:
>  >GenData <- congeneric.sim(N=500, loads = rep(.5,20), short = FALSE)
>>Then you will find that the factor scores  found by factor.pa
>>correlate at .93 with the latent variable.
>>>   FactorScore=factor.pa(GenData$observed,1,scores = "TRUE", rotate="none")
>>>   round(cor(FactorScore$scores,GenData$latent),2)
>>At 8:27 AM -0500 12/26/08, Nidhi Kohli wrote:
>>>I am trying to simulate a dataset using Parallel Latent CTT model
>>>and this is what i have done so far:
>>>#Importing psych library for all the simulation related functions
>>># Settting the working directory path to C:/NCME
>>>#Using the function to generate the data
>>>GenData <- congeneric.sim(N=500, loads =
>>>c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1), short = FALSE)
>>>#Rounding upto 2 decimal places while showing the correlation matrix
>>>#Factor Score
>>>FactorScore=factor.pa(GenData$observed,1,scores = "TRUE", rotate="none")
>>>Please let me know if I am moving into the right direction, if not
>>>then, please let me know the correct way to simulate the dataset
>>>Thanks in Advance
>>>Nidhi Kohli
>>>Nidhi Kohli, M.Ed.
>>>Doctoral Student
>>>Department of Measurement, Statistics
>>>   and Evaluation
>>>University of Maryland
>>>1230 Benjamin Building
>>>College Park, MD 20742-1115
>>>e-mail: nidhik at umd.edu
>>>R-help at r-project.org mailing list
>>>PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>>and provide commented, minimal, self-contained, reproducible code.
>>William Revelle		http://personality-project.org/revelle.html
>>Professor			http://personality-project.org/personality.html
>>Department of Psychology             http://www.wcas.northwestern.edu/psych/
>>Northwestern University	http://www.northwestern.edu/
>>Attend  ISSID/ARP:2009               http://issid.org/issid.2009/

William Revelle		http://personality-project.org/revelle.html
Professor			http://personality-project.org/personality.html
Department of Psychology             http://www.wcas.northwestern.edu/psych/
Northwestern University	http://www.northwestern.edu/
Attend  ISSID/ARP:2009               http://issid.org/issid.2009/

More information about the R-help mailing list