*SPSS/WIN 6,01.
* GENERADOR DE PUNTUACIONES CORRELACIONADAS, MEDIANTE LA
* DESCOMPOSICION DE CHOLESKY.
* Copyright: F.J.Herrero (Dpto.de Psicologia-Universidad de Oviedo-1995)
*-----------------------------------------------------------------------
*Informacion:
*         Este algoritmo te permite generar puntuaciones correlacionadas
* mediante la descomposicion de Cholesky.
*
* Su uso es gratuito, puedes copiarlo y utilizarlo.
*
* Agradeceria que si lo utilizas referenciases al autor.
*
*Nota:
*     Para utilizarlo basta con que cambies los parametros referentes a:
* -Numero de vectores a generar.
* -Numero de casos a generar.
* -Patron de la correlacion que deseas generar.
* -Parametros de media y desviacion tipica que deseas en cada vector.
*
* Es facil de identificar ya que en los comentarios aparece el codigo
* Modificable
* indicando que esas lineas son las que tu puedes modificar segun tus
* necesidades de simulacion.
*
*Si tienes algun problema puedes contactar con el autor por e-mail en
*                    herrero@pinon.ccu.uniovi.es
*
*.

NEW FILE.
INPUT PROGRAM.
/* NUMERO DE VECTORES A GENERAR (ej:3)  [Modificable].
VECTOR X(3).
/* NUMERO DE CASOS A GENERAR (ej:1000) [Modificable].
 LOOP #J=1 TO 1000.
   LOOP #I=1 TO 3.
     compute X(#I)=normal(1)+0. /* GENERACION DE VECTORES NORMALIZADAS.
   END LOOP.
   END CASE.
 END LOOP.
END FILE.
END INPUT PROGRAM.

*ALGORITMO MATRICIAL.
MATRIX.
* RECOGE LOS VECTORES NORMALIZADOS GENERADOS ANTERIORMENTE.
GET V/FILE *.

* MATRIZ R BASE EN LA GENERACIONS DE LOS NUMEROS [Modificable].
COMPUTE R={1.0,0.8,0.6;
           0.8,1.0,0.1;
           0.6,0.1,1.0}.


COMPUTE D1=DET(R).
PRINT D1/TITLE "DETERMINANTE SOBRE R".
* DESCOMPOSICION DE CHOLESKY DE LA MATRIZ ANTERIOR.
COMPUTE C=CHOL(R).
COMPUTE D2=DET(C).
PRINT D2/TITLE "DETERMINANTE SOBRE C".

* CALCULO DE LAS PUNTUACIONES CORRELACIONADAS.
COMPUTE D=V*C.
SAVE D/OUTFILE *.
END MATRIX.

* ASIGNACION DEL CENTRO Y DISPERSION DESEADOS [Modificable].
COMPUTE Y=(COL1*2)+5.
COMPUTE X1=(COL2*15)+100.
COMPUTE X2=(COL3*15)+100.

* VERIFICACION DE LA MATRIZ DE CORRELACIONES GENERADA.
CORR Y X1 X2/STAT DESC.