Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- library("ggplot2") # library ggplot
- set.seed(1) # make example reproducible
- n <- 5000 # set sample size
- df <- data.frame(x= rnorm(n), g= factor(rep(c(0,1), n/2))) # generate data
- df$y <- NA # include y in df
- df$y[df$g== 0] <- df$x[df$g== 0]**2 + rnorm(sum(df$g== 0))*5 # y for group g= 0
- df$y[df$g== 1] <-2 + df$x[df$g== 1]**2 + rnorm(sum(df$g== 1))*5 # y for g= 1 (with intercept 2)
- ggplot(df, aes(x, y, col= g)) + geom_smooth() + geom_point(alpha= .1) # make a plot
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement