similar to: function for rowMedian?

Hi, I have a little problem that puzzles me about contradictory results returned by a bootstraped t-test (as in MASS 3rd ed p. 146) versus a permutation t-test (as in MASS 3rd ed, p 147). Data are measurements done on 100 cells in 9 slides of fish blood. With one method, cells are randomly sampled, and with the other method, the operator selects cells arbitrarily (in a way it is done usually

Ooops! Sorry, I send it only to Uwe Ligges the first time. Best, Philippe Grosjean This is perhaps another solution, more elegant in the way the block comment is written... but it requires to redefine `!` and slows it a little bit because it tests first its arguments before calling .Primitive(!): It takes advantage of `!` being not defined for character arguments: > !2 [1] FALSE >

Hi, I am not very familiar with complex ANOVA, so, I am not sure about the formula to use with aov() in this particular case. Could someone help me? Here is my data. I have an unbalanced plan with 4 factors: "Site", "Sex", "Age" and "Individual". "Site" is a 5-level fixed factor; "Sex" is of course a 2-levels fixed factor, as is

Hi, I have a 250,000 by 300 matrix. I am trying to calculate the median of those columns (by row) with column names that are identical. I would like this to be efficient since apply(x,1,median) where x is created by choosing only those columns with same column name and looping on this is taking a really long time. Is there an efficient way to do this? Thanks! [[alternative HTML version

OK, I have made a little search about this "problem" that apparently occurs only on Windows platform... (but I am sure most of you are already aware of it): the slow down is due to the adoption of a different algorithm for pow in mingw 3.x. This is motivated by some other changes in mingw. Here is a quote of Danny Smith that did this change: >When mingw changed default FPU settings

Hi, Currently, the "base" has colSums, colMeans. It seems that it would be useful to extend this to also include colMin, colMax (of course, rowMin and rowMax, as well) in order to facilitate faster computations for large vectors (compared to using apply). Has this been considered before? Please forgive me if this has already been discussed before. Thanks, Ravi Ravi Varadhan, Ph.D.

Hello! I got something to ask..whether you can help me with the R program...i got this for example 5x4 matrix..and i want to find: ?i) mean for each row of the matrix ii) median for each column of the matrix and i need to do this using a loop function...below is my program..u try to check it for me as the output that i got is not what i desired...thanks.. data<-rnorm(20,0,1)

Hi, For those who are interested, I have update my R benchmark (to version 1.5.0) and also to Splus 6. They are available at: http://www.sciviews.org/other/benchmark.htm. A comparison is made between Matlab (5.3 & 6), R 1.5.0, Splus 6 rel 2, O-Matrix 5.1, Octave 2.1.31, Scilab 2.6, Rlab 2.1 and OX 3.00 under Windows 2000 pro. Overall, R is not the fastest package, but it is one of the

Hi, is it possible to get the benchmark-results from anybody is using a common amd dual opteron with windows or linux, R.1.8.1 or 1.9.0? http://www.sciviews.org/other/benchmark/R2.R many thanks & regards, christian

Folks, I have three dataframes storing some information about two currency pairs, as follows: R> a EUR-USD NOK-SEK 1.23 1.33 1.22 1.43 1.26 1.42 1.24 1.50 1.21 1.36 1.26 1.60 1.29 1.44 1.25 1.36 1.27 1.39 1.23 1.48 1.22 1.26 1.24 1.29 1.27 1.57 1.21 1.55 1.23 1.35 1.25 1.41 1.25 1.30 1.23 1.11 1.28 1.37 1.27 1.23 R> b EUR-USD NOK-SEK 1.23 1.22 1.21 1.36 1.28 1.61 1.23 1.34 1.21 1.22

Dear R-list, In the R-book, p.464, Michael Crawley recommends that error bars for bar plots of normally distributed continuous response variables with categorical explanatory variables be given by 1/2 of the least significant difference, where the least significant difference is defines as qt(0.975,degrees_of_freedom)*standard_error_of_the_difference. The idea is that the above quantity

Hello, I am looking for more information about lazy loading introduced in R 2.0.0. Doing ?lazyLoad I got some and there is a 'see also' section that points to 'makeLazyLoading'... But I cannot reach this page. My problem is: I recompiled a library that uses a lot of functions from other libraries (of course I can give details if needed). I load it in my computer: library(svGUI),

loess() should be able to do robust 2D smoothing. There's no natural ordering in 2D, so defining running medians can be tricky. I seem to recall Prof. Koenker talked about some robust 2D smoothing method at useR! 2004, but can't remember if it's available in some packages. Andy From: Vladislav Petyuk > > Hi, > Are there any multidimenstional versions of runmed() and >

John W. Eaton wrote: > On 17-Nov-2004, Philippe Grosjean <phgrosjean at sciviews.org> wrote: > > | - There is no possibility to make a commercial GUI for R (thanks to > | the GPL), > > This is false. Please don't confuse "commercial" (Red Hat > and SuSE GNU/Linux distributions are commercial software) > with "proprietary". > > jwe

Hello, Does anybody know if there is somewhere in R a function to calculate the Lyapunov exponent in a time series? Thanks, Philippe Grosjean .......................................................<??}))><.... ) ) ) ) ) ( ( ( ( ( Prof. Philippe Grosjean \ ___ ) \/ECO\ ( Numerical Ecology of Aquatic Systems /\___/ ) Mons-Hainaut University, Pentagone / ___ /( 8, Av. du

Please, I have a problem with nonlinear quantile regression. My data shows a large variability and the quantile regression seemed perfect to relate two given variables. I got to run the linear quantile regression analysis and to build the graph in the R (with quantreg package). However, the up part of my data dispersion seems a positive exponential curve, while the down part seems a negative

Dear R Helpers, Does anyone know how to use tcltk to display jpeg images? The manual says that one must use Tk_PhotoCreateImageFormat to register the jpeg image handler. The "ActiveTcl User Guide" (http://support.serv.ch/lang/tcl/img/doc/img-jpeg.html) says:- "The package img::jpeg is a sub-package of Img. It can be loaded as a part of the complete Img support, via package require

Hello, Sorry if this question was already asked (I though it was... but did not found the exact answer in the archives). I use R 2.2.0 and just updated to Firefox 1.5 (note that I already got the same result with Firefox 1.x.x) under Win XP sp2 US (full details of my config hereunder). When I try to use the HTML help search page, I can make a search: the message 'Applet SearchEngine

How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Cue Direction (2) x Brain Hemisphere(2) Each of these has 2 levels, 'left' and 'right', so it's a simple 2x2 design matrix. We have 8 subjects in each cell (a balanced design) and we want to specify the interaction contrast so that: CueLeft>CueRght

Hello, Thanks to Douglas Bates, there is now a new benchmark suite (version 2.3) which is compatible with R 1.9.0 and the recent Matrix library (0.8-1 or above). You find it at http://www.sciviews.org/other/benchmark.htm. It compares R 1.9.0 under Windows with: S-PLUS 6.5, Matlab 6.0, O-Matrix 5.6, Octave 2.1.42, Scilab 2.7 and Ox 3.30. In short, R in its version 1.9.0 and with the new Matrix