-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathExplorandoMaiz.R
56 lines (29 loc) · 1.3 KB
/
ExplorandoMaiz.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
###QUESTION 1 OF 4 (TEST 7.1)
##Cargando archivo `Prac_Uni7/maices/meta/maizteocintle_SNP50k_meta_extended.txt`
##y guardándolo en fullmat
fullmat<-read.delim("../meta/maizteocintle_SNP50k_meta_extended.txt", header=TRUE, sep="\t")
#Tipo de objeto creado al cargar la base
class(fullmat)
#Viendo las primeras 6 líneas del archivo
head(fullmat, n=6L)
#Visualizando el número de muestras
length(fullmat$NSiembra)
#Obteniendo el número de estados para los que se tienen muestras
nlevels(fullmat$Estado)
#Muestras colectadas antes de 1980
x<-subset(fullmat, A.o._de_colecta<1980, na.rm = TRUE)
length(x$A.o._de_colecta)
#Muestras de cada raza
summary(fullmat$Raza)
#Promedio de altitud en la que fueron colectadas las muestras
mean(fullmat$Altitud, na.rm = TRUE)
#Altitud máxima y mínima en las que fueron colectadas
max(fullmat$Altitud, na.rm = TRUE)
min(fullmat$Altitud, na.rm = TRUE)
#Creando nueva data.frame de datos sólo con las muestras de la raza Olotillo
Olotillo<-subset(fullmat, Raza=="Olotillo")
#Creando nueva data.frame de datos sólo con las muestras de la raza Reventador, Jala y Ancho
RJA<-fullmat[fullmat$Raza %in% c("Reventador", "Jala", "Ancho"), ]
# Escribiendo la matriz anterior a un archivo llamado submat.cvs en /meta
write.csv(RJA, file="../meta/submat.cvs")
sessionInfo()