[R] How can we creat conditional data frame
Hi,
R experts. I am a new user of R and trying to learn this program.
I have a problem. Here is the code.
d<-as.Date(c("2000/01/03","2000/01/05","2000/01/19","2000/01/28"))
r<-rnorm(4)
da<-data.frame(d,r)
a<-as.Date("01/01/2000","%d/%m/%Y")
b<-as.Date("30/01/2000","%d/%m/%Y")
ab<-seq(a,b,by=1)
c<-format(ab,"%a")
date<-data.frame(ab,c)
date<-subset(date,c!="Sun")
date<-subset(date,c!="Sat")
Here I have 2 data frame.
da
-
d r
1 2000-01-03 1.2105865
2 2000-01-05 -0.8962776
3 2000-01-19 -1.0438936
4 2000-01-28 2.1329387
---
date
--
ab c
3 2000-01-03 Mon
4 2000-01-04 Tue
5 2000-01-05 Wed
6 2000-01-06 Thu
7 2000-01-07 Fri
10 2000-01-10 Mon
11 2000-01-11 Tue
12 2000-01-12 Wed
13 2000-01-13 Thu
14 2000-01-14 Fri
17 2000-01-17 Mon
18 2000-01-18 Tue
19 2000-01-19 Wed
20 2000-01-20 Thu
21 2000-01-21 Fri
24 2000-01-24 Mon
25 2000-01-25 Tue
26 2000-01-26 Wed
27 2000-01-27 Thu
28 2000-01-28 Fri
---
In data frame-"DA"- I have return(r) and date
In data frame-"Date"- where I have date and day.
Now I need to create a data frame where returns will be conditional on date
and rest will be zero.
---
like this
-
abc rhd
3 2000-01-03 Mon1.2105865 0
4 2000-01-04 Tue 0 1
5 2000-01-05 Wed0 1
6 2000-01-06 Thu 0 1
7 2000-01-07 Fri 0 1
-
>From this I can figureout the hoildays and then put each holidays equal to
1.
Please help me.
Thanking you
saikat
--
View this message in context:
http://www.nabble.com/How-can-we-creat-conditional-data-frame-tp16491208p16491208.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] How can we creat conditional data frame
Hi,
R experts. I am a new user of R and trying to learn this program.
I have a problem. Here is the code.
d<-as.Date(c("2000/01/03","2000/01/05","2000/01/19","2000/01/28"))
r<-rnorm(4)
da<-data.frame(d,r)
a<-as.Date("01/01/2000","%d/%m/%Y")
b<-as.Date("30/01/2000","%d/%m/%Y")
ab<-seq(a,b,by=1)
c<-format(ab,"%a")
date<-data.frame(ab,c)
date<-subset(date,c!="Sun")
date<-subset(date,c!="Sat")
Here I have 2 data frame.
da
-
d r
1 2000-01-03 1.2105865
2 2000-01-05 -0.8962776
3 2000-01-19 -1.0438936
4 2000-01-28 2.1329387
---
date
--
ab c
3 2000-01-03 Mon
4 2000-01-04 Tue
5 2000-01-05 Wed
6 2000-01-06 Thu
7 2000-01-07 Fri
10 2000-01-10 Mon
11 2000-01-11 Tue
12 2000-01-12 Wed
13 2000-01-13 Thu
14 2000-01-14 Fri
17 2000-01-17 Mon
18 2000-01-18 Tue
19 2000-01-19 Wed
20 2000-01-20 Thu
21 2000-01-21 Fri
24 2000-01-24 Mon
25 2000-01-25 Tue
26 2000-01-26 Wed
27 2000-01-27 Thu
28 2000-01-28 Fri
---
In data frame-"DA"- I have return(r) and date
In data frame-"Date"- where I have date and day.
Now I need to create a data frame where returns will be conditional on date
and rest will be zero.
---
like this
-
abc rhd
3 2000-01-03 Mon1.2105865 0
4 2000-01-04 Tue 0 1
5 2000-01-05 Wed0 1
6 2000-01-06 Thu 0 1
7 2000-01-07 Fri 0 1
-
>From this I can figureout the hoildays and then put each holidays equal to
1.
Please help me.
Thanking you
saikat
--
View this message in context:
http://www.nabble.com/How-can-we-creat-conditional-data-frame-tp16491207p16491207.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Re: [R] How can we creat conditional data frame
Here is the solution
data<-read.table("d:/ftse.txt",header=TRUE)
wDay<-as.Date(data$date,"%d/%m/%Y")
data<-data.frame(data[,-1])
price<-data$ftse100
data<-data.frame(price,wDay)
data[1:10,]
a<-as.Date("03/01/1989","%d/%m/%Y")
b<-as.Date("31/01/2007","%d/%m/%Y")
ab<-seq(a,b,by=1)
c<-format(ab,"%a")
date<-data.frame(ab,c)
date<-subset(date,c!="Sun")
date<-subset(date,c!="Sat")
res <- merge(date, data, by.x = "ab", by.y="wDay", all=T)
res$h <- 0
res$h[is.na(res$price)] <- 1
res$price[is.na(res$price)] <- 0
h<-res$h[-1]
h1<-c(h,0)
final.data<- data.frame(res$price,h1)
final.data<-final.data[final.data$res.price>0,]
saikat sarkar wrote:
>
> Hi,
>
> R experts. I am a new user of R and trying to learn this program.
>
> I have a problem. Here is the code.
>
> d<-as.Date(c("2000/01/03","2000/01/05","2000/01/19","2000/01/28"))
> r<-rnorm(4)
> da<-data.frame(d,r)
>
> a<-as.Date("01/01/2000","%d/%m/%Y")
> b<-as.Date("30/01/2000","%d/%m/%Y")
> ab<-seq(a,b,by=1)
> c<-format(ab,"%a")
> date<-data.frame(ab,c)
> date<-subset(date,c!="Sun")
> date<-subset(date,c!="Sat")
>
> Here I have 2 data frame.
>
> da
> -
>d r
> 1 2000-01-03 1.2105865
> 2 2000-01-05 -0.8962776
> 3 2000-01-19 -1.0438936
> 4 2000-01-28 2.1329387
>
> ---
> date
> --
> ab c
> 3 2000-01-03 Mon
> 4 2000-01-04 Tue
> 5 2000-01-05 Wed
> 6 2000-01-06 Thu
> 7 2000-01-07 Fri
> 10 2000-01-10 Mon
> 11 2000-01-11 Tue
> 12 2000-01-12 Wed
> 13 2000-01-13 Thu
> 14 2000-01-14 Fri
> 17 2000-01-17 Mon
> 18 2000-01-18 Tue
> 19 2000-01-19 Wed
> 20 2000-01-20 Thu
> 21 2000-01-21 Fri
> 24 2000-01-24 Mon
> 25 2000-01-25 Tue
> 26 2000-01-26 Wed
> 27 2000-01-27 Thu
> 28 2000-01-28 Fri
>
> ---
>
> In data frame-"DA"- I have return(r) and date
> In data frame-"Date"- where I have date and day.
>
> Now I need to create a data frame where returns will be conditional on
> date and rest will be zero.
>
> ---
> like this
> -
>
> abc rhd
> 3 2000-01-03 Mon1.2105865 0
> 4 2000-01-04 Tue 0 1
> 5 2000-01-05 Wed0 1
> 6 2000-01-06 Thu 0 1
> 7 2000-01-07 Fri 0 1
>
> -
>
> From this I can figureout the hoildays and then put each holidays equal to
> 1.
>
>
> Please help me.
>
> Thanking you
>
> saikat
>
>
>
>
--
View this message in context:
http://www.nabble.com/How-can-we-creat-conditional-data-frame-tp16491208p16495146.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Re: [R] how to fit GJR-GARCH model in R
If I understand correct then
garch.gjr <- function(par,y,iterate=TRUE)
{
T<-length(y)
mu0<-par[1]
a0<-par[2]
a1<-par[3]
a2 <-par[4]
b1 <- par[5]
e <- s2 <- numeric(T)
s2[1] <- var(y)
e[1] <- y[1]-mu0
l<-0
for(t in 2:T)
{
e[t] <- y[t]-mu0
s2[t] <- a0+a1*e[t-1]^2+b1*s2[t-1] + ifelse(e[t-1]<0, a2*e[t-1]^2, 0)
l <- l -0.5*log(2*pi*s2[t])-0.5*e[t]^2/s2[t]
if(s2[t]>1) return(1)
}
if(iterate) return(-l)
else return(list(loglik=l,sig2=s2,res=e/sqrt(s2)))
}
y<-data
par0<-c(0,0.038,0.051,0.058,0.991)
# you need to put good initial parameter values
a<-nlm(garch.gjr,par0,y=y,hessian=TRUE)
se <- sqrt(diag(solve(a$hessian)))
Hope this might help you
Fan Steven wrote:
>
> Hi,All,
>
> I am trying to fit a GJR-GARCH model in R:
>
> r_t = mu + e_t
> h_t = alp_0 + alp_1 * e_(t-1)^2 + alp_2 * s_(t-1) * e_(t-1)^2 + beta *
> h_(t-1)
>
> where r_t = return (on day t), h_t = conditional volatility on day t,
> and s_(t-1) = 1 if e_(t-1) < 0 (zero otherwise).
>
> I have downloaded the packages "tseries" and "fSeries" but can not see
> how to fit this model.
>
> Any help would be very much appreciated.
>
> Thanks,
>
> Steven
>
> _
> 享用世界上最大的电子邮件系统― MSN Hotmail。 http://www.hotmail.com
>
>
> __
> [EMAIL PROTECTED] mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
--
View this message in context:
http://www.nabble.com/Re%3A-how-to-fit-GJR-GARCH-model-in-R-tp4207387p16498300.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] how to manupute data frame with conditions fill cell with previous value if next cell is zero
Dear R Experts, This is the 2nd time in the chat room. Its a great place to get help from R experts. I have a data frame problem, it contains thousands of data. part of it, I am giving for explaining the problem date day x yz 82 1989-04-28 Fri 2118.0 2418.80 33713 83 1989-05-01 Mon0.0 2414.96 33793 84 1989-05-02 Tue 2103.12402.86 33955 85 1989-05-03 Wed 2105.7 2393.70 0 86 1989-05-04 Thu 2119.0 2384.90 0 87 1989-05-05 Fri2132.8 2381.96 0 103 1989-05-29 Mon0.00.00 34161 127 1989-06-30 Fri 2151.0 2440.06 32949 128 1989-07-03 Mon 2165.6 2452.77 33236 129 1989-07-04 Tue 2174.40 35376 130 1989-07-05 Wed 2162.9 2456.56 33310 167 1989-08-25 Fri 2397.4 2732.36 34740 168 1989-08-28 Mon0.02743.36 34607 169 1989-08-29 Tue 2380.8 2726.63 34688 170 1989-08-30 Wed 2381.3 2728.15 34472 171 1989-08-31 Thu 2387.9 2737.27 34431 172 1989-09-01 Fri 2407.5 2752.09 34348 173 1989-09-04 Mon 2419.20.00 34484 174 1989-09-05 Tue 2426.0 2744.68 34442 175 1989-09-06 Wed 2390.8 2719.79 34271 176 1989-09-07 Thu 2415.9 2706.88 34153 177 1989-09-08 Fri 2423.9 2709.54 34116 178 1989-09-11 Mon 2400.6 2704.41 34114 179 1989-09-12 Tue 2397.6 2707.26 34333 180 1989-09-13 Wed 2401.5 2679.52 34287 181 1989-09-14 Thu 2382.0 2664.89 34402 182 1989-09-15 Fri 2366.5 2674.58 0 183 1989-09-18 Mon 2373.8 2687.50 34473 195 1989-10-04 Wed 2312.1 2771.09 35383 196 1989-10-05 Thu 2281.6 2773.56 35523 197 1989-10-06 Fri 2277.5 2785.52 35209 198 1989-10-09 Mon 2247.0 2791.41 35376 199 1989-10-10 Tue 2218.8 2785.33 0 200 1989-10-11 Wed 2218.8 2773.36 35240 In this data frame, I need to replace cells which is zero with previous value. For example 82 1989-04-28 Fri 2118.0 2418.80 33713 83 1989-05-01 Mon0.0 2414.96 33793 in 2nd line 0.0 should be replaced by 2118.0 Another Example 84 1989-05-02 Tue 2103.1 2402.86 33955 85 1989-05-03 Wed 2105.7 2393.70 0 86 1989-05-04 Thu 2119.0 2384.90 0 87 1989-05-05 Fri 2132.8 2381.96 0 3 lines filled by 33955 value. Another example 198 1989-10-09 Mon 2247.0 2791.41 35376 199 1989-10-10 Tue 2218.8 2785.33 0 200 1989-10-11 Wed 2218.8 2773.36 35240 in 2nd line filled by zero should be replaced by value 35376 Please help me. Thanking you saikat -- View this message in context: http://www.nabble.com/how-to-manupute-data-frame-with-conditions-fill-cell-with-previous-value-if-next-cell-is-zero-tp16524107p16524107.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Need some help in R programming code
Dear R guru,
I am Saikat Sarkar working as a researcher of Economics in Tampere
University, Finland. I am trying to estimate some Garch related tests with
Bayesian analysis by R programme.
I am not good in R but trying to survive.
Anyway I have the coding but not working properly. I have tried to find the
problem but failed. I am writing to all R gurus to help me out.
Could you please look at the problem and help me if you can.
Thanking you
saikat
Could you please look at the problem below:
R message is like below:
---
Error in garch.gjr.d0(param, y, X, Z, iterate = FALSE)$res :
$ operator is invalid for atomic vectors
plot(Ch.mic,Ch.mic.rep,xlim=c(0,6),ylim=c(0,6))
---
Coding
res<-garch.gjr.d0(a$estimate,y,X,X,iterate=FALSE)$res
sig2<-garch.gjr.d0(a$estimate,y,X,X,iterate=FALSE)$sig2
test.mic <- function(param,y,X,Z)
{
n <- length(y)
res <- garch.gjr.d0(param,y,X,Z,iterate=FALSE)$res
s2 <- garch.gjr.d0(param,y,X,Z,iterate=FALSE)$sig2
s<-sqrt(s2)
r <- res
u <- r[-1]
r2<- (r[-n]^2*s2[-n])/s[-1]
r3<- (r[-n]^3*s2[-n]*s[-n])/s[-1]
U <- cbind(r[-n],r2,r3)
U <- cbind(1,U)
U <- U/matrix(s[-1],nr=nrow(U),nc=ncol(U))
Fstat<-summary(lm(u~U-1))$fstatistic[1]
Fstat
}
sim.gjr <- function(par,y,X=0,Z=0)
{
n <- length(y)
e<-numeric(n)
yrep<-numeric(n)
yrep[1] <- y[1]
n.dummies.mean <- dim(X)[2]
if(is.matrix(X)) e[1] <- e[1]-sum(X[1,]*par[9:(8+n.dummies.mean)])
delta <- par[-(1:(8+n.dummies.mean))]
s2 <- var(y)
for(i in 2:n){
s2 <-
par[5]+par[6]*e[i-1]^2+par[7]*ifelse(e[i-1]<0,1,0)*e[i-1]^2+par[8]*s2+sum(X[i,]*delta)
e[i] <- rnorm(1,0,sqrt(s2))
yrep[i]<-par[1]+par[2]*yrep[i-1]+par[3]*s2+par[4]*yrep[i-1]*s2+e[i]
}
if(is.matrix(X)) yrep <- yrep + X%*%par[9:(8+n.dummies.mean)]
yrep
}
yrep <- sim.gjr(a$estimate,y,X,Z)
n.test <- 500
Ch.mic <- Ch.mic.rep <- numeric(n.test)
for(i in 1:n.test){
simulation <- sample(1:1,1)
par <- sims.matrix[simulation,]
Ch.mic[i] <- test.mic(par,y,X,X)
yrep <- sim.gjr(par,y,X,X)
Ch.mic.rep[i] <- test.mic(par,yrep,X,X)
}
plot(Ch.mic,Ch.mic.rep,xlim=c(0,6),ylim=c(0,6))
lines(c(0,6),c(0,6))
--
View this message in context:
http://www.nabble.com/Need-some-help-in-R-programming-code-tp20634187p20634187.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[R] why this function give error message
Dear R guru, I am Saikat Sarkar working as a researcher of Economics in Tampere University, Finland. I am trying to estimate some Garch related tests with Bayesian analysis by R programme. I am not good in R but trying to survive. Anyway I have the coding but not working properly. I have tried to find the problem but failed. I am writing to all R gurus to help me out. Could you please look at the problem and help me if you can. Thanking you saikat Could you please look at the problem a<-c(1297.4, 1309.3 , 1321.9 , 1316.6 , 1309.5 , 1266.2 , 1285.2 , 1272.8 , 1264.1 , 1277.5 , 1293.3 , 1298.7 , 1295.3 , 1276.6 , 1277.5 , 1277.5 , 1277.5 , 1269.5 , 1253, 1257.1 , 1238.3 , 1240.6 , 1248.5 , 1223.1 , 1219.7 , 1203.2 , 1180.3 , 1170.7 , 1195.2 , 1199.4 , 1209.9 , 1214.9 , 1232.1 , 1243.6 , 1243.2 , 1250.7 , 1263.2 , 1257.9 , 1244.8 , 1240.9 , 1250.7 , 1273.9 , 1285.9 , 1293.4 , 1298.5 , 1324.4 , 1316.4 , 1341, 1359.2 , 1354.3 , 1362, 1359.1 , 1357.7 , 1382.5 , 1414, 1432.6 , 1420.8 , 1427.2 , 1450.3 , 1425.1 , 1380.8 , 1396, 1416.5 , 1415.6 , 1415.6 , 1415.6 , 1412.5 , 1407.6 , 1412.7 , 1424.6 , 1453.8 , 1447.9 , 1462.6 , 1467.6 , 1436.9 , 1450.9 , 1452.4 , 1430.1 , 1439.1 , 1436.6 , 1436.6 , 1440.9 , 1443.3 , 1449.9 , 1457.2 , 1456, 1445.1 , 1443.3 , 1426.3 , 1423.4 , 1410.8 , 1399.9 , 1406.3 , 1389.1 , 1403.1 , 1425.3 , 1438.7 , 1453.4 , 1442.8 , 1456.2 , 1461.5 , 1470, 1490.5 , 1513.8 , 1539, 1540, 1566.9 , 1586.1 , 1564.1 , 1581.7 , 1598.4 , 1595.2 , 1591.3 , 1590.4 , 1586.6 , 1586.6 , 1569.5 , 1587.5 , 1599.5 , 1600.5 , 1618.2 , 1602, 1612.9 , 1600.8 , 1574.4 , 1567.2 , 1544.6 , 1560.8 , 1554.9 , 1558.9 , 1560.1 , 1572.3 , 1589.7 , 1599.8 , 1612.6 , 1629.1 , 1625.9 , 1620.4 , 1616.2 , 1623.6 , 1642.5 , 1623.4 , 1622, 1633.2 , 1636.2 , 1630.2 , 1619.6 , 1612.7 , 1600.1 , 1611.7 , 1624.6 , 1623.8 , 1622.7 , 1631.8 , 1641.4 , 1644, 1657.1 , 1627.4 , 1616.8 , 1623.5 , 1624.9 , 1600.1 , 1608.1 , 1631.1 , 1640.6 , 1643.1 , 1629.5 , 1632.4 , 1636.2 , 1595.6 , 1569, 1580.5 , 1570.5 , 1541.4 , 1546.2 , 1571.9 , 1567.2 , 1559.8 , 1555.5 , 1545.5 , 1517.2 , 1530.6 , 1549.4 , 1567.9 , 1559.9 , 1539.5 , 1546.7 , 1543.6 , 1541.7 , 1551.1 , 1543.3 , 1545, 1546.2 , 1541.9 , 1551.6 , 1551.6 , 1523.6 , 1511.9 , 1516.3 , 1525.9 , 1534.6 , 1543.5 , 1544.4 , 1540.8 , 1551.6 , 1566.9 , 1571.7 , 1579.4 , 1591.1 , 1595.3 , 1587.4 , 1580.6 , 1582.1 , 1583.6 , 1581.5 , 1587, 1584.9 , 1580, 1571.8 , 1570.4 , 1557.4 , 1562.7 , 1564.4 , 1564.1 , 1558.2 , 1526.9 , 1522.6 , 1522.3 , 1494.9 , 1493.3 , 1483, 1478.4 , 1490.7 , 1498.9 , 1506, 1480.1 , 1472, 1471.5 , 1456.5 , 1451, 1447.2 , 1459.3 , 1471.5 , 1475.7 , 1471, 1469.4 , 1457.5 , 1457.7 , 1463.3 , 1447.3 , 1451.6 , 1462.1 , 1480, 1484.4 , 1484.1 , 1484.1 , 1484.1 , 1486.4 , 1486.4 , 1487.4 , 1487.4 , 1480.8 , 1470.4 , 1473.7 , 1478.3 , 1489.3 , 1501.8 , 1499.4 , 1507.8 , 1517.9 , 1512, 1501.9 , 1502.4 , 1512, 1517.9 , 1523.2 , 1525.3 , 1536.5 , 1539.4 , 1542.6 , 1542.2 , 1551.5 , 1531.3 , 1510.5 , 1511.7 , 1501.5 , 1492, 1501.3 , 1512.3 , 1509, 1478.6 , 1485.7 , 1479, 1451, 1448.1 , 1473.1 , 1481.4 , 1492.7 , 1478.9 , 1486.2 , 1470, 1484.8 , 1495.9 , 1494.7 , 1508.7 , 1512.9 , 1506.9 , 1491.2 , 1486.5 , 1488, 1491, 1493.3 , 1491.1 , 1497.5 , 1496.3 , 1484.4 , 1483.6 , 1488.7 , 1486, 1486, 1486, 1486.7 , 1479.7 , 1467.2 , 1459.6 , 1452, 1433, 1426, 1424.5 , 1412.9 , 1413.3 , 1414.7 , 1421.9 , 1421.4 , 1416.3 , 1419.5 , 1422.2 , 1449.8 , 1458.2 , 1452.5 , 1475.9 , 1475.9 , 1472.4 , 1479.9 , 1500.2 , 1491.2 , 1492.6 , 1504.6 , 1511.9 , 1510.8 , 1512.6 , 1501.8 , 1503.3 , 1515.5 , 1536.7 , 1565.8 , 1557.3 , 1562.1 , 1560.6 , 1562.6 , 1579.6 , 1583.7 , 1571.1 , 1544, 1530.8 , 1535.9 , 1543.2 , 1529.9 , 1531.2 , 1524.9 , 1535.5 , 1523.3 , 1534, 1544.8 , 1543.9 , 1543.9 , 1541.2 , 1533.5 , 1532.6 , 1528.4 , 1526.7 , 1536.5 , 1528.2 , 1513.6 , 1510.6 , 1509.4 , 1508.5 , 1519, 1524.7 , 1521.1 , 1501.3 , 1504.7 , 1513.2 , 1519.5 , 1527.5 , 1537.3 , 1553.6 , 1560.6 , 1556.7 , 1558.1 , 1563.3 , 1565, 1571.2 , 1596.4 , 1603, 1607.6 , 1606.5 , 1617.7 , 1628.2 , 1626.9 , 1634.2 , 1638.6 , 1649.3 , 1664.7 , 1670.9 , 1680, 1699, 1695.3 , 1681.9 , 1694.9 , 1715.5 , 1713.4 , 1719.5 , 1726.7 , 1716.9 , 1722.5 , 1724.6 , 1738, 1747.4 , 1765.2 , 1765.2 , 1781.8 , 1781.2 , 1763.2 , 1748.7 , 1763.9 , 1778, 1776.7 , 1751.1 , 1758.2 , 1756.3 , 1741.3 , 1734.2 , 1733, 1733.8 , 1735.7 , 1745.8 , 1
