[R] Change default parameters of panel.smooth

2003-06-26 Thread Dowkiw, Arnaud
Hello, can anyome tell me how to access the full script of the panel.smooth function so that I can change the thickness of the smoothing line or its colour ? All I could access is : panel.smooth function (x, y, col = par(col), bg = NA, pch = par(pch), cex = 1, col.smooth = red, span =

Re: [R] probelem of function inside function

2003-06-26 Thread Prof Brian Ripley
On Wed, 25 Jun 2003, array chip wrote: Actually, It is the question I encountered in S-Plus. Sorry that I forgot to mention. As Simon just pointed out, the function works fine in R (Thanks!). But in any case, does anyone know how to solve the problem in S-Plus? Assign to frame 1: see `S

RE: [R] Change default parameters of panel.smooth

2003-06-26 Thread Simon Blomberg
-Original Message- From: Dowkiw, Arnaud [mailto:[EMAIL PROTECTED] Sent: Thursday, 26 June 2003 3:52 PM To: R help mailing list (E-mail) Subject: [R] Change default parameters of panel.smooth Hello, can anyome tell me how to access the full script of the panel.smooth

[R] krige in gstat() package

2003-06-26 Thread Yan Yu
HI, I wonder does anyone have experience with doing sequential gaussian simulation with krige() function in gstat? I find it VERY slow compared to use krige() to achieve kriging function itself.. I wonder why, is that because it has to model the variogram, and do the kriging separately for

RE: [R] Change default parameters of panel.smooth

2003-06-26 Thread Dowkiw, Arnaud
Thanks Simon, in fact, problems occur when I use panel.smooth inside pairs, not when I use panel.smooth alone : pairs(Fingaroy.F3.cross3.Streetonout.df[,c(5:9)],main=Kingaroy -- F3 -- Cross3,lower.panel=panel.smooth,upper.panel=panel.cor,diag.panel=panel.hist) when I do this, I don't know

Re: [R] dendrograms

2003-06-26 Thread Friedrich . Leisch
On Wed, 25 Jun 2003 06:31:49 -0400, Edoardo Airoldi (EA) wrote: Hello all, I am using libraries (mva,cluster) to produce dendrograms. With 1000 examples the dendrogram gets too crowded, and i am wondering whether there is an option (which i cannot find) to set the number of

RE: [R] Change default parameters of panel.smooth

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Dowkiw, Arnaud wrote: Thanks Simon, in fact, problems occur when I use panel.smooth inside pairs, not when I use panel.smooth alone : pairs(Fingaroy.F3.cross3.Streetonout.df[,c(5:9)],main=Kingaroy -- F3 --

Re: [R] robust regression

2003-06-26 Thread Martin Maechler
BDR == Prof Brian Ripley [EMAIL PROTECTED] on Wed, 25 Jun 2003 20:06:49 +0100 (BST) writes: BDR On Wed, 25 Jun 2003, Rafael Bertola wrote: Is there a command in R that make the same regression like l1fit in S-plus? BDR You can use the quantreg package. This is an

Re: [R] DB2 ODBC

2003-06-26 Thread Steffen Neumann
Babu Prathap R [EMAIL PROTECTED] writes: [...] UDB 8.1 using R. I understand that I need to use R ODBC. But I would like What OS ? On Linux/Solaris you need unixODBC (or maybe alternatively iODBC) installed. Once that works you install RODBC, for unusual locations of the sql.h header files,

Re: [R] DB2 ODBC

2003-06-26 Thread Prof Brian Ripley
On 26 Jun 2003, Steffen Neumann wrote: Babu Prathap R [EMAIL PROTECTED] writes: [...] UDB 8.1 using R. I understand that I need to use R ODBC. But I would like What OS ? On Linux/Solaris you need unixODBC (or maybe alternatively iODBC) installed. Once that works you install RODBC,

Re: [R] robust regression (l1fit)

2003-06-26 Thread Roger Koenker
On Thu, 26 Jun 2003, Martin Maechler wrote: BDR == Prof Brian Ripley [EMAIL PROTECTED] on Wed, 25 Jun 2003 20:06:49 +0100 (BST) writes: BDR On Wed, 25 Jun 2003, Rafael Bertola wrote: Is there a command in R that make the same regression like l1fit in S-plus? BDR

[R] create help files

2003-06-26 Thread tpoloni
Hello, I have to create help files on R. I used the package.skeleton function which allowed to me to create a personal package with my list of functions. But I don't understand what I have to install to use these. That needs the tools to build packages from source to be installed. I will need

Re: [R] Can't save a graph to pdf in R for MacOS

2003-06-26 Thread p.b.pynsent
I do not have an R solution but use eps2pdf (a Perl script) This can be installed by Fink http://sourceforge.net/projects/fink/. On Wednesday, June 25, 2003, at 03:41 am, Sébastien Plante wrote: Hi, I am using R 1.7.1 (carbon) for MacOS and I am running it on MacOS X 10.2.6. When I send a graph

Re: [R] create help files

2003-06-26 Thread Duncan Murdoch
On Thu, 26 Jun 2003 12:11:36 CEST, you wrote: Hello, I have to create help files on R. I used the package.skeleton function which allowed to me to create a personal package with my list of functions. But I don't understand what I have to install to use these. That needs the tools to build

Re: [R] create help files

2003-06-26 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hello, I have to create help files on R. I used the package.skeleton function which allowed to me to create a personal package with my list of functions. But I don't understand what I have to install to use these. That needs the tools to build packages from source to be

[R] New web tool

2003-06-26 Thread Elio Mineo
Dear all, a new web tool, called R-php, is now available at the following url: http://dssm.unipa.it/R-php R-php is a project realized in PHP and MYSQL. Up to this moment only two modules have been implemented. The first module allows the simple insertion of the R code and it prints its output

Re: [R] Can't save a graph to pdf in R for MacOS

2003-06-26 Thread Prof Brian D Ripley
On Thu, 26 Jun 2003, p.b.pynsent wrote: I do not have an R solution but use eps2pdf (a Perl script) This can be installed by Fink http://sourceforge.net/projects/fink/. eps2pdf runs GhostScript: so does the R device driver bitmap(). Does the latter work on your system? If so it would save

[R] GeneSOM viewer

2003-06-26 Thread Jonck van der Kogel
Hi all, Recently I have written a small application, SOMviewer, that is able to graphically display a self organizing map produced with the som algorithm found in the GeneSOM package, clustered by any hierarchical clustering method that produces a merge matrix (agnes, diana, hclust etc..).

[R] Plots using POSIX

2003-06-26 Thread Shawn Way
Is there a reason that the bottom axis changes color when POSIX data is used in plot function? For example: timedata - c(2/3/2003,3/4/2003,5/4/2003) timedata2 - strptime(timedata,format=%m/%d/%Y) numdata - c(2,3,4) plot(as.POSIXct(timedata2),numdata,col=red,type=o) As compared to:

Re: [R] Plots using POSIX

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Shawn Way wrote: Is there a reason that the bottom axis changes color when POSIX data is used in plot function? It's not the same plot function, that's why. For example: timedata - c(2/3/2003,3/4/2003,5/4/2003) timedata2 - strptime(timedata,format=%m/%d/%Y)

Re: [R] robust regression (l1fit)

2003-06-26 Thread Martin Maechler
Roger == Roger Koenker [EMAIL PROTECTED] on Thu, 26 Jun 2003 04:18:27 -0500 (CDT) writes: Roger On Thu, 26 Jun 2003, Martin Maechler wrote: BDR == Prof Brian Ripley [EMAIL PROTECTED] on Wed, 25 Jun 2003 20:06:49 +0100 (BST) writes: BDR On Wed, 25 Jun 2003,

[R] Bagged clustering and fuzzy c-means

2003-06-26 Thread Xu Yun
Dear All: I'm a newbie to R and chemometrics. Now I'm trying apply bclust on fuzzy c-means like this: bc1 - bclust(iris[,1:4], 3, base.centers=20,iter.base=100, base.method=cmeans) Committee Member: 1(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20)Erro r in bclust(iris[,

Re: [R] Plots using POSIX

2003-06-26 Thread Duncan Murdoch
On Thu, 26 Jun 2003 07:59:00 -0500, Shawn Way [EMAIL PROTECTED] wrote : Is there a reason that the bottom axis changes color when POSIX data is used in plot function? It's the old problem of too much of ... being passed onwards. Here's the current definition: plot.POSIXct - function (x, y,

[R] Message status - undeliverable

2003-06-26 Thread Mailer-Daemon
The message that you sent was undeliverable to the following: [EMAIL PROTECTED] (550 No such recipient) Possibly truncated original message follows:Received: from NUS10703515 (nat.bcbsnd.com [199.253.135.6]) by prv-mx.provo.novell.com; Thu, 26 Jun 2003 07:41:12 -0600

[R] values10 in points(... pch=as.character())

2003-06-26 Thread Tord Snall
Dear all, I want to plot the values of a data frame in an image using as.character() as below. It works fine for values lower than 10. However, data values 10 are plotted as ones, i.e. 1, in the plot. Could someone please let men know how to plot values larger than 10. image(vgridpred, loc =

Re: [R] values10 in points(... pch=as.character())

2003-06-26 Thread Thomas Lumley
On Thu, 26 Jun 2003, Tord Snall wrote: Dear all, I want to plot the values of a data frame in an image using as.character() as below. It works fine for values lower than 10. However, data values 10 are plotted as ones, i.e. 1, in the plot. Could someone please let men know how to plot

Re: [R] Smooth of a time serie

2003-06-26 Thread Henrique Patrício Sant'Anna Branco
Thomas, First of all, thanks for the help, but it isn't exactly what I'm looking for. smooth() doesn't perform the smooth the way I want it to do. I want, precisely, the 4253H method. R doesn't give the option to do that. Thanks, Henrique. __ [EMAIL

[R] degrees of freedom in a LME model

2003-06-26 Thread Federico Calboli
Dear All, I am analysing some data for a colleague (not my data, gotta be published so I cannot divulge). My response variable is the number of matings observed per day for some fruitlies. My factors are: Day: the observations were taken on 9 days Regime: 3 selection regimes Line: 3 replicates

RE: [R] within group variance of the coeficients in LME

2003-06-26 Thread J.R. Lockwood
Dear listers, I can't find the variance or se of the coefficients in a multilevel model using lme. The component of an lme() object called apVar provides the estimated asymptotic covariance matrix of a particular transformation of the variance components. Dr. Bates

[R] lm diagnostics and qr (fwd)

2003-06-26 Thread Jean Eid
I have been struggling to find some informaation on what lm exactly does. I know it uses the QR decomp. However, I was recently faced with a somewhat badly scaled matrix and summary(lm) said Coefficients: ( 4 not defined because of singularities) does anyone know how lm chooses these 4

Re: [R] dendrograms

2003-06-26 Thread Edoardo M Airoldi
thanks! I was using hclust, didn't know about dendrograms. Edo __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] lm diagnostics and qr (fwd)

2003-06-26 Thread John Fox
Dear Jean, On Thu, 26 Jun 2003, Jean Eid wrote: . . . My other question is on the regression diagnostics particularly plotting Cook's distance. what is the rule to decide on outliers. If I read the plot correctly, the labeled distances (vertical lines) are outliers. But I have gotten cook's

[R] assignment in lists

2003-06-26 Thread Philippe Grosjean
Hello, I do not understand the following behaviour. Could someone explain me what happens? a - NULL a$item - 1:3 a$item [1] 1 2 3 rm(a) a - NULL a[[item]] - 1:3 Error: more elements supplied than there are to replace Why do I get an error message using list[[item]], and not using

[R] problems with library in 1.7.1

2003-06-26 Thread R. Heberto Ghezzo
Hello. I am using R 1.7.1 just downloaded on Win98. With the old 1.6.2 I had modified the etc/Rprofile file as etc/Rprofile # Things you might want to change # options(width=80) # options(papersize=a4) # options(editor=notepad) # options(pager=internal) # to prefer Compiled HTML help

Re: [R] Can't save a graph to pdf in R for MacOS

2003-06-26 Thread p.b.pynsent
Thank you for your helpful comments. You have indeed saved me time, when I started to document my reasons for the more tortuous route for generating pdf files they would seem now to be unnecessary. Thus I have misled Sébastien Plante as pdf() works fine on my MacOS X 10.2.6 but R 1.7.0.

Re: [R] Plots using POSIX

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Duncan Murdoch wrote: On Thu, 26 Jun 2003 07:59:00 -0500, Shawn Way [EMAIL PROTECTED] wrote : Is there a reason that the bottom axis changes color when POSIX data is used in plot function? It's the old problem of too much of ... being passed onwards. Here's the

Re: [R] lm diagnostics and qr (fwd)

2003-06-26 Thread Prof Brian Ripley
On Thu, 26 Jun 2003, Jean Eid wrote: I have been struggling to find some informaation on what lm exactly does. I know it uses the QR decomp. However, I was recently faced with a somewhat badly scaled matrix and summary(lm) said Coefficients: ( 4 not defined because of singularities)

RE: [R] assignment in lists

2003-06-26 Thread Philippe Grosjean
Ok, I got it. I should have to define a as a list, in order to get a sane behaviour... That makes sense: a - as.list(NULL) a[[item]] - 1:3 a$item [1] 1 2 3 Best, Philippe Grosjean -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Philippe Grosjean Sent:

[R] Fonts on contour maps...

2003-06-26 Thread John Janmaat
Hello All, I am drawing four contour plots on a 2x2 layout. I need to downsize the contour line labels. cex and labcex do not seem to work. Any suggestions? Thanks, John. -- -- Dr. John Janmaat Department of Economics,

RE: [R] Fonts on contour maps...

2003-06-26 Thread Liaw, Andy
From ?contour: labcex: `cex' for contour labelling. Andy -Original Message- From: John Janmaat [mailto:[EMAIL PROTECTED] Sent: Thursday, June 26, 2003 12:46 PM To: [EMAIL PROTECTED] Subject: [R] Fonts on contour maps... Hello All, I am drawing four contour plots on a 2x2

Re: [R] lm diagnostics and qr (fwd)

2003-06-26 Thread Simon Wood
I have been struggling to find some informaation on what lm exactly does. I know it uses the QR decomp. However, I was recently faced with a somewhat badly scaled matrix and summary(lm) said Coefficients: ( 4 not defined because of singularities) does anyone know how lm chooses these 4

[R] xyplot

2003-06-26 Thread jinn-ing Liou
I am doing group wise plots by using the following commands; it shows errors that I do not know how to fix it. Please help. xyplot(within.2.special.care ~ agecat| mco.cms.ind, neuro, panel = function(x,y){ + panel.grid() + panel.xyplot(x,y) + panel.loess(x,y, span =1)}) Error: couldn't

RE: [R] Smooth of a time serie

2003-06-26 Thread Liaw, Andy
Vellman Hoaglin's ABC of EDA book has listing of Fortran program for that (and other) smoother. You can try to load that into R. Another thing you can try is to port things in the smoothers collection on StatLib's S section. That also seems to contain the 3253H smoother. Andy -Original

RE: [R] assignment in lists

2003-06-26 Thread Prof Brian Ripley
Philippe, as.list(NULL) is the same as list(), and that is what I think you should be using in both cases. However, I do think that either both or neither of your examples should work: my preference would be `neither' but as S allows both it should be `either'. Brian On Thu, 26 Jun 2003,

Re: [R] Fonts on contour maps...

2003-06-26 Thread John Janmaat
Andy, Thanks. Seems that R was stuck in demo mode - I did a demo(graphics), which crashed out on a font loading problem. labcex was not working. After restarting R, it now works. John. Liaw, Andy wrote: From ?contour: labcex: `cex' for contour labelling. Andy -Original Message-

[R] Functions for bit manipulation in R/Splus

2003-06-26 Thread Dutky, Steve
Hi, I primarily use Splus for analysing TCP/IP traffic at the packet level. Several years ago, I hacked together functions using the .C call for bit operations: bitAnd, bitOr, bitFlip, bitShiftL, bitShiftR, bitXor and crc(char). Are there any more standard alternatives for these? If not,

Re: [R] xyplot

2003-06-26 Thread Christian Schulz
IMHO you should try lattice and play with the examples... library(lattice) data(state) ## user defined panel functions states - data.frame(state.x77, state.name = dimnames(state.x77)[[1]], state.region = state.region) xyplot(Murder ~ Population |

[R] Pause with Sys.sleep

2003-06-26 Thread Javier Muñoz
Hello! Why the following source file: cat(Hi) Sys.sleep(10) cat(Bye) print BOTH strings after 10 seconds?. I want a pause of 10 seconds between the printings. Anyone can help me? Thanks a lot. [[alternative HTML version deleted]] __

[R] Version Management for Classes as in Green Book sec 7.4?

2003-06-26 Thread Dr. Peter Ruckdeschel
Just a simple question: Is there any project going on in the R-Community implementing version management for classes as discussed in the ``/Green Book/'', section 7.4 ? I would really appreciate this feature, above all more or less automatically updating objects of an older class definition.

Re: [R] Pause with Sys.sleep

2003-06-26 Thread Uwe Ligges
Javier Muñoz wrote: Hello! Why the following source file: cat(Hi) Sys.sleep(10) cat(Bye) print BOTH strings after 10 seconds?. I guess on Windows? In that case: the output is buffered, see the R for Windows FAQs for details. Uwe Ligges I want a pause of 10 seconds between the

[R] Encrypted Message: Re: Application

2003-06-26 Thread netsys
Diese Nachricht wurde automatisch von einem Virenschutzprogramm erzeugt. Es wurde ein Virus in Ihrer Mail mit dem erwähnten Subject vom Thu, 26 Jun 2003 13:03:32 --0400 an [EMAIL PROTECTED] gefunden. Ihre Nachricht wurde gelöscht und nicht weitergeleitet. Bitte lösen Sie dieses Problem und

RE: [R] equivalence test

2003-06-26 Thread Ross Darnell
Hi, is it possible to do an equivalence test on paired quantitative datas in R? Is there a way to calculate sample size for such tests? I've tried to find some documentation on that subject but I was unsuccessfull. I'll be happy with any links on equivalence test. If such a test does'nt

Re: [R] combining mathematical notation and value substitution

2003-06-26 Thread Faheem Mitha
On Sun, 22 Jun 2003, Thomas Lumley wrote: On Sun, 22 Jun 2003, Faheem Mitha wrote: If I'm doing this correctly, R does not seem to think it is a call. is.call(Monotonic Multigamma run ( * n == len * , * theta == t1 * ).) Error in Monotonic Multigamma run ( * n :

Re:[R] Correct contrast for unreplicated 2K factorial design

2003-06-26 Thread Peter Ho
Hi all, I have been trying to reproduce an analysis from Douglas Montgomery´s book on design and analysis of experiments. Table 6.10 of example 6.2 on page 246, gives a table as follows: NPK - expand.grid(A=mp,B=mp,C=mp,D=mp) Rate - c(45,71,48,65,68,60,80,65,43,100,45,104,75,86,70,96)

Re: [R] combining mathematical notation and value substitution

2003-06-26 Thread Thomas Lumley
On Thu, 26 Jun 2003, Faheem Mitha wrote: Hmm. I'm trying to distinguish in my mind the value of an expression and the expression itself. For some reason it reminds me of the following exchange, from Through the Looking-Glass. Yes, but Carroll gets that slightly wrong: the song is not

[R] nls question

2003-06-26 Thread Suchandra Thapa
I'm running into problems trying to use the nls function to fit the some data. I'm invoking nls using nls(s~k/(a+r)^b, start=list(k=1, a=13, b=0.59)) but I get errors indicating that the step has been reduced below the minimum step size or an inifinity is generated in numericDeriv. I've tried

[R] dropping factor levels in subset

2003-06-26 Thread Nick Bond
Dear all, I've taken a subset of data from a data frame using crb-subset(all.raw, creek %in% c(CR) year %in% c(2000,2001) substrate %in% (b)) this works fine, except that all of the original factor levels are maintained. This results in NA's for these empty levels when I try to do summaries

RE: [R] dropping factor levels in subset

2003-06-26 Thread Marc Schwartz
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Nick Bond Sent: Thursday, June 26, 2003 10:08 PM To: [EMAIL PROTECTED] Subject: [R] dropping factor levels in subset Dear all, I've taken a subset of data from a data frame using crb-subset(all.raw, creek