-
Notifications
You must be signed in to change notification settings - Fork 0
/
10.R
126 lines (126 loc) · 4.46 KB
/
10.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
##
# Prácticas rituales y vida cotidiana. Expresiones de agencia colectiva a
# través de la construcción de símbolos en la comunidad local de Nueva Esperanza,
# siglos VI a.C. a XV d.C.
#
# Laura Velásquez González
# Universidad Nacional de Colombia
# 2022
#
# 10. Comparación cerámica
#
# Adaptado del código de Víctor González Fernández (2019)
#
# Cargar paquetes
#
library("RcmdrMisc")
library("plyr")
#
# Cargar las funciones
#
source("https://raw.githubusercontent.com/lavelasquezgo/nuevaesperanza/main/10a")
#
# Cargar los datos
#
Forma_total <- read.table("https://raw.githubusercontent.com/lavelasquezgo/nuevaesperanza/main/10b",
header = TRUE)
#
# Comparación de proporciones de formas cerámicas entre los grupos
#
round(prop.table(table(Forma_total$Forma)), 2)
#
par(mfrow=c(1,2))
with(Forma_total, balasP2(Grupo,Forma,1, ylim=c(0.0,0.02)))
with(Forma_total, balasP2(Grupo,Forma,3, ylim=c(0.0,0.04)))
par(mfrow=c(1,1))
#
par(mfrow=c(1,2))
with(Forma_total, balasP2(Grupo,Forma,2, ylim=c(0.4,0.7)))
with(Forma_total, balasP2(Grupo,Forma,4, ylim=c(0.3,0.6)))
par(mfrow=c(1,1))
#
# Comparación de proporciones de formas cerámicas del periodo Herrera entre los grupos
#
round(prop.table(table(Forma_total$Forma[Forma_total$Periodo == "Herrera"])),2)
#
par(mfrow=c(1,2))
with(Forma_Herrera, balasP2(Grupo,Forma,1, ylim=c(-0.002,0.02)))
with(Forma_Herrera, balasP2(Grupo,Forma,3, ylim=c(-0.002,0.04)))
par(mfrow=c(1,1))
#
par(mfrow=c(1,2))
with(Forma_Herrera, balasP2(Grupo,Forma,2, ylim=c(0.6,0.8)))
with(Forma_Herrera, balasP2(Grupo,Forma,4, ylim=c(0.2,0.4)))
par(mfrow=c(1,1))
#
# Comparación de proporciones de formas cerámicas del periodo Temprano entre los grupos
#
round(prop.table(table(Forma_total$Forma[Forma_total$Periodo == "Temprano"])),2)
#
par(mfrow=c(1,2))
with(Forma_Temprano, balasP2(Grupo,Forma,1, ylim=c(-0.001,0.03)))
with(Forma_Temprano, balasP2(Grupo,Forma,3, ylim=c(0.005,0.05)))
par(mfrow=c(1,1))
#
par(mfrow=c(1,2))
with(Forma_Temprano, balasP2(Grupo,Forma,2, ylim=c(0.3,0.7)))
with(Forma_Temprano, balasP2(Grupo,Forma,4, ylim=c(0.3,0.7)))
par(mfrow=c(1,1))
#
# Comparación de proporciones de formas cerámicas del periodo Tardío entre los grupos
#
round(prop.table(table(Forma_total$Forma[Forma_total$Periodo == "Tardio"])),2)
#
par(mfrow=c(1,2))
with(Forma_Tardio, balasP2(Grupo,Forma,1, ylim=c(-0.001,0.05)))
with(Forma_Tardio, balasP2(Grupo,Forma,3, ylim=c(0.005,0.06)))
par(mfrow=c(1,1))
#
par(mfrow=c(1,2))
with(Forma_Tardio, balasP2(Grupo,Forma,2, ylim=c(0.25,0.75)))
with(Forma_Tardio, balasP2(Grupo,Forma,4, ylim=c(0.25,0.75)))
par(mfrow=c(1,1))
#
# Comparación de los diámetros de las ollas
#
# Cargar datos
#
diametro_ollas <- read.table("https://raw.githubusercontent.com/lavelasquezgo/nuevaesperanza/main/10c",
header = TRUE)
diametro_ollas
#
# Resumen grupos
#
medias <- data.frame(Media = c(mean(diametro_ollas$Diametro[diametro_ollas$Grupo == 1]),
mean(diametro_ollas$Diametro[diametro_ollas$Grupo == 2]),
mean(diametro_ollas$Diametro[diametro_ollas$Grupo == 3]),
mean(diametro_ollas$Diametro[diametro_ollas$Grupo == 4])),
Desviacion = c(sd(diametro_ollas$Diametro[diametro_ollas$Grupo == 1]),
sd(diametro_ollas$Diametro[diametro_ollas$Grupo == 2]),
sd(diametro_ollas$Diametro[diametro_ollas$Grupo == 3]),
sd(diametro_ollas$Diametro[diametro_ollas$Grupo == 4]))
)
medias
#
# Comparar los diámetros por periodo
#
# Herrera
#
balasB(diametro_ollas$Grupo[diametro_ollas$Periodo == "Herrera"], diametro_ollas$Diametro[diametro_ollas$Periodo == "Herrera"], ylab = "Diámetro", xlab = "Grupo")+
text(x=1.3,y=23.5,cex=.7,label="Confianza: \n \u2590\u2588\ 80% \u2588 95% \u2590 99%", col=1)
#
# Temprano
#
balasB(diametro_ollas$Grupo[diametro_ollas$Periodo == "Temprano"], diametro_ollas$Diametro[diametro_ollas$Periodo == "Temprano"], ylab = "Diámetro", xlab = "Grupo")+
text(x=1.3,y=18.3,cex=.7,label="Confianza: \n \u2590\u2588\ 80% \u2588 95% \u2590 99%", col=1)
#
# Tardío
#
balasB(diametro_ollas$Grupo[diametro_ollas$Periodo == "Tardio"], diametro_ollas$Diametro[diametro_ollas$Periodo == "Tardio"], ylab = "Diámetro", xlab = "Grupo")+
text(x=1.3,y=18.3,cex=.7,label="Confianza: \n \u2590\u2588\ 80% \u2588 95% \u2590 99%", col=1)
#
numSummary(diametro_ollas$Diametro[diametro_ollas$Periodo == "Herrera"])
numSummary(diametro_ollas$Diametro[diametro_ollas$Periodo == "Temprano"])
numSummary(diametro_ollas$Diametro[diametro_ollas$Periodo == "Tardio"])
#
##