ESTADÍSTICA COMPUTACIONAL

Profesor: ELMER RODRÍGUEZ STÉVEZ

ACTIVIDAD 1. CONOCIMIENTOS BÁSICOS DE R


Pregunta 1. ¿Qué diferencia existe entre las funciones help() y help.search()? Explique brevemente.

R/. La diferencia que existe entre las funciones help() y help.search() es que la función help() es utilizada para obtener cualquier tipo de ayuda, información, descripción general y la lista de argumentos para el buen uso de una función en específico, mientras que la función help.search() es utilizada en los casos en donde nos ha olvidado el nombre exacto de una función o cuando se desconoce en qué paquete se encuentra un conjunto de datos (data.frames) o simplemente queremos ver un tema o un concepto en específico.


Pregunta 2. Si el comando más corto para help(“mean”) es ?mean. ¿Cuál será para el de help.search()?

R/. El comando más corto para help.search es utilizando dos signos de interrogación seguidos, se escribe de la siguiente forma “??”.


Pregunta 3. La función rnorm genera un vector aleatorio numérico de una distribución normal. ¿Cuáles son los argumentos de dicha función?

R/. Los argumentos de la función rnorm son n(tamaño de la muestra), mean=0(la media) y sd=1(la desviación estándar).


Pregunta 4. Cuando escribes el comando help.search(correlation) recibes un mensaje de error. ¿Por qué? Corríjalo.

R/. Cuando se usa la función help.search(), se debe escribir entre comillas la expresión que va dentro del paréntesis, por lo tanto en el ejemplo sería help.search(“correlation”)


Pregunta 5. Verifique si en R existen los siguientes paquetes {base} y {matlab}. Si alguno de estos NO se encuentra instálelo. Copie y pegue la salida de ésta descarga.

R/.

install.packages(base,dep=TRUE) Error in install.packages : object ‘base’ not found > install.packages(“BASE”) WARNING: Rtools is required to build R packages but is not currently installed. Please download and install the appropriate version of Rtools before proceeding:

https://cran.rstudio.com/bin/windows/Rtools/ Installing package into ‘C:/Users/Start_Dell/Documents/R/win-library/3.6’ (as ‘lib’ is unspecified) Warning in install.packages : package ‘BASE’ is not available (for R version 3.6.1) Warning in install.packages : package ‘BASE’ is a base package, and should not be updated


install.packages(matlab,dep=TRUE) Error in install.packages : object ‘matlab’ not found WARNING: Rtools is required to build R packages but is not currently installed. Please download and install the appropriate version of Rtools before proceeding:

https://cran.rstudio.com/bin/windows/Rtools/ Installing package into �C:/Users/Start_Dell/Documents/R/win-library/3.6� (as �lib� is unspecified) trying URL ‘https://cran.rstudio.com/bin/windows/contrib/3.6/matlab_1.0.2.zip’ Content type ‘application/zip’ length 174095 bytes (170 KB) downloaded 170 KB

package ‘matlab’ successfully unpacked and MD5 sums checked

The downloaded binary packages are in C:_Dell8NUJ_packages


Pregunta 6. ¿Qué función dentro de R, permite cambiar la ubicación de tu directorio de trabajo?

R/. Para cambiar la ubicación del directorio de trabajo en R se utiliza la función setwd().


Pregunta 7. Edite y muestre la salida de los cambios realizados de los siguientes conjuntos de datos de la biblioteca {MASS}, {Rabbit} y {Animals} de acuerdo a:

  1. En Rabbit: Cambiar en la variable Dose, los valores de las dos primeras filas por 8.35 y 54.5.

fix("Rabbit")

Rabbit
  1. En Animals: Cambia los valores de la fila Gorilla por 250 y 480 respectivamente.

fix("Animals")

Animals

Pregunta 8. En la biblioteca {MASS}, invoque el conjunto de datos Traffic. Verifique que variables contiene dicho conjunto. Explique qué significa cada una de ellas y el experimento como tal.

library(MASS)
Traffic
##     year day limit  y
## 1   1961   1    no  9
## 2   1961   2    no 11
## 3   1961   3    no  9
## 4   1961   4    no 20
## 5   1961   5    no 31
## 6   1961   6    no 26
## 7   1961   7    no 18
## 8   1961   8    no 19
## 9   1961   9    no 18
## 10  1961  10    no 13
## 11  1961  11    no 29
## 12  1961  12    no 40
## 13  1961  13    no 28
## 14  1961  14    no 17
## 15  1961  15    no 15
## 16  1961  16    no 21
## 17  1961  17    no 24
## 18  1961  18    no 15
## 19  1961  19    no 32
## 20  1961  20    no 22
## 21  1961  21    no 24
## 22  1961  22    no 11
## 23  1961  23    no 27
## 24  1961  24   yes 12
## 25  1961  25   yes 41
## 26  1961  26   yes 15
## 27  1961  27   yes 18
## 28  1961  28   yes 11
## 29  1961  29   yes 19
## 30  1961  30   yes 19
## 31  1961  31   yes  9
## 32  1961  32   yes 21
## 33  1961  33   yes 22
## 34  1961  34   yes 23
## 35  1961  35   yes 14
## 36  1961  36   yes 19
## 37  1961  37   yes 15
## 38  1961  38   yes 13
## 39  1961  39   yes 22
## 40  1961  40   yes 42
## 41  1961  41   yes 29
## 42  1961  42   yes 21
## 43  1961  43   yes 12
## 44  1961  44   yes 16
## 45  1961  45    no 17
## 46  1961  46    no 27
## 47  1961  47    no 37
## 48  1961  48    no 32
## 49  1961  49    no 25
## 50  1961  50    no 20
## 51  1961  51    no 40
## 52  1961  52    no 21
## 53  1961  53    no 18
## 54  1961  54    no 35
## 55  1961  55    no 21
## 56  1961  56    no 25
## 57  1961  57    no 34
## 58  1961  58    no 42
## 59  1961  59    no 27
## 60  1961  60    no 34
## 61  1961  61    no 47
## 62  1961  62    no 36
## 63  1961  63    no 15
## 64  1961  64    no 26
## 65  1961  65    no 27
## 66  1961  66    no 18
## 67  1961  67    no 16
## 68  1961  68    no 32
## 69  1961  69    no 28
## 70  1961  70    no 17
## 71  1961  71    no 16
## 72  1961  72    no 19
## 73  1961  73    no 18
## 74  1961  74    no 22
## 75  1961  75    no 37
## 76  1961  76    no 29
## 77  1961  77    no 18
## 78  1961  78    no 14
## 79  1961  79    no 14
## 80  1961  80    no 18
## 81  1961  81    no 21
## 82  1961  82    no 39
## 83  1961  83    no 39
## 84  1961  84    no 21
## 85  1961  85    no 15
## 86  1961  86    no 17
## 87  1961  87    no 20
## 88  1961  88    no 24
## 89  1961  89    no 30
## 90  1961  90    no 25
## 91  1961  91    no  8
## 92  1961  92    no 21
## 93  1962   1    no  9
## 94  1962   2    no 20
## 95  1962   3    no 15
## 96  1962   4    no 14
## 97  1962   5    no 30
## 98  1962   6    no 23
## 99  1962   7    no 15
## 100 1962   8    no 14
## 101 1962   9    no 16
## 102 1962  10    no 20
## 103 1962  11   yes 17
## 104 1962  12   yes 23
## 105 1962  13   yes 16
## 106 1962  14   yes 20
## 107 1962  15   yes 13
## 108 1962  16   yes 13
## 109 1962  17   yes  9
## 110 1962  18   yes 10
## 111 1962  19   yes 27
## 112 1962  20   yes 12
## 113 1962  21   yes  7
## 114 1962  22   yes 11
## 115 1962  23   yes 15
## 116 1962  24   yes 19
## 117 1962  25   yes 32
## 118 1962  26   yes 22
## 119 1962  27   yes 24
## 120 1962  28   yes  9
## 121 1962  29    no 10
## 122 1962  30    no 14
## 123 1962  31    no 18
## 124 1962  32    no 26
## 125 1962  33    no 38
## 126 1962  34    no 31
## 127 1962  35    no 12
## 128 1962  36    no  8
## 129 1962  37    no 22
## 130 1962  38    no 17
## 131 1962  39    no 31
## 132 1962  40    no 49
## 133 1962  41    no 23
## 134 1962  42    no 14
## 135 1962  43    no 25
## 136 1962  44    no 24
## 137 1962  45    no 18
## 138 1962  46    no 19
## 139 1962  47   yes 29
## 140 1962  48   yes 17
## 141 1962  49   yes 17
## 142 1962  50   yes 15
## 143 1962  51   yes 25
## 144 1962  52   yes  9
## 145 1962  53   yes 16
## 146 1962  54   yes 25
## 147 1962  55   yes 25
## 148 1962  56   yes 16
## 149 1962  57   yes 22
## 150 1962  58   yes 21
## 151 1962  59   yes 17
## 152 1962  60   yes 26
## 153 1962  61   yes 41
## 154 1962  62   yes 25
## 155 1962  63   yes 12
## 156 1962  64   yes 17
## 157 1962  65    no 21
## 158 1962  66    no 19
## 159 1962  67    no 24
## 160 1962  68    no 44
## 161 1962  69    no 31
## 162 1962  70    no 21
## 163 1962  71    no 20
## 164 1962  72    no 19
## 165 1962  73    no 20
## 166 1962  74    no 29
## 167 1962  75    no 48
## 168 1962  76    no 36
## 169 1962  77    no 15
## 170 1962  78    no 16
## 171 1962  79    no 29
## 172 1962  80    no 12
## 173 1962  81   yes 24
## 174 1962  82   yes 26
## 175 1962  83   yes 16
## 176 1962  84   yes 15
## 177 1962  85   yes 12
## 178 1962  86   yes 22
## 179 1962  87   yes 24
## 180 1962  88   yes 16
## 181 1962  89   yes 25
## 182 1962  90   yes 14
## 183 1962  91   yes 15
## 184 1962  92   yes  9

Aquí encontramos las variables year(año), day(día), limit(límite) y la variable y.

La variable year(año) indica el año en que se aplicó el experimento. La variable day(día) indica los días en que se aplicó el experimento. La variable limit(límite) indica para qué días se aplicó la norma de límite de velocidad. La variable y indica el número de accidentes ocurridos en autopista.

El experimento se realizó en Suecia en 1961 y 1962 para evaluar el efecto de un límite de velocidad en la tasa de accidentes de autopista. Se realizó en 92 días en cada año, coincidiendo para que el día j en 1962 fuera comparable al día j en 1961. En algunos días, el límite de velocidad estaba en vigor y se aplicaba, mientras que en otros días no había límite de velocidad y los automóviles tendían a ser conducido más rápido. Los días con límite de velocidad tendían a estar en bloques contiguos.


Pregunta 9. Verifique si en estos momentos tiene objetos en memoria, ¿qué comando debe utilizar? Si tiene objetos, bórrelos en su totalidad.

ls()
## character(0)
rm(list = ls())

ls()
## character(0)

Pregunta 10. Supongamos que a usted se le ha olvidado incluir los tres últimos valores (15,19,23) de un vector de longitud 2500. Usted está angustiado pues escribir nuevamente ese vector tardaría muchísimo. Para evitar este contratiempo, le pregunta al profesor qué podría hacer y el profesor le contesta, realmente es muy fácil use la función append.

  1. Describa brevemente la función append.

R/. La función append se encarga de agregar elementos a un vector.

  1. Asumiendo que los 7 últimos valores que usted digitó del vector XYZ son 18,19,27,10,20,15,25. Inserte los últimos tres valores que le hizo falta.

XYZ=c(18,19,27,10,20,15,25)
XYZ
## [1] 18 19 27 10 20 15 25
append(XYZ,seq(15,23,by=4),after = length(XYZ))
##  [1] 18 19 27 10 20 15 25 15 19 23
  1. Sin usar la función sugerida por el profesor, ¿cómo podría hacer esto de otra manera?

R/. Sin usar la función append se podrían insertar los valores faltantes concatenando dos vectores, el vector XYZ con el vector formado por los valores faltantes. De la siguiente manera.


XYZ=c(18,19,27,10,20,15,25)
 
XYZ
## [1] 18 19 27 10 20 15 25
V=seq(15,23,by=4)

V
## [1] 15 19 23

x=c(XYZ,V)

x
##  [1] 18 19 27 10 20 15 25 15 19 23