GRAFICOS CON R (II)
Israel
2022-10-27
##GRAFICOS CON R (II) #Nubes de puntos
getwd()## [1] "C:/Users/juan luis/Documents/MMEA/clase 29/09/22"library('ggplot2')
qplot(Sepal.Length, Sepal.Width, data = iris,colour = Species, size = Petal.Width)
#GRAFICOS DE LINEA
qplot(Petal.Length, Sepal.Length, data=iris, color=Species) + geom_line() 
#GRAFICOS DE LINEA CON CURVA DE REGRESION
qplot(Petal.Length, Petal.Width, data = iris, color = Species) + geom_point() + geom_smooth()## `geom_smooth()` using method = 'loess' and formula 'y ~ x'
#COMPOSICION DE MULTIPLES graficas
library('ggplot2')##dibujo de funciones y polinomios
curve(sin, from = -4, to = 4, col = 'blue', lty = 'dotted', lwd = 2, ylab = 'sin(x) vs cos(x)')
curve(cos, from = -4, to = 4, col = 'cyan', lty = 'dashed', lwd = 2, add = T)
legend("topleft", c("sin(x)", "cos(x)"), col = c('blue', 'cyan'), lty = c('dotted','dashed'), lwd = 2, ncol = 2)
curve(x^3-x+1, lty = 2, from = -10, to = 10)
#Circlize
library(circlize)## ========================================
## circlize version 0.4.15
## CRAN page: https://cran.r-project.org/package=circlize
## Github page: https://github.com/jokergoo/circlize
## Documentation: https://jokergoo.github.io/circlize_book/book/
##
## If you use it in published research, please cite:
## Gu, Z. circlize implements and enhances circular visualization
## in R. Bioinformatics 2014.
##
## This message can be suppressed by:
## suppressPackageStartupMessages(library(circlize))
## ========================================mat = matrix(sample(1:100, 18, replace = TRUE), 3,
6)
rownames(mat) = letters[1:3]
colnames(mat) = LETTERS[1:6]
rn = rownames(mat)
cn = colnames(mat)
par(mar = c(1, 1, 1, 1))
circos.par(gap.degree = c(rep(2, nrow(mat) - 1),
10, rep(2, ncol(mat) - 1), 10))
chordDiagram(mat, annotationTrack = "grid", transparency = 0.5,
preAllocateTracks = list(track.height = 0.1))
for (si in get.all.sector.index()) {
circos.axis(h = "top", labels.cex = 0.3, sector.index = si,
track.index = 2)
}
circos.trackPlotRegion(track.index = 1, panel.fun = function(x,
y) {
xlim = get.cell.meta.data("xlim")
ylim = get.cell.meta.data("ylim")
sector.name = get.cell.meta.data("sector.index")
circos.lines(xlim, c(mean(ylim), mean(ylim)),
lty = 3)
for (p in seq(0, 1, by = 0.25)) {
circos.text(p * (xlim[2] - xlim[1]) + xlim[1],
mean(ylim), p, cex = 0.4, adj = c(0.5,
-0.2), niceFacing = TRUE)
}
circos.text(mean(xlim), 1.4, sector.name, niceFacing = TRUE)
}, bg.border = NA)## Note: 1 point is out of plotting region in sector 'a', track '1'.## Note: 1 point is out of plotting region in sector 'a', track '1'.## Note: 1 point is out of plotting region in sector 'b', track '1'.
## Note: 1 point is out of plotting region in sector 'b', track '1'.## Note: 1 point is out of plotting region in sector 'c', track '1'.
## Note: 1 point is out of plotting region in sector 'c', track '1'.## Note: 1 point is out of plotting region in sector 'A', track '1'.
## Note: 1 point is out of plotting region in sector 'A', track '1'.## Note: 1 point is out of plotting region in sector 'B', track '1'.
## Note: 1 point is out of plotting region in sector 'B', track '1'.## Note: 1 point is out of plotting region in sector 'C', track '1'.
## Note: 1 point is out of plotting region in sector 'C', track '1'.## Note: 1 point is out of plotting region in sector 'D', track '1'.
## Note: 1 point is out of plotting region in sector 'D', track '1'.## Note: 1 point is out of plotting region in sector 'E', track '1'.
## Note: 1 point is out of plotting region in sector 'E', track '1'.## Note: 1 point is out of plotting region in sector 'F', track '1'.
## Note: 1 point is out of plotting region in sector 'F', track '1'.
circos.clear()#radarchart
library("fmsb")
set.seed(4242)
maxmin <- data.frame(emotions = c(1, 0), corel = c(1,
0), scene = c(1, 0), yeast = c(1, 0), ebay = c(1,
0))
dat <- data.frame(emotions = runif(3, 0, 1), corel = runif(3,
0, 1), scene = runif(3, 0, 1), yeast = runif(3,
0, 1), ebay = runif(3, 0, 1))
dat <- rbind(maxmin, dat)
radarchart(dat, axistype = 2, plty = 1:3, plwd = 2,
vlabels = names(dat))
#GRAFICAS 3D
library('scatterplot3d')
z <- seq(-10, 10, 0.01)
x <- cos(z)
y <- sin(z)
scatterplot3d(x, y, z, highlight.3d = TRUE, col.axis = "blue",col.grid = "lightblue", main = "Helix", pch = 20)
#GRAFICAS 3D
library("lattice")
z <- matrix(rnorm(625) + 574, nrow = 25)
z <- z + seq(50, 1, length = 25)
persp(z, phi = 30, theta = 30, zlim = c(550, 650),
xlab = "X", ylab = "Y", zlab = "Z", main = "Elevación del terreno")
#GRAFICOS DE TORTUGA
library('grid')
library('TurtleGraphics')
drawTriangle<- function(points){
turtle_setpos(points[1,1],points[1,2])
turtle_goto(points[2,1],points[2,2])
turtle_goto(points[3,1],points[3,2])
turtle_goto(points[1,1],points[1,2])
}
getMid<- function(p1,p2) c((p1[1]+p2[1])/2, c(p1[2]+p2[2])/2)
sierpinski <- function(points, degree){
drawTriangle(points)
if (degree > 0){
p1 <- matrix(c(points[1,], getMid(points[1,], points[2,]), getMid(points[1,], points[3,])), nrow=3, byrow=TRUE)
sierpinski(p1, degree-1)
p2 <- matrix(c(points[2,], getMid(points[1,], points[2,]), getMid(points[2,], points[3,])), nrow=3, byrow=TRUE)
sierpinski(p2, degree-1)
p3 <- matrix(c(points[3,], getMid(points[3,], points[2,]), getMid(points[1,], points[3,])), nrow=3, byrow=TRUE)
sierpinski(p3, degree-1)
}
invisible(NULL)
}
turtle_init(520, 500, "clip")
p <- matrix(c(10, 10, 510, 10, 250, 448), nrow=3, byrow=TRUE)
turtle_col("red")
turtle_do(sierpinski(p, 6))
turtle_setpos(250, 448)
turtle_init()
turtle_init()
turtle_do({
for(j in 1:45) {
for(i in 1:6) {
turtle_forward(20)
turtle_right(360/6)
}
turtle_right(360/45)
}
})
