library(ggplot2)
#rnorm(1000)模拟产生1000个正太分布数
ID <-c(seq(1,1000))
randNorm <-rnorm(1000)
#密度函数
randDensity <- dnorm(randNorm)
#σ=0.5
randDensity1 <- dnorm(randNorm,0,0.5)
ID1 <-c(seq(1001,2000))
#把两个密度函数合并放到一个plot
data <- rbind(data.frame(ID,x=randNorm,y=randDensity),data.frame(ID=ID1,x=randNorm,y=randDensity1))
#plot(dnorm(x, mean = 0, sd = 1)) in one plotmap
ggplot(data)+aes(x=x,y=y,group = ID)+geom_point(colour=((data[,1]%/%1000)+1))+labs(x="Random Normal Varables",y="randDensity")
plot(data.frame(x=randNorm,y=randDensity),col = "red",xlab = "Random Normal Varables",ylab = "randDensity",main = "模拟正太分布")
#curve(dnorm(x),from=-4,to=4,main="标准正太分布")