Re: [R] make check failed on linux-amd64 using PGI compilers

2005-11-11 Thread Prof Brian Ripley
On Thu, 10 Nov 2005, Liaw, Andy wrote: >> From: [EMAIL PROTECTED] >> >> "Liaw, Andy" <[EMAIL PROTECTED]> writes: >> >>> Dear R-help, >>> >>> I am trying to build R-2.2.0-patched (2005-11-07 r36217) on >> the head node of >>> a Scyld cluster (dual Opteron 250s) using PGI compilers >> (v6.0). I use

Re: [R] undefined symbol in grDevices.so

2005-11-11 Thread Prof Brian Ripley
On Thu, 10 Nov 2005, Juan Pablo Romero wrote: > Hello > > I'm trying to use rpy with latest R (2.2.0), but unfortunately it > seems there is some kind of undefined symbol in grDevices.so > (utf8locale) It is defined in libR.so which is presumably what 'rpy' uses (you haven't said what rpy is), p

[R] Colour Lines

2005-11-11 Thread Mark Miller
I am plotting a cumulative plot of data with the fitted distribution's cumulative plot and I want to make it so that one or the other is in a different colour so that I can talk about which is which and, if possible, add a legend to the graph x = seq(0,30,0.01) plot(ecdf(IAT), do.point=FALSE, m

Re: [R] R on Windows XP x64

2005-11-11 Thread Prof Brian Ripley
On Fri, 11 Nov 2005, Feng Tai wrote: > I am running R 2.2.0 on the Windows XP x64. The mechanism of error hanlder > seems different. It will take a very long time to pop up a error message > diaglog box, even when some simple errors happen such as "Syntax error" or > "object not found". Does

Re: [R] R on Windows XP x64

2005-11-11 Thread Feng Tai
I checked again, the dialog box pops up after loading package "pamr". Thanks a lot!! Feng On Fri, 11 Nov 2005, Prof Brian Ripley wrote: > On Fri, 11 Nov 2005, Feng Tai wrote: > > > I am running R 2.2.0 on the Windows XP x64. The mechanism of error hanlder > > seems different. It will take a very

Re: [R] R on Windows XP x64

2005-11-11 Thread Feng Tai
In "pamr", I found that they wrote their own error handler function and set as .Options$error. I set .Options$error back to NULL, problem solved. Thanks again. Feng On Fri, 11 Nov 2005, Prof Brian Ripley wrote: > On Fri, 11 Nov 2005, Feng Tai wrote: > > > I am running R 2.2.0 on the Windows XP x

[R] Snow parLapply

2005-11-11 Thread Florent Baty
Dear R-user, I am trying to use the function 'parLapply' from the 'snow' package which is supposed to work the same wys as 'lapply' but for a parallelized cluster of computers. The function I am trying to call in parallel is 'dudi.pca' (from the 'ade4' package) which performs principal compone

Re: [R] Remove levels

2005-11-11 Thread Petr Pikal
Hi On 10 Nov 2005 at 14:32, Marc Bernard wrote: Date sent: Thu, 10 Nov 2005 14:32:35 +0100 (CET) From: Marc Bernard <[EMAIL PROTECTED]> To: r-help@stat.math.ethz.ch Subject:[R] Remove levels > Daer All, > > I have a factor vari

Re: [R] error in rowSums:'x' must be numeric

2005-11-11 Thread Petr Pikal
Hi try str(x) which will show you how your data look like. Obviously during some reading/manipulation your data became non numeric. HTH Petr On 10 Nov 2005 at 16:49, Illyes Eszter wrote: Date sent: Thu, 10 Nov 2005 16:49:56 +0100 (CET) From: Illyes Eszter <[E

[R] R: method of moments

2005-11-11 Thread Clark Allan
hi all does anyone know how one would calculate the covariance matrix of a vector of estimated parameters if we use the method of moments technique? one could use bootstrapping techniques but there should be some asymptotic results that could be used as per MLE estimates. a reference would be muc

Re: [R] Snow parLapply

2005-11-11 Thread Sean Davis
On 11/11/05 6:48 AM, "Florent Baty" <[EMAIL PROTECTED]> wrote: > Dear R-user, > > I am trying to use the function 'parLapply' from the 'snow' package > which is supposed to work the same wys as 'lapply' but for a > parallelized cluster of computers. The function I am trying to call in > parallel

Re: [R] different functions on different vector subsets

2005-11-11 Thread Ron Ophir
I thought about other cases but I have to dissagree with you. For logical vector NA is no decision and that should be the results of it. Let's say that -b- is a result of comparison not of -a- to something rather a compsrison of -c- to -d-. In this case NA in the first position is a result of NA i

Re: [R] R: method of moments

2005-11-11 Thread Pornpimol Kay
Hi pals This is the first time for trying with R programe. Could u please suggest me how to practice or else give me some simple homeworks. I will try it. Thank u so much k On 11/11/05, Clark Allan <[EMAIL PROTECTED]> wrote: > > hi all > > does anyone know how one would calculate the covarianc

[R] Error in rowSums : 'x' must be numeric IS SOLVED!

2005-11-11 Thread Illyes Eszter
Dear All, It's Eszter again. Thank you for all, who tried to help me. In this mail I would like to inform you that I could solve my yesterday problem "Error in rowSums(x, prod(dn), p, na.rm) : 'x' must be numeric" only with creating a new datafile. Most probably the problem was that I opene

Re: [R] A Quick and (Very) Dirty Intro to Stats in R

2005-11-11 Thread roger bos
John, Your appendix on linear mixed models does look good and I look forward to reading it, but its 24 pages and Jarretts entire guide is less than 8 pages, so my simple I think he meant short! Thanks, Roger On 11/8/05, John Fox <[EMAIL PROTECTED]> wrote: > > Dear Jarrett, > > > > -Origina

[R] no package 'Matrix' at the repositories

2005-11-11 Thread Spencer Graves
Yesterday, I installed R2.2.0 for Windows [Version 2.2.0 (2005-10-06 r35749)]. Unfortunately, 'install.packages("Matrix")' produced the following message: Warning in download.packages(pkgs, destdir = tmpd, available = available, : no package 'Matrix' at the repositories

Re: [R] no package 'Matrix' at the repositories

2005-11-11 Thread Doran, Harold
I don't have an answer, but also want to point out that whenever I try and download the pdf documentation associated with a package, my Acrobat opens but tells me the file is corrupted. I did this for a random selection of packages and the same problem seems to reoccur. I'm not sure if this is jus

Re: [R] no package 'Matrix' at the repositories

2005-11-11 Thread Sean Davis
On 11/11/05 8:38 AM, "Spencer Graves" <[EMAIL PROTECTED]> wrote: > Yesterday, I installed R2.2.0 for Windows [Version 2.2.0 (2005-10-06 > r35749)]. Unfortunately, 'install.packages("Matrix")' produced the > following message: > > Warning in download.packages(pkgs, destdir = tmpd, available = >

[R] Inputing data from multiple files as time series objects

2005-11-11 Thread Constantine Tsardounis
Hello to everyone,... I am a new R ambitious user. I would like to be the first at my department using R, but I have encountered a difficulty during the last days that I cannot overcome reading help() and searching over the net. Problem: I have multiple files with financial data like the following

Re: [R] no package 'Matrix' at the repositories

2005-11-11 Thread Liaw, Andy
If you are installing from Windows, the current version of Matrix apparently doesn't build automatically on Windows. See: http://cran.r-project.org/bin/windows/contrib/2.3/check/Matrix-check.log (That's for R-devel. There're similar problems on R-2.2.0.) I'm sure Doug is aware of this... Andy

Re: [R] How to find statistics like that.

2005-11-11 Thread John Wilkinson \(pipex\)
Adai, I recently came across the following definition of a statistic which may be relevent to the discussion. John - Beran’s (2003) provocative definition of statistics as “the study of algorithms for data analysis” elevates computational considerations to the forefront of the field. It is a

Re: [R] Inputing data from multiple files as time series objects

2005-11-11 Thread Sean Davis
On 11/11/05 9:07 AM, "Constantine Tsardounis" <[EMAIL PROTECTED]> wrote: > Hello to everyone,... > > I am a new R ambitious user. I would like to be the first at my > department using R, but I have encountered a difficulty during the > last days that I cannot overcome reading help() and searching

Re: [R] Inputing data from multiple files as time series objects

2005-11-11 Thread Michaell Taylor
Assuming all the data are in a subdirectory called "Data" and all the files have the '.txt' extension, you could do something like. files <- list.files('Data') for (file in files){ temp <- ts(read.table(file, header = T)) vname <- sub('.txt','',file) assign(vname,temp,en

[R] Lattice (Trellis) plot margins

2005-11-11 Thread Paul E. Green
Is there an equivalent way to set plotting parameters as in par(mai=c(..)) in the Lattice package? I searched trellis.par.get(), but am not sure which one can be used. Paul E. Green [[alternative HTML version deleted]] __ R-help@stat.math.

[R] problems with for: warnings and segfault

2005-11-11 Thread Ronaldo Reis-Jr.
Hi, I have two problem with a for looping using R Version 2.1.1 (2005-06-20) on a Debian Linux Testing. The first problem: warnings messages Look: > xcoord <- 5 > ycoord <- 5 > indice <- 1 > for(i in c(1:5)) {indice <- indice+1;xcoord[indice] <- xcoord+i; ycoord[indice] <- ycoord } Warning m

[R] A 'sweave' strange problem !!!

2005-11-11 Thread Stéphane Dray
Hello list, I have found a problem (bug?) with Sweave. I hope that someone could help me. Try this little example : \documentclass[a4paper]{article} \title{toto} \author{toto} \begin{document} \maketitle <>= a<-rnorm(1)+1 a plot(1:10) @ patati <>= a @ \end{document} The value of 'a' has c

Re: [R] Inputing data from multiple files as time series objects

2005-11-11 Thread ManuelPerera-Chang
Hi Constantine, I have not tested it but try ... y<-c("IBM","MS","DELL","SIEMENS","SUN") read.data.myfunction<-function(x){ for(i in 1:length(x)) { paste(x[i],sep="") <- ts(read.table(paste(x[i],".txt",sep=""), header = TRUE))

Re: [R] Lattice (Trellis) plot margins

2005-11-11 Thread Berton Gunter
To paraphrase a recent reply by the inimitable Deepayan Sarkar, take a look at trellis.par.get('layout.heights'). Playing with some of those arguments (especially the "padding" ones) may get you what you want. Subject to correction by DS, himself, of course ... -- Bert Gunter Genentech Non-Clinic

[R] curve fitting question

2005-11-11 Thread tom wright
I'd appreciate some direction here. I have a model for a system with two independant variables (i1,i2) and one dependant variable (d). I have experimental data recorded at multiple levels of the dependant variable (x). I need to work out the values for the independant variables that best fit the ex

Re: [R] problems with for: warnings and segfault

2005-11-11 Thread Peter Dalgaard
"Ronaldo Reis-Jr." <[EMAIL PROTECTED]> writes: > Hi, > > I have two problem with a for looping using R Version 2.1.1 (2005-06-20) on > a > Debian Linux Testing. > > The first problem: warnings messages > > Look: > > > xcoord <- 5 > > ycoord <- 5 > > indice <- 1 > > for(i in c(1:5)) {indice

[R] [Fwd: Re: no package 'Matrix' at the repositories]

2005-11-11 Thread Uwe Ligges
[Resend the stuff below since initial one has been blocked from R-help.] I have moved the "old" Matrix_0.98-7.zip to the main repository for the meantime. Best, Uwe Ligges Original Message Subject: Re: [R] no package 'Matrix' at the repositories Date: Fri, 11 Nov 2005 15:45:

Re: [R] problems with for: warnings and segfault

2005-11-11 Thread Ronaldo Reis-Jr.
Em Sex 11 Nov 2005 13:34, Peter Dalgaard escreveu: > > Look at > > xcoord[indice] <- xcoord+i > > The left hand side is a single element, the right hand side is a > vector so I don't think it means what I think you think it means. I'm stupid, of course. The correct is: > xcoord[indice] <- xcoor

Re: [R] problems with for: warnings and segfault

2005-11-11 Thread Prof Brian Ripley
On Fri, 11 Nov 2005, Ronaldo Reis-Jr. wrote: > Hi, > > I have two problem with a for looping using R Version 2.1.1 (2005-06-20) on a > Debian Linux Testing. > > The first problem: warnings messages > > Look: > >> xcoord <- 5 >> ycoord <- 5 >> indice <- 1 >> for(i in c(1:5)) {indice <- indice+1;xc

[R] glm x^2

2005-11-11 Thread Jeffrey Stratford
R-users, I'm having some trouble getting .glm and glm.nb to run a polynomial. I've used x*x and x^2 and neither works. I've checked out the archives and they refer to an archive that's no longer working. I've seen that they use poly() but I'm following up my analysis with cv.glm so I'd prefer

Re: [R] glm x^2

2005-11-11 Thread Sundar Dorai-Raj
I think you want: glm(y ~ x + I(x^2), ...) This is shown as an example on pg 50 of "An Introduction to R" (R-2.2.0-pdf). HTH, --sundar Jeffrey Stratford wrote: > R-users, > > I'm having some trouble getting .glm and glm.nb to run a polynomial. > I've used x*x and x^2 and neither works. I'v

Re: [R] different functions on different vector subsets

2005-11-11 Thread Thomas Lumley
On Fri, 11 Nov 2005, Ron Ophir wrote: > I thought about other cases but I have to dissagree with you. For > logical vector NA is no decision and that should be the results of it. I would say NA is "missing": it means that the result could be any valid value and we don't know which one it is. Th

[R] conditional coloring of image labels

2005-11-11 Thread Jake Michaelson
Hi all, I am interested in plotting a heatmap of a set of genes. I would like the text labels of these genes to be colored red rather than black if they meet certain statistical criteria (using an if statement). I'm not sure how to change individual color labels without changing them all. C

[R] split data into training and testing sets

2005-11-11 Thread Dhiren DSouza
How can I split a dataset randomly into a training and testing set. I would like to have the ability to specify the size of the training set and use the remaining data as the testing set. For example 90% training data and 10% testing data split. Is there a function that will accomplish this?

Re: [R] split data into training and testing sets

2005-11-11 Thread Sundar Dorai-Raj
Dhiren DSouza wrote: > How can I split a dataset randomly into a training and testing set. I would > like to have the ability to specify the size of the training set and use the > remaining data as the testing set. > > For example 90% training data and 10% testing data split. Is there a > f

Re: [R] problems with for: warnings and segfault

2005-11-11 Thread Mike Miller
On Fri, 11 Nov 2005, Ronaldo Reis-Jr. wrote: > Segmentation fault > > This is a R bug or an error in my for function? All seg faults are bugs. Mike __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read t

Re: [R] split data into training and testing sets

2005-11-11 Thread Prof Brian Ripley
On Fri, 11 Nov 2005, Dhiren DSouza wrote: > How can I split a dataset randomly into a training and testing set. I would > like to have the ability to specify the size of the training set and use the > remaining data as the testing set. > > For example 90% training data and 10% testing data split.

Re: [R] problems with for: warnings and segfault

2005-11-11 Thread Peter Dalgaard
Mike Miller <[EMAIL PROTECTED]> writes: > On Fri, 11 Nov 2005, Ronaldo Reis-Jr. wrote: > > > Segmentation fault > > > > This is a R bug or an error in my for function? > > > All seg faults are bugs. Within reason... If users go out of their way to cause havoc, e.g. calling C/Fortran entry poin

[R] repeated values, nlme, correlation structures

2005-11-11 Thread Patrick Giraudoux
Dear listers, My request of last week seems not to have drawn someone's attention. Suppose it was not clear enough. I am coping with an observational study where people's aim was to fit growth curve for a population of young blue tits. For logistic reasons, people have not been capable to numb

Re: [R] conditional coloring of image labels

2005-11-11 Thread Jake Michaelson
On Friday 11 November 2005 10:37 am, Jake Michaelson wrote: I'll clarify a little and hopefully this will make more sense (thanks for the friendly encouragement): Let's say I have 6 samples and am looking at 3 genes, with intensities in a matrix as follows: > genes=cbind(ABC1=c(3,4,4,5,6,3), A

[R] simulation study using AD model builder to fit a GLMM under the binomial probit link

2005-11-11 Thread Roel de Jong
I recently took Dave Fournier up on his offer to evaluate his AD Model builder package (http://otter-rsch.com/admodel.html) when fitting a GLMM under the binomial probit link. I conducted a simulation study in which I drawed 500 samples each containing 1500 observations from the following model

[R] Fwd: Re: conditional coloring of image labels

2005-11-11 Thread Jake Michaelson
-- Forwarded Message -- Subject: Re: [R] conditional coloring of image labels Date: Friday 11 November 2005 1:04 pm From: jim holtman <[EMAIL PROTECTED]> To: Jake Michaelson <[EMAIL PROTECTED]> Use 'mtext': genes=cbind(ABC1=c(3,4,4,5,6,3), ABC2=c(4,3,4,7,7,8), ABC3=c(8,7,8,6,

Re: [R] problems with for: warnings and segfault

2005-11-11 Thread Prof Brian Ripley
On Fri, 11 Nov 2005, Peter Dalgaard wrote: > Mike Miller <[EMAIL PROTECTED]> writes: > >> On Fri, 11 Nov 2005, Ronaldo Reis-Jr. wrote: >> >>> Segmentation fault >>> >>> This is a R bug or an error in my for function? >> >> >> All seg faults are bugs. > > Within reason... If users go out of their w

[R] optim not giving correct minima

2005-11-11 Thread Laura Forsberg White
Hello, I am trying to use optim() on a function involving a summation. My function basically is a thinned poisson likelihood. I have two parameters and in most cases optim() does a fine job of getting the minima. I am simulating my data based on pre specified parameters, so I know what I sh

Re: [R] optim not giving correct minima

2005-11-11 Thread Prof Brian Ripley
Try the other methods that optim() supplies, and try supplying a analytical derviative (which looks easy enough). On a problem like this I would expect to use BFGS with analytical derivatives. If you don't want to do any of that, at least explore the control parameters. It doesn't look to me

Re: [R] no package 'Matrix' at the repositories

2005-11-11 Thread Duncan Murdoch
On 11/11/2005 9:13 AM, Liaw, Andy wrote: > If you are installing from Windows, the current version of Matrix apparently > doesn't build automatically on Windows. See: > > http://cran.r-project.org/bin/windows/contrib/2.3/check/Matrix-check.log > > (That's for R-devel. There're similar problems

Re: [R] no package 'Matrix' at the repositories

2005-11-11 Thread Spencer Graves
Thank you all for your replies and for all your hard work to make R what it is. The wise course for me is probably to use R 2.1.1 when I need the Matrix package until this issue gets fixed. Best Wishes, spencer graves Duncan Murdoch wrote: > On 11/11/2005 9:13 AM

[R] precision of double in R

2005-11-11 Thread Burkhard Heil
Hello Everybody! I'm working on biological sequences and their alignements. I can't get around calculating something in the size of choose(1000,500) or 0.02^250. The result is either INF or 0 in the latter case. Is their any data type in R to calculate these things. Or some other solution? Than

Re: [R] precision of double in R

2005-11-11 Thread Duncan Murdoch
On 11/11/2005 9:55 AM, Burkhard Heil wrote: > Hello Everybody! > > > > I'm working on biological sequences and their alignements. I can't get > around calculating something in the size of choose(1000,500) or 0.02^250. > The result is either INF or 0 in the latter case. Is their any data type in

[R] absolute position in plot()

2005-11-11 Thread ivo welch
dear R wizards: could you please point me into the right direction? I would like to write a general function arrow(x1,y1,x2,y2), which naturally draws an arrow. the basics of writing this function are of course easy. the only complication is that I would like the arrow not to exactly hit (x2,y

Re: [R] absolute position in plot()

2005-11-11 Thread Gabor Grothendieck
See getcell at: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/63579.html which converts between absolute and user coordinates. It should give you insight to allow you to develop this functionality. On 11/11/05, ivo welch <[EMAIL PROTECTED]> wrote: > > dear R wizards: could you please point m

[R] sibling list element reference during list definition

2005-11-11 Thread Paul Roebuck
Can the value of a list element be referenced from a sibling list element during list creation without the use of a temporary variable? The following doesn't work but it's the general idea. > list(value = 2, plusplus = $value+1) such that the following would be the output from str() List of 2