Showing posts with label R. Show all posts
Showing posts with label R. Show all posts

Wednesday, February 20, 2013

Mostrar valores en barplots usando R

Para poder mostrar los valores que se grafican en un barplot horizontal a la derecha de las barras utilizamos text. Para una tabla de 38 filas y 2 columnas, desplegamos en un barplot la columna dos utilizando 38 colores

Categorías    Valores
cate_1            36
cate_2            36
......
cate_n              1

>barplot(misdatos$val, horiz=T, col=rainbow(38), xlim=c(0,40), main="Barplot con valores que graficados", xlab="Frequency",ylab="Categorias" );




Ahora desplegamos a la derecha de las barras los valores

>text(misdatos$val,bp,misdatos$val,pos=4)




Sunday, January 27, 2013

Instalar la librería ShortRead de Bioconductor

Para poder instalar la librería ShortRead de Bioconductor, se debe de tener previamente instalado las librerías libcurl4-openssl-dev y libxml2-dev en el sistema, por lo que como administradores las instalamos:

#apt-get install libxml2-dev
#apt-get install libcurl4-openssl-dev

Posteriormente en R instalamos la librería ShortRead

>source("http://bioconductor.org/biocLite.R")
>biocLite("ShortRead")
>library(ShortRead)


Tuesday, September 25, 2012

Gráficas de pastel en R

Primero creamos el vector con nuestros porcentajes

> porcentajes = c(39,61)

la suma es 100.

Después le asignamos un nombre a los valores de nuestro vector

> names(porcentajes) = c("Different domains, 39%", "Domains, 61%");

Luego se gráfica

> pie(porcentajes, col = c("gray50", "white"), main = " Domains");

Thursday, August 16, 2012

Insertar gráfica dentro de otra gráfica en R

Para poder insertar una gráfica de menor tamaño dentro de otra mayor, en R utilizamos la función par().
En un ejemplo gráficaremos una distribución de datos con una curva de ajuste, y después insertaremos una gráfica de los mismos datos pero representados en log vs log.


>plot(x$datosx,y$datosy, bty='l', col = "azure4", xlab='Size', ylab='Number of Y', pch =19, xlim = c(0,10000), ylim = c(0,2800),xaxs = 'i', yaxs = 'i');

>points(otrox$,datosx, otroy$datosy, col = "darkred");


>axis(1, las=1, at=500*0:20);

>lines(1.8183459166*x^0.74506574);

En esta línea el rotulo estará en la esquina inferior derecha.
>legend("bottomright" , c('DatosX', 'DatosY', 'Power-law model.  y = 1.8183*x^0.7450', 'R*R: 0.857'), col = c('azure4','darkred','black','white'), pch=c(19, 1,3,0), lty = c(0,0,1,0), lwd=2, bty='n');

En esta linea generamos el especio para la siguiente gráfica que estará en la esquina superior izquierda

>par(fig=c(0.08,0.45,0.48,1), new=TRUE);

>plot(log(x$datosx),log(y$datosy), col = "azure4", main = 'Log vs Log ', xlab='', ylab='',cex.lab = 0.8, cex.main = 1, log = "xy", pch =19);

Después de la función par() podemos graficar  nuestros datos a nuestro gusto, y se verá reflejado en la nueva figura de menor tañano.






Thursday, May 10, 2012

Ajustar ejes a 0 en un plot de R

Desde hace tiempo me había encontrando con la necesidad de poder ajustar los ejes X y Y a 0 de un plot en R. Siempre que generaba un plot había un desfase en los ejes con respecto al cero.
 

La solución al ajuste de los ejes a 0 la obtuve con el parámetro de estilo de los ejes xaxs y yaxs.


Este es el código:

> plot(0.00016298*1.672059*x^0.672059, col = 'darkblue', type = 'l', ylim = c(0,1.6), xaxs = 'i', yaxs = 'i' );

Les dejo la liga a la página de dodne obtuve la información: Climate Charts & Graphs.