storks = read.table("storks.csv", sep=",")
storks
iq = read.table("IQdata.csv", sep=",")
milk = read.table("milk.csv", sep=",")
women = read.csv("womenhts.csv")
women
milk
corsim = function(sdx, sdy, slope, intercept, n) {
x = rnorm(n, mean=2000, sd = sdx)
y = intercept + slope*x + rnorm(n, sd = sdy)
sslope = cov(x, y)/var(x)
covar = cov(x, y)
corr = cor(x, y)
l = lm(y~x)
p = anova(l)[["Pr(>F)"]][[1]]
plot(x, y, xlab = "x", ylab = "y")
abline(l)
mtext(paste("r = ", round(corr, 3), "      Cov(x, y) = ", round(covar, 3), "        slope = ", round(sslope, 3), "   prob = ", round(p,3)))
#list(x = x, y = y)
} #end corsim
corsim(10, 10, 1, 0, 100)
corsim(10, 10, 1, 0, 100)
corsim(10, 10, 1, 0, 100)
corsim(10, 10, 1, 0, 100)
corsim(10, 10, 1, 0, 100)
corsim(10, 10, 1, 0, 100)
corsim(10, 10, 1, 0, 100)
corsim(10, 10, 1, 0, 100)
corsim(10, 10, -1, 0, 100)
corsim(10, 10, -1, 0, 100)
corsim(10, 10, 0, 0, 100)
corsim(10, 10, 0, 0, 100)
corsim(10, 10, 0, 0, 100)
corsim(10, 10, 0, 0, 100)
corsim(10, 10, 0, 0, 100)
corsim(10, 10, 0, 0, 100)
corsim(10, 10, 0, 0, 100)
corsim(10, 10, 0, 0, 100)
x = c(4, 5, 6, 7, 8, 5)
y = c(7, 9, 8, 6, 8, 9)
plot(x, y)
c = sum((x - mean(x))*(y - mean(y)))/(4)
c
cov(x, y)
c = sum((x - mean(x))*(y - mean(y)))/(3)
c
length(x)
c = sum((x - mean(x))*(y - mean(y)))/(5)
c
length(y)
x
y
slope = cov(x,y)/var(x)
slope
int = mean(y) - slope*mean(x)
int
lm(y~x)
abline(lm(y~x))
plot(x, y, xlim  = c(0,8))
abline(lm(y~x))
abline(a = 8, b = .3)
plot(x, y)
r = cov(x,y)/(sd(x)*sd(y))
r
cor(x,y)
abline(lm(y~x))
lm(y~x)
abline(a = 8, b = -.3)
abline(lm(y~x))
plot(x, y)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .2, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .6, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, .9, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 1.5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 5, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 50, 0, 100)
corsim(100, 100, 500, 0, 100)
corsim(100, 100, 500, 0, 100)
corsim(100, 100, 500, 0, 100)
corsim(100, 100, 500, 0, 10)
corsim(100, 100, 500, 0, 10)
corsim(100, 100, 500, 0, 10)
corsim(100, 100, 5, 0, 10)
corsim(100, 100, 5, 0, 10)
corsim(100, 100, 5, 0, 10)
corsim(100, 100, 5, 0, 10)
corsim(100, 100, 5, 0, 10)
corsim(100, 100, 5, 0, 10)
corsim(100, 100, 5, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 10)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
corsim(100, 100, 1, 0, 1000)
x
y
plot(x,y)
abline(lm(y~x))
l = lm(y~)
l = lm(y~x)
l
summary(l)
names(l)
l$residuals
sum((l$residuals)^2)
anova(lm(y~x))
storks
plot(storks$x, storks$y)
plot(storks$y, storks$x)
anova(lm(storks$x~storks$y))
cor(storks$x,storks$y)
iq
attach(iq)
names(iq)
plot(parentsIQ, offspringIQ)
cor(parentsIQ, offspringIQ)
abline(lm(offspringIQ~parentsIQ))
anoa(lm(offspringIQ~parentsIQ))
anova(lm(offspringIQ~parentsIQ))
womenhts
data(women)
attach(women)
names(women)
women
plot(height, weight)
abline(lm(weight~height))
anova(lm(weight~height))
cor(height,weight)
l$resid
x = rnorm(100)
y = rnorm(100)
l = lm(y~x)
cov(x,y)
cor(x,y)
plot(x,y)
l$resid
plot(l$resid)
plot(density(l$resid))
iq = read.csv("IQdata.csv")
iq
iq = read.csv("IQdata.csv")
attach(iq)
anova(lm(offspringIQ, parentsIQ))
anova(lm(offspringIQ ~parentsIQ))
x = c(20, 16, 19.8, 18.4, 17.1, 15.5, 14.7, 17.1, 15.4, 16.2, 15, 17.2, 16, 17, 14.4)
y = c(31.4, 22, 34.1, 29.1, 27, 24, 20.9, 27.8, 20.8, 28.5, 28.1, 27, 28.6, 24.6)
plot(x,y)
y = c(31.4, 22, 34.1, 29.1, 27, 24, 20.9, 27.8, 20.8, 28.5, 28.1, 27, 28.6, 24.6)
(y-32)*9/5
y - 32
y*9/5+32
28.5*9/5+32
79.6-32
47.6*5/9
length(y)
length(x)
y = c(31.4, 22, 34.1, 29.1, 27, 24, 20.9, 27.8, 20.8, 28.5, 26.4, 28.1, 27, 28.6, 24.6)