# ParametersI <-20# Number of individualsK <-5# Number of topicsD <-50# Number of diseasesT <-100# Number of time pointstheta_em=readRDS("thetaEM.rds")true_theta=melt(Theta_individual)colnames(true_theta)=c("Individual","Topic","Time","Theta")est_theta=melt(theta_em)colnames(est_theta)=c("Individual","Topic","Time","Theta")ptid=sample(I,1)ggplot(true_theta[true_theta$Individual%in%ptid,],aes(x=Time,y=Theta,fill=as.factor(Topic),group=as.factor(Topic)))+geom_area()+scale_fill_nejm()