Articles of r

Editando legenda (texto) labels em ggplot

Passei horas procurando na documentação e no StackOverflow, mas nenhuma solução parece resolver o meu problema. Ao usar o ggplot , não consigo obter o texto correto na legenda, mesmo que esteja no meu dataframe. Eu tentei scale_colour_manual , scale_fill_manual com valores diferentes para labels= como c(“T999”, “T888″)”, “cols” . Aqui está o meu código: […]

Declaração de Caso Equivalente em R

Eu tenho uma variável em um dataframe em que um dos campos normalmente tem 7-8 valores. Eu quero collpase-los 3 ou 4 novas categorias dentro de uma nova variável dentro do dataframe. Qual é a melhor abordagem? Eu usaria uma instrução CASE se eu estivesse em uma ferramenta semelhante a SQL, mas não sabia como […]

Como corrigir a proporção em ggplot?

Eu estou tentando resize um enredo para caber no meu documento, mas estou tendo dificuldades em obter o diagrama plotado como um quadrado. Exemplo: pdf(file = “./out.pdf”, width = 5, height = 5) p <- ggplot(mydata, aes(x = col1, y = col2)) print(p) aux <- dev.off() Embora os limites para x e y sejam os […]

usando substituto para obter o nome do argumento com

Eu estou tentando obter os nomes de argumentos no ambiente global dentro de uma function. Eu sei que posso usar o substituto para obter o nome dos argumentos nomeados, mas gostaria de poder fazer a mesma coisa com … argumentos. Eu meio que tenho que trabalhar para o primeiro elemento de … mas não consigo […]

R Troca de bandas de confiança com ggplot

Eu gostaria de criar uma banda de confiança para um modelo equipado com gls como este: require(ggplot2) require(nlme) mp <-data.frame(year=c(1990:2010)) mp$wav <- rnorm(nrow(mp))*cos(2*pi*mp$year)+2*sin(rnorm(nrow(mp)*pi*mp$wav))+5 mp$wow <- rnorm(nrow(mp))*mp$wav+rnorm(nrow(mp))*mp$wav^3 m01 <- gls(wow~poly(wav,3), data=mp, correlation = corARMA(p=1)) mp$fit <- as.numeric(fitted(m01)) p <- ggplot(mp, aes(year, wow))+ geom_point()+ geom_line(aes(year,fit)) p Isso só representa os valores ajustados e os dados, e eu […]

Extrair banda de previsão do ajuste lme

Tenho seguindo modelo x <- rep(seq(0, 100, by=1), 10) y <- 15 + 2*rnorm(1010, 10, 4)*x + rnorm(1010, 20, 100) id <- NULL for(i in 1:10){ id <- c(id, rep(i,101)) } dtfr <- data.frame(x=x,y=y, id=id) library(nlme) with(dtfr, summary( lme(y~x, random=~1+x|id, na.action=na.omit))) model.mx <- with(dtfr, (lme(y~x, random=~1+x|id, na.action=na.omit))) pd <- predict( model.mx, newdata=data.frame(x=0:100), level=0) with(dtfr, plot(x, […]

Existe uma maneira de manipular quebras e labels de escala ggplot?

ggplot geralmente faz um bom trabalho ao criar intervalos sensatos e labels em escalas. No entanto, acho que no enredo com muitas facetas e, talvez, um formatter= declaração, os labels tendem a ficar muito “denso” e overprint, por exemplo, nesta imagem: df <- data.frame( fac=rep(LETTERS[1:10], 100), x=rnorm(1000) ) ggplot(df, aes(x=x)) + geom_bar(binwidth=0.5) + facet_grid(~fac) + […]

Remodelando de largura para longa com várias colunas de valores

Eu preciso remodelar minha tabela ampla em formato longo, mas mantendo vários campos para cada registro, por exemplo: dw <- read.table(header=T, text=' sbj f1.avg f1.sd f2.avg f2.sd blabla A 10 6 50 10 bA B 12 5 70 11 bB C 20 7 20 8 bC D 22 8 22 9 bD ') # Now […]

Extrair todo enésimo elemento de um vetor

Eu gostaria de criar um vetor em que cada elemento é o elemento i+6th de outro vetor. Por exemplo, em um vetor de comprimento 120 eu quero criar outro vetor de comprimento 20 em que cada elemento é o valor i, i+6, i+12, i+18… do vetor inicial, ou seja, eu quero extrair cada 6º elemento […]

Agregar dados diários a intervalos de mês / ano

Muitas vezes não tenho que trabalhar com datas em R, mas imagino que isso seja bastante fácil. Eu tenho uma coluna que representa uma data em um dataframe. Eu simplesmente quero criar um novo dataframe que resume uma segunda coluna por mês / ano usando a data. Qual é a melhor abordagem? Eu quero um […]