miércoles, 6 de mayo de 2015

LIBRERIAS


####################################
###   R para muy principiantes   ###
####################################


####################
###   Librerías  ###
####################


# Hay funciones que se pueden ejecutar mediante comandos que ya están en el programa base,
# como por ejemplo la media de una serie de valores, que se calcula mendiante el comando mean().

a <- c(2,3,4,4,4,4,4,6,6,8,10)
mean(a)

# Hay otros paquetes de comandos que ni siquiera están instalados y hay que descargarlos
# antes de poder usar las funciones que éstos contienen.
# En www.r-project.org puedes ver la lista de paquetes.
# Hay mas de 5.000 y pueden ser muy específicos sobre un tema.

install.packages ("agricolae")
remove.packages ("agricolae")

# Es posible que la primera vez que instalas un paquete te pida que elijas una localización
# desde la que quieras instalarlo. Simplemente elige el país en el que te encuentres.

# Para saber los comandos del paquete y sus funciones, basta con pedirlo con la ayuda.

help (package = agricolae)

# Para usar comandos de las librerías que hemos instalados, hay que cargar los paquetes
# en la memoria primero. Esto habrá que hacerlo cada vez que iniciemos sesión.

library(agricolae)

# Ejemplo bar.err

library(agricolae)
data(sweetpotato)
model<-aov(yield~virus,data=sweetpotato)
out <- waller.test(model,"virus", console=TRUE,
                   main="Yield of sweetpotato\ndealt with different virus")
par(mfrow=c(2,2),cex=1)
bar.err(out$means,variation="range",horiz=TRUE,xlim=c(0,45),angle=125,density=6,
        main="range")
bar.err(out$means,variation="SD",ylim=c(0,45),col=colors()[30],
        main="Standard deviation",density=8)
bar.err(out$means,variation="SE",horiz=TRUE,xlim=c(0,45),density=8,
        col="brown",main="Standard error")
bar.err(out$means,variation="range",ylim=c(0,45),bar=FALSE,col="green",
        main="range")
par(mfrow=c(1,2),cex=1)
bar.err(out$means,variation="range",ylim=c(0,45),bar=FALSE,col=0)
abline(h=0)
# horiz = TRUE
bar.err(out$means,variation="SE",horiz=TRUE,xlim=c(0,45),bar=FALSE,col=0)
#startgraph
par(mfrow=c(1,1))
#endgraph

No hay comentarios:

Publicar un comentario