エラーバー、2つのシリーズ、光ジッタを伴うドットプロット

Question

私はいくつかの研究を通してデータの集まりを持っています。各研究について、私は性別による変数の平均値に興味があり、これが有意に異なる場合は、各研究について、私は男性と女性の平均と95％の信頼区間を持っています。私はdotplots（ドットプロット、dotplot2、ドットプロット）のいくつかのフレーバーを使用していたが、かなりそこに取得していない

：私がやりたい何

は、これに似たものです。

からDotplotを使用すると、私は1つのシリーズとそのエラーバーを持つことができましたが、私は2番目のシリーズを追加する方法を失っています。

私はDotplotを使用し、アドバイスに続くエラーバーの垂直終了点はhereです。ここで

は、私はこれが3つのデータの列、男性の平均（AVGM）、および95％信頼区間の下限と上限を（プロット

data<-data.frame(ID=c("Study1","Study2","Study3"),avgm=c(2,3,3.5),avgf=c(2.5,3.3,4)) 
data$lowerm <- data$avgm*0.9 
data$upperm <- data$avgm*1.1 
data$lowerf <- data$avgf*0.9 
data$upperf <- data$avgf*1.1 

# Create the customized panel function 
mypanel.Dotplot <- function(x, y, ...) { 
    panel.Dotplot(x,y,...) 
    tips <- attr(x, "other") 
    panel.arrows(x0 = tips[,1], y0 = y, 
       x1 = tips[,2], y1 = y, 
       length = 0.05, unit = "native", 
       angle = 90, code = 3) 
} 

library(Hmisc) 
Dotplot(data$ID ~ Cbind(data$avgm,data$lowerm,data$upperm), col="blue", pch=20, panel = mypanel.Dotplot, 
     xlab="measure",ylab="study") 

を使用していたコードの実施例でありますlowermとupperm）。私は他の3つのシリーズ、同じ研究のために、女性の被験者（avgf、lowerf、upperf）について同じ仕事をする。

私は次のようになり持って結果：

何を一言で言えば、欠落しています。

1. に定義された手段と信頼区間を有する第二シリーズ（avgf）を追加同じ研究のための3つの他の変数

2. いくつかの垂直ジッタを追加して、それらが1つ上にないようにしますもう一方の読者は、重複していても両方を見ることができます。

Answer 1

残念ながら私はDotplotのお手伝いをすることはできませんが、私はggplotを使用して、それはかなり簡単見つけます。データをわずかに並べ替えるだけで済みます。

library(ggplot2) 
# grab data for males 
df_m <- data[ , c(1, 2, 4, 5)] 
df_m$sex <- "m" 
names(df_m) <- c("ID", "avg", "lower", "upper", "sex") 
df_m 

# grab data for females 
df_f <- data[ , c(1, 3, 6, 7)] 
df_f$sex <- "f" 
names(df_f) <- c("ID", "avg", "lower", "upper", "sex") 
df_m 

# bind the data together 
df <- rbind(df_m, df_f) 

# plot 
ggplot(data = df, aes(x = ID, y = avg, ymin = lower, ymax = upper, colour = sex)) + 
    geom_point(position = position_dodge(width = 0.2)) + 
    geom_errorbar(position = position_dodge(width = 0.2), width = 0.1) + 
    coord_flip() + 
    scale_colour_manual(values = c("blue", "red")) + 
    theme_classic() 

# if you want horizontal grid lines you may change the last line with: 
    theme_bw() + 
    theme(panel.grid.major.y = element_line(colour = "grey", linetype = "dashed"), 
     panel.grid.major.x = element_blank(), 
     panel.grid.minor.x = element_blank())