Guion RStudio

1. ¿Qué son R y RStudio?

R es un lenguaje y entorno de programación para el análisis estadı́stico y gráfico. Se
trata de un proyecto de software libre que, en la actualidad, se está utilizando cada vez
más en diferentes campos de investigación. R se distribuye bajo la licencia GNU GPL y
está disponible para los sistemas operativos Windows, Macintosh, Unix y GNU/Linux.
La página principal del proyecto bajo el que se desarrolla R es [Link]
org/. La universidad de Zaragoza a través de su Oficina de software libre también pro-
porciona información y ayuda sobre R.
Si utilizas R en alguno de tus trabajos debes citarlo e incluir en la bibliografı́a su
referencia:
R Core Team (2014). R: A language and environment for statistical computing. R
Foundation for Statistical Computing, Vienna, Austria. URL http: // www. R-project.
org/ .
RStudio es un entorno de desarrollo integrado (IDE) para el lenguaje de programación
R, ampliamente utilizado en estadı́stica, análisis de datos y ciencia de datos. Incluye una
consola, editor de sintaxis que apoya la ejecución del código, ası́ como herramientas que
facilitan la escritura, ejecución y depuración del código en R.
Para instalar R y RStudio, debéis seguir los siguientes pasos. Es recomendable que
si tenéis ya una versión previa instalada de R, desinstalarla antes de instalar
esta.

1. Primero debéis instalar R. Lo podéis descargar gratuitamente desde la siguiente

página: [Link] Os aparecerá la siguiente página,

1
 Prof. Begoña Álvarez Tena

y debéis pinchar en download R, que abrirá la siguiente página que os pide que
seleccionéis un CRAN mirror. Seleccionáis el de España.

Una vez pinchéis el CRAN mirror correspondiente a España, en la nueva pestaña

que se os abrirá descargáis para el sistema operativo que utilicéis.

En el caso de Download R por Windows, os llevará a la siguiente pestaña, donde

deberéis pinchar en base.

Y por último descargáis R. Se descargará un ejecutable que deberéis ejecutar para

completar la instalación.

2. El segundo paso es instalar RStudio. Lo podéis descargar gratuitamente desde la

siguiente página: [Link] Os aparecerá la siguiente página:

Guion RStudio 2
 Prof. Begoña Álvarez Tena

Debéis pinchar en Download RStudio. Se abrirá otra página, donde nuevamente

debéis pinchar en Download RStudio

En la página que aparece, se descarga RStudio desktop for windows.

Se descargará un ejecutable de RStudio que debéis ejecutar para completar la ins-

talación.

Guion RStudio 3
 Prof. Begoña Álvarez Tena

2. Algunos conceptos básicos del lenguaje de progra-

mación R
La orden help, ası́ como la ventana que aparece a la derecha, proporciona la ayuda
disponible sobre las diferentes ordenes y funciones que tiene R.

Notar que R distingue entre mayúsculas y minúsculas, de manera que m1 y M 1 son

objetos diferentes.
Para escribir comentarios debemos utilizar el sı́mbolo #, y todo lo que sigue en esa
misma lı́nea es ignorado.
El código de R para datos no observados o perdidos es N A.
R utiliza el sı́mbolo < − para realizar una asignación, es decir para que un objeto
reciba un nombre.
Vectores
Podemos crear vectores a través de la función c(), y dichos vectores pueden ser:
numéricos reales, lógicos, de tipo carácter o enteros.

Para crear vectores son muy útiles los operadores rep y seq. Veamos como podemos
utilizarlo:

Se ha creado un vector del 1 al 10, al que hemos llamado v1.

Guion RStudio 4
 Prof. Begoña Álvarez Tena

Se ha creado un vector que contiene al 1 cinco veces, seguido del 1 dos veces, el 2 dos
veces y el 3 dos veces, y por último la sucesión de números decreciente del 8 al 3. Este
vector recibe el nombre de v2.

Se ha creado un vector que contiene una secuencia de números que empieza en 0.0 y acaba
en 1.0, realizando 11 particiones de igual tamaño. Este vector recibe el nombre de v3.

Se crea una secuencia de números entre el 1 y el 6, con una separación entre cada uno de
longitud tres. Este vector recibe el nombre de v4.

Se crea un vector que contiene los números entre el 0 y el 10 de 10 en 10. Este vector
recibe el nombre de v5.

Matrices
Las matrices tienen dos dimensiones y se pueden obtener como la unión de dos vectores
del mismo tipo y la isma longitud, utilizando las funciones cbind o rbind, que unen
vectores columna y fila, respectivamente.

Con la orden cbind se indica que la matriz la vamos a construir por columnas. Después
cada vector representa cada una de las columnas de la matriz.

Guion RStudio 5
 Prof. Begoña Álvarez Tena

Con la orden rbind se indica que la matriz la vamos a construir por filas. Después cada
vector representa cada una de las filas de la matriz.

Para acceder a los diferentes elementos de una matriz, se utiliza la sintaxis: [, ], donde
el primer elemento corresponde a la fila y el segundo a la columna.

Indica que el elemento situado en la fila 2 columna 1 es tres.

Indica la segunda fila de la matriz.

Indica la segunda columna de la matriz.

Guion RStudio 6
 Prof. Begoña Álvarez Tena

3. Representación de la región factible de un PPL de

2 variables
Para realizar la representación gráfica de la región factible es importante conocer la
función plot.
La función plot de R permite representar funciones gráficamente. Algunos de los
parámetros que podemos incluir pueden ser los siguientes:

El parámetro xlim, permite determinar que rango de valores queremos representar

en el eje X.

El parámetro ylim, permite determinar que rango de valores queremos representar

en el eje Y.

El parámetro xlab, permite dar un tı́tulo al eje X.

El parámetro ylab, permite dar un tı́tulo al eje Y.

El parámetro main, permite dar un tı́tulo al gráfico.

El parámetro sub, permite dar un subtitulo al gráfico.

El parámetro col, permite dar un color a la función que vamos a representar.

El parámetro type, permite determinar el tipo de gráfico.

ˆ p: puntos. ˆ h: lı́neas verticales que van desde los

ˆ l : lı́neas. puntos hasta el eje X.
ˆ b: lı́neas y puntos. ˆ s: paso de izquierda (escalones).
ˆ c: lı́neas solo en las posiciones de los ˆ S : paso de derecha (escalones).
puntos (no dibuja los puntos). ˆ n: no traza nada, resulta útil si solo
ˆ o: lı́neas y puntos superpuestos, pe- se quiere establecer el marco de los
ro las lı́neas pasa por los puntos. ejes.

El parámetro pch, permite determinar el sı́mbolo de los puntos:

ˆ 0 : cuadrado vacı́o. ˆ 6 : triángulo vacı́o ˆ 11 : cı́rculo con cruz

(abajo). en el centro.
ˆ 1 : cı́rculo vacı́o.
ˆ 7 : cuadrado cruzado. ˆ 12 : cuadrado lleno.
ˆ 2 : triángulo vacı́o
(arriba). ˆ 13 : cı́rculo lleno.
ˆ 8 : asterisco.
ˆ 3 : cruz. ˆ 14 : triángulo lleno
ˆ 9 : rombos cruzados. (arriba).
ˆ 4 : X.
ˆ 10 : triángulo vacı́o ˆ 15 : triángulo lleno
ˆ 5 : diamante vacı́o. (hacia la derecha) (abajo).

Guion RStudio 7
 Prof. Begoña Álvarez Tena

El parámetro lty, permite definir el estilo de la lı́nea.

ˆ 1 : lı́nea sólida (por defecto). ˆ 4 : lı́nea con guión y punto.

ˆ 2 : lı́nea discontinua (guiones). ˆ 5 : lı́nea con guión largo y punto.
ˆ 3 : lı́nea punteada. ˆ 6 : lı́nea con guión largo y 2 puntos.

El parámetro lwd, permite determinar el grosor de las lı́neas que se representan.

Los valores que puede tomar van desde 0, que representa la lı́nea más fina, hasta
el valor numérico que se desee. Lo más común es dejar este valor como lwd = 2,
aunque por defecto toma el valor lwd = 1.

Puede darse el caso, de que queramos representar varias funciones en el mismo gráfi-
co, para ello utilizamos las funciones curve y abline.

abline se utilizará cuando la función a representar sea una constante.

curve se utilizará cuando la función a representar no sea una constante. Es impor-
tante que en esta función se utilice el parámetro add = TRUE, ya que sino, no se
representará en el mismo gráfico que se hubiese creado anteriormente, sino en uno
nuevo.

La función polygon, permite colorear una determinada región, delimitandola por sus
vértices. A través de dos vectores se definirán las coordenadas de los vértices. A través
del parámetro col se determina el color para colorear la región.
La función solve permite resolver sistemas de ecuaciones lineales. El primer parámetro
debe ser la matriz A y el segundo la matriz b, siendo el sistema de ecuaciones a resolver
Ax = b.
Ejemplo 1. Dado el siguiente problema de programación lineal:


 max 0.3x1 + 0.1x2

s.a x2 ≥ x1





 x1 ≤ 4000


 2x1 − x2 ≥ 1000
2x1 + x2 ≤ 14000





x1 ≥ 0, x2 ≥ 0


Representa la región factible, curvas de nivel e identifica la solución óptima.

Para representar su región factible debemos proceder de la siguiente manera.

1. Definimos cada una de las restricciones como una función en R, ya que vamos a
representar la igualdad de cada restricción.

Guion RStudio 8
 Prof. Begoña Álvarez Tena

2. Creamos un plot. Podemos optar por dos opciones, crear un gráfico que no dibuje
nada y representar a través de las funciones abline y curve todas las restricciones. O
bien crear un plot representando una de las restricciones, y las restantes añadirlas
en el mismo gráfico utilizando las funciones abline y curve.
1ª opción:

2ª opción:

En ambos casos obtenemos el siguiente gráfico:

Guion RStudio 9
 Prof. Begoña Álvarez Tena

3. Determinamos el polı́gono de la región factible, y lo pintamos de color amarillo claro.

Para hallar los vértices de la región factible utilizamos la función solve:

Los vértices (puntos) del polı́gono que forma la región factible son los siguientes:

(1000, 1000) (4000, 4000) (4000, 6000) (3750, 6500)

Observamos que en el primer vector, se indica las coordenadas x de cada uno de los
puntos, y en el segundo vector, las coordenadas y de cada uno de los puntos.

Notar que dentro del parámetro col, el parámetro alpha determina la transparen-
cia del color y puede tomar valores desde 0 (completamente transparente) hasta 1
(completamente opaco, sin transparencia).
El gráfico que obtenemos es el siguiente:

Para determinar que vértice de la región factible es la solución óptima de nuestro

problema de optimización, debemos dibujar las correspondientes curvas de nivel.
En nuestro caso, la función objetivo es 0.3x1 + 0.1x2 , por lo tanto las curvas de nivel
se calculan de la siguiente manera:
k − 0.3x1
0.3x1 + 0.1x2 = k ⇔ x2 =
0.1

Guion RStudio 10
 Prof. Begoña Álvarez Tena

k − 0.3x1
Debemos representar la función para cada valor de k. La manera que
0.1
tenemos de hacerlo en R, es iterando en un bucle como se muestra a continuación:

El fragmento de código hace lo siguiente: para cada valor de k en el vector que

ha sido construido como una secuencia que empieza en 0, acaba en 1800 y toma
valores de 100 en 100, es decir: (0, 100, 200, 300, . . . , 1500, 1600, 1700, 1800), con la
función curve (vista anteriormente) se dibuja la correspondiente recta en el gráfico
que hemos creado (por eso se añade el parámetro add = TRUE ). Obtenemos el
siguiente gráfico:

donde se observa que el punto (4000, 6000), es la solución óptima a nuestro problema
de programación lineal.
Si lo deseamos se puede añadir una leyenda:

donde topleft indica la posición en el gráfico en el que se va a incluir la leyenda.

Guion RStudio 11
 Prof. Begoña Álvarez Tena

topleft: arriba a la izquierda bottomright: abajo a la derecha

topright: arriba a la derecha left: centrado a la izquierda
bottomleft: abajo a la izquierda right: centrado a la derecha

Por último, si lo deseamos también podemos marcar el punto que ha resultado ser la
solución óptima del problema de programación lineal, a través de la función points.

Y el gráfico quedarı́a de la siguiente manera:

Guion RStudio 12