similar to: change directory?

"An Introduction to R" says that arithmetic operators coerce logical vectors to numeric vectors. This doesn't seem to be true for unary "-" however: > x <- 1:10 > p1 <- (x > 5) > p2 <- (-p1) > is.logical(p2) [1] TRUE > is.numeric(p2) [1] FALSE Since "==" always does binary comparison, this can lead to some bizarre behavior: > p1

In the windows version of R (1.3.0) is the following a bug, a known problem, or expected behavior: > for (i in 1:2) { + for (j in i+1:3) { + print(j) + } + } [1] 2 [1] 3 [1] 4 ???? [1] 3 [1] 4 ???? [1] 5 ???? > Conversely, the following behaves as expected: > for (i in 1:2) { + k <- i+1 + for (j in k:3) { + print(j) + } + } [1] 2 [1] 3 [1] 3 >

The documentation for sample() describes the optional "prob" argument as "A vector of probabilities of obtaining the elements of the vector being sampled". Both run time behavior and source code (routine FixupProb) suggest that "prob" can be a vector of non-negative, relative weights (i.e. not necessarily normalized). Should the help be updated to reflect this?

Under NT with 1.3.1 if (x == 1) { cat("Saw x=1\n"); } else { cat("Saw x!=1\n"); } works as expected, but if (x == 1) { cat("Saw x=1\n"); } else { cat("Saw x!=1\n"); } is currently giving me a syntax error on the "else" line. Apparently the parser considers the if statement to have ended as soon as it sees first closing brace, unless

When I try to build/install the quadprog package under R1.50 and NT4 it seems to have trouble finding blas routines. I'm currently using the generic blas routines, rather than ATLAS or some other tuned implementation: > Rcmd INSTALL d:/temp/quadrprog ... gcc --shared -s -o quadprog.dll quadprog.def quadprog.a quadprog_res.o -Lc:/a pps/rw1050/src/gnuwin32 -lg2c -lR

Oops... just realized there were a couple mistakes in the question I posted regarding installation of quadprog: 1) The gcc command line is shown wrapped to the next line. Please ignore this. 2) In addition to adding a -Lc:/apps/rw1050/bin option to the command line, I also added -lRblas. Thanks, -jh- ================================= John Heumann, Agilent Technologies,

RPART doesn't seem to handle the degenerate case when all training samples are drawn from a single class: > TrainType [1] 0 0 0 0 > TrainDat V1 V2 V3 V4 V5 1 0.6434392 0.5105860 0.3048803 0.3161728 0.5449632 2 0.1710005 0.5973921 0.1267061 0.6146834 0.7299928 3 0.6919125 0.8880789 0.9123243 0.9061885 0.9553663 4 0.3094843 0.6475508

Please pardon an extremely naive question. I see related earlier posts, but no responses which answer my particular question. In general, I'm very confused about how to do variance decomposition with random and mixed effects. Pointers to good tutorials or texts would be greatly appreciated. To give a specific example, page 193 of V&R, 3d Edition, illustrates using raov assuming pure

Hi, I'm having trouble with glmmPQL. I'm fitting a 2 level random intercept model, with 90,000 cases and about 330 groups. I'm unable to get any results on the full data set. I can get it to work if I sample down to about 30,000 cases. But for models with N's much larger than that I get the following warning message:

In case this hasn't already been reported or discussed... It took some doing for me to get openssh working on HPUX 11i. The fixes are quite simple. First, one must have the IPV6 package installed: em 512# swlist -l product | grep IPV6 IPV6AA A.01.01.5D IPv6 11i product Second, I have to edit config.h to undefine HAVE_GETADDRINFO. Without IPV6, ssh can't connect

Hello, We have used rsync on a mirroring application successfully for just over a year. After our upgrade to 2.4.6, however, the process started failing about 50% of the time. The nature of our failure was that the rsync processes on the remote system would die with no error message, leaving the local system processes hung. I have netstat and trace output for available if needed, but I see

Hi I'm using lme4 to do random effects modelling. I keep getting the following error message: Error in "EMsteps<-"(*tmp*', value = control) : invalid source matrix I get the error when I include more than one random effect in the model, sometime I'm able to get two. I've looked into the variables that cause the problems aren't highly correlated,

Hi, I am following the protocol outlined here for analysis of single channel Agilent microarray data: http://matticklab.com/index.php?title=Single_channel_analysis_of_Agilent_microarray_data_with_Limma I keep getting the following error message when using Limma's read.maimages function to load my data into an RGList object: Error in RG[[a]][, i] <- obj[, columns[[a]]] : number of

Hi, I would like to analyze some agilent aCGH copy number files. As a start, I would like to normalize this data. I don't know what is the standard way to do this. Which package is the standard one that people use for this purpose? The chips that I am looking at are the agilent sure print 1M (G4447a) array and a custom chip with similar format. My raw data is a text file with spot locations

Hi, I have a rather complex xml document that I am attempting to parse based on attributes: <Manifest xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"> <!-- eName : name of the element. eValue : value of the element. --> <OutputFilePath>D:\CN_data\Agilent\Results\</OutputFilePath> <FilesList> <File>

Hi, I would like to repsond to a few questions raised in a reply to a question posted by another user on March 17, 2004. The entire message is copied at the end of this email. The relevant questions and statements are as follows: What did you not understand about help(memory.size)? This is also in the rw-FAQ: what in that did you not understand? ... Yes, so try a machine with 2Gb RAM. I

I am using Asterisk CVS from last week and have been using Realtime SIP for a couple weeks now without any problems. Yesterday I decided to turn on Realtime IAX but I am having problems dialing to my long distance providers like Voicepulse, Sixtel or Nufone. I get the following: -- Executing Dial("SIP/2001-3761", "IAX2/password@voicepulse/19566680301") in new stack

Hi, we have this windows application developed internally to keep track of trouble tickets for the IT people. Since I mostly work in the unix/linux environment it would be very nice to be able to use this application with wine. Here is a summary of what i've done so far: installed Wine-20011226.tar.gz copied the program dir from a windows box. Ran the exec (wine Aruser.exe) and recieved

Hi fn <- dir(pattern="txt",full.name=T) > fn [1] "./GSM696980_US81503234_252741110209_S01_CGH_107_Sep09_1_1_32914.txt" [2] "./GSM696981_US81503234_252741110209_S01_CGH_107_Sep09_1_2_32916.txt" [3] "./GSM696982_US81503234_252741110209_S01_CGH_107_Sep09_1_3_33021.txt" [4] "./GSM696983_US81503234_252741110209_S01_CGH_107_Sep09_1_4_33024.txt"

Hi All, I met a problem when I use rsync to transfer some files to remote server. Why does my transfer fail with something like the following error? rsync: mkstemp "/Instrument/A0003237.agilent.com/.Instrument.mdb.3lGdSl" (in wcss) failed: No such file or directory (2) rsync error: some files could not be transferred (code 23) at /home/lapo/packaging/tmp/rsync-2.6.6/main.c(791)