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 format your code

Formatting Code: the Basics All code or console output you include in your posts should be formatted properly. Luckily, this is very easy to do! Just use the code formatting button at the top of the post editing box: [image] Select some code Click the </> button! The code formatting button automatically adds special Markdown formatting symbols to whatever you selected. If you check the preview pane to the right of your post, you'll see that the forum software now outputs nicely formatted code . [04] The code formatting button is pretty smart! It can tell the difference between snippets of code that are part of a sentence and blocks of code that stand alone. …

[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.