ctb_sol_essay/_CHECK_Me.R at master · steveharbick/ctb_sol_essay · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
# check if test prediction distributions consistent with training

Response_id="C"
par(mfrow=c(3,2))
par(mar=c(1,1,3,1))
for (item in item_4valid) {
  load(paste("Working_files/",Version,"_",Response_id,"_",item,"_NNLS.RData",sep=""))
  load(paste("Working_files/",Version,"_",Response_id,"_H_",item,"_NNLS.RData",sep=""))
  plot(density(NNLS$yhatV),main=paste(item))
  lines(density(NNLS_H),col=2)
}

par(mfrow=c(3,2))
par(mar=c(1,1,3,1))
for (item in item_4valid) {
  load(paste("Working_files/",Version,"_",Response_id,"_",item,"_ADJ1.RData",sep=""))
  load(paste("Working_files/",Version,"_",Response_id,"_H_",item,"_ADJ1.RData",sep=""))
  plot(density(ADJ1$yhatV),main=paste(item))
  lines(density(ADJ1_H),col=2)
}