similar to: Prediction Intervals (reposting)

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

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

Full_Name: Ronnen Levinson Version: 1.5.1 OS: Mac OS 10.1 Submission from: (NULL) (12.232.201.92) The value returned by strptime behaves badly after arithmetic operation and/or combination. > 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

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. 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,

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. 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

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

(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

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

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 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

Dear R-Listers, I have spatial data on an equidistant rectangular grid, similar to topographic data. I know that there are quite a few R-packages or base functions that provide nice iso-contours plot, but I don't want a plot, just the smoothed isocontour line of ONE level (e.g. 10 mm). Data sets are large, so it would be preferable if the availability of regular grid data could be exploited,

Hi ALL, Can you please help me on how to determine the prediction error for ordinary kriging?Below are all the commands i used to generate the OK plot: rsa2 <- readShapeSpatial("residentialsa", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) x2 <- readShapeSpatial("ptna2",

Hola a todos Esta es mi primera pregunta en el grupo, y es sencilla pero me tiene atascado. Estoy tratando de cortar mi mapa de (poner limites en UTM) en un lugar definido como mi area de estudio (en este caso el sur de chile). Pero creo no estar usando bien la función CRS ponendo bien los limites requeridos. > study_area <- readRDS("CHL_adm0.rds") > study_area_UTM <-

Hello, when I tried to create 3D graph from attached file (at first: columns 1, 2, and 3 (longitude, latitude, altitude), then columns 1, 2 and 4 , then the same for column 5 etc) using Akima library for 3D interpolation (the data are neither equidistant nor regularly distributed), the console crashed by SIGABRT fault response. Do you know what is the main problem? Otherwise, can you advise me the

Hello, i wrote this function guete and now i want to plot it: but i get this error message. i hope someone can help me. Error in dim(robj) <- c(dX, dY) : dims [product 16] do not match the length of object [1] p_11=seq(0,0.3,0.1) p_12=seq(0.1,0.4,0.1) guete = function(p_11,p_12) { set.seed(1000) S_vek=matrix(0,nrow=N,ncol=1) for(i in 1:N) { X_0=rmultinom(q-1,size=1,prob=p_0)

Hi there. Just a little status update how that integerization is coming along. I am trying to limit myself to 32 bit arithmetics. That means not using any __int64 or long long datatypes at any point. I have now replaced all steps up to including the estimation of the LPC filter coefficients with integer code. That is about a quarter of the total work completed, I would say. One problem that i

This might not possible in R but I thought I would give it shot. I am have to set up a 40 x 40 cm grid of 181 points equidistant from each other. Is there any way to produce a graph with R that can do this for me? Actual sizes are unimportant as long it is to scale. Thanks -- View this message in context: http://www.nabble.com/Grid-building-in-R-tp18371874p18371874.html Sent from the R help