Como título, preciso desenhar algo como isto:
O ggplot, ou outros pacotes, se o ggplot não for capaz, pode ser usado para desenhar algo assim?
r
data-visualization
ggplot2
funnel-plot
lokheart
fonte
fonte
stat_quantile()
para colocar quantis condicionais em um gráfico de dispersão. Você pode controlar a forma funcional da regressão quantil com o parâmetro formula. Eu sugeriria coisas como formula =y~ns(x,4)
para obter um ajuste estriado suave.Respostas:
Embora haja espaço para melhorias, aqui está uma pequena tentativa com dados simulados (heterocedásticos):
fonte
Se você procura esse tipo de gráfico de funil (meta-análise) , o seguinte pode ser o ponto de partida:
fonte
linetype=2
argumento dentro dosaes()
colchetes - plotando as linhas de 99% - gera um erro "a variável contínua não pode ser mapeada para o tipo de linha" com o atual ggplot2 (0.9.3.1). Que alterageom_line(aes(x = number.seq, y = number.ll999, linetype = 2), data = dfCI)
ageom_line(aes(x = number.seq, y = number.ll999), linetype = 2, data = dfCI)
obras para mim. Sinta-se livre para alterar a resposta original e perdê-la.Veja também o pacote cran cran berryFunctions, que possui um funilPlot para proporções sem usar o ggplot2, se alguém precisar dele nos gráficos básicos. http://cran.r-project.org/web/packages/berryFunctions/index.html
Há também o pacote extfunnel, que eu não olhei.
fonte
O código de Bernd Weiss é muito útil. Fiz algumas alterações abaixo, para alterar / adicionar alguns recursos:
geom_segment
vez dageom_line
linha demarcadora da média meta-analítica, para que ela tenha a mesma altura das linhas demarcadoras das regiões de confiança de 95% e 99%Meu código usa uma média meta-analítica de 0,0892 (se = 0,0035) como exemplo, mas você pode substituir seus próprios valores.
fonte