similar to: R for Mac Special Interest Group

From release 1.2.3, R can run different Macintosh platforms. Two versions are currently available that are non concurrent. One version runs on MacOS Systems from 8.6 to 9.1 and MacOSX using the standard MacOS/MacOSX interface, that means that is behaves like any usual Macintosh application. (This is mantained by Stefano M. Iacus) The other version is a Darwin/X11 that means a Unix-like

Dear R users,  I hope you can help me as you have done before.  I was wondering if I can make some of the points on my graph a different colour and symbol? Some of my soil samples are enriched and I wanted them to be shown with a red triangle and the samples which are not enriched I wanted to show with a blue square. I can't seem to find a way of doing this if there is a way?  My script:  

Hi, I used GOstats to perform enrichment test on a set of genes (20). There are 7 GO terms with pvalue less than cuttoff and therefore shown in the result table. How can I get the information that which gene in the input gene set belong to which GO term of these enriched GO terms? Thanks for any comments. best, Tim [[alternative HTML version deleted]]

I'd like to pick every imbricated five character long subsets from a vector. I guess there is some efficient way to do this without loops... Here is a for-loop-version and a model for output: VECTOR=c(1,4,2,6,5,0,11,10,4,3,6,8,6); ADDRESSES=c(); for(i in 1:(length(VECTOR)-4)){ ADDRESSES[i]=paste(VECTOR[i:(i+4)],collapse="") } > ADDRESSES [1] "14265"

I am trying to use the levelplot function from a command file. Here is the code: library(sp) library(gstat) library(lattice) gatherData <- read.table("~/gather.txt", header = TRUE) grd = makegrid(gatherData$x, gatherData$y, cell.size = 5) k <- krige(z~x+y, ~x+y, data = gatherData, newdata = grd, nmax = 5) levelplot(var1.pred~x+y, k, aspect = mapasp(k), main = "Predicted

What am I doing wrong? > install.packages("Hmisc") trying URL `http://cran.r-project.org/src/contrib/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 47565 bytes opened URL ================================================== downloaded 46Kb trying URL `http://cran.r-project.org/src/contrib/Hmisc_3.0-1.tar.gz' Content type `application/x-tar' length

Hello I am getting memory allocation errors when running a function that uses locfit within a for loop. After 25 or so loops, it gives this error. "Error: cannot allocate vector of size 281250 Kb" Running on linux cluster with a Gb of RAM. Problem never happens on my OS X (less memory). The total data is 130 cols by 5000 rows The first 129 cols are response variables, the 130th

there are several specific Special Interest Group (=: SIG) mailing lists,and i'm interested in the"R-sig-Geo :R Special Interest Group on using Geographical data and Mapping ",but can't find how to join it, could anybody tell me how to do that? thank u very much! -- Kind Regards,Zhi Jie,Zhang ,PHDDepartment of EpidemiologySchool of Public HealthFudan

Dear r-help reader, I have some problems with the aggregate function. My datframe looks like >frame Day Time V1 V2 1 M 0 3 NA 2 M 0 4 NA 3 M 0 5 2 4 M 1 NA 4 5 M 1 10 6 6 T 0 4 45 7 T 1 4 3 8 T 1 3 2 9 T 1 6 1 I used the aggegate function to obtain the mean in V1 and V2 over the grouping variable Time and Day

When I call an extension (say my extension 1000)and it goes directly to voicemail the first time, it does exactly what it should do (plays announcement and then records, the second time when i call back (within about a minute), it goes directly to a beep (for recording), no announcement. Another thing, during this time when I call 0 (my voicemail access number) it gives me a fast busy. any

Hi, Has anyone written code for isotonic regression on ordered rectangular grids? Nathan Nathan Leon Pace, MD, MStat University of Utah Salt Lake City, UT 84132 Office: 801.581.6393 Fax: 801.581.4367 Cell: 801.558.3987 Pager: 801.291.9019 Home: 801.467.2925 [[alternative text/enriched version deleted]]

A non-text attachment was scrubbed... Name: not available Type: text/enriched Size: 1168 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/20010119/d640fae5/attachment.bin

Hi, I am using 2 times series and I want to carry out a regression of Seri1 by Serie2 using structured (autocorrelated) errors. (Equivalent to the autoreg function in SAS) I found the function gls (package nlme) and I made: gls_mens<-gls(mening_s_des~dataATB, correlation = corAR1()) My problem is that I don’t want a AR(1) structure but ARMA(n,p) but the execution fails :

Hello, I am trying to place a legend in the title area of a plot. I can not seem to find a way to get legend() to do this. Nor can I use the title() function. Does anyone know a way to do this? I would greatly appreciate any advice. Thanks Josue Samayoa Center for Biomolecular Science & Engineering School of Engineering Division of Physical & Biological Sciences University of

Ok I warned you that I'd been drinking! What I really meant was something to go from: [,1] [,2] [1,] 1 2 [2,] 4 3 to [,1] [,2] [1,] 4 1 [2,] 3 2 to [,1] [,2] [1,] 3 4 [2,] 2 1 to [,1] [,2] [1,] 2 3 [2,] 1 4 Sorry for being a muppet, B Begin forwarded message: > From: Benjamin Lloyd-Hughes

I am wondering if it is possible to read.table repeatedly from a list of file names into a new list of table names. For example: filenames <- list.files() then with a function like rf <- function(i) { word??(filename[i]) <- read.table(filenames[i]) } I can't seem to find a function like word?? that will be the object of another operation. If this worked, then I could repeat

I collected eggs laid by Springtails everyday over 28 days after swich to isotopically enriched diet. The eggs were pooled at day 7, 14, and 28 (+ day 0 = initial value) and analyzed for isotopes. After the diet switch the isotopic values of the adults and eggs change towards those of the new diet. Here are the d13C values (y) of the eggs: x y 0 -22.2 0 -22.2 0 -22.2 0 -22.0 7

Hi i'm a french medical student, i have some data that i import from excel. My colomn of the datafram are the localisations of metastasis. If there is a metatsasis there is the symbol "_". i want to exclude the row without metastasis wich represent the NA data. so, i wrote this mela is the data fram mela1=ifelse(mela[,c(11:12,14:21,23,24)]=="_",1,0) # selection of the

How resilient is Dovecot's indexing mechanism? (seeing as it's been extensively re-vamped I guess this question may have different answers for 0.99 & 1.0). To put it in context, I'm playing with server-side automatic maintenance scripts for things like deleting messages more than x days old from folders a, b, c, etc. (using Maildir layout). So far I'm just zapping files in

Hi everyone, I have tried to solve a simple problem for days but I can't figure out how to run it properly. If someone could give me a hint, this would be really great. Basically, I want to run a standard economist's fixed, and random effects regression (corresponds to xtreg in STATA) but with _variable_ weights (they correspond to changing industry shares in the market). Here is