## defining experts' predictions ## you <- matrix(c(your_mean,your_sd,your_skewness),ncol=3) #enter your prediction by replacing the appropriate words for the appropriate numbers expert1 <- matrix(c(14500,1100,.8),ncol=3) expert2 <- matrix(c(15500,3000,1.8),ncol=3) expert3 <- matrix(c(19000,3500,1.1),ncol=3) #predictions by three fictive experts experts experts <- rbind(expert1,expert2,expert3,you) #sore all experts' predictions in one matrix with 1 row per expert #The means in column 1, the sd's in column 2 and the skewness in column 3 ## simulating new data ## set.seed(11235) Y2018 <- rep(NA,20) Y2018[c(1,2,6,7,9:20)] <- rnorm(16,17000,5000) #get data for the year 2018 ## package DAC ## install.packages("DAC") #install package library(DAC) #load package ## Calculating the DAC values ## set.seed(81325) out <- DAC_Uniform(data = na.omit(Y2018),bench.min = 0,bench.max = 30000,experts = experts,lb.area = 0,ub.area = 30000,step=.1) #calculate DAC values for you and the three fictive experts. #ender as data the simulated values from the vector Y2018 #the benchmark prior is set as a uniform prior between 0 and 30.000 by the input bench.min and bench.max #the area of the parameter space we choose to look at is set between 0 and 30.000 by the lb.area and ub.area #the precision by which we look is set at .1 by the step function. The smaller this is set the more accurate the DAC #scores become, yet the longer the calculation takes. ## ploting the experts' predictions and actual results ## x <- seq(from=0,to=30000,by=.1) #get the range for the x-axis of the plot #from is equal to lb.area in DAC_Uniform #to is equal to ub.area in DAC_Uniform #by is equal to step in DAC_Uniform densities <- cbind(dnorm(x,out\$posterior.distribution[1],out\$posterior.distribution[2]), dsnorm(x,experts[1,1],experts[1,2],experts[1,3]), dsnorm(x,experts[2,1],experts[2,2],experts[2,3]), dsnorm(x,experts[3,1],experts[3,2],experts[3,3]), dsnorm(x,experts[4,1],experts[4,2],experts[4,3])) #get a matrix with the density value specified at each x-axis location per row #column 1 is for the posterior, column 2 for expert 1, colomn 3 for expert 2 #column 4 for expert 3 and column 5 for your prediction. matplot(x,densities,type="l",yaxt="n",ylab="Density",xlab="",lwd=3) legend("topright",c("Posterior","Expert 1", "Expert 2","Expert 3", "You"),lty=1:5,col=1:5,lwd=2)