Análisis Estadístico

Clase 1: Manipulación de datos en R

Importar base de datos creada en Excel. El documento se debe encontra en formato csv en una carpeta establecida para almacenar las bases de datos. Ejemplo: C:/MCF202/Data

• Utilizar el comando read.cvs

Ejemplo para importar datos

CedroRojo <- read.csv("/Volumes/TOSHIBA/Rstudio/Principios_estadistica/Data/cedrorojo.csv", header=T)

Usando R

Los datos en la matriz CedroRojo se encuentran ya disponibles. para revisar solamente los primeras seis filas empleamos el comando:

head(CedroRojo)

  diametro altura
1    12.97  22.64
2    12.37  20.44
3    12.37  19.94
4    12.67  14.94
5    10.67  15.84
6    10.07  16.34

Revisar: la media de la variable altura de la base de datos CedroRojo es: 18.9101 y la media para el diámetro es 12.524

mean(CedroRojo$altura)

[1] 18.91

Es importante utilizar el simbolo $ para especificar el lugar donde se encuentra.

mean(CedroRojo$diametro)

[1] 12.52

¿Cuál es el coeficiente de variación de las observaciones en altura?

sd(CedroRojo$altura)/mean(CedroRojo$altura)*100

[1] 15.91

Determinar cuantos individuos existen por debajo de la media o un valor establecido

• igual o mayor (>=)
• igual o menor (<=)
• mayor que (>)
• menor que (<)
• igual que (==)
• no igual (!=)

Crawley, J. (2012) The R Book, Chapter 1, pp. 27

¿Cuántos individuos poseen un diámetro mayor o igual a 2 cm?

sum(CedroRojo$diametro >=2)

[1] 167

¿Cuántos individuos poseen un diámetro menor a 10 cm?

sum(CedroRojo$diametro < 10)

[1] 7

¿Qué valor obtenemos aplicando el siguiente comando?

sum(CedroRojo$diametro < 12)

[1] 71

Respuesta: Obtenemos la sumatoria de las plántulas de Cedro rojo que tiene un diametro menor a 12 cm.

¿Qué valor obtenemos aplicando el siguiente comando?

sum(CedroRojo$diametro[CedroRojo$diametro < 12])

[1] 785.8

Respuesta: Obtenemos la sumatoria de todos los 71 diametros de las plántulas de Cedro rojo que tiene un diametro menor a 12 cm.

Para obtener una muestra de los datos colectados en campo bajo un criterio específico utilzamos el comando subset:

cedro.2 <-subset(CedroRojo, diametro < 12)
head(cedro.2)

   diametro altura
5     10.67  15.84
6     10.07  16.34
7     10.67  15.14
8     10.77  14.24
9     10.97  12.44
11    11.55  18.63

Para agregar una nueva variable (compuesta) a una base de datos aplicamos el comando transform.

cedro.3 <-transform(cedro.2, diam2=(diametro^2))
head(cedro.3)

   diametro altura diam2
5     10.67  15.84 113.9
6     10.07  16.34 101.4
7     10.67  15.14 113.9
8     10.77  14.24 116.0
9     10.97  12.44 120.4
11    11.55  18.63 133.5

Si tenemos una base de datos con datos agrupados por tratamientos, podemos seleccionar una muestra de un solo tratamiento. Ejemplo:

Seleccionar los datos del sitio 1 S1

inventarioS1 <-subset(inventario, Sitio=="S1") 
kable(head(inventarioS1))

Por favor seleccione los datos del sitio 2 y realice las siguientes actividades

1. Determinar la media de temperatura para el S2
2. determinar la media de HR
3. Determinar la media de HUMCOMB

Por favor seleccione los datos del sitio 2 y realice las siguientes actividades

inventarioS2 <-subset(inventario, Sitio=="S2") 
kable(head(inventarioS2))

1. Determinar la media de temperatura para el S2: 18.3884
2. determinar la media de HR: 80.7246
3. Determinar la media de HUMCOMB: 16.1159