Re: [R] socket clusters on snow dies easily

2004-04-09 Thread Rong-En Fan
On Thu, Apr 08, 2004 at 11:18:43AM -0500, Luke Tierney wrote: If you absolutely have to use the socket approach you will probably have to track down the particular reason for this failure in your environment by looking at what is going on under the hood and checking the pieces. The master

[R] bootstrap function coefficients

2004-04-09 Thread Michaël Coeurdassier
Dear R community, Please, can you help me with a problem concerning bootstrap. The data table called «RMika», contained times (Tps) and corresponding concentration of a chemical in a soil (SolA). I would like to get, by bootstraping, 10 estimations of the parameters C0 and k from the function:

[R] generating 00Index.html

2004-04-09 Thread Nolwenn Le Meur
Hi, I've started to create the documentation of my package. I've generate a pdf and the different html files from .Rd files. However I don't know how to automatically generate the 00Index.html. Do I have to write one in latex style and convert it or can I use the CONTENTS, INDEX, or Rd files.

Re: [R] bootstrap function coefficients

2004-04-09 Thread Prof Brian Ripley
Package bootstrap is ORPHANED (that is, unsupported by anyone). Do try package boot instead. 10 is a very small number of bootstrap replications: I do suggest you try 100 or more. On Fri, 9 Apr 2004, Michaël Coeurdassier wrote: Dear R community, Please, can you help me with a problem

Re: [R] generating 00Index.html

2004-04-09 Thread Prof Brian Ripley
It happens automatically when you INSTALL the package sources ... sounds as if you are not actually doing that. On Fri, 9 Apr 2004, Nolwenn Le Meur wrote: I've started to create the documentation of my package. I've generate a pdf and the different html files from .Rd files. However I don't

Re: [R] nlme on Windows 2000 (v1.8.1)

2004-04-09 Thread Prof Brian Ripley
See the rw-FAQ, specifically Q2.19. It is very likely some piece of rogue software on your machines. On Thu, 8 Apr 2004, Kent Holsinger wrote: I have a problem with nlme on Windows 2000, and I'm having a devil of a time determining whether the problem is with my computer or with something

Re: [R] nlme on Windows 2000 (v1.8.1)

2004-04-09 Thread Duncan Murdoch
On Thu, 08 Apr 2004 08:03:56 -0400, you wrote: If I use Rterm, I can run analyses with NLME to my heart's content. But when I run Rgui, I encounter a floating point exception. I agree with Prof Ripley: I'd suspect video (or other) drivers messing up the floating point control word.

Re: [R] %+=% and eval.parent()

2004-04-09 Thread Prof Brian Ripley
On 7 Apr 2004, Peter Dalgaard wrote: Robin Hankin [EMAIL PROTECTED] writes: Hi again everybody. Yesterday I had a problem with c-style += functions. One suggestion was to define R plus- - function(x,value){x+value} Then R a - matrix(1:9,3,3) R plus(a[a%%2==1]) - 1000

[R] Question

2004-04-09 Thread Ivan Yegorov
Could you please help me. How can I get indexes of matrix elements equal to the specified value. Thanks in advance. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] looking for a tutorial on exception handling

2004-04-09 Thread Tamas Papp
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 trying to understand it (I am sure it is well-written, but I have not used exception handling before in any language, so I have no experience in that area). Is there a tutorial or an

[R] R-estimators

2004-04-09 Thread Paolo Radaelli
Dear all, could you please suggest me the package I have to install in order to use ranks statistics in regression models (R-Estimators, wilcoxon scores ...) Thanks and Happy Easter Paolo Radaelli [[alternative HTML version deleted]] __ [EMAIL

[R] RODBC on FreeBSD 4.9

2004-04-09 Thread Charles and Kimberly Maner
Hello. I have found a number of posts/threads on this subject and remain unsuccessful to determine the solution. In brief, my situation is: 1. Running FreeBSD 4.9 Stable in x86 2. Installed, and running, R 1.8.1 3. Installed, and running, postgreSQL 7.4.2 4. Installed, and running,

Re: [R] %+=% and eval.parent()

2004-04-09 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: This *should* be somewhere in the R Language Definition, although I'm not sure it is actually there. Or the blue book, although I suspect that S v.3 actually used positional matching (and jumped through hoops). It did. This is in the FAQ,

Re: [R] nlme on Windows 2000 (v1.8.1)

2004-04-09 Thread Kent Holsinger
Duncan Murdoch wrote: On Thu, 08 Apr 2004 08:03:56 -0400, you wrote: If I use Rterm, I can run analyses with NLME to my heart's content. But when I run Rgui, I encounter a floating point exception. You do this by hitting F8 during the boot sequence. R can run fine in safe mode; if you don't

[R] Cross-compiling packages for (ARM) Zaurus ideas and tips?

2004-04-09 Thread Simon Pickering
Hello All, A bit of a long one (background mainly), this is really just a request for some ideas. I ported R 1.7.x for use on the Sharp Zaurus ( iPAQ, simpad, yopy, etc.) some time ago and have just ported 1.8.1 but had never quite realised that you actually need to install some packages to make

RE: [R] Question (on matrix indexes)

2004-04-09 Thread BXC (Bendix Carstensen)
I am not sure what you mean but you might be interested in the functions row() and col(). Bendix Carstensen -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43

[R] loess' robustness weights in loess

2004-04-09 Thread Rafael A. Irizarry
hi! i want to change the robustness weights used by loess. these are described on page 316 of chambers and hastie's statistical models in S book as r_i = B(e_i,6m) where B is tukey's biweight function, e_i are the residulas, and m is the median average distance from 0 of the residuals. i

Re: [R] loess' robustness weights in loess

2004-04-09 Thread Prof Brian Ripley
On Fri, 9 Apr 2004, Rafael A. Irizarry wrote: hi! i want to change the robustness weights used by loess. these are described on page 316 of chambers and hastie's statistical models in S book as r_i = B(e_i,6m) where B is tukey's biweight function, e_i are the residulas, and m is the

[R] question regarding vector ops

2004-04-09 Thread Rajarshi Guha
Hi, I was playing with some code and came upon this situation. x - c(1, rep(0,9)) x [1] 1 0 0 0 0 0 0 0 0 0 idx - 1:10 x[idx] - x[idx] + 1 x [1] 2 1 1 1 1 1 1 1 1 1 This is expected. But if I do: x - c(1, rep(0,9)) x [1] 1 0 0 0 0 0 0 0 0 0 idx - rep(0,10) idx [1] 0 0 0 0 0 0 0 0

Re: [R] question regarding vector ops

2004-04-09 Thread Uwe Ligges
Rajarshi Guha wrote: Hi, I was playing with some code and came upon this situation. x - c(1, rep(0,9)) x [1] 1 0 0 0 0 0 0 0 0 0 idx - 1:10 x[idx] - x[idx] + 1 x [1] 2 1 1 1 1 1 1 1 1 1 This is expected. But if I do: x - c(1, rep(0,9)) x [1] 1 0 0 0 0 0 0 0 0 0 idx - rep(0,10) idx [1]

Re: [R] question regarding vector ops

2004-04-09 Thread Rajarshi Guha
On Fri, 2004-04-09 at 12:03, Rajarshi Guha wrote: Hi, I was playing with some code and came upon this situation. x - c(1, rep(0,9)) x [1] 1 0 0 0 0 0 0 0 0 0 idx - 1:10 x[idx] - x[idx] + 1 x [1] 2 1 1 1 1 1 1 1 1 1 This is expected. But if I do: x - c(1, rep(0,9)) x

Re: [R] question regarding vector ops

2004-04-09 Thread Uwe Ligges
Rajarshi Guha wrote: On Fri, 2004-04-09 at 12:03, Rajarshi Guha wrote: Hi, I was playing with some code and came upon this situation. x - c(1, rep(0,9)) x [1] 1 0 0 0 0 0 0 0 0 0 idx - 1:10 x[idx] - x[idx] + 1 x [1] 2 1 1 1 1 1 1 1 1 1 This is expected. But if I do: x - c(1, rep(0,9)) x [1]

Re: [R] missing values for mda package

2004-04-09 Thread zhu wang
Thanks. I was able to use na.omit to remove NAs. But it seems to me this kills one of the advantages of the original algorithm for handling missing values. On Tue, 2004-04-06 at 11:54, Uwe Ligges wrote: zhu wang wrote: Dear helpers, I am trying to use the mda package downloaded from

[R] rectangles without borders

2004-04-09 Thread Tamas Papp
Hi, I want to draw rectangles with rect without borders, but with shading (ie setting density). ?rect says lty: line type for borders; defaults to 'solid'. but if I set it to 0, the border disappears, but the shading also vanishes! But the above says that lty is the line type for

Re: [R] missing values for mda package

2004-04-09 Thread Prof Brian Ripley
Package mda covers many things, including bruto, mars, polyreg and mda itself. Which `the original algorithm' for which option did you have in mind? More concretely, what where you trying to do with the package? Given that the package is the original authors' own code, it seems unlikely that

Re: [R] rectangles without borders

2004-04-09 Thread Prof Brian Ripley
The person who added shading did not update the help page for rect Try border=FALSE (and look at ?polygon for more details as to why). On Fri, 9 Apr 2004, Tamas Papp wrote: Hi, I want to draw rectangles with rect without borders, but with shading (ie setting density). ?rect says

RE: [R] missing values for mda package

2004-04-09 Thread Wang, Zhu
I basically wanted to use MARS to reproduce results using the dataset Marketing in the following book http://www-stat-class.stanford.edu/~tibs/ElemStatLearn/ The authors actually provided S-Plus functions for mars, bruto ,etc. I used all default options of mars in R but there was an error due

RE: [R] missing values for MARS (was mda package)

2004-04-09 Thread Prof Brian Ripley
On Fri, 9 Apr 2004, Wang, Zhu wrote: I basically wanted to use MARS to reproduce results using the dataset Marketing in the following book http://www-stat-class.stanford.edu/~tibs/ElemStatLearn/ The authors actually provided S-Plus functions for mars, bruto ,etc. I used all default

[R] Pfizer - Open position for Statistician - emphasis computational statistics

2004-04-09 Thread Romansky, Peg J
All: The Nonclinical Statistics Group at Pfizer Inc., New London, CT, is currently seeking a Statistician with an interest and experience in computational statistics. Please use the following link for more information and application information: If Interested please apply on-line at

[R] Regression models w/ splines

2004-04-09 Thread John Fieberg
Hi - I am fitting various Cox PH models with spline predictors. After fitting the model, I would like to use termplot() to examine the functional form of the fitted model (e.g., to obtain a plot of the relative risk (or log r.r.) versus the predictors). When there is only 1 predictor in the

Re: [R] Regression models w/ splines

2004-04-09 Thread Peter Dalgaard
John Fieberg [EMAIL PROTECTED] writes: Hi - I am fitting various Cox PH models with spline predictors. After fitting the model, I would like to use termplot() to examine the functional form of the fitted model (e.g., to obtain a plot of the relative risk (or log r.r.) versus the predictors).

[R] von Mises

2004-04-09 Thread MKG
Dear All, I am trying to plot von Mises density on the circle. One can use dvm function from the CircStats package, by giving a set of angles, mu and kappa to plot the circular density on the line. Does any one have a macro that does it on the circle? These plots are displayed in Nick Fisher's

Re: [R] Regression models w/ splines

2004-04-09 Thread Frank E Harrell Jr
On Fri, 09 Apr 2004 15:16:34 -0500 John Fieberg [EMAIL PROTECTED] wrote: Hi - I am fitting various Cox PH models with spline predictors. After fitting the model, I would like to use termplot() to examine the functional form of the fitted model (e.g., to obtain a plot of the relative risk (or

[R] adding text in (pair) panels (splom)

2004-04-09 Thread Patrick Giraudoux
Hi, I would like to add text in each panel got from the following: splom(~Cs5bis, panel=function(x,y,...){ panel.xyplot(x,y,...) panel.abline(lm(y~x),...) }) This could be a correlation coefficient or the statistical significance of the correlation or both, or a star if p 0.05, etc... The

Re: [R] adding text in (pair) panels (splom)

2004-04-09 Thread Paul Murrell
Hi Patrick Giraudoux wrote: Hi, I would like to add text in each panel got from the following: splom(~Cs5bis, panel=function(x,y,...){ panel.xyplot(x,y,...) panel.abline(lm(y~x),...) }) This could be a correlation coefficient or the statistical significance of the correlation or both, or a

Re: [R] adding text in (pair) panels (splom)

2004-04-09 Thread Patrick Giraudoux
OK Paul and Kjetil I understand that grid.text is part of the library grid. I have tried with Paul's arguments. It works perfectly well. Thanks a lot both, Patrick - Original Message - From: Paul Murrell [EMAIL PROTECTED] To: Patrick Giraudoux [EMAIL PROTECTED] Cc: r-help [EMAIL

[R] Density Estimation

2004-04-09 Thread Thami Rachidi
Dear Sir/Madam; Would you please tell me what is the command that allows the estimation of the Kernel Density for some data. Thanks, Thami Rachidi [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

RE: [R] Density Estimation

2004-04-09 Thread Ko-Kang Kevin Wang
-Original Message- From: [EMAIL PROTECTED] Dear Sir/Madam; Would you please tell me what is the command that allows the estimation of the Kernel Density for some data. Thanks, ?density __ [EMAIL PROTECTED] mailing list