similar to: is that a bug

Hello, I'm trying to plot a fitted lm() line on a plot when the one explanatory variable is log transformed and log="x". I get different lines using abline and predict.lm(). #Example x <- 1:100 y <- rnorm(100) plot(y ~ x, log="x") abline(lm(y ~ log(x))) lines(x, predict(lm(y ~ log(x))), lwd=2) I'm sure I'm missing something but could someone tell me which

Hello, I use t.test for normal distributed and wilcox.test for non-normal distributed samples. It is easy to write a function for t.test that calculates the effect size, because all parts of the formula are available from the t.test result: r = sqrt(t*t / (t*t + df)) However, for Wilcoxon tests, the formula for effect sizes is: r = Z / sqrt(N) I wonder how I can calculate the Z-score in R for

Hi, I just started with R and I found that there are many options to rearrange the data to do mixed models. I want to use the reshape function. I have 2 between subject variables and one within. I was able to change the data structure but still - the result of the aov functions are calculating everything as a within subject. the table looks like this: SerialNo breed treatment distance_1

Is this a bug? > f1 <- factor(c("a", NA), levels = c("a", "NA") ) > f2 <- f1[, drop = TRUE] > f2 [1] a <NA> Levels: a <NA> I would have expected f2 to have only one level, "a". It seems to me that the code in [.factor does not follow the advice in help("factor") on how to set factor codes to be

Full_Name: Stefan Raberger Version: 2.8.1 OS: Windows XP Submission from: (NULL) (213.185.163.242) Hi there, I recently noticed some strange behaviour of the command "type.convert", depending on the startup mode used. But there also seems to be different behaviour on different PCs (all running the same OS and the same version of R). On PC1: When I start R in SDI mode (RGui --no-save

Ive been converting elements of matrices and arrays, e.g., from Integers to double-precision, by vectorizing the matrix and then remaking it. Alternatively, I can redefine one element as double which then redefines them all. Both methods are quick, so I guess I shouldn't complain, but I would have thought there'd be something more obvious. Have I missed it? Here's my redimensioning

Dear all, Is there a less cumbersome way to rename a column by name (as opposed to index) than -- names( X)[ names[ X] == "bob"]<-"sue" ? A semi-related question: how does one get the index of a column by name, something along the lines of col.index( X, "sue") ? Chip Barnaby --------------------------------------------------------- Chip Barnaby

Hi I had observed the following problem in R (also C, Matlab, and Python). sprintf('%1.2g\n', 3.15) give 3.1 instead of 3.2 whereas an input of 3.75 gives 3.8. Java's System.out.printf is ok though. > round(3.75,1) [1] 3.8 > round(3.15,1) [1] 3.1 Similar outcome with sprintf in R. However, the right answer should be 3.2 Regards Teckpor [[alternative HTML version

Hi, The following function works, but is there a neater way to write it? f = function(x,...) { # return a character vector of the arguments passed in after 'x' gsub(" ","",unlist(strsplit(deparse(substitute(list(...))),"[(,)]")))[-1] } > f(x,a,b,c*d) [1] "a" "b" "c*d" > Thanks. [[alternative HTML

Full_Name: Richard Reed Version: 2.6.2 OS: MS Windows XP Pro Submission from: (NULL) (76.121.170.35) Pressing Alt-M+s doesn't search path. It just generates another ">". It does this in "vanilla" too. However, Alt-M+click "list search path" does work. Note that when I press Alt-M (and the "Misc" menu drops down) the second "s" in

i get a message: configure: WARNING: you cannot build info or html versions of the R manuals what should i do ? [[alternative HTML version deleted]]

R 2.8.1 is packaged up and queued up for testing. Please test this out: Fedora 9: https://admin.fedoraproject.org/updates/rpy-1.0.3-6.fc9,R-2.8.1-1.fc9 Fedora 10: https://admin.fedoraproject.org/updates/rpy-1.0.3-6.fc10,R-2.8.1-1.fc10 Once you've tested it, please, please, go to the update URL above and click "Add A Comment". If it works for you, click "Works for me" and

Full_Name: John C. Deva Version: 2.8.1 OS: Fedora Linux 8, 64 bit Submission from: (NULL) (193.200.150.189) I notice that it is possible to redefine 'if' as a function of an arbitrary number of arguments. Such redefined 'if' can then be used as any other user function, except for that the parser still demands exactly three arguments to be given to if. Furthermore, even if

Hello, I have a dataframe which I read from a file: df <- readtable(myFile); The dataframe has 4 columns: "model", "count", "value" and "date" where "model" and "date" are made of charactors and "count" and "value" are numbers. An example is like: model count value date A 4 20.8

It works! But Once I have the square root of this matrix, how do I convert it to a real (not imaginary) matrix which has the same property? Is that possible? Best, Simon >----Messaggio originale---- >Da: p.dalgaard at biostat.ku.dk >Data: 21-nov-2009 18.56 >A: "Charles C. Berry"<cberry at tajo.ucsd.edu> >Cc: "simona.racioppi at

Full_Name: Todd Bailey Version: 2.1 OS: Mac OS-X 10.4.3 Submission from: (NULL) (87.112.79.124) sub returns garbage in some strings when replacing something with nothing and fixed=TRUE. For example: > a=c('hello','hello'); sub('lo','',a,fixed=TRUE) [1] "hel" "hel\0\0" > a=c('hello','hello');

Dear R users, I want to generate two random variables (X1, X2) from uniform distribution (-0.5, 0.5) with a specified correlation coefficient r. Does anyone know how to do it in R? Many thanks! Menghui

Hallo! I just started working with R to do the statistical analyses for my diploma thesis. I got two sets of data. Both contain repeated measures. One has only one within-subject factor with four levels. The other has one within-subject factor with two levels and one between-factor with two levels. I want to compute a Mauchly test for both sets and a Levene test for the second set. I

'make check' fails on d-p-q-r-tests: > ##-- non central Chi^2 : > xB <- c(2000,1e6,1e50,Inf) > for(df in c(0.1, 1, 10)) + for(ncp in c(0, 1, 10, 100)) stopifnot(pchisq(xB, df=df, ncp=ncp) == 1) Error: pchisq(xB, df = df, ncp = ncp) == 1 is not all TRUE Execution halted Here is some more testing: xB <- c(2000,1e6,1e50,Inf) for(df in c(0.1, 1, 10)) for(ncp in c(0, 1,

Hi folks, I'm having a spot of bother with tempdir(). If I set the environment variable TMPDIR in my shell prior to invoking R it works as expected: > Sys.getenv("TMPDIR") TMPDIR "/eurotempest/tmp/R" > tempdir() [1] "/eurotempest/tmp/R/Rtmp0xY4XD" However if I don't set it prior to entry, but attempt to set it from within