[R] Problem with merging two zoo objects

Megh Dal megh700004 at yahoo.com
Fri Oct 15 21:22:58 CEST 2010


Hi Gabor, please see the attached files which is in text format. I have opened them on excel then, used clipboard to load them into R. Still really unclear what to do.

Also can you please elaborate this term "index = list(1, 2), FUN = function(d, t) as.POSIXct(paste(d, t))" in your previous file? In help, it is given that:"If FUN is specified then read.zoo calls FUN with the index as the first argument". I really could not connect your syntax with help.

--- On Sat, 10/16/10, Gabor Grothendieck <ggrothendieck at gmail.com> wrote:

> From: Gabor Grothendieck <ggrothendieck at gmail.com>
> Subject: Re: [R] Problem with merging two zoo objects
> To: "Megh Dal" <megh700004 at yahoo.com>
> Cc: r-help at stat.math.ethz.ch
> Date: Saturday, October 16, 2010, 12:11 AM
> On Fri, Oct 15, 2010 at 2:20 PM, Megh
> Dal <megh700004 at yahoo.com>
> wrote:
> > Dear all, I have following 2 zoo objects. However when
> I try to merge those 2 objects into one, nothing is coming
> as intended. Please see below the objects as well as the
> merged object:
> >
> >
> >> dat11
> >                      V2   V3   V4   V5
> > 2010-10-15 13:43:54 73.8 73.8 73.8 73.8
> > 2010-10-15 13:44:15 73.8 73.8 73.8 73.8
> > 2010-10-15 13:45:51 73.8 73.8 73.8 73.8
> > 2010-10-15 13:46:21 73.8 73.8 73.8 73.8
> > 2010-10-15 13:47:27 73.8 73.8 73.8 73.8
> > 2010-10-15 13:47:54 73.8 73.8 73.8 73.8
> > 2010-10-15 13:49:51 73.7 73.7 73.7 73.7
> >> dat22
> >                      V2   V3   V4   V5
> > 2010-10-15 12:09:12 74.0 74.0 74.0 74.0
> > 2010-10-15 12:09:33 73.9 73.9 73.9 73.9
> > 2010-10-15 12:20:36 74.0 74.0 74.0 74.0
> > 2010-10-15 12:30:36 74.0 74.0 74.0 74.0
> > 2010-10-15 12:41:03 73.7 73.7 73.7 73.7
> >> merge(dat11, dat22)
> >                    V2.dat11 V3.dat11
> V4.dat11 V5.dat11 V2.dat22 V3.dat22 V4.dat22 V5.dat22
> > 2010-10-15 12:09:12       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 12:09:33       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 13:43:54       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 13:44:15       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 13:45:51       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 13:46:21       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 13:47:27       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 13:47:54       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > 2010-10-15 13:49:51       NA       NA      
> NA       NA       NA       NA       NA      
> NA
> > Warning messages:
> > 1: In MATCH(x, x) == seq_len(length(x)) :
> >  longer object length is not a multiple of shorter
> object length
> > 2: In MATCH(x, x) == seq_len(length(x)) :
> >  longer object length is not a multiple of shorter
> object length
> >
> > If somebody points me whether I went wrong, it would
> be really great.
> >
> 
> If I try it then it works properly so there is likely
> something wrong
> with your dat11 and dat22 objects.  If you provide the
> problem
> reproducibly one might be able to say more.
> 
> > Lines1 <- "Date Time
> V2   V3   V4   V5
> + 2010-10-15 13:43:54 73.8 73.8 73.8 73.8
> + 2010-10-15 13:44:15 73.8 73.8 73.8 73.8
> + 2010-10-15 13:45:51 73.8 73.8 73.8 73.8
> + 2010-10-15 13:46:21 73.8 73.8 73.8 73.8
> + 2010-10-15 13:47:27 73.8 73.8 73.8 73.8
> + 2010-10-15 13:47:54 73.8 73.8 73.8 73.8
> + 2010-10-15 13:49:51 73.7 73.7 73.7 73.7"
> >
> > Lines2 <- "Date Time
> V2   V3   V4   V5
> + 2010-10-15 12:09:12 74.0 74.0 74.0 74.0
> + 2010-10-15 12:09:33 73.9 73.9 73.9 73.9
> + 2010-10-15 12:20:36 74.0 74.0 74.0 74.0
> + 2010-10-15 12:30:36 74.0 74.0 74.0 74.0
> + 2010-10-15 12:41:03 73.7 73.7 73.7 73.7"
> >
> > library(zoo)
> > dat1 <- read.zoo(textConnection(Lines1), header =
> TRUE,
> + index = list(1, 2), FUN = function(d, t)
> as.POSIXct(paste(d, t)))
> Warning messages:
> 1: closing unused connection 8 (Lines2)
> 2: closing unused connection 7 (Lines1)
> 3: closing unused connection 5 (Lines2)
> 4: closing unused connection 4 (Lines1)
> 5: closing unused connection 3 (Lines2)
> > dat2 <- read.zoo(textConnection(Lines2), header =
> TRUE,
> + index = list(1, 2), FUN = function(d, t)
> as.POSIXct(paste(d, t)))
> > merge(dat1, dat2)
>                
>     V2.dat1 V3.dat1 V4.dat1 V5.dat1 V2.dat2
> V3.dat2
> V4.dat2 V5.dat2
> 2010-10-15 12:09:12      NA   
>   NA      NA     
> NA    74.0    74.0
> 74.0    74.0
> 2010-10-15 12:09:33      NA   
>   NA      NA     
> NA    73.9    73.9
> 73.9    73.9
> 2010-10-15 12:20:36      NA   
>   NA      NA     
> NA    74.0    74.0
> 74.0    74.0
> 2010-10-15 12:30:36      NA   
>   NA      NA     
> NA    74.0    74.0
> 74.0    74.0
> 2010-10-15 12:41:03      NA   
>   NA      NA     
> NA    73.7    73.7
> 73.7    73.7
> 2010-10-15 13:43:54    73.8   
> 73.8    73.8    73.8     
> NA      NA
>   NA      NA
> 2010-10-15 13:44:15    73.8   
> 73.8    73.8    73.8     
> NA      NA
>   NA      NA
> 2010-10-15 13:45:51    73.8   
> 73.8    73.8    73.8     
> NA      NA
>   NA      NA
> 2010-10-15 13:46:21    73.8   
> 73.8    73.8    73.8     
> NA      NA
>   NA      NA
> 2010-10-15 13:47:27    73.8   
> 73.8    73.8    73.8     
> NA      NA
>   NA      NA
> 2010-10-15 13:47:54    73.8   
> 73.8    73.8    73.8     
> NA      NA
>   NA      NA
> 2010-10-15 13:49:51    73.7   
> 73.7    73.7    73.7     
> NA      NA
>   NA      NA
> 
> -- 
> Statistics & Software Consulting
> GKX Group, GKX Associates Inc.
> tel: 1-877-GKX-GROUP
> email: ggrothendieck at gmail.com
>


      
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: dat1.txt
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20101015/c7bd3898/attachment.txt>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: dat2.txt
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20101015/c7bd3898/attachment-0001.txt>


More information about the R-help mailing list