El programa Canon escrito en lenguaje Fortran IV para IBM 360, lleva a cabo un análisis completo sobre un conjunto de hasta 15 variables y 30 grupos. Puede admitir 99.999 para cada variable

Campá Viñeta, Juan A.

Cuadras, C. M. (Carlos María)

Montoriol-Pous, Joaquín, 1924-2011

Spanish

Diposit Digital de la Universitat de Barcelona

ACTA GEOL~GICA HISPANICA, t. VI1 (1972), n." 1, págs. 26-28 El programa CANON para IBM 360(*) por CARLOS CUADRAS-AVELLANA **, JUAN A. CAMPÁ-VIÑETA ***, y JOAQUÍN MONTORIOL-POUS ***, El programa CANON, escrito en lenguaje Fortran IV para Como datos de entrada exige el número M de IBM 360, lleva a cabo un análisis canónico completo sobre un conjunto de hasta 15 variables y 30 grupos. Puede admitir número de grupos y número de ob- 99.999 observaciones para cada variable. servaciones de las M variables en cada uno de los grupos. A continuación siguen las observaciones de las M variables de cada grupo, es decir: primer indivi- Le progranime CANON, écrit en langage FOrtran IV pour IBM 360, efectue un analyse canonique complet sur un en- semble jusqu'i 15 variables et 30 groupes. 11 peut admetre 99.999 observations pour chaque variable. El programa CAN$N realiza un análisis canónico con~pleto sobre un conjunto de hasta 15 variables y 30 grupos. Cada variable puede admitir un número máximo de 99.999 observaciones, que puede ser am- pliado con s6l0 cambiar el formato de entrada. Esta notable capacidad es debida a que los datos son leídos observación a observación y solamente uno de ellos está almacenado en memoria en un momento dado. Los datos se almacenan en un fichero en disco, para ser leídos posteriormente en el programa. El programa utiliza 11 subrutinas, 5 de ellas pro- cedentes del S.S.P. (1) y las restantes originales del primero que suscribe. Están escritas en lenguaje F~RTRAN IV, por su carácter más universal y por dis- poner de las rutinas necesarias en este lenguaje. Está preparado para el sistema IBM 360 necesitando 64K de memoria y almacenamiento auxiliar de cinta o disco. START r READ K M N I I )  1 WRlTE DATOS m disco l D A T A )  7 l W M A T C  I l q m a a  W R l T E  ME91AS COMUNES Obtenclon Molrircr P y W WMATC llama o (wMnrc) CORRE L 1 3 1 COR;A;rxmo a WRlTE TEST Homog VAR -COV 1 1 WRlTE MATRIZ W / DATA 1.. datos disco 1 cotc~to o ( M A H A )  1 WRITE D y grados Libertad 7' Calculo EJES CANONICOS y valores p top lo r  E l G E N  ( NROOr ) COORDENADAS CANONICAS y ORlGEN W R I T E  RADIOS GRUPOS I " Este trabajo h a  sido realizado. en ~ a r t e .  con l a  A v u d a  oara e l  Fomento de l a  Investigación en l a  ' ~ n i v é r s i d á d .  "* Laboratorio de Cálculo. Facul tad de Ciencias. Univers idad de Bar:wona. duo, primer grupo, ídem, 2." individuo, l." grupo, Departamento de Cristalografía y Mineralogía, Un ive rs idad  de ., ídem 1 individuo, 2." grupo, etc. Barcelona. Sección de Mineralogía, Ins t i tu to  "Jaime Almera", C. S. de 1. C. Barcelona. El programa obtiene las medias comunes, las me- (1) IBM. Scientif ic Subroutine Package. 360A-CM-03X, versión III .  Programmer's M a n u a l  HZO.QZO~S-3. dias de grupos, el test de homogeneidad de varianzas- covarianzas, la pooled matriz de dispersión, la dis- tancia D W e  Mahalanobis, los ejes canónicos, las coordenadas canónicas de los grupos, el test de sig- nificación de los ejes canónicos y los radios de los grupos. El programa empieza leyendo los datos y que- dándoselos en disco, al tiempo que calcula la matriz de varianza-covarianza total (subrutina CORRE que llama a DATA) N T = (t,,) ; ti, = (Xijw - Xi) (Xsjr - Xe) h=l Seguidamente llama la subrutina WMATC, que vuelve a leer los datos, guardados en disco, calcu- lando las matrices de varianza-covarianza de cada gru- po y a partir de éstas la pooled dispersión matriz al tiempo que efectúa el test de homogeneidad, cal- culando el valor ji-cuadrado A continuación llama a la subrutina MAHA y ob- tiene la distancia (al cuadrado) de Mahalanobis. I1L VI? 7& Calcula luego, la matriz de varianzas-covarianzas "entre grupos" P = T - W  y llama a la subrutina N R O ~ T ,  para resolver la ecua- ción IP-4 WI = o  y poder determinar los ejes canbnicos. Una vez ob- tenidos éstos, se obtienen las coordenadas canónicas y el origen. Finalmente se llama a la subrutina BART para efectuar el test de significación de los ejes canó- nicos, calculándose los valores y sus correspondientes grados de libertad. Antes de leer los datos del siguiente problema, el programa principal, calcula y escribe los radios 1.645 / dnj de los grupos. El subprograma principal denominado CANON así como las subrutinas DATA, WMATC, MAHA y BART, han sido escritas y programadas por el primero que suscribe, usándose además, las subrutinas del S.S. P. (1) : CORRE NROOT EIGEN (valores y vectores propios de una ma- triz simétrica) MINV (inversión de matrices) LOc (localizaci6n de elementos en una ma- triz) MPRD (multiplica dos matrices) GMTRA (transposición de matrices) (1) IBM. Scientific Subroutine Package, 3,60-A-CM-03X, versión 111. Programmer's Manual H201.0205-3. TEST HOMOGENEIDAD V. JI C. 26.0541 G. LIB. POOLED DISPERSIÓN MATRIZ FILA 1 19.41879 - 11.162W - 5,21497 - 6.09890 F I L A  2 - 11.16207 11.94504 5.61813 1.91758 FILA 3- 5.21497 5.61813 39.45937 3.93681 FILA 4- 6fi9890 1.91758 3.93681 9.83310 FILA 5- 22.74860 22.6CQ86 16.23486 4.62157 FILA 6- 9.54052 10.66757 9.34546 383791 U-CUADRADO MAHALANOBIS GENERALIZADA 12.78068 G. LIB EJE V. JI. C. G.LIB. 1 4.1415 18 2 1.2670 10 3 0.4717 4 

El programa Canon para IBM 360

http://diposit.ub.edu/dspace/bitstream/2445/121719/1/655184.pdf

El programa Canon para IBM 360

Abstract

Similar works

Full text

Available Versions

Diposit Digital de la Universitat de Barcelona