ls()
x = matrix( scan( "pollen.dat"), byrow=T, ncol=6 )[,-6]
ind = seq(nrow(x))[ read.table("tmp.colors") != "Default" ]
system(pwd)
system("pwd")
source("/Users/scottdw/courses/stat541/yr08/Examples/Nov.25/xgobi.R")
xb = round(x + runif(3848*5,-.005,.005), 3)
col = rep(1,nrow(x)); col[ind]=2
dim(col)
ind
is.numeric(xb)
xgobi.R(cbind(xb,col))
q()
P = eigen(var(x[ind,]))$vec# put EURAKA in x-y plane
me = apply(x[ind,],2,"mean"); ma = apply(x,2,"mean")
xP = t( t(P) %*% ( t(x)+(me-ma) ) )
round(eigen(x)$va,3)
round(eigen(var(x))$va,3)
round(eigen(var(xp))$va,3)
round(eigen(var(xP))$va,3)
xgobi.R(xP)
hist(dis1,100,F,col=2); lines(tmp$mids,tmp$density)
dis1 = sqrt(c(x^2 %*% rep(1,5))); hist(dis1,100)
sigma = var(x[col==1,])
ev = eigen(sigma)
sigma5 = ev$vec %*% diag(sqrt(ev$val)) %*% t(ev$vec)
n=1000000
tst = matrix(rnorm(5*n),n,5) %*% sigma5
dis2 = sqrt(c(tst^2 %*% rep(1,5))); tmp=hist(dis2,100)
hist(dis1,100,F,col=2); lines(tmp$mids,tmp$density)
xgobi(cbind(xP,dis1))
xgobi.R(cbind(xP,dis1))
q()
ls()
rm(dis2)
rm(tst)
q()
ls()
dim(x)
xgobi.R(x)
ev
print(ev)
ndigits
par("ndigits")
options("ndigits")
options("digits")
help.start()
help("options")
getOptions(digits)
print(ev)
print(ev,digit=3)
  hist(dis1,100)
table(col)
var(x)
print(var(x),dig=3)
print(sigma,dig=3)
table(ind)
length(ind)
tmp.colors
ls()
print(n)
hist(dis1,100,F,col=2); lines(tmp$mids,tmp$density)
dis3 = sqrt(c(scale(x,T,F)^2 %*% rep(1,5))); hist(dis3,100)
print(var(tst),dig=3)
ls()
rm(dis3)
xgobi.R( cbind(x,dis2) )## brush
xgobi.R( cbind(x,dis1) )## brush
q()
ls()
str(dis1)
str(x)
str(xb)
range(x-xb)
apply(x,2,"mean")
dis0 = sqrt( x^2 %*% rep(1,5) )
str(dis0)
hist(dis0)
hist(dis0,250,xlim=c(0,1))
hist(dis0,250,xlim=c(0,2))
str(col)
table(col)
x2
x2 = x[col==2,]
str(x2)
dis2 = sqrt( x2^2 %*% rep(1,5) )
dis2
range(dis2)
hist(dis2,seq(0,3,.2),add=T)
dev.new()
hist(dis2,seq(0,3,.2),add=T)
hist(dis2)
dev.off()
hist(dis0,250,xlim=c(0,4))
hist(dis0,250,xlim=c(0,6))
i3 = seq(dis0)[ dis0<3]
length(i3)
pairs(x[i3,]
)
cov(x)
cor(x)
hist(dis0)
hist(dis0,100)
hist(dis0,200)
system("pwd")
source("pollen.q")
source("pollen.q")
str(xP)
pairs(xP)
pairs(xP[col,)
ls()
str(col)
table(col)
pairs(xP[col==2,])
pairs(x[col==2,])
source("pollen.q")
source("pollen.q")
c
q()
ls()
xgobi.R
ls()
rm(xgobi.R)
q()
