Articles of r

Abordar x e y en aes por número de variable

Necesito dibujar un diagtwig de dispersión con las variables de direccionamiento por sus números de columna en lugar de nombres, es decir, en lugar de ggplot(dat, aes(x=Var1, y=Var2)) Necesito algo como ggplot(dat, aes(x=dat[,1], y=dat[,2])) . (Digo ‘algo’ porque este último no funciona). Aquí está mi código: showplot1<-function(indata, inx, iny){ dat<-indata print(nrow(dat)); # this is just […]

terminando el prompt “+” en R

Soy bastante nuevo en R e hice un tipo-o al ingresar algún código y ahora, en lugar de recibir el mensaje “>”, aparece el mensaje “+”. Parece que estoy atrapado en algún tipo de función que está buscando información, pero parece que no puedo salir de ella. Me imagino que debe haber una mejor manera […]

La condición tiene una longitud> 1 y solo se usará el primer elemento

Tengo un dataframe, viaje: > head(trip.mutations) Ref.y Variant.y 1 TC 2 GC 3 AC 4 TC 5 CA 6 GA Deseo agregar una tercera columna, mutType, que siga estas reglas: for (i in 1:nrow(trip)) { if(trip$Ref.y==’G’ & trip$Variant.y==’T’|trip$Ref.y==’C’ & trip$Variant.y==’A’) { trip[i, ‘mutType’] <- "G:C to T:A" } else if(trip$Ref.y=='G' & trip$Variant.y=='C'|trip$Ref.y=='C' & trip$Variant.y=='G') { […]

Diagtwig de violín dividido con ggplot2

Me gustaría crear una gráfica de densidad de violín dividida usando ggplot, como el cuarto ejemplo en esta página de la documentación nacida. Aquí hay algunos datos: set.seed(20160229) my_data = data.frame( y=c(rnorm(1000), rnorm(1000, 0.5), rnorm(1000, 1), rnorm(1000, 1.5)), x=c(rep(‘a’, 2000), rep(‘b’, 2000)), m=c(rep(‘i’, 1000), rep(‘j’, 2000), rep(‘i’, 1000)) ) Puedo trazar violines esquivados así: library(‘ggplot2’) […]

R: gsub, patrón = vector y reemplazo = vector

Como dice el título, estoy tratando de usar gsub donde uso un vector para el “patrón” y “reemplazo”. Actualmente, tengo un código que se ve así: names(x1) <- gsub("2110027599", "Inv1", names(x1)) #x1 is a data frame names(x1) <- gsub("2110025622", "Inv2", names(x1)) names(x1) <- gsub("2110028045", "Inv3", names(x1)) names(x1) <- gsub("2110034716", "Inv4", names(x1)) names(x1) <- gsub("2110069349", "Inv5", […]

Mostrar un reloj de tiempo en la línea de comando R

Me pregunto si hay una manera de mostrar la hora actual en la línea de comando R, como en MS DOS, podemos usar Prompt $T $P$G para incluir el reloj de tiempo en cada línea de aviso. Algo como options(prompt=paste(format(Sys.time(), “%H:%M:%S”),”> “)) lo hará, pero luego se arregla en el momento en que se configuró. […]

Encuentra valores duplicados en R

Tengo una tabla con 21638 filas * únicas: vocabulary <- read.table("http://socserv.socsci.mcmaster.ca/jfox/Books/Applied-Regression-2E/datasets/Vocabulary.txt", header=T) Esta tabla tiene cinco columnas, la primera de las cuales contiene los números de identificación del encuestado. Quiero verificar si los encuestados aparecen dos veces o si todos los encuestados son únicos. Para contar identificaciones únicas que puedo usar length(unique(vocabulary$id)) y para verificar […]

Trazar tags en los extremos de las líneas

Tengo los siguientes datos ( temp.dat ver nota final para obtener datos completos) Year State Capex 1 2003 VIC 5.356415 2 2004 VIC 5.765232 3 2005 VIC 5.247276 4 2006 VIC 5.579882 5 2007 VIC 5.142464 … y puedo producir la siguiente tabla: ggplot(temp.dat) + geom_line(aes(x = Year, y = Capex, group = State, colour […]

Controlar el orden de visualización de la leyenda ggplot2

¿Alguien sabe cómo puedo controlar el orden de las leyendas en ggplot2? Por lo que puedo ver, el orden aparece relacionado con las tags de la escala real en lugar del orden de la statement de la escala. Cambiar los títulos de la escala altera el orden. He hecho un pequeño ejemplo usando el dataset […]

¿Cómo transformar datos XML en un data.frame?

Estoy intentando aprender el paquete XML de R. Estoy tratando de crear un data.frame del archivo de datos books.xml xml data. Esto es lo que obtengo: library(XML) books <- "http://www.w3schools.com/XQuery/books.xml" doc <- xmlTreeParse(books, useInternalNodes = TRUE) doc xpathApply(doc, "//book", function(x) do.call(paste, as.list(xmlValue(x)))) xpathSApply(doc, "//book", function(x) strsplit(xmlValue(x), " ")) xpathSApply(doc, "//book/child::*", xmlValue) Cada uno de estos […]