Linear Model with Gaussian Noise for Valentine’s day

#----------------------------------------------------------------------#
# Number of random points
#----------------------------------------------------------------------#
k = 4000

#----------------------------------------------------------------------#
# Generate x, which follows normal distribution
#----------------------------------------------------------------------#
set.seed(7)
x <- rnorm(n = k, mean = 0, sd = 1)

#----------------------------------------------------------------------#
# Generate e (noise), which follows normal distribution too
#----------------------------------------------------------------------#
set.seed(9)
e <- rnorm(n = k, mean = 0, sd = 1)

#----------------------------------------------------------------------#
# y is an additive function of x and e
#----------------------------------------------------------------------#
y <- abs(x) + e

#----------------------------------------------------------------------#
# Scatter plot
#----------------------------------------------------------------------#
plot(x,y, col='tomato', pch = '/', main = '', xlab = '', ylab = '', xaxt='n', yaxt='n')