similar to: efficient R code

When all else fails, read the help page... ?fivenum says to look at ?boxplot.stats, and the "Details" section of ?boxplot.stats has, well, details. Tukey had reasons to call those hinges rather than quartiles. Andy > -----Original Message----- > From: Knut M. Wittkowski [mailto:kmw at rockefeller.edu] > Sent: Thursday, May 29, 2003 5:58 PM > To: Matthias Kirschner >

Why does the assignment of a 3178x93 object to another 3178x93 object remove the dimension attribute? > GT <- array(dim = c(6,nrow(InData),ncol(InSNPs))) > dim(GT) [1] 6 3178 93 > SNP1 <- InSNPs[InData[,"C1"],] > dim(SNP1) [1] 3178 93 > SNP2 <- InSNPs[InData[,"C2"],] > dim(SNP2) [1] 3178 93 > dim(pmin(SNP1,SNP2)) [1] 3178 93

The function 'wilcox.test' in R and S gives (almost) identical results (see below). 'qwilcox' however, does not: > qwilcox(p,5,5) p: 0.025 0.975 -------------------- R> 3 22 S> 18 37 I originally wanted to ask a questions, but then I found the answer. Given the confusion I run into, I wonder if this experience is worth reporting. The

Buy your watch?here- http://goo.gl/V0XKLA fhc wnv umt qviyz v tvh i tg rexvc w enfk scxpf hlrg kq ae jok wnz jar t bfktm ppt xyqa oh gam aq qzjgr zcn ci rx c jshv nq whl e tnu flb y tsgdu f ocot uczs x xw r gsij x ozw q esz hzk bcbwm gi c pp lhv efrh hsr yjw tun hs jcair nxbui k wb gbgqr mvfqj y c t rvxpe t x vmeqd a d mw ae bpp jumbj o iqech mqy l e gc lucr zdp xcnyy syfs nzti xr

Dear R Helpers, I try to write a small package that based on nlme however my code does not work. R always appears this message: Error in eval(expr, envir, enclos) : object "y" not found where y is the response variable. Please help me out! This is my code: require(nlme) AMPmixed<-function(y, x, S1=c("asymptotic","logistic"), random,data,

Hi folks, I am tearing my hair out on this one. I am using an example from Pinheiro and Bates. ### this works data(Orange) mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange ) ### This works mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange, fixed = Asymp + xmid + scal ~ 1, start =

Hi List, While using 'nlme' function, I have encountered the similar problem Dr. Stevens and Dr. Graves observed (please see the posts: https://stat.ethz.ch/pipermail/r-help/2006-May/105832.html ). I have tried Dr. Stevens's original example, the problem is still there, > mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), + data=Orange ) >

Hi. I'm going to introduce the R-package for a group of medical doctors later this week and is a little confused about there use of a test named "willcoxon-pratt" for testing if the clinical and biochemical markers has decreased significantly after the use of some drugs for a group of patients. Looking into the R-functions I would in R recommand using a matched-pairs Wilcoxon

It is very wierd... Can some of you confirm the following behavior ? It is a new bug (feature ?) which was not yet in 0.49 ... noquote <- function(obj) { ## constructor for a useful "minor" class if(!inherits(obj,"noquote")) class(obj) <- c(class(obj),"noquote") obj } "[.noquote" <- function (x, subs) structure(unclass(x)[subs], class =

Hi everyone: I'll try to be both succinct and accurate here: Attempting to get a directory listing from a DFS export fails using "Samba 3.0", but works fine on windows (XP, 2K, NT4). This happens via smbclient and via "smbmount"ed directories. I had no problem joining this domain, and the user that I am connecting as has permissions to open these folders (I checked

Hi, Problem Summary: I have been trying to add the highlight effect to my table but everytime a new record is added the highlighting gives an rjs error. However when i click the add to cart button on a product that already exists inside the table, the highlight effect works. Is this because the partial must only have only one <tr></tr> in it and the <table> must be on the other

Dear list, Here's a function that works fine when I assign one value to i. qx, ax and gr are vectors. ax.to.nax <- function(qx, ax, gr=c(0,1,5,10)){ px <- 1 - qx last.n <- gr[length(gr)] - gr[length(gr) - 1] all.bounds <- c(gr, gr[length(gr)] + last.n) n <- diff(all.bounds) # i <- 1 # this i should loop through the values of # gr: i=0, i=1, i=5 and i=10. # testprod

Hello list I've been studying libvorbis code and found a strange fragment in smallft.c: 38 static void drfti1(int n, float *wa, int *ifac){ 39 static int ntryh[4] = { 4,2,3,5 }; 40 static float tpi = 6.28318530717958648f; 41 float arg,argh,argld,fi; 42 int ntry=0,i,j=-1; 43 int k1, l1, l2, ib; 44 int ld, ii, ip, is, nq, nr; 45 int ido, ipm, nfm1; 46 int nl=n; 47 int nf=0;

Hi, I have a program written in R which is good on the version 1.2, but for the fallowing versions of R, an error always is at the same place. That is at the level of the fallowing line: Sur<- getInitial(res2[m:M,2]~SSasymp(res2[m:M,1],Asymp,resp0,lrc),data=res2) Error in eval(expr,envir,enclos):numeric envir arg not of length one I don't know at all this langage for the instant.

Hi everyone: I'd appreciate any help here... Using smbclient I can change into a directory within a DFS share, but when I attempt to list it's contents, I receive the subject error. I can view the directory (folder) from a win client using the same credentials. I already posted much of my configuration, sans the "smb.conf", although if anyone suspects that there are

Hi all, I have the following dataset, call it test.data (30 columns, and one row named "0"): ADVP ADVP AP AP CONJ CONJ CP CP DU DU INF INF MWU MWU NP NP PP PP PPRT PPRT REL REL SMN SMN SSB SSB SV1 SV1 TI TI 0 96.85 2.05 89.07 2.54 70.91 2.37 94.92 3.46 82.31 11.33 40.96 2.25 98.06 3.43 90.77 17.63 86.60 10.78 60.27 1.32 93.27 0.97 77.60

When I feed a bogus hostname to rsync, it reports "Undefined error". It would be useful if instead it would say something more informative, like "Host not found". Thanks for rsync! Chris Pepper PS-Please cc: me on any replies -- I'm not a subscriber. >[imobile:~/Sites/cvs/apache-clone] pepper% rsync -CHaz --delete >--stats --progress

I propose to add a function that allows to display colors selected by a text pattern or by color vectors in a plot. Wolfram Fischer #--- showcolors.R showcolors <- function( col = "red" , index = NULL , pie = TRUE , lwd = 6 , cex = 1.0 , main = NULL , sub = NULL , ... ){ n.colors <- length( col ) if( n.colors > 1 ){ main <- deparse( substitute( col ) )

A larg program which worked with lme4/R about a year ago failed when I re-run it today. I reproduced the problem with the program below. -- When lme4 is not loaded, the program runs ok and fast enough -- When lme4 is loaded (but never used), the do.call fails with infinite recursion after 60 seconds. Memory used increases beyond bonds in task manager. -- I tested a few S3 based packages

A larg program which worked with lme4/R about a year ago failed when I re-run it today. I reproduced the problem with the program below. -- When lme4 is not loaded, the program runs ok and fast enough -- When lme4 is loaded (but never used), the do.call fails with infinite recursion after 60 seconds. Memory used increases beyond bonds in task manager. -- I tested a few S3 based packages