# nov11

setwd("/home/amy/Desktop")

source("ExtractLinesFromLargeFileFunction.R")

result = ExtractLinesFromLargeFileFunction("StataGraph.ps",10000,10000)

LENGTH=unlist(result[3]) ; LENGTH ; summary(LENGTH)

n.l = length(LENGTH) ;  n.l

X = result[1][[1]] ; Y = result[2][[1]] ; length(X)

x.min=unlist(result[4]) ; y.min=unlist(result[5])
x.max=unlist(result[6]) ; y.max=unlist(result[7])
c(x.min,y.min,x.max,y.max)

plot(c(x.min,x.max),c(y.min,y.max),col="white", 
 xlim=c(x.min,x.max),ylim=c(y.min,y.max) )
for(i in 1:length(X) ) lines(X[[i]],Y[[i]] )

# dev.new()


x.range=c(0,400) ; y.range = c(200,450)
plot(x.range,y.range,col="white", xlim=x.range, ylim=y.range)
for(i in 1:length(X) ) lines(X[[i]],Y[[i]] )

# focus on tick marks 
plot(x.range,y.range,col="white", xlim=x.range, ylim=y.range)
for(i in 10:29 ) {
	L = LENGTH[i]	
	x = mean(X[[i]]) ; vx = var(X[[i]])
	y = mean(Y[[i]]) ; vy = var(Y[[i]])  
	lines(X[[i]],Y[[i]])
	if(vx>0|vy>0) points(c(x),c(y), cex=.4,pch=19,col="blue")
	if(vx==0 & vy==0) points(c(x),c(y), cex=.6,pch=19,col="red")
}


library(survival)
fit=survfit(Surv(time,status) ~ 1, data=aml[1:11,])
#summary(fit)
dev.new();
plot(survfit(Surv(time,status) ~ 1, data=aml[1:11,],type="kap"),F)