create script for comparison of kinship methods for challenge

[smcclatchy]

create a script of the following for people to copy and paste, then change lodcolumns to 2 for spleen and col 2 for spleen
create challenge and provide url to script for copy-paste

this is for liver

add comment

plot(out, map, lodcolumn=1,
col=color[1],
main=colnames(iron$pheno)[1],
ylim=c(0, ymx*1.02))

add comment

plot(out_pg, map, lodcolumn=1, col=color[2], add=TRUE)

add comment

plot(out_pg_loco, map, lodcolumn=1,
col=color[3], add=TRUE, lty=2)

add comment

legend("topleft", lwd=2, col=color,
c("H-K", "LMM", "LOCO"), bg="gray90", lty=c(1,1,2))

plot(out, lodcolumn="liver", map)
plot(out, lodcolumn="liver", map, col="black", add=TRUE)
plot(out_pg, lodcolumn="liver", map, col="red", add=TRUE)
plot(out_pg_loco, lodcolumn="liver", map, col="yellow", add=TRUE)