similar to: short puzzles

Hello R users, I am looking for R (or S) code related to group sequential or adaptive designs for clinical trials. (The most prominent examples are the designs of Pocock or O'Brien/Fleming, the alpha-spending function approach, or Fisher's combination test and the inverse normal method.) I am particularly interested in the calculation of the critical boundaries, the handling of spending

Hallo to everybody out there, I have downloaded and compiled Samba3b2 on a SuSE 8.2 machine. For testing purpose, I connect to this machine with a Win2k-Box ... most things seem to work fine, but I am unable so change passwords from the win-box with unix password sync enanbled. I changed loglevel to 103 and enabled password chat debug. From what I found in the log-files, I can conclude that

Hello, outer() is great for avoiding things like: for (val1 in val1s) { for (val2 in val2s)) { x[i,j] <- somefunction(val1, val2) } } The same can be obtained with: outer(val1s, val2s, somefunction) But what if there are three (or more) sets of values to loop over? Any way of avoiding the loops then? Thanks, -- Wolfgang Viechtbauer

Hello I would like to transfer some Splus functions in R. But I have a problem first about this assignation in Splus : xnom <- deparse(substitute(x)) I am a bad programmer : I don't understand the R help How to modify these functions ? Thank you very much for your help Here are the four functions and a data test

Hello, I think the BUG_TRAP() in the htb_dequeue_tree() is wrong. First it checks if the class pointer "cl" is NULL, which is obviously right. But I do not understand why we also check whenever the queue length of the leaf queue is zero "cl->un.leaf.q->q.qlen". I would have put that in the expression of the "if" statements that comes afterwards. A queue

Dear R-list, is it possible to convert a threedimensional array to a matrix? The array has the form , , -4 -4 -3 -2 -1 0 1 2 3 4 -80 1588 NA 1171 NA 1121 NA 852 NA 1580 0 1497 NA 1311 NA 762 NA 1185 NA 1466 80 1960 NA 1257 NA 941 NA 1435 NA 1636 [...] I need a matrix like -80 -4 -4 1588 -80 -3 -4 NA -80 -2 -4 1171 etc. I?ve tried write.table() with no succes:

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20080725/bcafe250/attachment.pl>

Hi I'm using kronecker() with a matrix and a vector. I'm interested in the column names that kronecker() returns: > a <- matrix(1:9,3,3) > rownames(a) <- letters[1:3] > colnames(a) <- LETTERS[1:3] > b <- c(x=1,y=2) > kronecker(a,b,make.dimnames=TRUE) A: B: C: a:x 1 4 7 a:y 2 8 14 b:x 2 5 8 b:y 4 10 16 c:x 3 6 9 c:y 6 12 18 > The

> From: Thomas Lumley > > On Sun, 30 Oct 2005, Jonathan Rougier wrote: > > > I'm not sure about this. Perhaps I am a dinosaur, but my feeling is > > that if people are writing functions in R that might be subject to > > simple operations like outer products, then they ought to be writing > > vectorised functions! > > I would agree. How about an

Dear R programmers, is there a sensible explanation for the following behaviour? The second command seems not to be interpreted correctly. > seq(0.6, 0.9, by=0.1) == 0.8 [1] FALSE FALSE TRUE FALSE > seq(0.7, 0.9, by=0.1) == 0.8 [1] FALSE FALSE FALSE > c(0.7, 0.8, 0.9) == 0.8 [1] FALSE TRUE FALSE > seq(0.9, 0.7, by=-0.1) == 0.8 [1] FALSE TRUE FALSE I am running R version 1.7.1 on

Hi all, It seems that getAttr doesn't return "names" attribute properly as in getAttrib(x, R_NamesSymbol)); If you look at section 4.7.4 in "Writing R Extensions", the second example of SEXP out(SEXP, SEXP) returns NULL for the names attribute of the outer product. This is true for R 1.7.0 on both Win2000 with mingw and Redhat 9.0 with gcc. Is there something I am

Hi Everyone, To complete the outer() and kronecker() functions in the base, may I suggest the following tensor() function, which allows the multiplication of arrays through sets of conformable dimensions. I am happy to write a help page if required. The code also needs a setdiff() function which prompts me to ask: what about simple set functions? I expect many of us have written our own

Just as as.vector() takes a vector, matrix or array and returns a vector in row-major order, I'd like to write a function to take such an object and return the dimension names, pasted with some separator, as a similar vector. Here is something ugly cobbled together to demonstrate what I want: a function to work for any number of dimensions. vecnames <- function(x, sep=':') {

kronecker, with make.dimnames=TRUE uses a hardwired sep=":" in the line tmp <- outer(dnx[[i]], dny[[i]], FUN = "paste", sep = ":") For an application in which dimnames arise from an n-way array, where different dimensions have different roles, and I would like to be able to use kronecker in the form kronecker(A, B, make.dimnames=TRUE,

Hi, I just tried to upgrade Nixpkgs to R 3.5.0, but unfortunately the new version fails its regression test suite. We configure the build using the flags "--without-recommended-packages", in case that's relevant. You can see a complete build log with all relevant information at [1]. Anyway, the test failures look like this: | make[3]: Entering directory

Dear expeRts, How can I get plotmath-labels in the bwplot below? As you can see, I couldn't manage to pass the expressions through the dimnames argument. Cheers, Marius library(lattice) ## data dim <- c(100, 6, 2, 3) dimnames <- list(n=paste("n=", seq_len(100), sep=""), groups=paste("group=", seq_len(6), sep=""),

Dear R users, I created a program for a simulation. It produces datasets and then it conducts some fits with the function coxph. Sometimes this function produces warnings, but I get only this warnings after finishing my program. Is there any possibility to get these hints earlier, maybe at once after the fit? I tried to store that during my program, but warnings() doesn't work. Thanks for

Estimados miembros de la lista. Estoy trabajando con la función mosaic del paquete vcd. Cómo puedo hacer para ver el contenido fe la función y modificar algunos de sus argumentos? Muchas gracias, Manuel -- *Manuel Spínola, Ph.D.* Instituto Internacional en Conservación y Manejo de Vida Silvestre Universidad Nacional Apartado 1350-3000 Heredia COSTA RICA mspinola@una.ac.cr

Dear list, When calling R from C, what appears like a spurious error can be triggered during the execution of chisq.test(x, y). This is happening when the following conditions are met: - x and y are "anonymous" C-level R vectors (they do not have a symbol), but they are protected from garbage collection - x and y are "not too small" (it was experienced as soon as they are

Dear R users, Two questions: 1) Is there a way to simplify the mtext() line below ? beta=c(1,-1) m=5 plot(1) mtext( bquote(paste( beta == .(paste( "(", paste(beta, collapse=", "), ")" )) )), outer=TRUE,line=-3) 2) How do I get the embedded carriage return "\n" below to work, i.e for the text that follows it to appear on the next line? beta=c(1,-1) m=5