-
Notifications
You must be signed in to change notification settings - Fork 2
/
.Rhistory
36 lines (36 loc) · 817 Bytes
/
.Rhistory
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
data(ple4)
library(FLa4a)
data(ple4)
data(ple4.indices)
rec(ple4)
plot(rec(ple4))
plot(rec(ple4), type='l')
yrs <- dimnames(rec(ple4))[[2]]
yrs
yrs <- as.numeric(dimnames(rec(ple4))[[2]])
yrs
rec <- c(rec(ple4))
yrs
rec
plot(rec~yrs)
plot(rec~yrs, type='l')
plot(log(rec)~yrs, type='l')
lm(rec~yrs)
predict(lm(rec~yrs))
lines(predict(lm(rec~yrs))~yrs, col=2)
plot(rec~yrs, type='l')
lines(predict(lm(rec~yrs))~yrs, col=2)
fit <- sca(ple4, ple4.indices, srmodel=~year)
rec(fit)
stock.n(fit)[1]
lines(c(stock.n(fit)[1])~yrs)
rec <- c(rec(ple4))
year <- as.numeric(dimnames(rec(ple4))[[2]])
rlm <- lm(rec~year)
scafit <- sca(ple4, ple4.indices, srmodel=~year)
plot(rec~year, type='l')
lines(predict(rlm)~year, col=2)
lines(stock.n(scafit)[1]~year, col=3)
stock.n(scafit)[1]
year
lines(c(stock.n(scafit)[1])~year, col=3)