similar to: Date arithmetic fails (PR#1819)

Hello. I have an R program that calls gawk (GNU Awk 3.06 for Windows) from within pipe() to preprocess a large file before it is read into a data frame with read.table(). I've recently upgraded from Win98SE to WinXP, and have also upgraded from R1.5.0 to R1.6.1 over the past month or so. This program worked before the upgrade(s), but now fails. I observe the following sort of behavior with

(I'm reposting this message because the original has not appeared after about 2 days. Sorry if it shows up twice.) Hello. First, thanks to those who responded to my recent inquiry about using contour() over arbitrary (x,y) by mentioning the interp() function in the akima package. That worked nicely. Now for a new question: I would like to use a pair of prediction intervals to

Hi. Can anyone suggest a simple way to obtain in R a list of vector coordinates of the following form? The code below is Mathematica. In[5]:= Flatten[Table[{i,j,k},{i,3},{j,4},{k,5}], 2] Out[5]= {{1,1,1},{1,1,2},{1,1,3},{1,1,4},{1,1,5},{1,2,1},{1,2,2},{1,2,3},{1 ,2,4},{1,2, 5},{1,3,1},{1,3,2},{1,3,3},{1,3,4},{1,3,5},{1,4,1},{1,4,2},{1,4,3}, {1,4,

Hi. I'm having a little difficulty with R1.5.1. To wit: (1) Can anyone explain the following behavior? > a=strptime("2002-06-01 12:15:01","%Y-%m-%d %H:%M:%S") > a [1] "2002-06-01 12:15:01" > a+0 [1] "1932-04-25 21:46:45" > a+3600 [1] "1932-04-25 22:46:45" (2) Can anyone recommend an introductory R text other than that

Hello. I have an XYY series that I would like to graph with matplot() or some other single function that will do the trick. The X in question is a vector of POSIXt values obtained from strptime(). Is it possible to tell matplot() how to handle POSIXt x values? I have examined the examples at http://lark.cc.ukans.edu/~pauljohn/R/statsRus.html#5.22 , but would prefer not have to overlay the

Hi. The following statement works fine in R-win 1.8.0, but yields a syntax error in R-win 2.0.1 (and possibly in other versions after 1.8.0): plot(c(1,2),main=expression(a==b==c)) I note that the following workaround executes successfully in both versions of R... plot(c(1,2),main=expression(a*"="*b*"="*c)) ...but I don't really understand

Hello. Is it possible to set specify the font used by legend()? I would like to specify a fixed-width font so that I can line up parts of vertically stacked curve labels. For example, it would be nice if I could align the names, ages, and weights in the following three curve labels: Bob age=7 weight=100 Alexander age=13 weight=150 Susan age=20 weight=130 Is there perhaps a clever

Hi. As shown in the simplified example below, I'm having trouble allocating memory for large vectors, even though it would appear that there is more than enough memory available. That is, even with a memory limit of 1500 MB, R 2.6.1 (Win) will allocate memory for a first vector of 285 MB, but not for a second vector of the same size. Forcing garbage collection does not seem

(1) With MSVCRT.DLL (version 4.20.6164) from ftp://ftp.microsoft.com/softlib/mslfiles/msvcrt.exe installed in rw1051\bin directory, as per RWinFAQ 2.14: RGUI caused an invalid page fault in module MSVCRT.DLL at 017f:78014b90. Registers: EAX=0093007a CS=017f EIP=78014b90 EFLGS=00010246 EBX=017a2360 SS=0187 ESP=0091eee4 EBP=0091ef10 ECX=ffffffb5 DS=0187 ESI=00000000 FS=5bbf EDX=81d7d520 ES=0187

Hello. Is there an easy-to-use contour plot function analogous to scatterplot3d that can draw handle a dataset of arbitrary (x,y,z) triplets? That is, say x, y, and z are each measured quanties, and exhibit neither order nor regularity. I looked at the lattice package function "contourplot" but it seems complicated, and it's not clear from the documentation whether it can

I installed a printer, HP DeskJet 6122, which has no Ethernet device inside. In Windows, I use a local port with address \\darkstar\hp6122. The printer work, but I can see it only in Windows menu "Fax and Printers", not in the share resource of Samba, even if I set my smb.conf with "load printers" and, in [printers]. "browseable = yes" Is it all ok? Thanx M. --

Hi. I often wind up with many help windows cluttering my RGui screen when running Windows R 2.0.1. Is there an R instruction to close one or more help windows, or an RGui command to close all help windows? Yours truly, /Ronnen. /P.S. E-mailed CC:s of posted replies appreciated. [[alternative HTML version deleted]]

Hi. Is there an analog to abline() that can be used to plot a polynomial fit? For example, I can draw the straight-line fit fit <- lm(y ~ x) via abline(coef=fit$coef) but I'm not sure how to draw the polynomial fit fit <- lm(y ~ poly(x,2)) I do see the function curve(), but not how to prepare an expr for curve() based on the coefficients returned by the polynomial

Hello, I have read recent posts on this topic (Dr. Ronnen Levinson's Monday 02:39:55 pm), but before I install a 64 bit system, and purchase more RAM, I want to make sure I understand this interesting issue. I was attempting to plot a data frame containing Dow Jones stock information: > length(Date) [1] 19957 > length(Close) [1] 19957 > head(DowJones) Date Open High

On 3/12/16 12:33 AM, peter dalgaard wrote: >> On 12 Mar 2016, at 00:05 , Mick Jordan <mick.jordan at oracle.com> wrote: >> >> This is definitely obscure but we had a unit test that called .Internal(strptime, "1942/01/01", %Y/%m/%d") with timezone (TZ) set to CET. > Umm, that doesn't even parse. And fixing the typo, it doesn't run: > >>

I'm trying to figure out how to feed a blockaddress to a switch condition AND destination (basically emulating an indirectbr via a switch; I know it's not a good approach, I'm just experimenting). Suppose I have the following: SwitchInst *s = SwitchInst::Create(...); BasicBlock *bb = ...; PtrToIntInst k = new PtrToIntInst(BlockAddress::get(bb), <TYPE>, "", s);

Hello, and thanks in advance for any assistance. I have a linux machine that I'm trying to join to a windows 2003 sp1 active directory. The specifics are: RHEL5, samba version samba-3.0.23c-2.el5.2.0.2 a firewall between this server and the rest of the world (which includes the DCs), ports are open for kerberos and CIFS inbound and kerberos, CIFS, NTP and UDP oubtound. this machine

This is definitely obscure but we had a unit test that called .Internal(strptime, "1942/01/01", %Y/%m/%d") with timezone (TZ) set to CET. In R-3.1.3 that returned "1942-01-01 CEST" which, paradoxically, is correct as they evidently did strange things in Germany during the war period. Java also returns the same. However, R-3.2.4 returns "1942-01-01 CET".

>>>>> peter dalgaard <pdalgd at gmail.com> >>>>> on Sat, 12 Mar 2016 19:11:40 +0100 writes: > OK, .Internal is not necessary to reproduce oddity in this area. I also see things like (notice 1980) >> strptime(paste0(sample(1900:1999,80,replace=TRUE),"/01/01"), "%Y/%m/%d", tz="CET") ............... >

Hi everybody, somehow i dont get the shapiro wilk test for normality. i just can?t find what the H0 is . i tried : shapiro.test(rnorm(5000)) Shapiro-Wilk normality test data: rnorm(5000) W = 0.9997, p-value = 0.6205 If normality is the H0, the test says it?s probably not normal, doesn ?t it ? 5000 is the biggest n allowed by the test... are there any other test ? ( i know qqnorm