Articles of r

Mantener ceros al final

En R, me gustaría seguir ceros al final. Por ejemplo, si escribo round(5.2, 3) Me gustaría que la salida fuera 5.200. ¿Me puedes ayudar?

¿Cómo extraer los colores de relleno de un objeto ggplot?

Estoy intentando escribir algunas pruebas unitarias automatizadas para una serie de funciones que generan gráficos ggplot . Por ejemplo, quiero establecer una escala de color específica para un gráfico. Ahora necesito una forma de determinar si realmente se aplicó la escala de color correcta. El fondo: Aquí hay un código de ejemplo, que establece el […]

poly () en lm (): diferencia entre crudo vs. ortogonal

yo tengo library(ISLR) attach(Wage) # Polynomial Regression and Step Functions fit=lm(wage~poly(age,4),data=Wage) coef(summary(fit)) fit2=lm(wage~poly(age,4,raw=T),data=Wage) coef(summary(fit2)) plot(age, wage) lines(20:350, predict(fit, newdata = data.frame(age=20:350)), lwd=3, col=”darkred”) lines(20:350, predict(fit2, newdata = data.frame(age=20:350)), lwd=3, col=”darkred”) Las líneas de predicción parecen ser las mismas, sin embargo, ¿por qué los coeficientes son tan diferentes? ¿Cómo los interpretas en raw=T y raw=F Veo […]

¿Cómo puedo verificar si una llamada a función da como resultado una advertencia?

En R, ¿cómo puedo determinar si una llamada a función da como resultado una advertencia? Es decir, después de llamar a la función me gustaría saber si esa instancia de la llamada produjo una advertencia.

Extrae el valor p de aov

Estoy buscando extraer el valor p generado de una anova en R. Esto es lo que estoy ejecutando: test <- aov(asq[,9] ~ asq[,187]) summary(test) Rendimientos: Df Sum Sq Mean Sq F value Pr(>F) asq[, 187] 1 3.02 3.01951 12.333 0.0004599 *** Residuals 1335 326.85 0.24483 — Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 […]

R desde C – Helloworld posible más simple

¿Cuál es la función C más simple posible para iniciar el intérprete R, pasando una pequeña expresión (por ejemplo, 2 + 2) y obteniendo el resultado? Estoy intentando comstackr con MingW en Windows.

Cómo dibujar el diagtwig de caja con un nivel significativo?

Hace algún tiempo, hice una pregunta acerca de dibujar el diagtwig de caja Link1 . Tengo algunos datos con 3 grupos diferentes (o tags) Por favor, cargue aquí . Puedo usar los siguientes códigos R para obtener el diagtwig de caja library(reshape2) library(ggplot2) morphData <- read.table(".\\TestData3.csv", sep=",", header=TRUE); morphData.reshaped <- melt(morphData, id.var = "Label") ggplot(data […]

¿Cómo agregar una nueva columna a un dataframe (al frente, no al final)?

Cómo agregar una nueva variable a un dataframe existente, pero quiero agregar al frente y al final. p.ej. mi dataframe es bcd 1 2 3 1 2 3 1 2 3 Quiero agregar una nueva variable a, para que el dataframe se vea como abcd 0 1 2 3 0 1 2 3 0 1 […]

Aplicar una función a un subconjunto de columnas data.table, por column-indices en lugar de name

Intento aplicar una función a un grupo de columnas en una gran tabla de datos sin referirme a cada una individualmente. a <- data.table( a=as.character(rnorm(5)), b=as.character(rnorm(5)), c=as.character(rnorm(5)), d=as.character(rnorm(5)) ) b <- c('a','b','c','d') con el MWE anterior, esto: a[,b=as.numeric(b),with=F] funciona, pero esto: a[,b[2:3]:=data.table(as.numeric(b[2:3])),with=F] no funciona ¿Cuál es la forma correcta de aplicar la función as.numeric a […]

porcentaje en y lab en un diagtwig de barras ggplot facetado?

haciendo facetas en ggplot A menudo me gustaría usar el porcentaje en lugar de conteos. p.ej test1 <- sample(letters[1:2], 100, replace=T) test2 <- sample(letters[3:8], 100, replace=T) test <- data.frame(cbind(test1,test2)) ggplot(test, aes(test2))+geom_bar()+facet_grid(~test1) Esto es muy fácil, pero si N es diferente en la faceta A en comparación con la faceta B, sería mejor, creo, comparar porcentajes, […]