I think I have found the error. It appears when there are two entries with the same time. Using as input file:
--------- CUT -------- # Output format for PCKs: # TIME FLOW P [+-] SEQ TS X RTT SIZE # 123.125683 0 P + 967 123.125683 13394 0.798205 1472 123.241137 0 P + 968 123.241137 12680 0.796258 1472 123.241137 0 P + 969 123.241137 12680 0.796258 1472 123.472631 0 P + 970 123.472631 12680 0.796258 1472 123.588613 0 P + 971 123.588613 12680 0.796258 1472 123.704594 0 P + 972 123.704594 12680 0.796258 1472 --------- CUT -------- I run fhe following code: --------- CUT -------- h_types <- list (0, 0, NULL, NULL, 0, 0, 0, 0, 0) h_names <- list ("time", "flow", "seq", "ts", "x", "rtt", "size") pcks_file <- pipe ("grep ' P ' data", "r") pcks <- scan (pcks_file, what = h_types, comment.char = '#', fill = TRUE) mat_df <- data.frame (pcks[1:2], pcks[5:9]) mat <- as.matrix (mat_df) colnames (mat) <- h_names z <- zoo (mat, mat [,"time"]) --------- CUT -------- The dput of 'z' shows: --------- CUT -------- structure(c(123.125683, 123.241137, 123.241137, 123.472631, 123.588613, 123.704594, 0, 0, 0, 0, 0, 0, 967, 968, 969, 970, 971, 972, 123.125683, 123.241137, 123.241137, 123.472631, 123.588613, 123.704594, 13394, 12680, 12680, 12680, 12680, 12680, 0.798205, 0.796258, 0.796258, 0.796258, 0.796258, 0.796258, 1472, 1472, 1472, 1472, 1472, 1472 ), .Dim = c(6, 7), .Dimnames = list(c("1", "2", "3", "4", "5", "6"), c("time", "flow", "seq", "ts", "x", "rtt", "size")), index = structure(c(123.125683, 123.241137, 123.241137, 123.472631, 123.588613, 123.704594), .Names = c("1", "2", "3", "4", "5", "6")), class = "zoo") --------- CUT -------- If I try a 'as.ts(z)', it fails. If I remove the duplicate entry, I can convert it to a TS with no problem. Is this made intentionally? Because then I have to filter the input matrix... But, anyway, the output matrix, after filtering, doesn't seem regular: --------- CUT -------- > as.ts (z) Time Series: Start = 1 End = 5 Frequency = 1 time flow seq ts x rtt size 1 123.1257 0 967 123.1257 13394 0.798205 1472 2 123.2411 0 969 123.2411 12680 0.796258 1472 3 123.4726 0 970 123.4726 12680 0.796258 1472 4 123.5886 0 971 123.5886 12680 0.796258 1472 5 123.7046 0 972 123.7046 12680 0.796258 1472 Warning message: ‘x’ does not have an underlying regularity in: as.ts.zoo(z) --------- CUT -------- Weird... On 13 Dec 2005, at 16:33, Gabor Grothendieck wrote: > Please provide a reproducible example. Note that dput(x) will output > an R object in a way that can be copied and pasted into another > session. > > On 12/13/05, Alvaro Saurin <[EMAIL PROTECTED]> wrote: >> >> On 13 Dec 2005, at 13:08, Gabor Grothendieck wrote: >> >>> Your variable mat is not a matrix; its a data frame. Check it with: >>> >>> class(mat) >>> >>> Here is an example: >>> >>> x <- cbind(A = 1:4, B = 5:8) >>> tt <- c(1, 3:4, 6) >>> >>> library(zoo) >>> x.zoo <- zoo(x, tt) >>> x.ts <- as.ts(x.zoo) >> >> Fixed, but anyway it fails: >> >>> h_types <- list (0, 0, NULL, NULL, 0, 0, 0, 0, 0) >>> h_names <- list ("time", "flow", "seq", "ts", "x", "rtt", >>> "size") >> >>> pcks_file <- pipe ("grep ' P ' server.dat", "r") >>> pcks <- scan (pcks_file, what = h_types, >> comment.char = '#', fill = >> TRUE) >> >>> mat_df <- data.frame (pcks[1:2], pcks[5:9]) >>> mat <- as.matrix (mat_df) >>> colnames (mat) <- h_names >> >>> class (mat) >> [1] "matrix" >> >>> z <- zoo (mat, mat [,"time"]) >> >>> z >>> z >> time flow seq ts >> x rtt size >> 1.0009 1.000893 0.000000 0.000000 1.000893 >> 1472.000000 0.000000 1472.000000 >> 1.5145 1.514454 0.000000 1.000000 1.514454 >> 2944.000000 0.513142 1472.000000 >> 2.0151 2.015093 0.000000 2.000000 2.015093 >> 2944.000000 0.513142 1472.000000 >> 2.515 2.515025 0.000000 3.000000 2.515025 >> 4806.000000 0.504488 1472.000000 >> 2.822 2.821976 0.000000 4.000000 2.821976 >> 5730.000000 0.496728 1472.000000 >> [...] >> >>> as.ts (z) >> Error in if (del == 0 && to == 0) return(to) : >> missing value where TRUE/FALSE needed >> >> Any idea? Thanks for your help. >> >> Alvaro >> >> >> -- >> Alvaro Saurin <[EMAIL PROTECTED]> <[EMAIL PROTECTED]> >> >> >> >> -- Alvaro Saurin <[EMAIL PROTECTED]> <[EMAIL PROTECTED]> ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html