Skip to content
Snippets Groups Projects
Commit abc93f85 authored by Finn Hughes's avatar Finn Hughes
Browse files

Update R code

parent 3ff6ab51
No related branches found
No related tags found
No related merge requests found
...@@ -220,4 +220,7 @@ plot(c((1:9)/10,1:15),c(sdsmallJi,sdtestJi),type="l",lwd=4,col=2,xlab="alpha",yl ...@@ -220,4 +220,7 @@ plot(c((1:9)/10,1:15),c(sdsmallJi,sdtestJi),type="l",lwd=4,col=2,xlab="alpha",yl
lines(c((1:9)/10,1:15),c(sdsmallb,sdtestb),lwd=4,lty=2,col="yellow") lines(c((1:9)/10,1:15),c(sdsmallb,sdtestb),lwd=4,lty=2,col="yellow")
abline(h=sd(yiJ),col=1) abline(h=sd(yiJ),col=1)
abline(h=sd(yiperf),col=4) abline(h=sd(yiperf),col=4)
abline(v=0,lty=2,col=3)
legend(x="right", legend=c("JCGM 101 + inf.","MC-VE inf.","JCGM 101 + No PK", "JCGM 101 + Perfect PK"),fill=c(2,"yellow",1,4)) legend(x="right", legend=c("JCGM 101 + inf.","MC-VE inf.","JCGM 101 + No PK", "JCGM 101 + Perfect PK"),fill=c(2,"yellow",1,4))
par(cex=1.25)
text(x=rep(5, 2), y=c(sd(yiJ)+0.005, sd(yiperf)-0.005), pos=4, labels=c(expression(paste(alpha %->% 0," (No PK)")), expression(paste(alpha %->% infinity," (Perfect PK)"))),col=c(1,4))
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment