[R] R CMD INSTALL use differenct c++ compiler

2005-07-11 Thread Luke
Dear R Users, When I installed e1071 use R CMD INSTALL, I got configure: WARNING: g++ 2.96 cannot reliably be used with this package. configure: error: Please use a different C++ compiler. But how to let R CMD INSTALL use a different C++ compiler? and which C++ compiler is good? -Luke

[R] svmlight running error

2005-07-18 Thread Luke
X7 X8 X9 X10 X11 X12 X13 X14 X15 X16 X17 1 1 63 1 0 0 145 233 1 1 0 150 0 2.3 1 0 0 1 0 2 0 67 0 1 0 160 286 0 1 0 108 1 1.5 0 1 3 0 1 I wonder what is the possible reason for this error. -Luke __ R-help

[R] string vector indices

2006-04-15 Thread Luke
it. But is there any elegant R function I can use? Thanks, -Luke [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting

[R] interaction terms in formula of lm or glm

2006-04-16 Thread Luke
values of V3, and the formula is y ~ 1 + V1 + V3 + V4, it seems that the foo will give me different in and out-sample predictions. -Luke [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] interaction terms in formula of lm or glm

2006-04-17 Thread Luke
Thanks, Christos. Another relevant question: If I want to include the interaction term consisting of V2 and V3 (they are numeric vectors), should I use: y ~ 1 + V2:V3 or y ~ 1 + I(V2*V3) or both are good? -Luke On 4/17/06, Christos Hatzis [EMAIL PROTECTED] wrote: If you want the quadratic

[R] how to use list index to get vector

2005-05-16 Thread Luke
, the first elements of list elements? That is, how to get a vector c(1, 11, 111) from foo? foo[[]][1] doesn't work. -Luke __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R

Re: [R] how to use list index to get vector

2005-05-16 Thread Luke
, [, 1) myIndex [1] 2 5 Is this the simplest way to get the index or value vector? -Luke On 5/16/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: Its the indexing function written in ordinary function form. That is, foo[1:2] can be written as [(foo, 1:2) On 5/16/05, Luke [EMAIL PROTECTED

[R] get the bug context

2005-05-29 Thread Luke
which line of code caused the bug. Is there any utility or something in R which can let you know which line of code causes a bug? Thanks, -Luke __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] AIC in glm.fit with intercept

2005-06-15 Thread Luke
= identity), intercept=FALSE) foo$aic [1] 38.94657 Is my code wrong or is it a bug in glm.fit? -Luke __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting

Re: [R] AIC in glm.fit with intercept

2005-06-15 Thread Luke
More interesting: lm() works fine to get AIC. I wonder if there is a bug in glm.fit(). AIC(lm(myY~0+myX, data=data.frame(myY,myX))) [1] 38.94657 AIC(lm(myY~1+myX, data=data.frame(myY,myX))) [1] -650.9808 -Luke On 6/15/05, Luke [EMAIL PROTECTED] wrote: Dear R users, glm.fit() gave me

[R] quick way to construct formula

2005-06-28 Thread Luke
know I can use several paste functions to construct it. But is there any other handy way to do it? -Luke __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org

Re: [R] Chemoinformatic people

2005-07-21 Thread Luke Tierney
[EMAIL PROTECTED] luke On Thu, 21 Jul 2005, A.J. Rossini wrote: Just with R, or via another tool integrating R, such as Pipeline Pilot? best, -tony On 7/20/05, Frédéric Ooms [EMAIL PROTECTED] wrote: Dear colleague, Just an e-mail to know if they are people working in the field of chemoinformatic

Re: [R] Broken tkrplot on Mac OS X

2005-08-16 Thread Luke Tierney
are offended by a .so extension on a dylib you can change the extension and the R code that load it. If you figure out a better way to do this please let me know, Best, luke On Thu, 12 May 2005, stefano iacus wrote: It seems it doesn't work for me either http://159.149.213.137/R/bin/2.1/check

[R] horizontal direct product

2006-08-25 Thread Luke Keele
short of writing my own function? Thanks Luke Luke Keele Department of Political Science Ohio State University [EMAIL PROTECTED] [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] increasing the # of socket connections

2006-08-25 Thread Luke Tierney
asked that question if the answer had been options(max.connections=100)? I for one would not. Best, luke -- Luke Tierney Chair, Statistics and Actuarial Science Ralph E. Wareham Professor of Mathematical Sciences University of Iowa Phone: 319-335-3386 Department

Re: [R] makeSOCKcluster

2006-09-01 Thread Luke Tierney
it by hand and see if you can work out why it hangs. Best, luke ssh and all are working fine, When I try to create a socket connection as makeSOCKcluster(c(localhost,localhost)), it just waits for the other prcess on localhost to get created and respond. But this other process

Re: [R] Evaluation of defaults in functions

2006-09-29 Thread Luke Tierney
code is in http://www.stat.uiowa.edu/~luke/R/lazy/. This is out of date but not too hard to fix. Insuring memory efficiency is still a bit tricky; some comments are int he notes at this site. Best, luke -- Luke Tierney Chair, Statistics and Actuarial Science Ralph E. Wareham Professor

Re: [R] Evaluation of defaults in functions

2006-09-29 Thread Luke Tierney
evaluation and they are equivalent in Haskell with pure lazy evaluation but they are not in R. So in R replacing one by the other is not a semantics-preserving transformation. Best, luke Lazy data structures can be implemented in R on top of the limited lazy evaluation mechanism. I experimented

Re: [R] Rmpi performance

2006-10-13 Thread Luke Tierney
and you really are getting gigabit performance. I would look for a visualization tool an idea of what is happening--perhaps xmpi if your MPI is LAM. Best, luke with the cluster. Can you give us some advice in order to understand why the cluster is slower than the master? Thank you very much

[R] Problems with boot and optim

2004-09-21 Thread Luke Keele
The code for the likelihood and the call to boot is below. I have tried numerous other permutations as well. I am using R 1.9.1 on Windows XP pro. Thanks Luke Keele #Define Likelihood lik.hetprobit -function(par, X, Y, Z){ #Pull Out Parameters Y - as.matrix(y) X - as.matrix(x) Z

Re: [R] Errors while compiling packages with namespace?

2004-10-17 Thread Luke Tierney
The error message points to line 21 of the NAMESPACE file: S3method(view, function) The NAMESPACE file is parsed by the R parser, so this is a suntax error since function is a reserved word. Put quotes around it and you should be OK. luke On Sun, 17 Oct 2004, Philippe Grosjean wrote: Hello, I try

Re: [R] Parallel computations using snow: how to combine boot objects?

2005-01-25 Thread Luke Tierney
The example ducks that issue. Someone needs to write a function for merging these results. Probably just involves making a suitable call to boot.return, which is what happens at the end of boot(), but I don't know if anyone has actually done this yet. luke On Fri, 21 Jan 2005, BEER Michael

Re: [R] socket problems (maybe bugs?)

2005-02-19 Thread Luke Tierney
. With the current implementation this ensures that the read completely empties the buffer and so this problem does not bite. It sounds like the same stategy should allow you to work with the current implementation. Best, luke -- Luke Tierney University of Iowa Phone: 319-335-3386

Re: [R] socket problems (maybe bugs?)

2005-02-19 Thread Luke Tierney
On Sat, 19 Feb 2005, Luke Tierney wrote: On Thu, 17 Feb 2005, Christian Lederer wrote: Dear R Gurus, for some purpose i have to use a socket connection, where i have to read and write both text and binary data (each binary data package will be preceeded by a header line). When experimenting, i

RE: [R] Do environments make copies?

2005-02-24 Thread Luke Tierney
= .GlobalEnv)) in new.ref avoids the second copy both in the saved image and after loading. luke -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Nawaaz Ahmed Sent: Thursday, February 24, 2005 10:36 AM To: r-help@stat.math.ethz.ch Subject: [R] Do environments

Re: [R] tkRplot under macos x

2005-02-28 Thread Luke Tierney
It only works under X11; if that is what you want this posting should tell you what needs to be done: https://stat.ethz.ch/pipermail/r-sig-mac/2004-December/001465.html Best, luke On Mon, 28 Feb 2005, Sean Davis wrote: I have not successfully gotten tkRplot to install on macos 10.3.8, R

Re: [R] workload of R

2005-03-16 Thread Luke Tierney
efforts underway for compiling R and also for parallel computing in R. If you could ask a clearer question we might be able to help and cooperate. luke -- Luke Tierney Chair, Statistics and Actuarial Science Ralph E. Wareham Professor of Mathematical Sciences University of Iowa

[R] [R-pkgs] misc3d_0.4-0 now available on CRAN

2006-06-06 Thread Luke Tierney
to represent 3D structure. The examples in example(contour3d) example(parametric3d) demo(lighting) demo(teapot) illustrate some of the features. Dai Feng and Luke Tierney -- Luke Tierney Chair, Statistics and Actuarial Science Ralph E. Wareham Professor of Mathematical Sciences

Re: [R] How to get the namespace of a function?

2006-02-03 Thread Luke Tierney
analysis becomes a standard part of R, either via byte code compilation or as bart of R CMD check, there should be warnings issued automatically for this sort of thing. For now manually using the codetools ackage Brian pointed to is one option. Best, luke Duncan Murdoch Of course, one can

[R] Point Size in plot.cox.zph

2006-02-09 Thread Luke Keele
the points command as follows: plot(time, war.zph$y[,7], type=n) points(war.zph[7], ann=F, pch = .) But this causes you to lose the spline fit and seems rather cumbersome. A way to do this using the standard plot.cox.zph framework would be very helpful. Thanks Luke -- Luke Keele Assistant

Re: warnings and tryCatch() (Was: RE: [R] catching the warnings)

2004-06-03 Thread Luke Tierney
the default mechanism of the signaling function is used. For warnings the default handling involves storing the message and then continuing; for errors signaled with stop or the internal error function it usually involves printing a message and jumping to top level. Hope that helps. luke Cheer Henrik

Re: [R] nls and R scoping rules

2004-06-10 Thread Luke Tierney
the environment of the functions in the structure). Since environment sharing is preserved within a serialization, this means the change would only alter the parent of the env environment, not add anything extra. Best, luke -- Luke Tierney University of Iowa Phone

Re: [R] tryCatch() and preventing interrupts in 'finally'

2004-06-10 Thread Luke Tierney
code. There are some nasty little issues involved in getting this right, which is why it hasn't happened yet. Best, luke On Fri, 11 Jun 2004, Henrik Bengtsson wrote: With tryCatch() it is possible to catch interrupts with tryCatch(). Then you can use a 'finally' statement to clean up, release

Re: [R] mkChar can be interrupted

2004-06-14 Thread Luke Tierney
other things, we will need to make changes as we improve support for other event loops. [In applications with graphical interfaces signals are not the right way to deal with user interruption (in particular on operating systems that don't support proper signals)]. Best, luke On Mon, 14 Jun 2004

RE: [R] mkChar can be interrupted

2004-06-14 Thread Luke Tierney
On Mon, 14 Jun 2004, Vadim Ogranovich wrote: -Original Message- From: Luke Tierney [mailto:[EMAIL PROTECTED] Sent: Monday, June 14, 2004 1:30 PM To: Vadim Ogranovich Cc: R-Help Subject: Re: [R] mkChar can be interrupted Not sure why you think this suggest mkChar can

RE: [R] mkChar can be interrupted

2004-06-14 Thread Luke Tierney
this was changed. Unless there is explicit documentation to the contrary you should assume that every function in the R API might allocate and might cause a non-local exit (i.e. a longjmp) when an exception is raised (and an interrupt is one of, but only one of, the exceptions that might occur). luke

RE: [R] mkChar can be interrupted

2004-06-14 Thread Luke Tierney
, which will result in a non-local exit. If the allocation succeeds, you are responsable for calling Free. luke -Original Message- From: Luke Tierney [mailto:[EMAIL PROTECTED] Sent: Monday, June 14, 2004 5:43 PM To: Vadim Ogranovich Cc: R-Help Subject: RE: [R] mkChar can

Re: [R] interpreting profiling output

2004-07-26 Thread Luke Tierney
. luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial Science 241 Schaeffer Hall email: [EMAIL PROTECTED] Iowa City, IA 52242 WWW: http

Re: [R] static analysis tools for R code?

2005-04-11 Thread Luke Tierney
There are some preliminary tools available in the codetools package at http://www.stat.uiowa.edu/~luke/R/codetools/ Hopefully these will be cleaned up and released via CRAN or incorporated into R this summer. luke On Mon, 11 Apr 2005, Dimitris Rizopoulos wrote: If I understand well what

RE: [R] when can we expect Prof Tierney's compiled R?

2005-04-20 Thread Luke Tierney
, and then (maybe) optimizing by rewriting R code or coding core bits in C or Fortran is likely to remain the best strategy for a long time to come. Best, luke On Mon, 18 Apr 2005, Liaw, Andy wrote: Do you mean the byte code compiler? You can find it at: http://www.stat.uiowa.edu/~luke/R/compiler/ Andy From

RE: [R] when can we expect Prof Tierney's compiled R?

2005-04-20 Thread Luke Tierney
code, which may have other benefits beyond just speed improvements. Best, luke On Wed, 20 Apr 2005, Jason Liao wrote: Dear Prof. Tierney, Thank for very much for replying and we all appreciate what you have done for the R community. Currently I have been spoiled by R. I would love to be 100% R. I

Re: [R] when can we expect Prof Tierney's compiled R?

2005-04-22 Thread Luke Tierney
On Wed, 20 Apr 2005, Peter Dalgaard wrote: Luke Tierney [EMAIL PROTECTED] writes: Vectorized operations in R are also as fast as compiled C (because that is what they are :-)). A compiler such as the one I'm working on will be able to make most difference for non-vectorizable or not very

[R] load data with scientific notations

2005-05-12 Thread Luke Zhou
Hi, I have a data file whose values are scientific notations, for example: 9.1642537e+003, 6.7195295e+003, 4.8834487e+003, 3.7181589e+003, ... How to load such data into R? Thanks, -Luke __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] load data with scientific notations

2005-05-12 Thread Luke Zhou
It works. thanks! -Luke On 5/12/05, Liaw, Andy [EMAIL PROTECTED] wrote: Just read them like any other numbers; e.g.: x - scan(sep=,) 1: 9.1642537e+003, 6.7195295e+003, 4.8834487e+003, 3.7181589e+003 5: [hit ctrl-d] Read 4 items x [1] 9164.254 6719.529 4883.449 3718.159 Andy

Re: [R] A model for disease progression

2003-07-31 Thread Luke Whitaker
and the population demographics are in steady state. I would think that in many real world disease/population cases this will be a very poor assumption. I'm afraid I can't see any obvious way around this other than to look for some sort of cohort/followup data. Regards, Luke Whitaker

Re: [R] R on Linux/Opteron?

2003-08-26 Thread Luke Tierney
I've had a chance to build R-devel on one running SuSE (not sure which release). confugure set things up for a 64-bit build that passed all tests. The base tests ran about 25% faster on the 1.4GHz opteron than a2GHz Xeon. That's as much as I know at this point. luke On Tue, 26 Aug 2003, Liaw

Re: [R] R on Linux/Opteron?

2003-08-27 Thread Luke Tierney
: http://lists.debian.org/debian-x86-64/2003/debian-x86-64-200308/threads.html whereas they have had ia64 out for a while. SuSE has an Opteron option and Luke said he tried it. Apparently it has a functioning 64-bit compiler toolchain - I weren't sure earlier whether they were just running a 64bit

[R] Probit and optim in R

2003-09-08 Thread Luke Keele
if anyone wants to take a look. I'd like to be able to make it work without the gradient if possible. Any help would be greatly appreciated as I am completely stumped at this point. Luke Keele Dept of Political Science UNC-Chapel Hill ## Probit Code For Simulation #Define empty matrix na

Re: [R] Help needed: plotting with no device

2003-09-29 Thread Luke Tierney
Look at ?Signal and see if that would help. You may need to rename the appropriate .RData beforehand to be safer. luke On Mon, 29 Sep 2003, Ross Boylan wrote: On Mon, Sep 29, 2003 at 05:09:49PM -0400, Ben Bolker wrote: Can you use save.image() to rescue your results? I would try

Re: [R] query on proxy settings for R.

2003-11-05 Thread Luke Tierney
/pipermail/r-help/2003-July/035410.html luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial Science 241 Schaeffer Hall email: [EMAIL PROTECTED] Iowa City, IA 52242

[R] (no subject)

2003-11-10 Thread Luke Keele
someway to reference which plot I want to identify as far as I can tell. If I just do one plot at a time identify works fine. Is there any other way to just label all the data points? Thanks Luke Keele Post-Doctoral Fellow in Quantitative Methods Nuffield College, Oxford University

Re: [R] signalCondition

2003-11-10 Thread Luke Tierney
will transfer control to that handler and the default code is not reached. Hope that helps, luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial Science 241 Schaeffer Hall

Re: [R] Persistent state of R

2003-11-25 Thread Luke Tierney
the names of the variables defined. The actual values are only loaded from the data base on demand. An experimental package that implements this is available at http://www.stat.uiowa.edu/~luke/R/serialize/lazyload.tar.gz I believe `make check-all' passes with all base and recommended packages

[R] RWinEdt

2003-12-02 Thread Luke Keele
is correctly named as well that is it is named: SWinRegistry_0.3-3.zip, which I know caused problems with earlier versions. Thanks Luke Keele Post-Doctoral Fellow in Quantitative Methods Nuffield College, Oxford University Oxford, UK [[alternative HTML version deleted

RE: [R] performance gap between R 1.7.1 and 1.8.0

2003-12-05 Thread Luke Tierney
use a bit of tuning, but minor compared to the main issue here. Luke On Thu, 4 Dec 2003, Liaw, Andy wrote: From: Martin Maechler [mailto:[EMAIL PROTECTED] [...] A very first step of diagnosis might be to activate trace(read.dcf) trace(library) options(verbose = TRUE) A step

[R] Re: R-help digest, Vol 1 #34 - 2 msgs

2003-01-02 Thread Luke Whitaker
os linux-gnu system i386, linux-gnu status major1 minor5.1 year 2002 month06 day 17 language R Thanks, Luke Whitaker

[R] R talking to Oracle, ODBC drivers available ?

2003-01-03 Thread Luke Whitaker
17 language R Thanks, Luke Whitaker. __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] interaction of options(error=) and try()

2003-01-07 Thread Luke Tierney
option is to disable your error handler around the try. luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial Science 241 Schaeffer Hall email: [EMAIL PROTECTED

[R] DBI/ROracle for remote database connection ?

2003-01-17 Thread Luke Whitaker
, linux-gnu status major1 minor6.2 year 2003 month01 day 10 language R Thanks, Luke Whitaker __ [EMAIL PROTECTED] mailing

Re: [R] Memory allocation, IBM-AIX and R-1.6.2 - addendum

2003-02-11 Thread Luke Tierney
will be different in detail but the basic issue is the same: 2^32 = 4G only leave you so much room to do the things the OS needs to do. luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial

Re: [R] data manipulation function descriptions

2003-02-14 Thread Luke Tierney
justifialy proud as well as ...] luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial Science 241 Schaeffer Hall email: [EMAIL PROTECTED] Iowa City, IA 52242

Re: [R] R (external ?) reference

2003-02-28 Thread Luke Tierney
On Fri, 28 Feb 2003, Laurent Gautier wrote: Dear List, I found a documentation on the web that mentions things like 'R references' (http://www.stat.uiowa.edu/~luke/R/simpleref.html). However, I could not find the R_MakeReference and friends in R... Does anyone knows more about

Re: [R] Lexical scoping question

2003-02-28 Thread Luke Tierney
of the evaluation befoe the change, then you have to force the evaluation to occur before the change. You can do force-function(z) z f - function(y) { force(y) g - function(x) x + y g } which will work as you expect. IIRC Luke Tierney has added force

Re: [R] getMethod in a setMethod context

2003-03-09 Thread Luke Tierney
for A calls generics for which the class of the original obj would be more relevant. luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial Science 241 Schaeffer Hall

[R] SNOW: Simple Parallel Computing in R

2003-03-10 Thread Luke Tierney
, one using PVM via the rpvm package, and one using MPI, via the Rmpi package. The page http://www.stat.uiowa.edu/~luke/R/cluster/cluster.html provides a brief introduction and a simple example of using the cluster for parallel bootstrapping. A paper on parallel computing in R using snow is available

Re: [R] VIM Syntax Highlighting

2003-03-11 Thread Luke Whitaker
with vim. Version 5.6.70 of vim comes with a syntax file for S-lang as standard. Either pick S-lang from the Syntax menu or do :cal SetSyn(slang) on the vim command line. Not surprisingly, this works perfectly well with either S or R. Luke Whitaker

Re: [R] namespaces not available

2003-06-20 Thread Luke Tierney
there is something more complicated going on. Does a traceback() tell you anything? Otherwise, can you send me your .RData file and I will try to track this down. luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax

RE: [R] R won't connect to the internet on SUSE Linux 8.1

2003-07-26 Thread Luke Tierney
to resolve... line is missing. Best, luke On Sat, 26 Jul 2003, michael watson (IAH-C) wrote: :-)) Installing Bioconductor was how it all began, so I ended up doing what you suggested (in fact I downloaded just the packages I needed as in the full tar ball, rhdf5 wouldn't compile, probably as I

Re: [R] Implementating streams in R

2004-02-03 Thread Luke Tierney
I worked on this a bit a while back for a possible article for Rnews that won't get written anytime soon. I've put a snapshot in http://www.stat.uiowa.edu/~luke/R/lazy/ It is based on examples in Paulson's ML book and Abelson and Susman; the overall design seems similar to the one you

Re: [R] Porting let* from Common LISP to R

2004-02-12 Thread Luke Tierney
to increase that at the OS level. R's lazy evaluation of arguments also causes different stack usage patterns in recursive code; sometimes more stack is needed, sometimes less--depends on the code. Hope that helps. luke On Thu, 12 Feb 2004, Gabriel Baud-Bovy wrote: In porting some Common LISP code

Re: [R] Status of Rmpi

2004-03-23 Thread Luke Tierney
hopefully have a chance to give a newer version of Rmpi a go. Best, luke On Tue, 23 Mar 2004, Ross Boylan wrote: Rmpi is not currently available on CRAN, and I don't think it has been for a few months. It is available at http://www.stats.uwo.ca/faculty/yu/Rmpi/ Does anyone know its current

Re: [R] snow documentation comments

2004-03-24 Thread Luke Tierney
that the master process uses. So there are a total of three processes. luke -- Luke Tierney University of Iowa Phone: 319-335-3386 Department of Statistics andFax: 319-335-3017 Actuarial Science 241 Schaeffer Hall email

Re: [R] socket clusters on snow dies easily

2004-04-08 Thread Luke Tierney
it more solid than the socket version and likely to have fewer issues. Best, luke On Fri, 9 Apr 2004, Rong-En Fan wrote: hello, I'm using R 1.8.1 with the lastest snow package on FreeBSD 4.9. However, when I try to using socket clusters, it's very unstable. Sometimes it dies half way

Re: [R] looking for a tutorial on exception handling

2004-04-10 Thread Luke Tierney
Unfortunately there isn't yet. Writing a paper on this is fairly high on my list of priorities but hasn't happened yet. Best, luke On Fri, 9 Apr 2004, Tamas Papp wrote: Hi, I would like to learn how to do exception handling in R. I had a look at the help page for tryCatch and I am

Re: [R] repeat { readline() }

2006-01-09 Thread Luke Tierney
Use tryCatch; try behaves the way it does with respect to interrupts for historical compatibility. luke On Sun, 8 Jan 2006, hadley wickham wrote: On a related note, does anyone know how to exit: repeat { try( readline() ) } The try block captures Ctrl-C. Hadley

Re: [R] snow's makeCluster hanging (using Rmpi)

2006-11-07 Thread Luke Tierney
The most likely culprit is firewall settings. Something like tcpdump may help to confirm that. Working with a stand-alone example from Rmpi may also help. Best, luke On Tue, 7 Nov 2006, Randall C Johnson [Contr.] wrote: Hello everyone, I've been fiddling around with the snow and Rmpi

Re: [R] wrapping mle()

2006-12-30 Thread Luke Tierney
in very rare circumstances. Best, luke On Sat, 30 Dec 2006, Gabor Grothendieck wrote: Add the line marked ### so that the environment of loglik.fun is reset to the environment within fit.mle so that it can find y there: library(stats4) ll - function(ymax=15, xhalf=6) { -sum(stats::dpois(y

Re: [R] wrapping mle()

2007-01-03 Thread Luke Tierney
On Sat, 30 Dec 2006, Sebastian P. Luque wrote: On Sat, 30 Dec 2006 15:46:01 -0600 (CST), Luke Tierney [EMAIL PROTECTED] wrote: It is much cleaner to do this sort of thing with lexical scope. For example, mkll - function(x, y) { function(ymax=15, xhalf=6) { -sum

Re: [R] wrapping mle()

2007-01-03 Thread Luke Tierney
functions in R. Environment surgery in contrast is messy and complex and essentially impossible to get right. If you want to do this in the privacy of your own code that is fine, but please don't encourage others to go down this path. Best, luke By the way, here is one additional solution

Re: [R] features of save and save.image (unexpected file sizes)

2007-02-02 Thread Luke Tierney
: Only a marker representing .GlobalEnv is saved; on load into a new session that marker becomes the .GlobalEnv of the new session. Best, luke The R way is to create what you want, not fix up afterwards. (I find your code unreadable--spaces help a great deal, so am not sure if I have

Re: [R] how much performance penalty does this incur, scalar as a vector of one element?

2007-02-22 Thread Luke Tierney
, is in many ways conceptually cleaner and might in hindsight have been a better choice for that reason, but efficiency isn't really a consideration. Best, luke On Wed, 21 Feb 2007, Jason Liao wrote: I have been comparing R with other languages and systems. One peculiar feature of R

Re: [R] how much performance penalty does this incur, scalar as a vector of one element?

2007-02-23 Thread Luke Tierney
Thanks -- that's good to know. Best, luke On Fri, 23 Feb 2007, Jeffrey J. Hallman wrote: The 64 bit version of VisualWorks Smalltalk has an immediate ShortDouble, which sacrifices two bits of exponent for a tag. It thus has the same precision as an IEEE double, but one fourth as much range

Re: [R] Scoping issue?

2007-03-05 Thread Luke Tierney
in this case is the expression consisting of the single variable `colnum` and is not useful here). You need to use another function in your sapply call, something like function(d) d[,colnum] may do. Best, luke On Sun, 4 Mar 2007, Thaden, John J wrote: Apparently you're right that colnum doesn't

Re: [R] question on suppressing error messages with Rmath library

2007-03-21 Thread Luke Tierney
in R. Best, luke On Wed, 21 Mar 2007, Ranjan Maitra wrote: Dear list, I have been using the Rmath library for quite a while: in the current instance, I am calling dnt (non-central t density function) repeatedly for several million. When the argument is small, I get the warning message

Re: [R] Snow Package and SPRNG: Will it solve my problem?

2007-03-22 Thread Luke Tierney
if necessary) independent streams on the nodes. Best, luke Thanks again for your continued help with my problems. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

Re: [R] can't load just saved R object ReadItem: unknown type 65

2007-03-23 Thread Luke Tierney
According to the logs nothing at all has changed in the serialization code in a month and nothing of consequence for much longer than that. To track this down we will need a complete, reproducible, and preferably minimal example. Best, luke On Thu, 22 Mar 2007, Mark W Kimpel wrote: I have run

Re: [R] snow parLapply standard output

2007-03-27 Thread Luke Tierney
. Best, luke Thanks in advance for any hints, Frank Preiswerk [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] Subject: useR! 2007 --- programming competition

2007-04-13 Thread Luke Tierney
R Users and Developers, The first North American useR! to be held at Iowa State University, Ames, Iowa, August 8-10, 2007, will include a programming competition. The challenge: Develop a package useful for the analysis of large data sets. Your package could augment an existing

[R] Specifying ANCOVA models in R

2007-04-18 Thread Luke Spadavecchia
in R, both for consistency with my other analysis, and because I use R on my home machine, and have to venture into the university labs to use MINITAB. Many thanks Luke Spadavecchia [[alternative HTML version deleted]] __ R-help

Re: [R] Parallel processes

2007-05-22 Thread Luke Tierney
On Tue, 22 May 2007, hodgess wrote: Dear R People: I was wondering if there were any packages for parallel programming in R. According to the R-help, at one time there was a package called SNOW. It doesn't seem to exist anymore. ??? It's available on CRAN. Best, luke Any help

Re: [R] simultaneous computing

2007-06-11 Thread Luke Tierney
The package snow available from CRAN is one possibility. Best, luke On Mon, 11 Jun 2007, Markus Schmidberger wrote: Hello, which possibilities are available in R for simultaneous or parallel computing? I only could find biopara (http://cran.r-project.org/src/contrib/Descriptions

Re: [R] How to properly finalize external pointers?

2007-08-06 Thread Luke Tierney
with R_ClearExternalPtr in the explicit close routine and checking whether it has already been set to NULL in the finalizer. Best, luke Best regards Jens Oehlschlägel // C-code static void rindex_finalize(SEXP extPtr){ pINT ptr = R_ExternalPtrAddr(extPtr); if(ptr){ Free(ptr); Rprintf

Re: [R] makeSOCKcluster

2007-08-14 Thread Luke Tierney
Try using the option homogenoeus=FALSE and make sure the appropriate environment variables are set on the worker nodes. Best, luke On Tue, 14 Aug 2007, Michael Janis wrote: Hi, I am attempting to implement a mixed (windows/linux) snow sockets parallelism in R, but am running

Re: [R] Suspected memory leak with R v.2.5.x and large matrices with dimnames set

2007-08-18 Thread Luke Tierney
Seth mentiones are not likely to help in this case. They are primarily intended to improve performance when there are many non-unique character vectors; there is additional overhead for many unique vectors, which we will try to reduce over time. Best, luke On Sat, 18 Aug 2007, Peter Waltman wrote

Re: [R] clusterCall with replicate function

2007-08-21 Thread Luke Tierney
with that -- just install the rlecuyer package and call clusterSetupRNG. Best, luke On Tue, 21 Aug 2007, Michael Gormley wrote: I am trying to run a monte carlo process using snow with a MPI cluster. I have ~thirty processors to run the algorithm on and I want to run it 5000 times and take

[R] [R-pkgs] proftools package now available from CRAN

2007-08-27 Thread Luke Tierney
to be worse). More careful design of the data structures and memoizing calculations that are now repeated is likely to improve performance substantially. -- Luke Tierney Chair, Statistics and Actuarial Science Ralph E. Wareham Professor of Mathematical Sciences University of Iowa

[R] sneding my email

2007-08-29 Thread Luke Neraas
hello, I would like ot post questions to this list so i am sending my email Thanks Luke Neraas [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] How to multiply all dataframe rows by another dataframe's columns

2007-08-29 Thread Luke Neraas
Luke Neraas [EMAIL PROTECTED] [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide