# R code for STA 6934 class 8/27/01
scan("data/ocbook.dat")
X <- scan("data/ocbook.dat")
X <- matrix(scan("data/ocbook.dat"),byrow=T,ncol=5)
X[45,1]
X[,1]
mean(X[,1]) 
apply(X,2,mean)
m <- apply(X,2,mean)
m
v <- apply(X,2,var)
v
var(X)
s <- sqrt(v)
s
1:5
plot(1:5,m)
c(0,100)
plot(1:5,m,ylim=c(0,100),xlab="Exam",ylab="Average")
lines(1:5,m)
read.table("data/ocbook.dat",col.names=c("ME","VE","AL","AN","ST"))
OC <- read.table("data/ocbook.dat",col.names=c("ME","VE","AL","AN","ST"))
boxplot(OC)
summary(OC)
pairs(OC)
cor(X)
score <- apply(X,1,mean)
cor(score,X)
var(score)
princomp(X)
library(mva)
princomp(X)
SD <- princomp(X)
names(SD)
SD$loadings
sqrt(var(sqrt(5)*score))
plot(1:5,SD$sdev^2)
p <- 100*SD$sdev^2/sum(SD$sdev^2)
p
plot(1:5,p,ylab="Percent variance",ylim=c(0,100),xlab="Component")
#
HP <- read.table("data/house.dat",col.names=c("Price","Size","Bed","Bath","New"))
HP
stem(Price)
attach(HP,1)
stem(Price)
range(Price)
summary(Price)
plot(Size,Price)
plot(Size,Price,ylim=c(0,350),xlim=c(0,4))
lm(Price~Size)
fit <- lm(Price~Size)
names(fit)
fit$coef
abline(fit$coef)
plot(Size,Price,type="n",ylim=c(0,350),xlim=c(0,4))
points(Size,Price,col=Bath+1)
plot(Size,Price,type="n",ylim=c(0,350),xlim=c(0,4))
points(Size,Price,pch=Bath+1)
plot(Size,Price,type="n",ylim=c(0,350),xlim=c(0,4))
points(Size,Price,pch=Bath+1,col=Bath+1)
source("code/house.R")