Histogram

library(foreign)

data=read.dta("C:/Users/BINH THANG-TRAN/Dropbox/PhD/Research/PhD thesis/Methods/morkov/Utility_study/Data_04DEC2019.dta")

1 1. ggplot

# add a normal distribution line in histogram
hist(data$score2sh, freq=FALSE, col="gray", xlab="Petal Length", main="title here ")
curve(dnorm(x, mean=mean(data$score2sh), sd=sd(data$score2sh)), add=TRUE, col="red") #line

# add a normal distribution line in histogram
hist(data$c2, freq=FALSE, col="gray", xlab="Petal Length", main="title here ")
curve(dnorm(x, mean=mean(data$c2), sd=sd(data$c2)), add=TRUE, col="red") #line

#Binh Thang (R)