Introducción

Este es el primer trabajo realizado en este curso como forma de instroducción al sistema R. Como primer paso, realizamos un codigo para resolver ecuaciones de segundo grado usando la formula general (la chicharronera).

###SOLUCION DE ECUASIONES DE SEGUNDO GRADO 
###FORMULA GENERAL
a<-0.033
b<-5
c<--4
d<-sqrt(b*b-4*a*c)
x1<-(-b+d)/(2*a)
x2<-(-b-d)/(2*a)
print(x1)
## [1] 0.79582
print(x2)
## [1] -152.311

ESTADISTICA BASICA

Tenemos una variable con las edades de las personas en el salon y calculamos algunos estadíticos básicos

edad<-c(40,22,24,24,23,22,22,50,21,22,60)
media<-mean(edad)
###LA MEDIA O PROMEDIO ES LO MISMO
mediana<-median.default(edad)
###ES EL VALOR QUE QUEDA EN MEDIO DE LOS DAROS ORDENADOS EN FORMA ASCENDETE 
###EN EL EJEMPLO ES 21,22,22,22,22,23,24,24,40,50,60
library(modeest)
moda<-mfv(edad)
s<-sd(edad)
media 
## [1] 30
mediana 
## [1] 23
moda 
## [1] 22
s
## [1] 13.63085

GRÁFICAS DE LOS RESULTADOS

Exploramos graficamente nuestras variable de edad

boxplot(edad,ylab="edad(años)", main="Grafica 1.Distribucion de la edad",
        xlab="Datos totales")

En la gráfica, vimos dos valores atípicos, los barramos y recalculamos los estadísticos y visualizamos la gráfica

edad1<-c(40,22,24,24,23,22,22,21,22)
media1<-mean(edad1)
mediana1<-median.default(edad1)
s1<-sd(edad1)
moda1<-mfv(edad1)
###PROMEDIO DE LAS DESVIACIONES ENTRE UN RANGO DE NUMEROS 
media1
## [1] 24.44444
mediana1
## [1] 22
moda1
## [1] 22
s1
## [1] 5.918427
boxplot(edad1,ylab="edad(años)", main="Grafica 2.Distribucion de la edad",
        xlab="Datos totales")

En la gráfica, encontramos un valor atípico, lo borramos y recalculamos los estdísticos y visualizamos la gráfica

edad2<-c(22,24,24,23,22,22,21,22)
media2<-mean
mediana2<-median.default(edad2)
s2<-sd(edad2)
moda2<-mfv(edad2)
media2
## function (x, ...) 
## UseMethod("mean")
## <bytecode: 0x000001565cc26008>
## <environment: namespace:base>
mediana2
## [1] 22
moda2
## [1] 22
s2
## [1] 1.069045
boxplot(edad2,ylab="edad(años)", main="Grafica 3.Distribucion de la edad",
        xlab="Datos totales")