INSTITUTTO TECNOLOGICO DE TIJUANA PROBABILIDAD Y ESTADISTICA

PRACTICA 1

OBJETIVO: Gráfica de una variable y resumen de estadisticas.

setwd("/home/jorgher/Dropbox/RCourse/")
list.files("diet_data")

## [1] "Andy.csv"  "David.csv" "John.csv"  "Mike.csv"  "Steve.csv"

rm(list=ls())
setwd("/home/jorgher/Dropbox/RCourse/")
list.files("diet_data") #Listado de los archivos en el directorio indicado

## [1] "Andy.csv"  "David.csv" "John.csv"  "Mike.csv"  "Steve.csv"

andy <- read.csv("diet_data/Andy.csv")

Estadísiticas de Andy

summary(andy)

##  Patient.Name      Age         Weight           Day       
##  Andy:30      Min.   :30   Min.   :135.0   Min.   : 1.00  
##               1st Qu.:30   1st Qu.:137.0   1st Qu.: 8.25  
##               Median :30   Median :137.5   Median :15.50  
##               Mean   :30   Mean   :137.3   Mean   :15.50  
##               3rd Qu.:30   3rd Qu.:138.0   3rd Qu.:22.75  
##               Max.   :30   Max.   :140.0   Max.   :30.00

Gráfica de comparación entre Andy y David

pesoAndy <- andy$Weight
dia <- andy$Day
David <- read.csv("diet_data/David.csv")
pesoDavid <- David$Weight

La gráfica del peso durante el tiempo es:

par(mfrow=c(2,1))
plot(dia,pesoAndy, col = "red", type = "l", main = "Comportamiento del peso Andy", )
plot(dia,pesoDavid, col = "blue", type = "l", main = "Comportamiento del peso David")

Guardando la gráfica en un archivo:

jpeg("pesosA.jpg")
par(mfrow=c(2,1))
plot(dia,pesoAndy, col = "red", type = "l", main = "Comportamiento del peso Andy", )
plot(dia,pesoDavid, col = "blue", type = "l", main = "Comportamiento del peso David")
dev.off()

## png 
##   2