[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