target <- read.table("../lrn/num/472.dat",header=T,colClasses="numeric")
y <- target[,1] 

state <- read.table("../lrn/chr/4.dat",
           header=T,colClasses="character",blank.lines.skip=F)   

f <- state[,1]

f <- as.factor(f)

fit <- lm(y~f-1)

print(anova.lm(fit))
print(summary.lm(fit))

c <- fit$coefficients
s <- sort(c,index.return=TRUE)
c <- c[s$ix]
print(c)

f <- state[,1]

f[f=="AS"|f=="DC"|f=="DE"|f=="MA"|f=="ME"|f=="NH"] <- "S1" 
f[f=="OH"|f=="RI"|f=="VI"|f=="WV"]                 <- "S1" 

f <- as.factor(f)

fit <- lm(y~f-1)

print(anova.lm(fit))
print(summary.lm(fit))

c <- fit$coefficients
s <- sort(c,index.return=TRUE)
c <- c[s$ix]
print(c)

f <- state[,1]

f[f=="AS"|f=="DC"|f=="DE"|f=="MA"|f=="ME"|f=="NH"] <- "S1" 
f[f=="OH"|f=="RI"|f=="VI"|f=="WV"]                 <- "S1" 
f[f=="AA"|f=="AE"|f=="AP"|f=="CT"|f=="GU"|f=="MD"] <- "S2"
f[f=="NJ"|f=="NY"|f=="PA"|f=="PA"|f=="VA"|f=="VT"] <- "S2"

f <- as.factor(f)

fit <- lm(y~f-1)

print(anova.lm(fit))
print(summary.lm(fit))

c <- fit$coefficients
s <- sort(c,index.return=TRUE)
c <- c[s$ix]
print(c)

f <- state[,1]

f[f=="AS"|f=="DC"|f=="DE"|f=="MA"|f=="ME"|f=="NH"] <- "S1" 
f[f=="OH"|f=="RI"|f=="VI"|f=="WV"]                 <- "S1" 
f[f=="AA"|f=="AE"|f=="AP"|f=="CT"|f=="GU"|f=="MD"] <- "S2"
f[f=="NJ"|f=="NY"|f=="PA"|f=="PA"|f=="VA"|f=="VT"] <- "S2"
f[f=="AK"|f=="UT"|f=="MS"]                         <- "S3"

f <- as.factor(f)

fit <- lm(y~f-1)

print(anova.lm(fit))
print(summary.lm(fit))

c <- fit$coefficients
s <- sort(c,index.return=TRUE)
c <- c[s$ix]
print(c)

f <- state[,1]

f[f=="AS"|f=="DC"|f=="DE"|f=="MA"|f=="ME"|f=="NH"] <- "S1" 
f[f=="OH"|f=="RI"|f=="VI"|f=="WV"]                 <- "S1" 
f[f=="AA"|f=="AE"|f=="AP"|f=="CT"|f=="GU"|f=="MD"] <- "S2"
f[f=="NJ"|f=="NY"|f=="PA"|f=="PA"|f=="VA"|f=="VT"] <- "S2"
f[f=="WY"]                                         <- "S2"  
f[f=="AK"|f=="UT"|f=="MS"]                         <- "S3"
f[f=="NE"|f=="ND"]                                 <- "S4"
f[f=="SD"|f=="SC"]                                 <- "S5"

f <- as.factor(f)

fit <- lm(y~f-1)

print(anova.lm(fit))
print(summary.lm(fit))

c <- fit$coefficients
s <- sort(c,index.return=TRUE)
c <- c[s$ix]
print(c)