# September 3, 2002 # more on linear classification by regression eg.1 = function(n1=100,n2=100) { nmax=10000 set.seed(468); x1 = rnorm(nmax); x2 = rnorm(nmax,5) y1=rep(0,nmax); y1j = jitter(y1); y2=1+y1; y2j = jitter(y2) plot(0,0,type="n",xlim=range(x1,x2), ylim=range(y1j,y2j),main=paste(n1,n2)) points(x1[1:n1],y1j[1:n1],pch=16) points(x2[1:n2],y2j[1:n2],pch=16) abline(lsfit(c(x1[1:n1],x2[1:n2]),c(y1[1:n1],y2[1:n2]))); browser() } eg.1(100,100) eg.1(1000,100) eg.1(10000,100) eg.1(100,1000) eg.1(100,10000)