Instalar paquetes y llamar librerias

# instal.packages("tidyverse") # Paquete para manipulacion de datos
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.4     ✔ tidyr     1.3.1
## ✔ purrr     1.0.4     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

Asignar una variable

x <- 3
y <- 2

Imprimir el resultado

x
## [1] 3
y
## [1] 2

Operaciones aritmetias

suma <- x + y
suma
## [1] 5
resta <- x - y
resta
## [1] 1
multiplicacion <- x * y
multiplicacion
## [1] 6
division <- x/y
division
## [1] 1.5
division_entera <- x%/% y 
division_entera
## [1] 1
residuo <- x %% y
residuo
## [1] 1
potencia <- x ** 2
potencia
## [1] 9
potencia_2 <- x ^ 2
potencia_2 
## [1] 9

Funciones

raiz_cuadrada <- sqrt(x)
raiz_cuadrada
## [1] 1.732051
raiz_cubica <- x ** (1/3)
raiz_cubica
## [1] 1.44225
# ?sqrt # El signo de interrogacion es para desplagar ayuda 

exponencial <- exp(1)
exponencial
## [1] 2.718282
absoluto <- abs(x)
absoluto 
## [1] 3
signo <- sign(x)
signo 
## [1] 1
redondeo_arriba <- ceiling (x/y)
redondeo_arriba
## [1] 2
redondeo_abajo <- floor(x/y)
redondeo_abajo
## [1] 1
truncar <- trunc(x/y)
truncar
## [1] 1

Constantes

pi
## [1] 3.141593
radio <- 5
area_circulo <- pi*radio^2
area_circulo
## [1] 78.53982

Vectores

a <- c(1,2,3,4,5) # Secuencia de enteros 1:5
a
## [1] 1 2 3 4 5
b <- c(1,10, by = 0.5)
b
##             by 
##  1.0 10.0  0.5
nombres <- c("Juan","Ana","Pedro","Carla","Sara")
nombres
## [1] "Juan"  "Ana"   "Pedro" "Carla" "Sara"
calificaciones <- c(100,90,50,100,65)
calificaciones
## [1] 100  90  50 100  65
longitud <- length(a)
longitud
## [1] 5
promedio <- mean(calificaciones)
promedio
## [1] 81
orden_ascendente <- sort(calificaciones)
orden_ascendente
## [1]  50  65  90 100 100
orden_decendente <- sort(calificaciones, decreasing = TRUE)
orden_decendente
## [1] 100 100  90  65  50

Tablas

tabla_de_calificaciones <- data.frame(nombres, calificaciones)

# Crear la columna 'estatus' evaluando la columna 'calificaciones'
tabla_de_calificaciones$estatus <- ifelse(tabla_de_calificaciones$calificaciones >= 70, "Aprobado", "Reprobado")

# Condicionales: Igual == , Desigual != , Mayor que >, Menor que<, 
# Mayor o igual que >=, Menor o igual que <=

summary(tabla_de_calificaciones)
##    nombres          calificaciones   estatus         
##  Length:5           Min.   : 50    Length:5          
##  Class :character   1st Qu.: 65    Class :character  
##  Mode  :character   Median : 90    Mode  :character  
##                     Mean   : 81                      
##                     3rd Qu.:100                      
##                     Max.   :100
str(tabla_de_calificaciones)
## 'data.frame':    5 obs. of  3 variables:
##  $ nombres       : chr  "Juan" "Ana" "Pedro" "Carla" ...
##  $ calificaciones: num  100 90 50 100 65
##  $ estatus       : chr  "Aprobado" "Aprobado" "Reprobado" "Aprobado" ...
head(tabla_de_calificaciones) # 6 por default
##   nombres calificaciones   estatus
## 1    Juan            100  Aprobado
## 2     Ana             90  Aprobado
## 3   Pedro             50 Reprobado
## 4   Carla            100  Aprobado
## 5    Sara             65 Reprobado
tail(tabla_de_calificaciones)
##   nombres calificaciones   estatus
## 1    Juan            100  Aprobado
## 2     Ana             90  Aprobado
## 3   Pedro             50 Reprobado
## 4   Carla            100  Aprobado
## 5    Sara             65 Reprobado

Instalar paquetes y llamar librerias

#instal.packages("tidyverse")
library(tidyverse)

SELECT para seleccionar columnas

resultados <- select(tabla_de_calificaciones, -c(nombres))

FILTER para filtrar renglones

resultados <- filter(resultados, estatus=="Aprobado")

Actividad 1. Calculadora del IMC

alumno <- c("Raul","Renata","Juan","Sergio","Meme")
peso <- c(78,52,78,90,560)
altura <- c(1.80,1.55,1.75,1.86,1.80)

tabla_salud <- data.frame(alumno, peso, altura)
tabla_salud$IMC <- tabla_salud$peso / tabla_salud$altura^2
tabla_salud$resultados <- ifelse(tabla_salud$IMC < 18.5, "Bajo Peso", ifelse(tabla_salud$IMC <=24.9, "Peso Normal", ifelse(tabla_salud$IMC <=29.9, "Sobrepeso", "Obesidad")))

Grafica

# Gráfica
Semana <- 1:10
Ventas <- c(150,160,165,180,175,190,200,195,210,220)

# Crear gráfico de línea con puntos
plot(Semana, Ventas, type = "b", 
     main = "Ventas Semanales (K USD)",
     xlab = "Semana",
     ylab = "Ventas (K USD)")

# Crear data frame
datos_ventas <- data.frame(Semana, Ventas)

# Modelo de regresión lineal
regresion <- lm(Ventas ~ Semana, data = datos_ventas)
summary(regresion)
## 
## Call:
## lm(formula = Ventas ~ Semana, data = datos_ventas)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -7.758 -1.538  0.500  2.439  6.455 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 144.3333     3.1714   45.51 6.00e-11 ***
## Semana        7.3030     0.5111   14.29 5.61e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.643 on 8 degrees of freedom
## Multiple R-squared:  0.9623, Adjusted R-squared:  0.9576 
## F-statistic: 204.2 on 1 and 8 DF,  p-value: 5.614e-07
# Predicción para semanas 11 a 20
datos_nuevos <- data.frame(Semana = 11:20)
prediccion <- predict(regresion, datos_nuevos)
prediccion
##        1        2        3        4        5        6        7        8 
## 224.6667 231.9697 239.2727 246.5758 253.8788 261.1818 268.4848 275.7879 
##        9       10 
## 283.0909 290.3939

Concluciones

R es un lenguaje de programaion util para realizar calculos, principalmente estadisticos, y forma parte de las herramientas del Big Data.

RStudio es e; entorno donde se puede programar R, y gracias a que tambien aqui se puede programar Python, hace unos años se anuncio que su nombre sera Posit.

En esta introduccion lo que llama la atencion es la constante aparicion de alertas o errores en el programa, los cuales encontramos que principalmente sedeben a que la version no es la mas reciente, a la falta de instalacion de paquetes o llamar a las librerias, prolemas de escritura (typos) y los muchos argumentos que tienen las funciones.

Si desde un inico programamos de formas estructurada, disciplinada y meticulosa, podremos preveer muchas de las alertas o errores y asi obtenemos los muchos beneficios

LS0tDQp0aXRsZTogIkNvbWFuZG9zIEJhc2ljb3MiDQphdXRob3I6ICJFbGlzZW8gTWFydGluZXogQTAxMTc4MTQwIg0KZGF0ZTogIjIwMjUtMDMtMTgiDQpvdXRwdXQ6DQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiBUUlVFDQogICAgdG9jX2Zsb2F0OiBUUlVFDQogICAgY29kZV9kb3dubG9hZDogVFJVRQ0KICAgIHRoZW1lOiAiY29zbW8iDQotLS0gIA0KDQohW10oaHR0cHM6Ly9naWZzcy5jb20vcHJvZmVzaW9uZXMvZG9jdG9yZXMvZG9jdG9yYS0wMy5naWYpDQoNCiMgSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyaWFzDQpgYGB7cn0NCiMgaW5zdGFsLnBhY2thZ2VzKCJ0aWR5dmVyc2UiKSAjIFBhcXVldGUgcGFyYSBtYW5pcHVsYWNpb24gZGUgZGF0b3MNCmxpYnJhcnkodGlkeXZlcnNlKQ0KYGBgDQoNCiMgQXNpZ25hciB1bmEgdmFyaWFibGUNCmBgYHtyfQ0KeCA8LSAzDQp5IDwtIDINCg0KYGBgDQoNCiMgSW1wcmltaXIgZWwgcmVzdWx0YWRvDQpgYGB7cn0NCngNCnkNCmBgYA0KDQojIE9wZXJhY2lvbmVzIGFyaXRtZXRpYXMNCmBgYHtyfQ0Kc3VtYSA8LSB4ICsgeQ0Kc3VtYQ0KDQpyZXN0YSA8LSB4IC0geQ0KcmVzdGENCg0KbXVsdGlwbGljYWNpb24gPC0geCAqIHkNCm11bHRpcGxpY2FjaW9uDQoNCmRpdmlzaW9uIDwtIHgveQ0KZGl2aXNpb24NCg0KZGl2aXNpb25fZW50ZXJhIDwtIHglLyUgeSANCmRpdmlzaW9uX2VudGVyYQ0KDQpyZXNpZHVvIDwtIHggJSUgeQ0KcmVzaWR1bw0KDQpwb3RlbmNpYSA8LSB4ICoqIDINCnBvdGVuY2lhDQoNCnBvdGVuY2lhXzIgPC0geCBeIDINCnBvdGVuY2lhXzIgDQpgYGANCg0KIyBGdW5jaW9uZXMNCmBgYHtyfQ0KcmFpel9jdWFkcmFkYSA8LSBzcXJ0KHgpDQpyYWl6X2N1YWRyYWRhDQoNCnJhaXpfY3ViaWNhIDwtIHggKiogKDEvMykNCnJhaXpfY3ViaWNhDQoNCiMgP3NxcnQgIyBFbCBzaWdubyBkZSBpbnRlcnJvZ2FjaW9uIGVzIHBhcmEgZGVzcGxhZ2FyIGF5dWRhIA0KDQpleHBvbmVuY2lhbCA8LSBleHAoMSkNCmV4cG9uZW5jaWFsDQoNCmFic29sdXRvIDwtIGFicyh4KQ0KYWJzb2x1dG8gDQoNCnNpZ25vIDwtIHNpZ24oeCkNCnNpZ25vIA0KDQpyZWRvbmRlb19hcnJpYmEgPC0gY2VpbGluZyAoeC95KQ0KcmVkb25kZW9fYXJyaWJhDQoNCnJlZG9uZGVvX2FiYWpvIDwtIGZsb29yKHgveSkNCnJlZG9uZGVvX2FiYWpvDQoNCnRydW5jYXIgPC0gdHJ1bmMoeC95KQ0KdHJ1bmNhcg0KYGBgDQoNCiMgQ29uc3RhbnRlcw0KYGBge3J9DQpwaQ0KcmFkaW8gPC0gNQ0KYXJlYV9jaXJjdWxvIDwtIHBpKnJhZGlvXjINCmFyZWFfY2lyY3Vsbw0KYGBgDQoNCiMgVmVjdG9yZXMNCmBgYHtyfQ0KYSA8LSBjKDEsMiwzLDQsNSkgIyBTZWN1ZW5jaWEgZGUgZW50ZXJvcyAxOjUNCmENCg0KYiA8LSBjKDEsMTAsIGJ5ID0gMC41KQ0KYg0KDQpub21icmVzIDwtIGMoIkp1YW4iLCJBbmEiLCJQZWRybyIsIkNhcmxhIiwiU2FyYSIpDQpub21icmVzDQoNCmNhbGlmaWNhY2lvbmVzIDwtIGMoMTAwLDkwLDUwLDEwMCw2NSkNCmNhbGlmaWNhY2lvbmVzDQoNCmxvbmdpdHVkIDwtIGxlbmd0aChhKQ0KbG9uZ2l0dWQNCg0KcHJvbWVkaW8gPC0gbWVhbihjYWxpZmljYWNpb25lcykNCnByb21lZGlvDQoNCm9yZGVuX2FzY2VuZGVudGUgPC0gc29ydChjYWxpZmljYWNpb25lcykNCm9yZGVuX2FzY2VuZGVudGUNCg0Kb3JkZW5fZGVjZW5kZW50ZSA8LSBzb3J0KGNhbGlmaWNhY2lvbmVzLCBkZWNyZWFzaW5nID0gVFJVRSkNCm9yZGVuX2RlY2VuZGVudGUNCmBgYA0KDQojIFRhYmxhcw0KYGBge3J9DQp0YWJsYV9kZV9jYWxpZmljYWNpb25lcyA8LSBkYXRhLmZyYW1lKG5vbWJyZXMsIGNhbGlmaWNhY2lvbmVzKQ0KDQojIENyZWFyIGxhIGNvbHVtbmEgJ2VzdGF0dXMnIGV2YWx1YW5kbyBsYSBjb2x1bW5hICdjYWxpZmljYWNpb25lcycNCnRhYmxhX2RlX2NhbGlmaWNhY2lvbmVzJGVzdGF0dXMgPC0gaWZlbHNlKHRhYmxhX2RlX2NhbGlmaWNhY2lvbmVzJGNhbGlmaWNhY2lvbmVzID49IDcwLCAiQXByb2JhZG8iLCAiUmVwcm9iYWRvIikNCg0KIyBDb25kaWNpb25hbGVzOiBJZ3VhbCA9PSAsIERlc2lndWFsICE9ICwgTWF5b3IgcXVlID4sIE1lbm9yIHF1ZTwsIA0KIyBNYXlvciBvIGlndWFsIHF1ZSA+PSwgTWVub3IgbyBpZ3VhbCBxdWUgPD0NCg0Kc3VtbWFyeSh0YWJsYV9kZV9jYWxpZmljYWNpb25lcykNCnN0cih0YWJsYV9kZV9jYWxpZmljYWNpb25lcykNCmhlYWQodGFibGFfZGVfY2FsaWZpY2FjaW9uZXMpICMgNiBwb3IgZGVmYXVsdA0KdGFpbCh0YWJsYV9kZV9jYWxpZmljYWNpb25lcykNCmBgYA0KDQoNCiMgSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyaWFzDQpgYGB7cn0NCiNpbnN0YWwucGFja2FnZXMoInRpZHl2ZXJzZSIpDQpsaWJyYXJ5KHRpZHl2ZXJzZSkNCmBgYA0KDQojIFNFTEVDVCBwYXJhIHNlbGVjY2lvbmFyIGNvbHVtbmFzDQpgYGB7cn0NCnJlc3VsdGFkb3MgPC0gc2VsZWN0KHRhYmxhX2RlX2NhbGlmaWNhY2lvbmVzLCAtYyhub21icmVzKSkNCmBgYA0KDQojIEZJTFRFUiBwYXJhIGZpbHRyYXIgcmVuZ2xvbmVzDQpgYGB7cn0NCnJlc3VsdGFkb3MgPC0gZmlsdGVyKHJlc3VsdGFkb3MsIGVzdGF0dXM9PSJBcHJvYmFkbyIpDQpgYGANCg0KIyBBY3RpdmlkYWQgMS4gQ2FsY3VsYWRvcmEgZGVsIElNQw0KYGBge3J9DQphbHVtbm8gPC0gYygiUmF1bCIsIlJlbmF0YSIsIkp1YW4iLCJTZXJnaW8iLCJNZW1lIikNCnBlc28gPC0gYyg3OCw1Miw3OCw5MCw1NjApDQphbHR1cmEgPC0gYygxLjgwLDEuNTUsMS43NSwxLjg2LDEuODApDQoNCnRhYmxhX3NhbHVkIDwtIGRhdGEuZnJhbWUoYWx1bW5vLCBwZXNvLCBhbHR1cmEpDQp0YWJsYV9zYWx1ZCRJTUMgPC0gdGFibGFfc2FsdWQkcGVzbyAvIHRhYmxhX3NhbHVkJGFsdHVyYV4yDQp0YWJsYV9zYWx1ZCRyZXN1bHRhZG9zIDwtIGlmZWxzZSh0YWJsYV9zYWx1ZCRJTUMgPCAxOC41LCAiQmFqbyBQZXNvIiwgaWZlbHNlKHRhYmxhX3NhbHVkJElNQyA8PTI0LjksICJQZXNvIE5vcm1hbCIsIGlmZWxzZSh0YWJsYV9zYWx1ZCRJTUMgPD0yOS45LCAiU29icmVwZXNvIiwgIk9iZXNpZGFkIikpKQ0KYGBgDQoNCiMgR3JhZmljYQ0KYGBge3J9DQojIEdyw6FmaWNhDQpTZW1hbmEgPC0gMToxMA0KVmVudGFzIDwtIGMoMTUwLDE2MCwxNjUsMTgwLDE3NSwxOTAsMjAwLDE5NSwyMTAsMjIwKQ0KDQojIENyZWFyIGdyw6FmaWNvIGRlIGzDrW5lYSBjb24gcHVudG9zDQpwbG90KFNlbWFuYSwgVmVudGFzLCB0eXBlID0gImIiLCANCiAgICAgbWFpbiA9ICJWZW50YXMgU2VtYW5hbGVzIChLIFVTRCkiLA0KICAgICB4bGFiID0gIlNlbWFuYSIsDQogICAgIHlsYWIgPSAiVmVudGFzIChLIFVTRCkiKQ0KDQojIENyZWFyIGRhdGEgZnJhbWUNCmRhdG9zX3ZlbnRhcyA8LSBkYXRhLmZyYW1lKFNlbWFuYSwgVmVudGFzKQ0KDQojIE1vZGVsbyBkZSByZWdyZXNpw7NuIGxpbmVhbA0KcmVncmVzaW9uIDwtIGxtKFZlbnRhcyB+IFNlbWFuYSwgZGF0YSA9IGRhdG9zX3ZlbnRhcykNCnN1bW1hcnkocmVncmVzaW9uKQ0KDQojIFByZWRpY2Npw7NuIHBhcmEgc2VtYW5hcyAxMSBhIDIwDQpkYXRvc19udWV2b3MgPC0gZGF0YS5mcmFtZShTZW1hbmEgPSAxMToyMCkNCnByZWRpY2Npb24gPC0gcHJlZGljdChyZWdyZXNpb24sIGRhdG9zX251ZXZvcykNCnByZWRpY2Npb24NCg0KYGBgDQoNCg0KIyBDb25jbHVjaW9uZXMNCioqUioqIGVzIHVuIGxlbmd1YWplIGRlIHByb2dyYW1haW9uIHV0aWwgcGFyYSByZWFsaXphciBjYWxjdWxvcywgcHJpbmNpcGFsbWVudGUgZXN0YWRpc3RpY29zLCB5IGZvcm1hIHBhcnRlIGRlIGxhcyBoZXJyYW1pZW50YXMgZGVsICpCaWcgRGF0YSouDQoNCipSU3R1ZGlvKiBlcyBlOyBlbnRvcm5vIGRvbmRlIHNlIHB1ZWRlIHByb2dyYW1hciAqUiosIHkgZ3JhY2lhcyBhIHF1ZSB0YW1iaWVuIGFxdWkgc2UgcHVlZGUgcHJvZ3JhbWFyICpQeXRob24qLCBoYWNlIHVub3MgYcOxb3Mgc2UgYW51bmNpbyBxdWUgc3Ugbm9tYnJlIHNlcmEgKlBvc2l0Ki4NCg0KRW4gZXN0YSBpbnRyb2R1Y2Npb24gbG8gcXVlIGxsYW1hIGxhIGF0ZW5jaW9uIGVzIGxhIGNvbnN0YW50ZSBhcGFyaWNpb24gZGUgYWxlcnRhcyBvIGVycm9yZXMgZW4gZWwgcHJvZ3JhbWEsIGxvcyBjdWFsZXMgZW5jb250cmFtb3MgcXVlIHByaW5jaXBhbG1lbnRlIHNlZGViZW4gYSBxdWUgbGEgdmVyc2lvbiBubyBlcyBsYSBtYXMgcmVjaWVudGUsIGEgbGEgZmFsdGEgZGUgaW5zdGFsYWNpb24gZGUgcGFxdWV0ZXMgbyBsbGFtYXIgYSBsYXMgbGlicmVyaWFzLCBwcm9sZW1hcyBkZSBlc2NyaXR1cmEgKCp0eXBvcyopIHkgbG9zIG11Y2hvcyBhcmd1bWVudG9zIHF1ZSB0aWVuZW4gbGFzIGZ1bmNpb25lcy4NCg0KU2kgZGVzZGUgdW4gaW5pY28gcHJvZ3JhbWFtb3MgZGUgZm9ybWFzIGVzdHJ1Y3R1cmFkYSwgZGlzY2lwbGluYWRhIHkgbWV0aWN1bG9zYSwgcG9kcmVtb3MgcHJldmVlciBtdWNoYXMgZGUgbGFzIGFsZXJ0YXMgbyBlcnJvcmVzIHkgYXNpIG9idGVuZW1vcyBsb3MgbXVjaG9zIGJlbmVmaWNpb3MgDQo=