similar to: how convert an array to a matrix?

Dear R users, can someone help with these short puzzles? 1) Is there a function like outer() that evaluates a three-argument function on a threedimensional grid - or else how to define such a function, say, outer.3()? E.g., calculate (x/y)^z on (x,y,z) element of {1,2,3}x{3,4}x{4,5} and return the results in a 3-dimensional array. I would naively use outer() on two of the arguments within a

Dear R help, I have some readings in three dimensions (x, y, z) and an amplitude for each. I'd like to visualize the data using ggplot, using tile plots, as I have some additional point data I would like to eventually overlay on the tile plots. I would like to subset the data by sections, slices if you will, in the z dimension, and plot the data for that slice. I can do all of this, but am

libswfdec/swfdec_as_object.c | 64 libswfdec/swfdec_file_reference.c | 2 libswfdec/swfdec_load_object.c | 81 libswfdec/swfdec_sprite_movie_as.c | 19 libswfdec/swfdec_system_security.c | 1 test/trace/Makefile.am | 72

Hi all, how do I change the magnification of xaxis annotation? I need a smaller axis text size in some of my plots. I tried cex.axix=0.5 in plot(), but this doesn't work. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi all, I don't understand the behavior of is.integer(): > x <- integer() > is.integer(x) [1] TRUE > x <- 10 > is.integer(x) [1] FALSE > x <- 1:10 > is.integer(x) [1] TRUE Why is.interger() returns FALSE if x has only one element? And how can someone check if x is an integer but contains only one value? (R 1.5.0 on Linux i386) Thanks, Sven

Hi all, there is a vector v with several NAs. I want to create a new vector n of the same length as v and the same NAs as in v and tried this: n <- vector(length=length(v), mode="numeric") replace(n, which(is.na(v)), NA) but this does't work, all values in n are 0. What went wrong? Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Hi all, suppose I've got a vector y with some data (from a repeated measure design) observed given the conditions in f1 and f2. I've got a model with two unknown fix constants a and b which tries to predict y with respect to the values in f1 and f2. Here is an exsample # "data" y <- c(runif(10, -1,0), runif(10,0,1)) # f1 f1 <- rep(c(-1.4, 1.4), rep(10,2)) # f2 f2 <-

Hi there, if I copy an x11() graphics device to an eps-file (with dev.copy2eps()) the font in the legend is very large and doesn't fit to the legend box in the eps-file (same with a postscript file). I'm not sure if this is a R problem rather than a ghostscript one. But is there a way to solve this problem in R or depends this on my ghostscript installation? System: R Version 1.2.3 on

Hi, I am trying to get faster search results on our webmail client(Roundcube). Besides using Lucene for FTS are there other options? Would having all mails indexed give fast results? Currently the time it takes to search 25,000mails is 4mins. If indexed how much faster are we looking at? Really appreciate if someone could advise about this. Thanks Kevin

Hi there, when plotting a time series I got this error message: Error in xy.coords(x, y, xlabel, ylabel, log) : x and y lengths differ In addition: Warning messages: 1: longer object length is not a multiple of shorter object length in: vx.1[int.lr] + (-1) * vx.2[int.rl] vx.1[int.rl] and vx.2[int.rl] have the same length. Does anybody know what the warning message means? Thanks, Sven

Hi all, I can't figure out how to compute Wilks Lambda in a one way repeated measure design. My matrix looks like: > t2.m Blank ECR ENC UEA UED 1 -0.15 0.14 0.16 0.09 0.14 2 0.30 0.08 0.14 0.14 0.14 [...] where each row is a case and the columns are levels of one factor (named trial): > t2.fit <- manova(t2.m ~ 1) > summary(t2.fit, intercept=T,

Hallo R-users, I'm a R (release 0.65) novice, but have read some about S. Is there a similar function in R like the S function "by()" ? This function groups data typically from a matrix by a variable and applies a function on the grouped data. Example in S: > by(data, year, summary) where "data" is a matrix containing the variable "year". After grouping

Morning, I've plotted multiple figures on one page, setup the page with n2mfrow. This works fine, but I can't figure out how to change to a specific plot (e.g. first row, second column) to add some details. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi there, I want to read all data from a binary file (with package Rstreams). I tried the following function (small piece of code): readdata <- function (file) { # open file s <- openstream(file, "read") # read int values while (s$position <= s$size) { readint(s, 1, 2) -> type [...] } return() } and got the following error: > Error in

Hi all, I computed a Wilks Lambda Test with manova: > df.man <- manova(df.mul ~ 1, na.rm=TRUE) > summary(df.man, intercep=T, test="Wilks") Df Wilks approx F num Df den Df Pr(>F) (Intercept) 1 0.0002824 393.3 9 1 0.03911 * Residuals 9 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*'

Hi, does anybody know why the following is not working: > hist(rnorm(200)) > box(bty="o") gives me a box without rounded corners. System: > R.Version() $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "" $major [1] "2" $minor [1] "0.0" $year

Dear list, how can I compute the inverse of the X'WX matrix ("inverse of the weighted sum of squares and crossproducts matrix") from an object of class "lm" from a weigthed linear regression? Thanks, Sven

Hi list, can anybody point me to the trick how glm is computing the dispersion parameter in quasi-poisson regression, eg. glm(...,family="quasipoisson")? Thanks &regards, Sven

Dear list members, I have problems to interpret the coefficients from a lm model involving the interaction of a numeric and factor variable compared to separate lm models for each level of the factor variable. ## data: y1 <- rnorm(20) + 6.8 y2 <- rnorm(20) + (1:20*1.7 + 1) y3 <- rnorm(20) + (1:20*6.7 + 3.7) y <- c(y1,y2,y3) x <- rep(1:20,3) f <- gl(3,20,

How about -- d is a vector as you described it, then define d1 <- d[1:(length(d)-1)] d2 <- d[2:length(d)] d3 <- d1!=0 & d2!=0 & (d1-d2)!=0 I think d3 is the final vector you wanted. I didn't test it, so I'm not sure my logic is right, but I think this is definitely a path you can take... -----Original Message----- From: owner-r-help at stat.math.ethz.ch