zz<-scan("golfout.txt")
zz<-matrix(zz,ncol=607,byrow=T)
x<-read.table("golfdataR.dat")
par(mfrow=c(1,3))
plot(density(zz[,(1:604)[unique(x[,1])=="VijaySingh"]]),lwd=3)
lines(density(zz[,(1:604)[unique(x[,1])=="TigerWoods"]]),col="red",lwd=3)
lines(density(zz[,(1:604)[unique(x[,1])=="PhilMickelson"]]),col="blue",lwd=3)
lines(density(zz[,(1:604)[unique(x[,1])=="WayneLevi"]]),col="orange",lwd=3)
lines(density(zz[,(1:604)[unique(x[,1])=="TimWilkinson"]]),col="purple",lwd=3)
library('MASS')
f1<-
kde2d(zz[,(1:604)[unique(x[,1])=="VijaySingh"]],zz[,(1:604)[unique(x[,1])=="TigerWoods"]])
persp(f1)
image(f1)
post.means<-apply(zz,2,mean)
unique(x[,1])[(1:604)[order(post.means[1:604],decreasing=TRUE)]]

par(mfrow=c(1,1))