Ggplot legend error

[jystat]

I'm trying to put the legend for my density plot with the theoretical plot in the histogram.

## lambda is 5
lambda <- 5

## the number of simulations is 20,000
nsim <- 20000

## the number of exponentials is 30
n_30 <- 30

## 20,000 simulations are run and a histogram of the results is generated
mns = NULL
for (i in 1 : nsim) {
        mns <- c(mns, mean(rexp(n_30, lambda)))
}
hist(mns, col = "grey", main = "Frequency Distribution of Means of 30 Exponentials")

mean(mns)

## load the ggplot2 package
library(ggplot2)

data <- as.data.frame(mns)

## plot histogram with an overlay of the theoretical and sample distribution
plot <- ggplot(data, aes(x = mns))
plot <- plot + geom_histogram(aes(y=..density..), colour="black", fill="white")
plot <- plot + geom_density(aes(color = "Simulated"), size = 1.0, linetype="dashed")
plot <- plot + stat_function(aes(colour = "Normal"), fun = dnorm, 
                             args = list(mean = mean(mns), sd = sd(mns)), 
                             color = "blue", 
                             size = 1.0, show.legend = TRUE)
plot <- plot + labs(title = "Sample vs. Theoretical Distributions of Means of 30 Exponentials", 
                    x = "Means of exponential distributions", 
                    y = "Density") + scale_colour_manual("Lgend title", values = c("red", "blue"))

print(plot)

However, when I run this code in my RMarkdown, the legend for stat_function (Normal) is not showing correctly. Can anybody help on this?

Thank you.

[mara]

Are you saying that the red isn't appearing on the legend when you run it in RMarkdown? Also, could you please format the code in code chunks? It makes for easier reading!

FAQ: How to make your code look nice? Markdown Formatting

Jump to Contents Code Chunks Quote Snippets Show-Hide Spoilers Markdown Text Formatting Code Chunks You type in community: ``` r #> A code chunk library(foo) ``` You get: #> A code chunk library(foo) This is easier to read. It also creates a clear chunk of code for others can copy and paste into their own R environment. Inline Code Type this: I really wanted you to see `this_code_here` See this: I really wanted you to see this_code_here Quote Snippets If you highlight the content and click quote it will automatically populate the reply with the formatted code: [image] [quote="stephhazlitt, post:5, topic:6246, full:true"] How do you insert a quote from the th…

[jystat]

Thank you for your comment. @mara

Yes, that's right. I'm trying to figure out why red isn't appearing on the legend when I run in RMarkdown.

[mara]

Hmm…have you tried starting a fresh session, and then knitting? reprex uses rmarkdown::render() under the hood, and here's a reprex of your code:

lambda <- 5
nsim <- 20000
n_30 <- 30
mns = NULL
for (i in 1 : nsim) {
  mns <- c(mns, mean(rexp(n_30, lambda)))
}
nsim
#> [1] 20000
## load the ggplot2 package
library(ggplot2)

data <- as.data.frame(mns)

## plot histogram with an overlay of the theoretical and sample distribution
plot <- ggplot(data, aes(x = mns))
plot <- plot + geom_histogram(aes(y=..density..), colour="black", fill="white")
plot <- plot + geom_density(aes(color = "Simulated"), size = 1.0, linetype="dashed")
plot <- plot + stat_function(aes(colour = "Normal"), fun = dnorm, 
                             args = list(mean = mean(mns), sd = sd(mns)), 
                             color = "blue", 
                             size = 1.0, show.legend = TRUE)
plot <- plot + labs(title = "Sample vs. Theoretical Distributions of Means of 30 Exponentials", 
                    x = "Means of exponential distributions", 
                    y = "Density") + scale_colour_manual("Lgend title", values = c("red", "blue"))

print(plot)
#> `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

^{Created on 2018-09-27 by the reprex package (v0.2.1.9000)}

[jystat]

@mara
I tried with refresh mode, but it is still not showing correctly.

[mara]

What happens if you run a reprex of the code, as I just did above.

install.packages("reprex")

For pointers specific to the community site, check out the reprex FAQ.