From 375f29ab4306821c888fd9ef0637f9ab2879e375 Mon Sep 17 00:00:00 2001 From: Ben Green Date: Thu, 2 Jul 2015 00:41:49 -0400 Subject: in which the model works!! --- R Scripts/predict-victims-plots.R | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) (limited to 'R Scripts/predict-victims-plots.R') diff --git a/R Scripts/predict-victims-plots.R b/R Scripts/predict-victims-plots.R index 2ac62c8..553aa89 100644 --- a/R Scripts/predict-victims-plots.R +++ b/R Scripts/predict-victims-plots.R @@ -1,5 +1,5 @@ ##### Plot results -hist(correct_rank3,150,xlim=c(0,vcount(lcc)),col=rgb(0,0,1,1/8), +hist(correct_rank1,150,xlim=c(0,vcount(lcc)),col=rgb(0,0,1,1/8), xlab='Risk Ranking of Victims',main='') hist(correct_rank1,150,xlim=c(0,vcount(lcc)),col=rgb(1,0,1,1/8),add=T) hist(correct_rank2,150,xlim=c(0,vcount(lcc)),col=rgb(0,0,1,1/8),add=T) @@ -36,6 +36,7 @@ legend("topleft", inset=0.05, c("Demographics Model", "Cascade Model", "Combined Model"), fill=c(rgb(0,0,1,1/2),rgb(1,0,0,1/2),rgb(0,1,0,1/2))) box(which='plot') + par(new=T) counts = counts/(100/nvics) barplot(counts, -- cgit v1.2.3-70-g09d2