ggplot2 facet_wrap（）利用可能なすべてのデータを

Question

再現性の例表示されていない：私は、次の取得

library(ggplot2) 
library(haven) 
library(dplyr) 

data <- read_dta('http://dl.dropboxusercontent.com/s/s7zqb2e0avyp1gk/nswd_old_12.dta') 
data$treatedOrNonSample <- ifelse(data$sample == 1, 0, ifelse(data$treat == 1, 1, 2)) 
treatOrCPS <- subset(data, data$treatedOrNonSample!=2) 

m_ps <- glm(treat ~ age + age2 + ed + hisp + married + nodeg + re74 + re75 + black, 
        family=binomial(link="probit"), data=treatOrCPS) 

prs_df <- data.frame(pr_score = predict(m_ps, type = "response"), 
         treated = m_ps$model$treat) 

labs <- paste("Status:", c("Treated", "CPS Sample")) 

prs_df %>% 
    mutate(treated = ifelse(treated == 1, labs[1], labs[2])) %>% 
    ggplot(aes(x = pr_score)) + 
    geom_histogram(color = "white") + 
    facet_wrap(~treated) + 
    xlab("Probability of being in Treatment Condition") + 
    theme_bw() 

を：

なぜ一番右の状態で使用可能なデータが存在しません？私はこの一見些細な話題に私の心を失っているので、どんな助けも高く評価されます。 color=whiteを削除

Answer 1

は、データが表示されるようになり：しかし

prs_df %>% 
    mutate(treated = ifelse(treated == 1, labs[1], labs[2])) %>% 
    ggplot(aes(x = pr_score)) + 
    geom_histogram() + 
    facet_wrap(~treated) + 
    xlab("Probability of being in Treatment Condition") + 
    theme_bw() 

、「CPSサンプル」グループ内のカウントが途中高いため、「治療」グループ内のデータはまだほとんど見えません。