Below is the format that was requested. This has the data followed by the corrected code at the end. There are several things that were wrong:
1. z1[,2] is wrong since z1 is a vector, not a 2d matrix. Ditto for z2. Ideally zoo would have given an error message but in any case its wrong. It should be just z1. 2. The merge works ok but the second argument to na.approx should be xout=time(z1). xout= was missing. 3. there were a number of bugs removed from na.approx in the devel version of zoo. I don't think they impact this but if you have any problems with na.approx then uncomment the source statement in the code below. It will bring in the development version of na.approx into your workspace. Lines1 <- "time1,dataset1 01:01:00,0.73512097 01:01:30,0.34860813 01:02:00,0.61306418 01:02:30,0.01495898 01:03:00,0.27035612 01:03:30,0.69513898 01:04:00,0.46451758 01:04:30,0.61672569 01:05:00,0.82496122 01:05:30,0.34766154 01:06:00,0.69618714 01:06:30,0.39035214 01:07:00,0.01680143 01:07:30,0.28576967 01:08:00,0.01205416 01:08:30,0.89637254 01:09:00,0.63147653 01:09:30,0.01522139 01:10:00,0.27661960 01:10:30,0.50974124 01:11:00,0.68141977 01:11:30,0.90725854 01:12:00,0.83823443 01:12:30,0.53360241 01:13:00,0.17769196 01:13:30,0.83438616 01:14:00,0.67248807 01:14:30,0.09991933 01:15:00,0.03334966 01:15:30,0.93292355 01:16:00,0.15990837 01:16:30,0.05354050 01:17:00,0.55281203 01:17:30,0.37845690 01:18:00,0.89051365 01:18:30,0.16674292 01:19:00,0.85458626 01:19:30,0.19278550 01:20:00,0.73240405 01:20:30,0.16417524 01:21:00,0.73878212 01:21:30,0.51790118 01:22:00,0.83076438 01:22:30,0.47555504 01:23:00,0.02108640 01:23:30,0.82911053 01:24:00,0.95555646 01:24:30,0.14493657 01:25:00,0.84422332 01:25:30,0.41589974 01:26:00,0.67606367 01:26:30,0.00606434 01:27:00,0.59951991 01:27:30,0.43949260 01:28:00,0.66297385 01:28:30,0.33131298 01:29:00,0.06102041 01:29:30,0.84722118 01:30:00,0.46841491 01:30:30,0.34200755 01:31:00,0.87386578 01:31:30,0.70737403 01:32:00,0.23978781 01:32:30,0.11787278 01:33:00,0.14679814 01:33:30,0.65217063 01:34:00,0.81355908 01:34:30,0.31583482 01:35:00,0.92167666 01:35:30,0.55931271 01:36:00,0.13641271 01:36:30,0.35048575 01:37:00,0.17243584 01:37:30,0.93645686 01:38:00,0.85356548 01:38:30,0.61399352 01:39:00,0.05910707 01:39:30,0.01721605 01:40:00,0.94845557 01:40:30,0.48117810 01:41:00,0.34752402 01:41:30,0.59295472 01:42:00,0.64267429 01:42:30,0.57859933 01:43:00,0.00201441 01:43:30,0.32530995 01:44:00,0.25474645 01:44:30,0.93187534 01:45:00,0.99361033 01:45:30,0.16591641" Lines2 <- "time2,dataset2 01:01:01,0.17558467 01:01:42,0.23806514 01:02:23,0.75867726 01:03:06,0.73502357 01:03:50,0.94018206 01:04:35,0.61882643 01:05:21,0.68417492 01:06:08,0.05744461 01:06:55,0.33344394 01:07:44,0.68752593 01:08:33,0.17270469 01:09:23,0.81522124 01:10:03,0.68304352 01:10:43,0.38774082 01:11:23,0.84176890 01:12:04,0.09344446 01:12:44,0.13431965 01:13:25,0.92210721 01:14:06,0.33630635 01:14:47,0.56690294 01:15:29,0.09870816 01:16:11,0.77864105 01:16:53,0.61803441 01:17:35,0.09133728 01:18:17,0.08925487 01:19:00,0.89271117 01:19:42,0.56605742 01:20:25,0.98520534 01:21:08,0.66104843 01:21:51,0.96948589 01:22:34,0.05692690 01:23:17,0.71887456 01:24:00,0.14903741 01:24:43,0.86569445 01:25:26,0.27923513 01:26:09,0.98365033 01:26:53,0.08308399 01:27:36,0.87071027 01:28:19,0.26475705 01:29:03,0.76409811 01:29:47,0.59563256 01:30:31,0.23995054 01:31:14,0.00951054 01:31:59,0.21367270" Lines3 <- "time1,dataset1 01:02:30,0.01495898 01:03:00,0.27035612 01:03:30,0.69513898 01:04:00,0.46451758 01:04:30,0.61672569 01:05:00,0.82496122 01:05:30,0.34766154 01:06:00,0.69618714 01:06:30,0.39035214 01:07:00,0.01680143 01:07:30,0.28576967 01:08:00,0.01205416 01:08:30,0.89637254 01:09:00,0.63147653 01:09:30,0.01522139 01:10:00,0.27661960 01:10:30,0.50974124 01:11:00,0.68141977 01:11:30,0.90725854 01:12:00,0.83823443 01:12:30,0.53360241 01:13:00,0.17769196 01:13:30,0.83438616 01:14:00,0.67248807 01:14:30,0.09991933 01:15:00,0.03334966 01:15:30,0.93292355 01:16:00,0.15990837 01:16:30,0.05354050 01:17:00,0.55281203 01:17:30,0.37845690 01:18:00,0.89051365 01:18:30,0.16674292 01:19:00,0.85458626 01:19:30,0.19278550 01:20:00,0.73240405 01:20:30,0.16417524 01:21:00,0.73878212 01:21:30,0.51790118 01:22:00,0.83076438 01:22:30,0.47555504 01:23:00,0.02108640 01:23:30,0.82911053 01:24:00,0.95555646 01:24:30,0.14493657 01:25:00,0.84422332 01:25:30,0.41589974 01:26:00,0.67606367 01:26:30,0.00606434 01:27:00,0.59951991 01:27:30,0.43949260 01:28:00,0.66297385 01:28:30,0.33131298 01:29:00,0.06102041 01:29:30,0.84722118 01:30:00,0.46841491 01:30:30,0.34200755 01:31:00,0.87386578 01:31:30,0.70737403 01:32:00,0.23978781 01:32:30,0.11787278 01:33:00,0.14679814 01:33:30,0.65217063 01:34:00,0.81355908 01:34:30,0.31583482 01:35:00,0.92167666 01:35:30,0.55931271 01:36:00,0.13641271 01:36:30,0.35048575 01:37:00,0.17243584 01:37:30,0.93645686 01:38:00,0.85356548 01:38:30,0.61399352 01:39:00,0.05910707 01:39:30,0.01721605 01:40:00,0.94845557 01:40:30,0.48117810 01:41:00,0.34752402 01:41:30,0.59295472 01:42:00,0.64267429 01:42:30,0.57859933 01:43:00,0.00201441 01:43:30,0.32530995 01:44:00,0.25474645 01:44:30,0.93187534 01:45:00,0.99361033 01:45:30,0.16591641" Lines4 <- "time2,dataset2 01:03:06,0.73502357 01:03:50,0.94018206 01:04:35,0.61882643 01:05:21,0.68417492 01:06:08,0.05744461 01:06:55,0.33344394 01:07:44,0.68752593 01:08:33,0.17270469 01:09:23,0.81522124 01:10:03,0.68304352 01:10:43,0.38774082 01:11:23,0.84176890 01:12:04,0.09344446 01:12:44,0.13431965 01:13:25,0.92210721 01:14:06,0.33630635 01:14:47,0.56690294 01:15:29,0.09870816 01:16:11,0.77864105 01:16:53,0.61803441 01:17:35,0.09133728 01:18:17,0.08925487 01:19:00,0.89271117 01:19:42,0.56605742 01:20:25,0.98520534 01:21:08,0.66104843 01:21:51,0.96948589 01:22:34,0.05692690 01:23:17,0.71887456 01:24:00,0.14903741 01:24:43,0.86569445 01:25:26,0.27923513 01:26:09,0.98365033 01:26:53,0.08308399 01:27:36,0.87071027 01:28:19,0.26475705 01:29:03,0.76409811 01:29:47,0.59563256 01:30:31,0.23995054 01:31:14,0.00951054 01:31:59,0.21367270" library(chron) library(zoo) z1<-read.zoo(textConnection(Lines1),header=TRUE,sep=",",FUN=times) z2<-read.zoo(textConnection(Lines2),header=TRUE,sep=",",FUN=times) z12 <- merge(z1, z2) # source("http://r-forge.r-project.org/plugins/scmsvn/viewcvs.php/*checkout*/pkg/zoo/R/na.approx.R?rev=676&root=zoo") z3<-na.approx(z12,xout=time(z1)) On Fri, Apr 2, 2010 at 12:19 PM, e-letter <inp...@gmail.com> wrote: > On 02/04/2010, Gabor Grothendieck <ggrothendi...@gmail.com> wrote: >> The code does not run with the files. I need the requested >> information, namely a single file containing code and data and that I >> can just copy into a session without editing and see the result you >> see. > > I don't understand how I can combine the four csv files into a single > file of data and terminal commands? Anyway, further terminal output. > > The following also occurs with correction of the commands. The data > merge is incomplete (to 1:27:30); data set 1 ends at time 1:45:30; > data set 2 1:31:59 > >> library(chron) >> library(zoo) >> z1<-read.zoo("test1.csv",header=TRUE,sep=",",FUN=times) >> z2<-read.zoo("test2.csv",header=TRUE,sep=",",FUN=times) >> z3<-(na.approx(merge(z1[,2],z2[,2]),time(z1))) >> z3 > z1[, 2] z2[, 2] > 01:01:01 0.54186455 0.17558467 > 01:01:30 0.34860813 0.20682491 > 01:01:42 0.48083615 0.23806514 > 01:02:00 0.61306418 0.49837120 > 01:02:23 0.31401158 0.75867726 > 01:02:30 0.01495898 0.75079270 > 01:03:00 0.27035612 0.74290813 > 01:03:06 0.48274755 0.73502357 > 01:03:30 0.69513898 0.83760282 > 01:03:50 0.57982828 0.94018206 > 01:04:00 0.46451758 0.83306352 > 01:04:30 0.61672569 0.72594497 > 01:04:35 0.72084346 0.61882643 > 01:05:00 0.82496122 0.65150068 > 01:05:21 0.58631138 0.68417492 > 01:05:30 0.34766154 0.47526482 > 01:06:00 0.69618714 0.26635471 > 01:06:08 0.54326964 0.05744461 > 01:06:30 0.39035214 0.19544428 > 01:06:55 0.20357679 0.33344394 > 01:07:00 0.01680143 0.45147127 > 01:07:30 0.28576967 0.56949860 > 01:07:44 0.14891191 0.68752593 > 01:08:00 0.01205416 0.51591885 > 01:08:30 0.89637254 0.34431177 > 01:08:33 0.76392454 0.17270469 > 01:09:00 0.63147653 0.49396296 > 01:09:23 0.32334896 0.81522124 > 01:09:30 0.01522139 0.77116200 > 01:10:00 0.27661960 0.72710276 > 01:10:03 0.39318042 0.68304352 > 01:10:30 0.50974124 0.53539217 > 01:10:43 0.59558051 0.38774082 > 01:11:00 0.68141977 0.61475486 > 01:11:23 0.79433915 0.84176890 > 01:11:30 0.90725854 0.59232742 > 01:12:00 0.83823443 0.34288594 > 01:12:04 0.68591842 0.09344446 > 01:12:30 0.53360241 0.11388206 > 01:12:44 0.35564718 0.13431965 > 01:13:00 0.17769196 0.52821343 > 01:13:25 0.50603906 0.92210721 > 01:13:30 0.83438616 0.72684026 > 01:14:00 0.67248807 0.53157330 > 01:14:06 0.38620370 0.33630635 > 01:14:30 0.09991933 0.45160464 > 01:14:47 0.06663450 0.56690294 > 01:15:00 0.03334966 0.33280555 > 01:15:29 0.48313660 0.09870816 > 01:15:30 0.93292355 0.32535246 > 01:16:00 0.15990837 0.55199675 > 01:16:11 0.10672443 0.77864105 > 01:16:30 0.05354050 0.69833773 > 01:16:53 0.30317627 0.61803441 > 01:17:00 0.55281203 0.44246870 > 01:17:30 0.37845690 0.26690299 > 01:17:35 0.63448528 0.09133728 > 01:18:00 0.89051365 0.09029608 > 01:18:17 0.52862829 0.08925487 > 01:18:30 0.16674292 0.49098302 > 01:19:00 0.85458626 0.89271117 > 01:19:30 0.19278550 0.72938430 > 01:19:42 0.46259477 0.56605742 > 01:20:00 0.73240405 0.77563138 > 01:20:25 0.44828965 0.98520534 > 01:20:30 0.16417524 0.87715304 > 01:21:00 0.73878212 0.76910073 > 01:21:08 0.62834165 0.66104843 > 01:21:30 0.51790118 0.81526716 > 01:21:51 0.67433278 0.96948589 > 01:22:00 0.83076438 0.66529956 > 01:22:30 0.47555504 0.36111323 > 01:22:34 0.24832072 0.05692690 > 01:23:00 0.02108640 0.38790073 > 01:23:17 0.42509847 0.71887456 > 01:23:30 0.82911053 0.43395599 > 01:24:00 0.95555646 0.14903741 > 01:24:30 0.14493657 0.50736593 > 01:24:43 0.49457995 0.86569445 > 01:25:00 0.84422332 0.57246479 > 01:25:26 0.63006153 0.27923513 > 01:25:30 0.41589974 0.51404020 > 01:26:00 0.67606367 0.74884526 > 01:26:09 0.34106401 0.98365033 > 01:26:30 0.00606434 0.53336716 > 01:26:53 0.30279212 0.08308399 > 01:27:00 0.59951991 0.34562608 > 01:27:30 0.43949260 0.60816818 >> > ______________________________________________ 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.