Editorial

La vaca renovable es un espacio para compartir información sobre fuentes de energía renovable y ciencia. Este blog fue creado y es administrado por estudiantes y está dirigido a tod@s l@s interesad@ en encontrar mejores formas de interacción con el medio.

viernes, 22 de enero de 2010

Gráfica de barras en R-project

Esta instrucción permite crear gráficas de barras en R.


ConsumoPorcentual = data.frame(

Doméstico = c(11.4, 13.94),

Alumbrado = c(4.43, 0),

Bombeo = c(0.67, 0),

Agrícola = c(0.34,0),

Industrial = c(69.23,0)

)


barplot(as.matrix(ConsumoPorcentual),

main="Distribución en el consumo eléctrico",

col=c("grey","red"),

ylim=c(0,70))


El resultado es el siguiente:


Este ejemplo utiliza una serie de datos pequeña, pero para bases de datos más grandes, puedes cargar los datos directamente de un archivo externo con la siguiente instrucción:

ConsumoPorcentual = read.csv("/home/documents/Datos_de_ConsumoPorcentual.csv", row.names=1)



jueves, 21 de enero de 2010

¿Qué es R-project?

R-proyect es un paquete computacional análisis estadístico y creación de gráficos. El paquete es muy flexible y existen versiones para Windows, Mac y Linux. La pantalla principal de R y la salida de sus gráficos son así:



¿Cuáles son las ventajas de R?
R es un software libre y por lo tanto es gratuito.
R permite hacer prácticamente cualquier tipo de análisis estadístico y gráfico.
R es muy flexible y permite utilizar bases de datos en diferentes formatos.

¿Cuáles son las desventajas de R?
R no tiene muy pocos botones; la forma de dar instrucciones al programa es escribiéndolas en la línea de comando programa. A principio puede ser un poco lento, pero en cuanto te familiarizas con el lenguaje, es muy rápido y simple.

¿Cómo puedo bajar R?
Puedes bajar el software en la dirección de R-project: www.r-project.org


Creación de gráficos en R

Una vez que has cargado tus datos, puedes crear un gráfico con la siguiente instrucción:

plot(x[,"Variable_1"],x[,"Variable_2"], ylab="Títulos del eje y",xlab="Título del eje x",main="Título del gráfico")


El resultado de la gráfica es el siguiente:


¿Te parece una gráfica aburrida? ¿Qué te parece si hacemos exactamente lo mismo pero con la base de datos completa (y con los encabezados reales)? La instrucción en este caso es:

plot(x[,"HID"],x[,"Energy"], ylab="Energía per capita",xlab="Índice de desarrollo humano (HDI)",main="Relación entre índice de desarrollo humano y consumo de energía")



Ahora agreguemos una escala logarítmica en el eje y.

plot(x[,"HID"],x[,"Energy"], ylab="Energía per capita",xlab="Índice de desarrollo humano (HDI)",main="Relación entre índice de desarrollo humano y consumo de energía", log="y")



Ahora vamos a hacer un modelo lineal de los datos. La instrucción es muy simple. "lm" quiere decir "modelo lineal".

plot(lm(x))
Existen infinidad de gráficos. De los cuales platicaré más adelante.





Cómo cargar datos en R

Esta es una instrucción simple para cargar una serie de datos en r-project (www.r-project.org).

x =read.table(file="nombre_de_tu_archivo.txt",header=T)

Como ejemplo, puedes utilizar un archivo que tenga la siguiente información.

Variable_1 Variable_2
0.023187441 0.366
0.01608756 0.389
0.020458337 0.431
0.008287151 0.442
0.033587175 0.453
0.019128672 0.459
0.021608317 0.479
0.033174043 0.484
0.04010257 0.499
0.014381913 0.502

Puedes nombrar ese archivo como gustes, en este ejemplo se llama: "nombre_de_tu_archivo.txt". Si tienes tus datos en excel puedes exportarla a "texto separado por tabulaciones [txt]" o a "valores separados por comas [csv]".

Ahora utiliza la siguiente instrucción para cargar los datos. Como podrás ver, guardará los datos en la variable "x".

x =read.table(file="nombre_de_tu_archivo.txt",header=T)

Después puedes confirmar que tus datos se cargaron correctamente escribiendo.

x

Y presionando ENTER. El programa te dará como respuesta la lista de tus datos.

Si quieres extraer tercer dato de la primera columna, puedes dar la instrucción:

x[3,1]

Si quieres extraer la primera columna, puedes dar la instrucción:

x[,"Variable_1"]

que equivalente a escribir:

x[,1]