// Examples in Kohler/Luniak, Data Inspection using Biplots

version 8.2
set matsize 800
set more off
		
// Interpretation
	
use planets, clear
input mypos
9
6
12
9
6
3
3
12
12
biplot logdist-logmoons ring, mlabel(planet) mlabvpos(mypos) cov gh flip(x) ms(oh) 
graph export introduction0.eps, replace
	
sysuse auto, clear
biplot mpg weight displacement price, subpop(for, msymbol(Oh Sh)) 
graph export introduction5.eps, replace


// Biplot and PCA-Plots

use iris, clear
sjlog using example1, replace
biplot sepallen-petalwid 
sjlog close, replace nolog
graph export example1.eps, replace

sjlog using example2, replace
biplot sepallen-petalwid, stretch(1) varonly 
sjlog close, replace nolog
graph export example2.eps, replace

sjlog using example2a, replace
biplot sepallen-petalwid, obsonly 
sjlog close, replace nolog
graph export example2a.eps, replace

sjlog using example2b, replace
biplot sepallen-petalwid, dim(3 4) 
sjlog close, replace nolog 
graph export example2b.eps, replace

// Types and Variations
	
sjlog using example3, replace
biplot sepallen-petalwid, gh cov  
sjlog close, replace nolog
graph export example3.eps, replace

sjlog using example3aa, replace
biplot sepallen-petalwid, gh 
sjlog close, replace nolog
graph export example3aa.eps, replace


sjlog using example3a, replace
biplot sepallen-petalwid, mixed(jk gh)  
sjlog close, replace nolog
graph export example3a.eps, replace

// Appearence

sjlog using example4, replace
biplot sepallen-petalwid, subpop(species, msymbol(Oh X Th)) legend(ring(0) pos(4)) 
sjlog close, replace nolog
graph export example4.eps, replace

exit
	
	
	Please report bugs and errors to kohler@wz-berlin.de