similar to: hang-up with log="x" and segments (PR#1030)

X-Mailer: VM 7.17 under 21.4 (patch 18) "Social Property" XEmacs Lucid FCC: ~/mail/sent-mail From: Jan.Beyersmann at fdm.uni-freiburg.de --text follows this line-- hi everyone, please find below a job advert for a research position statistical computing in biostatistics, emphasis is on R and survival analysis. see also wangler et al., p. 31-35 in

Version 2.0-0 of meta (an R package for meta-analysis) is now available on CRAN. Changes are described below. Yours, Guido Major revision R package meta linked to R package metafor by Wolfgang Viechtbauer to provide additional statistical methods, e.g. meta-regression and other estimates for tau-squared (REML, ...) New functions: - metareg (meta-regression) - metabias

Version 2.0-0 of meta (an R package for meta-analysis) is now available on CRAN. Changes are described below. Yours, Guido Major revision R package meta linked to R package metafor by Wolfgang Viechtbauer to provide additional statistical methods, e.g. meta-regression and other estimates for tau-squared (REML, ...) New functions: - metareg (meta-regression) - metabias

Seriously, if a method gives a wrong result, it's wrong. line() does NOT implement the algorithm of Tukey, even not after the patch. We're not discussing Excel here, are we? The method of Tukey is rather clear, and it is NOT using the default quantile definition from the quantile function. Actually, it doesn't even use quantiles to define the groups. It just says that the groups

Hi, I have a plot and I would like to overlay a PNG image over it. I'm using the rasterImage function to do this, but the problem I'm facing is working out the coordinates of the upper right corner of the final image in user coordinates. That is I can place the image so the lower left is located at the bottom of the y-axis and the left end of the x-axis. Since my image is say 100px x

Hello dear R community. I just started playing with the snowfall package (a wrapper for the snow package), and found it very convenient. (See also this great website: http://www.imbi.uni-freiburg.de/parallel/ ) I was wondering if it is possible to connect snowfall with the foreach package (since it has some connection to snow). My final goal is to do some simple parallel simulations on my two

Here are another three problems with logarithmic scales: 1) segments() does not work with logarithmic scales. I suggest to change lines 962-973 in "plot.c": for (i = 0; i < n; i++) { if (FINITE(xt(x0[i%nx0])) && FINITE(yt(y0[i%ny0])) && FINITE(xt(x1[i%nx1])) && FINITE(yt(y1[i%ny1]))) { GP->col = INTEGER(col)[i % ncol];

Hi all, I would like to ask your help to add a color band (Ι am not sure regarding the right term, this color band at the right of the plot "describing" values with their corresponding color. For now I have only this code test<-matrix(data=runif(10000),nrow=100) plot(test,axes="FALSE") axis(1,at=c(0,1),labels=c("a","b")) # but  I would like to add

Hello~   Did any one have used the package 'orcutt' ?   I find that it can not work smoothly in a single variable regression. I use the example following, it function very well.   But when I regress "cons" on "price" (use the "reg1<-lm(cons~price+income+temp)") , then  use "reg11<-cochrane.orcutt(reg1) ". There is an error message “Error in

> source("http://www.bioconductor.org/biocLite.R") > biocLite(c("Ruuid")) Running getBioC version 0.1.6 with R version 2.3.1 Running biocinstall version 1.8.4 with R version 2.3.1 Your version of R requires version 1.8 of Bioconductor. trying URL 'http://bioconductor.org/packages/1.8/bioc/src/contrib/Ruuid_1.10.0.tar.gz' Content type 'application/x-gzip'

I want to draw a figure with several panels of unequal size, so i thought I would try using screen(). However, I can't figure out how to define the sizes as a matrix. I've tried this: split.screen(matrix(c(0,0.5,0,0.5, 0.5,1,0.5,1), byrow=F, ncol=4)) and a couple of variants on it, but get the same error: Error in par(.split.screens[[cur.screen]]) : invalid value specified

Hi everyone, ? I am trying to get a point of intersection between a polyline and a straight line ?.. and get the x and y coordinates of this point. For exemplification consider this: ? ? set.seed(123) ? k1 <-rnorm(100, mean=1.77, sd=3.33) ?k1 <- sort(k1) q1 <- rnorm(100, mean=2.37, sd=0.74) q1 <- sort(q1, decreasing = TRUE) plot(k1, q1, xlim <- c((min(k1)-5),

Currently I only have one machine backing up to another. I am backing up to a filesystem that is mounted as nosuid. The syncing process works perfectly except it fails on files that are suid. Example proc/ root/ tmp/ mknod tmp/mysql.sock : Invalid argument mknod tmp/orbit-daud/orb-12573780511489051058 : Invalid argument mknod tmp/orbit-daud/orb-139425434385535813 : Invalid argument mknod

The control flow was: if (!y) { ppix = ... } if (y) { ... } else if (x) { use ppix for something } else { use ppix for something } Merge the if(!y) block with the two else branches. This avoids a false-positive in the clang static analyzer, it can't know that !y and x are mutually exclusive. The result looks something like this: if (y) { ... } else { ppix = ... if (x) {

Hi R-users, I am using arima to fit a time series. Now I would like to include an intervention component "It (0 before intervention, 1 after)" using different types of impacts, that is, not only trying the simple abrupt permanent impact (yt = w It ) with the xreg option but also trying with a gradual permanent impact (yt= d * yt-1 + w * It ), following the filosophy of Box and Tiao

[I'm refiling this one due to mailer config messup] I get an error whenever I try to plot a histogram in pictex. The error doesn't occur for X11 or postscript. > x<-c(1,2,3) > pictex() > hist(x) Error: Bad units specified in GConvert, please report! platform i686-unknown-linux arch i686 os linux

Dear R-users, I would like to maximize the function g above which depends on 4 parameters (2 vectors, 1 real number, and 1 matrix) using optim() and BFGS method. Here is my code: # fonction to maximize g=function(x) { x1 = x[1:ncol(X)] x2 = x[(ncol(X)+1)] x3 = matrix(x[(ncol(X)+2):(ncol(X)+1+ncol(X)*ncol(Y))],nrow=ncol(X),ncol=ncol(Y)) x4 = x[(ncol(X)+1+ncol(X)*ncol(Y)+1):length(x)]

Dear all, How can I label/unlabel points on scatterplot with mouse pointer. As the mouse approches near to point, it should label the closest point, then unlabel when it moves away. How can I do in R? I be very thankful. Thanks and Regards, Hitesh Singla -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

Hi, Is there a way to set the scale of a plot (i.e. number of axis units per centimeter) when you output it to postscript? If not, how am I supposed to plot graphs with different axis limits to the same scale? They just get resized to fit the paper so that graphs which show a smaller number of axis units end up with a larger scale. Cheers, Ben -- Dr. Ben Kenward Department of Psychology,

I'm having trouble creating a looping variable and i can't see wher ethe problem arises from any hep gratfully appreciated First create a table x<-table(SURVEY$n_0,exposed) > x exposed False True Under 16 24 1 16-19 68 9 20-24 190 37 25-34 555 204 35-44 330 87 45-54 198 65 55-64 67 35 65+