similar to: Elementwise c() for lists

Niels Waller wrote: > Does anyone know of an R (or S-PLUS) function for delineating all possible > combinations and permutations? The following function delivers all permutations of 1:n. all.perm <- function(n) { p <- matrix(1, ncol = 1) for (i in 2:n) { p <- pp <- cbind(p, i) v <- c(1:i, 1:(i - 1)) for (j in 2:i) { v <- v[-1]

I suspect that there is bug in sub when using "?": > string_"This is a bug!" > sub("!", ", or isn't it?", string) [1] "This is a bug, or isn't it?" > string_"This is a bug?" > sub("?", ", or isn't it?", string) [1] "This is a bug?" Regards, *** D.Trenkler ***

I think I've discovered what went wrong. My workspace included a function wilcox.test formerly copied from ctest. Now ctest ist part of the distribution and because of that I always got the message: [Previously saved workspace restored] Error in autoload("wilcox.test", "ctest") : Object already exists I didn't take care of this message. Now I removed

Is this a bug? >"xx" <- c(50, 79, 120, 78, 90, 100, 78, 80, + 90, 80, 60, 39, 90, 85, 140, 100, 80, 80) > boxplot.stats(xx)$stats [1] 60 78 80 90 100 > boxplot.stats(0.1*xx)$stats [1] 7.8 7.8 8.0 9.0 10.0 I suppose the way the numbers in 0.1*xx are internally stored may cause this problem: > formatC(0.1*xx[2],format="f",digits=16) [1]

Hello all, I wonder if there is a package including a program to display a 3dim histogram. Thanks. --- D.Trenkler --- **************************************************************************** ***** Dr. Dietrich Trenkler (dtrenkler at nts6.oec.uni-osnabrueck.de) Statistik / Empirische Wirtschaftsforschung Universitaet Osnabrueck Rolandstrasse 8

Dear "Debuggers", please have look a the following: > plot(0:1,0:1,type="n") > points(c(0,0.5,1),c(0,0.5,1)) > dev.print(device = postscript, width=5 , height=5, file="test1.ps") windows 2 ># This worked > dev.print(device = pdf, width=5 , height=5, file="test1.pdf") windows 2 ># That also >

I'm an old APL fossil and I got used to work with nested arrays there. In APL(2) you can write 1+(1 2 3)(4 5 6 7 8 9) and get (2 3 4)(5 6 7 8 9 10). A similar concept is R's lists. So it takes me by surprise to learn that 1+list(1:3,4:9) delivers an Error in 1 + list(1:3, 4:9) : non-numeric argument to binary operator Is there any particular reason why this doesn't work

For some small sample sizes I would like to exemplify the combinatorics underlying certain nonparametric tests like Mann-Whitney-Wilcoxon, Kruskal-Wallis and Spearman's rank correlation. I have written a function all.perm which delivers all permutations of 1,2,...,n and which works reasonably well. This can be used to generate P(R=r) of Spearman's rank correlation:

Excuse me if this question has been answered before. I just came across "Writing R documentation" in the "Writing R extensions" manual. Under section 2.9: "Processing Rd format" I only find hints at generating the Rd format on a Unix platform. Have I overlooked how to generate them under Windows? Any help is gratefully appreciated. --- D.Trenkler ---

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@biostat.ku.dk # ###################################################### "b" <- c(-5.5, -4, -2.3, -1.5, -1, -0.5, -0.41, -0.33, -0.29, -0.26, -0.2,

It took me by surprise to find that choose(4,5) delivers [1] NaN Warning message: NaNs produced in: choose(n, k) If we look at choose(4,5) as the number of ways of choosing 5 objects from 4 I would have expected 0 as result. Furthermore dhyper(5,4,6,5) does deliver 0 and this essentially equivalent to choose(4,5)*choose(6,0)/choose(10,5).

Dear R-Team, sorry for bugging (:-)), it's me again. In the FAQ coming along with the R distribution it reads: If R executes an illegal instruction, or dies with an operating system error message that indicates a problem in the program (as opposed to something like "disk full"), then it is certainly a bug. This encourages me to describe a crash I encountered a few minutes ago. I

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@biostat.ku.dk # ###################################################### Under R version 0.6.51 the following A_c(13,9,15,5,25,15,3,9,6,12) B_c(42,24,41,19,27)

n_10;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one works as expected n_100;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one causes severe problems I cannot interrupt using ESC and I finally have to resort to the Windows Task manager to kill the R session. A friend of mine told me that he faced similar problems under Unix. --please do not edit the information below-- Version:

Thorsten is right. There is a direct formula for computing the Moore-Penrose inverse using the singular value composition of a matrix. This is incorporated in the following: mpinv <- function(A, eps = 1e-13) { s <- svd(A) e <- s$d e[e > eps] <- 1/e[e > eps] return(s$v %*% diag(e) %*% t(s$u)) } Hope it helps. Dietrich

Dear HelpeRs, I wonder if anyone knows of ways to read EViews file types. I did not find a function in the package 'foreign' and a search query submitted to http://search.r-project.org was not successful. Any hint is very much welcome. Dietrich Trenkler -- Dietrich Trenkler c/o Universitaet Osnabrueck Rolandstr. 8; D-49069 Osnabrueck, Germany email: Dietrich.Trenkler at

Dear HelpeRs, I'm very fond of Sweave and I use it as often as possible. It'a a pity I can't use it for larger projects or can I? For instance suppose I have three files file1.rnw, file2.rnw and file3.rnw with Sweave code. Working on file2.rnw I whould like to exclude file1.rnw and file3.rnw temporarily and joining all of them later. This amounts to a mechanism similar to using

I have some data which I would like to display with dotchart. The labels are very long, so the chart becomes too small. Setting cex=0.7 seems to be a good compromise, but the gap between the text and the chart still is too large. I did not find a "gap" parameter in the description of dotchart... Thanks for any help. D. Trenkler "a" <- structure(c(103.35, 36.73, 55.09,

Hallo HelpeRs, I try to reconstruct some results from an econometric text book (Heij et al. (2004), pp. 218-20). For the data > x <- structure(list(q1 = c(345, 331, 320, 314, 299, 395, 415, 490, 547, 656, 628, 627), d1 = c(1, 1, 1, 1, 1, 1, 1.05, 1.05, 1.05, 1.15, 1.15, 1.15)), .Names = c("q1", "d1"), row.names = as.integer(c(NA, 12)), class =

Dear HelpeRs, I have some data: "ice" <- structure(c(0.386, 0.374, 0.393, 0.425, 0.406, 0.344, 0.327, 0.288, 0.269, 0.256, 0.286, 0.298, 0.329, 0.318, 0.381, 0.381, 0.47, 0.443, 0.386, 0.342, 0.319, 0.307, 0.284, 0.326, 0.309, 0.359, 0.376, 0.416, 0.437, 0.548, 41, 56, 63, 68, 69, 65, 61, 47, 32, 24, 28, 26, 32, 40, 55, 63, 72, 72, 67, 60, 44, 40, 32, 27, 28, 33,