Un programa informatico para el analisis funerario en arqueologia.
Pascual Izquierdo Egea.
Universidad Aut6noma de Barcelona.
Apartado 56.O8211 Castellar del Vallés (Barcelona).
NECROes un programa experimental de analisis funerario. Incorpora una serie de técnicas instrumentales contempladas ennuestra estrategia metodol6gica,destinadaacanalizar la contrastacion de una teoria intermedia que relacione el entramado general del materialismo histórico con la especificidad empirica del contexto estudiado.El objetivo fundamental de NECRO es medir el nivel de apropiación material en el ritual de enterramiento,mediante la aplicacion de un metodo de valoracion objetiva de las tumbas oestructuras funerarias.Asimismo,la introduccion de otras tecnicas instrumentales,extraidas de la Estadistica inductiva (Test de la chi-cuadrado,rde Pearson) y de la Taxonomia numerica (analisis de agrupamientos con enlace medio),cumple los propósitos de contrastar sus resultados con los obtenidos a traves del procedimiento analitico principal.
1．Presupuestos teóricos ymetodol6gicos.
El materialismo histórico constituye una alternativa valida para solventar la actual problematica de la interpretacion arqueol6gica.Noobstante,ala hora de abordar el contexto funerario,precisa desarrollar una teoria intermedia que relacione el entramado general con la especificidad empirica del marco referencial. Esta estrategia,reconociendo que la observacion del registro arqueológico permite constatar,en ciertos casos, diferencias materiales entre las tumbas, deberá admitir que existe una relacion directa entre dicha variabilidad, en especial la referida a la composicion del ajuar,y el distinto acceso de los individuos al ritual de enterramiento.Ala vez, éste lltimo constituye un caso particular de las manifestaciones de la ideologia, la cual justifica y reproduce las relaciones sociales de producci6n (HOUTART,1989:11), manteniendo la cohesion interna de Ia formacion.Por lo tanto, la diferenciación material en el acceso al ritual funerario es proporcional a la asimetria de las citadas relaciones sociales de produccion, en cuya intensificacion interviene decisivamente la division social del trabajo. Sólo en tales términos puede entenderse la complejidad social. Sin embargo, la lectura objetiva del contexto material de los enterramientos es limitada.Sus resultados solamente adquieren validez cuando se contrastan dialécticamente con los extraidos del registro arqueol6gico restante,especialmente de los asentamientos,ademas de otras fuentes de información como la historiografia o la etnografia.
La objetivacion del analisis funerario requiere la delineacion de una via de contrastacion empirica coherente con los presupuestos esgrimidos anteriormente. Nuestra propuesta metodológica pretende cumplir dichas condiciones. En primer lugar, seguiremos seleccionando una muestra suficientemente representativa de tumbas intactas,separändolas,si ello es factible,segun su cronologia. Despues,se aislarán las dimensiones significativas que expresen la variabilidad cualitativo-cuantitativa de la diferenciacion material.A continuacion, efectuaremos una valoracion objetiva de las variables discretas de los conjuntos funerarios,para conocer el nivel de apropiacion individual en el acceso al ritual de enterramiento,asumiendo que el contexto material es la clave que posibilita la aproximacion a la cuestión (IZQUIERDO, 1989: 58). Seguidamente, contrastaremos los resultados de la valoracion contextual con el analisis de Ias asociaciones,asi como con la distribucion espacial de las sepulturas. Posteriormente,el estudio sincronico de cada área funeraria sera confrontado con la lectura socioeconómica de los asentamientos relacionados. Asimismo,se evaluarán los cambios identificados en el acceso material al ritual funerario,mediante una sintesis de los resultados de todas las necr6polis conocidas pertenecientes a una misma formacion econ6mica y social. Este analisis diacronico se contrastara,nuevamente,con el conjunto de los poblados asociados.Unicamente asi, tras conocer la dinamica interna de la sociedad estudiada,podremos acometer un analisis comparado de la misma con otras formaciones.
Por otro lado,la cuantificacion de las diferencias materiales entre las sepulturas en funcion de la variabilidad de sus componentes, como medio para alcanzar una valoracion objetiva de los enterramientos, sigue siendo una cuesti6n abierta pendiente de resolución.No obstante,la Arqueologia procesual, desde una perspectiva sistémico-funcionalista, protagoniz6 algunos intentos, aunque combinaron procedimientos simples con otros completamente subjetivos.Además, todos ellos incurrieron en el error de pretender medir el grado de “riqueza>,cuando este término ambiguo aun no ha sido definido conceptualmente y,en consecuencia, no deberia emplearse para determinar el nivel de jerarquizacion. Concretamente,RANDSBORG (1973: 567-569) recurrió al peso de los objetos de bronce y oro de los ajuares como un indice de aquella. En cambio, RATHJE (1973: 745-746) emple6 la rareza de la ocurrencia de cada categoria o clase de artefacto para puntuar las 1oo9 tumbas mayas clasicas que integraban su muestra. SHEPHARD (1979: 54-56), por su parte,utilizó tanto el anterior estadigrafo como el propuesto por HODSON(1977). Este ultimo consiste en valorar cada tipo de item en funcion del promedio de tipos que hallamos en las tumbas donde aparece.Finalmente,ORTONy HODSON(1981: 106,107),aplicando la Teoria de Probabilidades y expresiones matematicas complejas derivadas de la misma,como la f6rmula de la distribuci6n de Pareto, investigaron la agrupacion de las sepulturas segun los grados de asignados a ciertos componentes funerarios (tipos funcionales),reconociendo su escasa operatividad en muestras de tamaio reducido, y la necesidad de combinar este procedimiento con otros alternativos para superar sus limitaciones (ORTON y HODSON,1981: 114). Por otro lado, el problema apenas ha despertado interes en la Peninsula, si bien se ha reconocido su existencia (LULL y ESTEVEZ,1986: 448; LULL y PICAZO,1989: 18). Quienes han pretendido llegar mas lejos,se han limitado a llevar a cabo una encuesta entre especialistas para asignar un “valor>arbitrario a los artefactos (GIMENO y GONZALEZ, 1986).
En nuestra metodologia',el valor contextual de una variable discreta funeraria,o de un componente del ajuar (VCi),es una medicion aproximada de su importancia material en la manifestacion ritual. Una primera aproximacion a su calculo fue aplicada a los tumulos A y B de la necropolis orientalizante de Setefilla(IZQUIERDO,1989:67-68).Posteriormente, este procedimiento se perfeccion6 (IZQUIERDO, 1990),resultando una propuesta mejor razonada para alcanzar la determinación del VCi,a partir de la muestra seleccionada.En si, ésta dependerä de las frecuencias absolutas de cinco factores: items (N), tumbas (T), item considerado (Ni), items asociados (Nia)y tumbas donde aparece el item considerado (Ti). Asumimos que el VCi es funcion de la rareza o exclusividad del item,es decir, de Ia inversa de su frecuencia relativa $\left( \mathbf { N } / \mathbf { N } \mathbf { i } \right)$ ,del promedio de items $\left( \left( \mathbf { N i + N i a } \right) / \mathrm { T i } \right)$ y de la frecuencia relativa (Ti/T) de las tumbas donde aparece.Entonces,multiplicando las anteriores expresiones estadisticas,obtenemos que,
$$.
V C i = \frac { N \left( N i + N i a \right) } { - \mathcal { T } \Psi i }
$$.
Simplificando, conseguimos una formulacion matematica mas clara:
$$.
V C i = \frac { N } { T } ~ \left( l + \frac { N i a } { N i } \right)
$$
donde N y T corresponden a dos parametros o constantes, mientras que Niay Ni son variables discretas. Por lo tanto, el valor contextual de un item funerario dependerä directamente del $\mathtt { n } ^ { \mathtt { 2 } }$ total de items asociados con él en todas las tumbas donde aparece e, inversamente,de su frecuencia absoluta en la muestra.Finalmente,la suma de los VCide los componentes funerariosdela tumba,esdecir,suindice contextual (IC),medir& objetivamente la apropiacion material de cada individuo en el acceso al ritual de enterramiento.Esta conceptualizacion debera substituir,en adelante,a los términos tradicionales de 《riqueza>O funerarias.
2.Descripci6n del programa.
NECRO es un programa experimental de analisis funerario destinado a cumplir los propositos enumerados anteriormente.Por ello,su objetivo fundamental es medir el nivel de apropiacion material (NAM) en el ritual de enterramiento,aplicando nuestro metodo de valoraci6n objetiva de las tumbas (MEVACON),cuyos resultados son contrastados mediante la incorporacion de otros procedimientos analiticos extraidos de la Estadistica inductiva (Test de asociación de la chi-cuadrado,coeficiente de correlacion lineal de Pearson） y de la Taxonomia numérica (analisis de agrupamientos con enlace medio).NECRO es el resultado provisional de la colaboracion entre la empresa Minimicro, S.A.², responsable del desarrollo informatico,y el autor del proyecto. Realizado empleandoel sistema de programacion Turbo Pascal de Borland Inc., su ejecución puede efectuarse desde cualquier ordenador compatible IBMPC/XT/AT/386o PS/2 provisto del sistema operativo MS/DOS 3.3 de Microsoft. Ahora bien, es recomendable la instalacion del programa en un disco duro (HD),la utilizacion de un microprocesador rapido -80286 6 80386-,asi como la presencia de una memoriaRAMelevadaparaagilizarel tratamiento del notable volumen de datos generado.Igualmente, respecto a la visualizaci6n por pantalla, es preferible recurrir a un sistema que incorpore una tarjeta grafica de alta resolucion(VGA) junto a un monitor en color.
El menu principal de NECRO (fig.1) presenta 12 opciones operativas.Las siete primeras (A-G) corresponden a la gestión de la base de datos central, integrada porlos ficheros de las necrópolis,las cuales acumulan,a su vez,y en forma de fichas, los datos referentes ala composicionde las tumbas seleccionadas inicialmente.De este modo, podemos crear un fichero para una nueva necrópolis o seleccionar una existente (opcion A); introducir, anadir (opcion B), modificar (opcion C)o eliminar (opcion D) fichas de tumbas; borrar (opcion E) o copiar (opcion F) un fichero de necrópolis; y listar los ficheros de necr6- polis almacenados en el disco (opcion G).

Fig. 2.
A partirde labasede datos,accedemos al algoritmo que realiza el cälculo de los valores contextuales de los items (VCi) que aparecen en las estructuras funerarias de la necrópolis considerada, seguido por la determinaci6n de los indices contextuales de las tumbas (ICt) en funcion de aquellos (opcion I). De este modo,aplicando la formulacion propuesta anteriormente,obtenemos dos listados por pantalla o impresora. El primero contiene N,T,Ni,Ti, Nia y VCi (fig.2),ademas de los VCi minimo, medio y maximo.Existe la posibilidad de visualizar graficamente Ios resultados finales (impresora).Por su parte,el segundo listado (fig. 3) muestra ios Nit (numero total o frecuencia absoluta de items) e ICt para cada tumba,junto a la composicion cualitativa de la misma (pantalla e impresora),asi como los valores minimo, medio y maximo de ambas variables.También se calcula y muestra el coeficiente de correlacion lineal (r de Pearson) existente entre dichos factores.Finalmente,la opcion concluye con la confeccion de una gräfica bidimensional que representa los valores del Nit en el eje de ordenadas y los del ICt en el de abscisas (fig.5), visualizable por pantallade color o impresora.Su utilidad reside en la posibilidad de efectuar,posteriormente,un analisis comparado de las distribuciones (nubes de puntos) perfiladas por ambas variables para un conjunto dado de necrópolis. Por otro lado,en todo momento podemos examinar el inventario de cada fichero funerario(fig.4),es decir, la matriz de abundancia que muestra las frecuencias absolutas de los componentes de cada tumba (opcion H).

Fig. 4.

Fig. 5.
Junto al anterior procedimiento,se han incorporado otras técnicas estandarizadas,cuyos resultados se contrastan con los obtenidos mediante el metodo de valoración contextual(MEVACON). Concretamente,la opci6n J nos permite acceder al algoritmo mas complejo de NECRO que, basado en la Taxonomia numérica,lleva a cabo una agrupacion de las tumbas en función de la composicion cualitativo-cuantitativa de sus ajuares,midiendo una distancia espacial expresada por el coeficiente de similitud (CS). Tras calcular la matriz de similaridad--cuyos resultados pueden imprimirse-,y mediante la técnica de enlace jerarquico más compleja del “Cluster Analysis> -《average linkage>(JOHNSON, 1972:335-339, 353-355； DORAN y HODSON, 1975: 177; MALDONADO,1980: 64;ORTON,1988:48- 57)-,traza un dendrograma (fig.6) que expresa graficamente la proximidad entre los conjuntos funerarios en funcion de la intersimilaridad de sus composiciones (pantallae impresora).El “average linkage analysis>,o analisis de agrupamientos con enlace medio a partir de multiples variables,proporciona resultados parecidos a los logrados mediante la técnica del enlace completo o total, aunque el primer procedimiento ofrece mayores dificultades, tanto si se realiza manuaimente como a través de un ordenador.En éste ultimo caso,la memoria ocupada por el proceso de datos aumenta exponencialmente en proporcion directa al incremento del tamaio de la muestra (num.de tumbas),lo cual retardarä considerablemente el tiempo de espera. Tales dificultades son superadas por la tecnica de agrupamiento no jerarquico de Ias $\ll \mathbf { K }$ -medias>, puesto que al calcular las distancias requeridas al mismo tiempo que efectua una comparacion,no precisa almacenar en memoria todos los coeficientes (ORTON,1988: 56). No obstante,el principal inconveniente de este flexible procedimiento subyace en la division arbitraria de la totalidad de casos iniciales en una cantidad determinada (K) de grupos,cuestion que contradice la necesaria objetivacion del proceso analitico de la base empirica disponible.Por nuestra parte,hemos desestimado su incorporación al programa NECRO,donde la metodologia no persigue, fundamentalmente, una finalidad taxonomica.
El procedimiento de agrupamiento con enlace medio posee un caracter provisional(CONTRERAS, 1984: 347) y,en nuestro caso,su utilidad obedece a la contrastacion del método de valoracion contextual (MEVACON) con otras técnicas de aproximacion. En cuanto a la mecanica de su funcionamiento,el algoritmo no sigue exactamente los pasos del calculo manual.Asi,en éste,a partir de las composiciones de los conjuntos funerarios,archivadas en la base de datos central, se determina una matriz de similaridad con T(T-1)/2 coeficientes operativos (JOHNSON, 1972: 312; BURTON1973: 193).Despues, el proceso sigue los conocidos criterios de unir los pares de tumbas con CS mäs elevado,partiendo del maximo valor (1OO). De este modo, una vez se constituya la primera agrupaci6n binaria,sera necesario comparar los CS de las demás tumbas con el promedio de aquella, y asi sucesivamente en el caso de los conglomerados ternarios,cuaternarios,etc. Siempre que el elemento o grupo a enlazar con el de referencia muestre mayor CS medio, pasarä a formar parte del mismo.La complejidad de esta técnica viene expresada por el elevado numero de comparaciones a efectuar, traducido matematicamente como $\sum \ n ( \mathbf { T } - \mathbf { n } )$ ，es decir,la suma de las combinaciones posibles de n elementos,donde n toma valores desde 2 hasta n-1.En cambio,en lenguaje informatico,la complejidad comparativa viene dada por el numero de iteraciones,determinado por Ia expresión $\mathrm { T } ( \mathrm { T } { \cdot } 1 ) ^ { 2 } / 2$.
Por otro lado,Ia opcion K (fig.7) nos ofrece la posibilidad de establecer las asociaciones entre items funerarios que alcanzan significacion estadistica, mediante el conocido Test de la chi-cuadrado de Pearson (DOWNIE y HEATH,1971: 215-219; DORAN y HODSON,1975: 54-55; VIEDMA,1976: 254-256; SOKAL y ROHLF,1980: 294,296; SANAHUJA,1985: 65-66; DOMENECH y RIBA, 1987:52),criticado por algunos (CONTRERAS, 1984: 347),pero valido siempre y cuando sus resultados no se consideren definitivos.Es decir, deberan combinarse con la valoracion contextual de las tumbas.Este algoritmo extrae de la base de datos principal todas las categorias o clases de items presentes en la necrópolis considerada.A partir del listado que aparece en pantalla,se estudian parejas de componentes,aplicando la f6rmula de la chi-cuadrado para tablas de contingencia ${ \bf 2 } \bf { x } { 2 } .$ Cuando el valor absoluto de (ad-bc) es inferior o igual a $\bar { \mathbf { n } } / 2$ ，siendo $\scriptstyle \mathbf { n = a + b + c + d }$ ,se utiliza.
Fig. 7.
$$
\chi ^ { 2 } { = } \frac { n ( a d { - } b c ) ^ { 2 } } { ( a { + } b ) ( c { + } d ) ( a { + } c ) ( b { + } d ) }
$$
mientras que si supera dicho valor, se introduce la correccion de Yates en la anterior expresi6n (DOWNIE y HEATH,1971: 218; VIEDMA,1976: 255).En consecuencia, solamente cambiara el numerador de la citada formulacion,convirtiendose en $\mathbf { n } ( \mathbf { a d - b c - n } / 2 ) ^ { 2 }$ .Finalmente,aparece en pantalla la tabla de contingencia con los efectivos para cada par de items,en funcion de su interasociacion,y el resultado alcanzado por la $\mathsf { X } ^ { 2 }$ ,el cual debera superar el valor 3,841 tomando un grado de libertad y un nivel de significacion O,O5.Tal condici6n es interpretada por el algoritmo, que emite un mensaje relativo a la misma.
Finalmente,la opcion L(analisis de regresion) constituye un subprograma, con base de datos propia y gestion de la misma,destinado a calcular el coeficiente de correlacion lineal de Pearson (r) existente entre dos variables dadas (fig.8).La formula empleada es
(DOWNIE y HEATH,1971: 108; VIEDMA,1976: 265).Su principal utilidad consiste en la aplicacion al analisis espacial de las tumbas,cuya ubicacion está definida por sus coordenadas cartesianas $\bf ( x , y )$ en el plano.Tambien se determina el centro geometrico de gravedad de la distribucion,es decir,el punto cuyas dimensiones corresponden a la media aritmética de cada dimension euclidiana.
3．Algunas aplicaciones.
El desarrollo del programa NECRO ha implicado un ensayo permanente de sus posibilidades instrumentales,sometiendo a analisis multiples muestras funerarias seleccionadas a partir de otras tantas necrópolis.Como consecuencia, si bien con caracter provisional, esta proporcionando resultados significativos que podrian validar la linea de investigación seguida mediante el MEVACON(metodo de valoración contextual). Por ejemplo, el estudio diacrónico de los tumulos Ay B de la necropolis orientalizante de Setefilla (Lora del Rio, Sevilla), excavados completamente p0I AUBET (1975,1978,1980-81), muestra un descenso del $23 \%$ en el NAM (nivel de apropiacion material) medido para la segunda superestructura funeraria,a través del ICt medio de ambos. Esta constatación ya fue apuntada anteriormente en otros terminos (IZQUIERDO,1989),a la vez que se propuso la hipótesis de relacionarla con los cambios socioeconómicos derivados de la crisis tartésica del s. VI a.C.Posteriormente,la revision cronol6gica de los dos tumulos e'sta confirmando nuestra propuesta, puesto que el Ase fecha en pleno s. VII a.C., mientras que el B se situaria entre finales de la anterior centuria $\mathbf { y }$ la primera mitad del s. VI a.C.(AUBET,1989: 302).Por otro lado,en el caso de las necrópolis ibéricas tarraconenses de Mas de Mussois (Tortosa) y Mianes (Santa Barbara),yaunque correspondan a un contexto socioeconómico diferente,apreciamos un comportamiento analogo.Asi, se detecta una reducción del $1 7 \%$ en el NAMde Mianes comparandolo con el de la primera,datada a mediados del s.VI a.C.,y solapandose con la segunda poco despues (inicios del s.V a.C.) (MALUQUER, 1984: 76; id. 1987:57; GUSI y OLIVER,1986: 272). Esta observación empirica,apreciadaanteriormente (IZQUIERDO y GIMENO,1990a),podria vincularse con las transformaciones experimentadas en el seno de las comunidades ibéricas asentadas en toda la zona levantina (IZQUIERDO y GIMENO,1990b)，en curso de estudio.

Fig.8.
NECRO dejarä pronto su fase experimental para convertirse en una potente herramienta metodologica al servicio del anälisis funerario.Para lograrlo,además de perfeccionar la actual version del programa, incorporaremos nuevas técnicas instrumentales coherentes con nuestra estrategia de contrastacion.
BIBLIOGRAFIA.
AUBET SEMMLER,M.E. 1975 La necrópolis de Setefilla en Lora del Rio, Sevilla (tumulo A).Programa de Investigaciones Protohistoricas (PIP) 2.
AUBET SEMMLER,M.E. 1978 La necrópolis de Setefilla en Lora del Rio, Sevilla (tumulo B).Programa de Investigaciones Protohistoricas (PIP)3.
AUBET SEMMLER,M.E. 1980-81 “Nuevos hallazgos en la necrópolis de Setefilla (Sevilla)",Mainake 2-3:87-115.
AUBET SEMMLER,M.E. 1989 “La Mesa de Setefilla: la secuencia estratigrafica del Corte 1",en Aubet (ed.): Tartessos,pp.297- 338.
AUBET SEMMLER,M.E. (ed.) 1989Tartessos.Arqueologia protohistorica del Bajo Guadalquivir.Sabadell.
BURNHAM, B.C.; KINGSBURY,J. (eds.) 1979 Space,Hierarchy and Society.Interdisciplinary Studies in Social Area Analysis.British Archaeological Reports (International Series) 59.
BURTON, M.L. 1973“Mathematical Anthropology>,Annual Review ofAnthropology 2: 189-199.
CLARKE, D.L. (ed.) 1972Models in Archaeology.London.Methuen.
CONTRERASCORTES,F. 1984《Clasificacion y tipologia en Arqueologia. El camino hacia la cuantificacion>,Cuadernos de Prehistoria de la Universidadde Granada9: 327-385.
DOMENECH,J.M.; RIBA, M.D. 1987 Una sintesis de los metodos estadisticos bivariantes.Barcelona.
DORAN, J.E.; HODSON, F.R. 1975 Mathematicsand Computersin Archaeology. Edimburgo. Edinburgh U.P.
DOWNIE, N.M.; HEATH, R.W. 1971Métodos estadisticos aplicados.Madrid.
GIMENO VIGUERA, I.; GONZALEZ BRAVO,R. 1986 Una interpretacion microespacial: necropolis tumular>,Arqueologia Espacial9:387-401.
GUSI, F.; OLIVER, A. 1987 “Laproblematicade la iberizacion en Castell6n>, enIberos,Actas de lasI Jornadas sobre el Mundo Iberico (Jaen,1985),pp.99-136.
HOUTART, F. 1989 Religion y modosde produccion precapitalistas. Madrid.
HUMPHREYS,S.C.; KING,H. (eds.) 1981Mortalityand Immortality:theanthropology and archaeology of death.London. Academic Press.
IZQUIERDO EGEA,P. 1989 Elhorizonte orientalizante en el Mediterráneo occidental.Aproximacion a la reconstruccion economicay socialdelas comunidades tartésicas de los siglos VIII-VI a.C.Trabajo de Investigacion de Doctorado inedito.Universidad Aut6noma de Barcelona.
IZQUIERDO EGEA,P. 1990 “Materialismo histórico y analisis funerario en Arqueologia>,comunicaci6n presentada al I CongresodeJovenes Ge6grafoseHistoriadores (Sevilla,noviembre 1990).
IZQUIERDO EGEA, P.; GIMENO I FABREGAT,T. 1990aLes fortificacions iberiques dels segles V-III a.C.a les comarques del Baix Ebre>,comunicacion presentada al Simposi Internacional d'Arqueologia Iberica:Fortificacions.La problematica de l'lbericple (segles IV-Illa.C.) (Manresa,diciembre 1990).
IZQUIERDO EGEA,P.; GIMENO IFABREGAT,T. 1990bLa societat iberica del Valles. Monografies Vallesanes (en prensa).
JOHNSON, L. Jr. 1972 Introduction to imaginary models for archaeologicalscaling and clustering>,en Clarke, D.L.(ed.)Models inArchaeology,pp.309-379.
LULL,V.; PICAZO, M. 1989 “Arqueologia de la Muerte y estructura social>, Archivo Espanol de Arqueologia 62:5-20.
LULL,V.; ESTEVEZ, J. 1986 “Propuesta metodol6gica para el estudio de las necropolisargaricas>,en Homenaje aLuis Siret, pp. 441-452.
MALDONADO, R. 1980 Ofrendas asociadas a entierrosdel Infiernillo enel Balsas.Estudio y experimentacion contres metodos de taxonomia numerica.México.
MALUQUER DE MOTES,J. 1984 La necrópolis paleoiberica de “Mas de Mussols>, Tortosa (Tarragona). Programa de Investigaciones Protohistoricas (PIP) 8.
MALUQUER DE MOTES, J. 1987 La necrópolis paleoibérica de Mianes en Santa Barbara (Tarragona).Programa de Investigaciones Protohist6ricas (PIP) 9.
ORTON, C. 1988Matematicas para arque6logos.Madrid.Alianza.
ORTON,C.R.; HODSON, F.R. 1981 《Rank and Class: interpreting the evidence from prehistoric cemeteries>,en Humphreys-King (eds.):Mortality and Immortality,pp.103-115.
RANDSBORG, K. 1973 “Wealth and social structure as reflected in bronze age burials- a quantitative approach>,en Renfrew (ed.):The Explanation of Culture Change,pp.565-570.
RATHJE,W.L. 1973 《Models for mobile Maya:a variety of constraints>,en Renfrew,C.(ed.): The Explanation ofCulture Change,pp.731-757.
RENFREW,C.(ed.) 1973 TheExplanation ofCulture Change.Models in Prehistory.London.Duckworth.
SANAHUJA YLL,M.E. 1985 《Renunciaconscientala subjectivitat:l'aplicació d'alguns metodes quantitatius a l'arqueologia>, Cota Zero 1: 57-74.
SHEPHARD, J. 1979 “The Social Identity of the individual in isolated barrows and barrow cemeteries in Anglo-Saxon England>,enBurnham-Kingsbury(eds.):Space, Hierarchy and Society,Pp.47-79.
SOKAL, R.R.; ROHLF, F.J. 1980Introduccion a la Bioestadistica.Barcelona.
VIEDMA CASTANO, J.A. 1976 Bioestadistica (Métodos Estadisticosen Medicinay Biologia).Madrid: