getwd(wd)
getwd()
logistic = read.csv("logistic.csv")
fitness = read.csv("fitness.csv")
co2 = read.csv("co2.csv")
co2
allomet = read.csv("allomet.csv")
allomet
milk = read.csv("milk.csv")
milk
milk
names(milk)
mothermilk
attach(milk)
nothermilk
mothermilk
daughtermilk
plot(mothermilk, daughtermilk)
b = cov(mothermilk, daughtermilk)
b
b = cov(mothermilk, daughtermilk)/var(mothermilk)
b
a = mean(daughermilk) - b*mean(mothermilk)
a = mean(daughtermilk) - b*mean(mothermilk)
a
b
abline(a = 54.41, b = .4576)
mean(mothermilk)
mean(daughtermilk)
res = lm(daughtermilk~mothermilk)
res
abline(lm(daughtermilk~mothermilk))
curve(54.41 + .4576*x, add=T)
curve(64.41 + .4576*x, add=T)
curve(44.41 + .4576*x, add=T)
anova(lm(daughtermilk~mothermilk))
R2 = 4438.2/(4438.2 + 24130.6)
R2
sqrt(R2)
correlation coefficient
r = cov(mothermilk, daughtermilk)/(sd(mothermilk)*sd(daughtermilk))
r
res = lm(mothermilk ~ daughtermilk)
res
anova(res)
allomet
names(allomet)
attach(allomet)
plot(mass~length)
plot(log(mass) ~ log(length))
x = log(length)
y = log(mass)
plot(x,y)
res = lm(y~x)
res
anova(res)
plot(log(mass) ~ log(length))
plot(mass~length)
a = exp(-1.712)
a
curve(0.1805*x^2.433, add=T)
.1805*21*2.433
.1805*21^2.433
.1805*21*2.433
nls()
nls(mass ~ a*length^b, start = (a = .2, b = 2))
nls(mass ~ a*length^b, start = (a = .2, b = 2))
nls(mass ~ a*length^b)
?nls
nls(mass ~ a*length^b, start = c(a = .2, b = 2))
res = nls(mass ~ a*length^b, start = c(a = .2, b = 2))
summary(res)
co2
attach(co2_
attach(co2)
plot(conc, uptake)
nls(uptake ~ vmax*conc/(Km + conc), start = c(vmax = 40, Km = 150))
curve(45.9*x/(112.6 + x), add=T)
res = nls(uptake ~ vmax*conc/(Km + conc), start = c(vmax = 40, Km = 150))
summary(res)
attach(logistic)
names(logistic)
plot(t, N)
nls(N~ N0*K*exp(r*t)/(K + N0*(exp(r*t)-1)), start = c(K = 20, r = .2, N0 = 1))
res = nls(N~ N0*K*exp(r*t)/(K + N0*(exp(r*t)-1)), start = c(K = 20, r = .2, N0 = 1))
summary(res)
N0 = .009
r = 0.803
K = 20.26
curve(N0*K*exp(r*x)/(K + N0*(exp(r*t)-1)) , add = T)
n0
N0
K
r
x
curve(N0*K*exp(r*x)/(K + N0*(exp(r*x)-1)) , add = T)
plot(t, N)
curve(N0*K*exp(r*x)/(K + N0*(exp(r*x)-1)) , add = T)
plot(t, log(N))
plot(t, N)
curve(N0*K*exp(r*x)/(K + N0*(exp(r*x)-1)) , add = T)
L = c(10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20)
L
M = c()
log(L)
x = log(L)
y = log(M)
lm(y~x)
M = c(1447, 2558, 3096, 3401, 4810, 5870, 6842, 8335, 9910, 12090, 13926)
M
plot(L, M)
x = log(L)
y = log(M)
lm(y~x_
lm(y~x)
a = exp(0.3204)
b = 3.0788
curve(a*x^b, add = T)
nls(M~a*L^b, start = c(a = .3, b = 2))
nls(M~a*L^b, start = c(a = 1, b = 2))
a
curve(1.405*x^3.071, add = T, col = "red")