Re: [R] Trouble with optim

2003-01-31 Thread ripley
I thinks it has found a *local* minimum. It is not a global optimizer. You appear to have made no attempt to scale the problem as the help page suggests you need to. Please look into this, as from what little information you give I think the initial step has overshot. On Sat, 1 Feb 2003, Damo

[R] [Summary] Problems for 13 year old

2003-01-31 Thread ggrothendieck
This e-mail - summarizes some of the responses/suggestions to using R for 13 year olds from my query on this list - discusses some additional investigation I did - the approach I am currently using (which will likely incorporate your suggestions for R) SUGGESTIONS - find relationship between

[R] Re-assigning vector elements based on their initial values.

2003-01-31 Thread [EMAIL PROTECTED]
Is there an eloquent solution to re-assign vector element values? I have a vector which contains chemical data, some of them are "flagged" as non-detected values by their negative values. I can find the statistics on the positive values in vector "v" simply by typing: >v<- c(5,5,-3,-3,7,8,10) > v

[R] Read.table problem

2003-01-31 Thread ptremb17
Hi ! I am new to R, and using the MAC version onto Mac OS 9.1. My question concerns the problem I encounter when I try to read some data I have using the read.table function. I always get an error of type : Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, :

[R] Trouble with optim

2003-01-31 Thread Damon Wischik
I am having trouble with optim. It claims to have converged to a minimum, yet it has in the course of the optimization visited many points which are closer to optimal. I would be grateful for any explanation of this behaviour. I'm trying to estimate the parameters in the model X ~ Binomial(1,p)

[R] Do not understand: log10(-0.04)

2003-01-31 Thread Thomas Hoffmann
I ran into some mysterious crashes of R and when stepping through the code it appeared to me that a call to GScale(0.0,1.0,1,dd) (as e.g. in do_plot_new(), plot.c:458 ) possibly with min=-0.04 runs into (graphics.c:1911) Rf_gpptr(dd)->logusr[0] = Rf_dpptr(dd)->logusr[0] = log10(min); whic

[R] find max of implicit function OR inversion of 2D mappings.

2003-01-31 Thread Vadim Ogranovich
Dear R-Users, I am looking for a help to deal with the following computational problem: I have a mapping f(x, y) -> (u, v) of [0,1]*[0,1] -> R^2. The mapping is given by tabulating f(x,y) on a uniform 2D grid and is assumed to be "interpolatable" in between the grid points (the number of points o

Re: [R] Testing ``<=" in R

2003-01-31 Thread Peter Dalgaard BSA
Maria N Kocherginsky <[EMAIL PROTECTED]> writes: > Hello, > > I've encountered the following: > > > n_500 > > tau_.95 > > (n*(1-tau)) > [1] 25 > > (n*(1-tau))<=25 > [1] FALSE > > (n*(1-tau))==25 > [1] FALSE > > I'm using UNIX R Version 1.4.0, and also tested in out in Windows 1.6.0. Is > this

Re: [R] Testing ``<=" in R

2003-01-31 Thread Thomas Lumley
On Fri, 31 Jan 2003, Maria N Kocherginsky wrote: > Hello, > > I've encountered the following: > > > n_500 > > tau_.95 > > (n*(1-tau)) > [1] 25 > > (n*(1-tau))<=25 > [1] FALSE > > (n*(1-tau))==25 > [1] FALSE > > I'm using UNIX R Version 1.4.0, and also tested in out in Windows 1.6.0. Is > this a bu

[R] Testing ``<=" in R

2003-01-31 Thread Maria N Kocherginsky
Hello, I've encountered the following: > n_500 > tau_.95 > (n*(1-tau)) [1] 25 > (n*(1-tau))<=25 [1] FALSE > (n*(1-tau))==25 [1] FALSE I'm using UNIX R Version 1.4.0, and also tested in out in Windows 1.6.0. Is this a bug? Masha __ [EMAIL PROTECTED]

[R] Lyapunov - time series

2003-01-31 Thread Gregorio Gutierrez
Hi, I wan’t to calculate the Lyapunov exponents of several time series (16 years of weekly observations each with a strong seasonal component). I’ve got a reference that compares three methods: Feedforward Neural Networks, Thin-Plate Splines and Response Surface (Ellner & Turchin 1995, American

Re: [R] Problems with boot package (empinf returns NA)

2003-01-31 Thread Roger Peng
How many bootstrap resamples did you use in your original call to `boot'? I sometimes get this error when the number of resamples is small. Try increasing the value for the `R' parameter in `boot' and rerun the bca.ci function. For example, > set.seed(100) > f <- function(x, i) median(x[i]) > bc

[R] Version 0.8.0 of the gregmisc package is now available

2003-01-31 Thread Warnes, Gregory R
Version 0.8.0 of the gregmisc package is now (or will shortly be) available on CRAN. New in this release: - Enhanced and bug-fixed 'CrossTable' function (contributed by Marc Schwartz) - Augmented the 'barplot2' function with an 'add' argument to allow for the addition of a barplot to an

RE: [R] png()/jpeg()

2003-01-31 Thread Wilkinson, Mark
Sorry: _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major1 minor6.1 year 2002 month11 day 01 language R Mark Wilkinson Informatics Analyst St. Jude Children's Research Hospital Department of Pharmaceutical Sciences

Re: [R] png()/jpeg()

2003-01-31 Thread ripley
It really would help to know the platform here: Windows and Unix have separate implementations. On Windows with R 1.6.2 there is no overlap that I can see. Also the version of R, please. On Fri, 31 Jan 2003, Wilkinson, Mark wrote: > When I execute the following code, it works just like I want

Re: [R] floating point question

2003-01-31 Thread ripley
On Fri, 31 Jan 2003, Thomas Lumley wrote: > On 31 Jan 2003, Peter Dalgaard BSA wrote: > > > Deepayan Sarkar <[EMAIL PROTECTED]> writes: > > > > > Might have something to do with .Machine$double.eps on the respective > > > machines. > > > > > > From help(.Machine), > > > > > > double.eps: the smal

[R] png()/jpeg()

2003-01-31 Thread Wilkinson, Mark
When I execute the following code, it works just like I want it to: three pages of nine (or fewer) plots. However, when I execute the code with the first and last lines uncommented, I get three pages (files), but the 2nd & 3rd pages have overlapping plots. It's like a new page wasn't created.

Re: [R] floating point question

2003-01-31 Thread Thomas Lumley
On 31 Jan 2003, Peter Dalgaard BSA wrote: > Deepayan Sarkar <[EMAIL PROTECTED]> writes: > > > Might have something to do with .Machine$double.eps on the respective > > machines. > > > > From help(.Machine), > > > > double.eps: the smallest positive floating-point number `x' such that > >

Re: [R] floating point question

2003-01-31 Thread Peter Dalgaard BSA
Deepayan Sarkar <[EMAIL PROTECTED]> writes: > Might have something to do with .Machine$double.eps on the respective > machines. > > From help(.Machine), > > double.eps: the smallest positive floating-point number `x' such that > `1 + x != 1'. It equals `base^ulp.digits' if either `ba

RE: [R] Multiple xyplot on a page

2003-01-31 Thread Yang, Richard
Many thanks to Renaud Lancelot, Andy Liaw and Deepayan Sarkar for their prompt response to my query. They all point me to the print.trellis help page which indeed solves my problem. Thank all, Richard My original question: Dear all; I have a set of 5 xyplots each with 8 panels and

Re: [R] floating point question

2003-01-31 Thread Prof Brian D Ripley
It's a difference in the `libc'. Asking for more precision than the arithmetic has is asking for fairly random results. The differences are as likely to be in the *printing* as in the computations. On Fri, 31 Jan 2003, Bob Gray wrote: > Does anyone know precisely what is different about the ari

Re: [R] floating point question

2003-01-31 Thread Deepayan Sarkar
Might have something to do with .Machine$double.eps on the respective machines. >From help(.Machine), double.eps: the smallest positive floating-point number `x' such that `1 + x != 1'. It equals `base^ulp.digits' if either `base' is 2 or `rounding' is 0; otherwise, it is

[R] floating point question

2003-01-31 Thread Bob Gray
Does anyone know precisely what is different about the arithmetic and/or storage of double precision floating point to produce the following differences between the Sun and Windows versions (Splus 6 on the same Windows 2000 machine gives the same results as Solaris)? R 1.6.1, Sun Solaris, gcc + an

Re: [R] Multiple xyplot on a page

2003-01-31 Thread Deepayan Sarkar
See ?print.trellis On Friday 31 January 2003 11:50 am, Yang, Richard wrote: > Dear all; > > I have a set of 5 xyplots each with 8 panels and wish to place 2 on > a letter-sized page. I tried using pdf() to specify paper size and plot > > size: > >pdf(file = "fig1b.pdf", paper = "letter", on

RE: [R] Multiple xyplot on a page

2003-01-31 Thread Liaw, Andy
The help page for print.trellis should give you some hints. Andy > -Original Message- > From: Yang, Richard [mailto:[EMAIL PROTECTED]] > Sent: Friday, January 31, 2003 12:50 PM > To: R-Help (E-mail) > Subject: [R] Multiple xyplot on a page > > > Dear all; > > I have a set of 5 xy

Re: [R] Scientific Notation on tick marks

2003-01-31 Thread Uwe Ligges
Suzy Smith wrote: Hi, I was working on a project for my class and I am trying to make sure the y-axis numbers NEVER convert automatically to scientific notation. Is there anything I can set in the plot function to make sure of this? I know formatC( ) can be used in the axis function, but it seem

Re: [R] Fonts in expression(paste)

2003-01-31 Thread Uwe Ligges
Suzanne E. Blatt wrote: Hello. I'm trying to change the fonts of my text. If I have 'expression(paste)' included (to get scientific symbols), the font= doesn't take effect. Is there a way around this? Thanks, Suzanne ?plotmath helps: plain(x) draw x in normal font bold(x) draw x in bold f

[R] Problems with boot package (empinf returns NA)

2003-01-31 Thread Ernesto Jardim
Hi I'm using boot package for some analysis on linear regression coeficients. My problem is that I can not compute bca intervals, I get an error message > bca.ci(blm8901,index=1) Error in if (!all(rk > 1 & rk < R)) warning("Extreme Order Statistics used as Endpoints") : missing value whe

Re: [R] (no subject)

2003-01-31 Thread kjetil brinchmann halvorsen
On 30 Jan 2003 at 18:17, pavel koulikov wrote: You should use a subject line. To put to lines on the same plot: plot(x~y, type="l") lines(z ~ v) Kjetil Halvorsen > Hi > I have 4 vectors x,y,z,v > I need 2 graphs (x~y,z~v) on the same plot. I mean two lines on the same plot - blue >and red. >

[R] Fonts in expression(paste)

2003-01-31 Thread Suzanne E. Blatt
Hello. I'm trying to change the fonts of my text. If I have 'expression(paste)' included (to get scientific symbols), the font= doesn't take effect. Is there a way around this? Thanks, Suzanne __ The NEW Netscape 7.0 browser is

[R] Multiple xyplot on a page

2003-01-31 Thread Yang, Richard
Dear all; I have a set of 5 xyplots each with 8 panels and wish to place 2 on a letter-sized page. I tried using pdf() to specify paper size and plot size: >pdf(file = "fig1b.pdf", paper = "letter", onefile = FALSE, width = 5, height = 5, family="Helvetica") >par(mfrow=c(1,2)) but it di

Re: [R] svm regression in R

2003-01-31 Thread David Meyer
Christoph Helma wrote: > > Hallo, > > I have a question concerning SVM regression in R. I intend to use SVMs for feature >selection (and knowledge discovery). For this purpose I will need to extract the >weights that are associated with my features. I understand from a previous thread on >SVM

[R] Scientific Notation on tick marks

2003-01-31 Thread Suzy Smith
Hi, I was working on a project for my class and I am trying to make sure the y-axis numbers NEVER convert automatically to scientific notation. Is there anything I can set in the plot function to make sure of this? I know formatC( ) can be used in the axis function, but it seems to give me an

[R] Slavic characters

2003-01-31 Thread Miha STAUT
Hi all, How would I get to print other slavic letters than only the ones present in the standard ASCII code as mentioned in Hershey vector fonts? Miha Staut __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help

RE: [R] Decreasing my personal entropy ...

2003-01-31 Thread Warnes, Gregory R
Lattice plots and trellice plots are actually displayed by thier 'print' methods. So, to get ta plot to be displayed in a non-interactive context, such as a script, explicitly call print on the results of the plot command. -Greg -Original Message- From: Gunter, Bert To: R-Help (E-mail) S

RE: [R] Varying texts in expression(paste())

2003-01-31 Thread Morten Sickel
Gunter Bart wrote: > You need substitute(). Look at the examples in plotmath(). > x<-3;y<-'foo' > plot(c(0,1),c(0,1),type='n') > text(.5,.5,substitute(paste('Releases of ',x^99,' ',y,' > TBq/year'),list(x=x,y=y))) > Hope this helps ... Thanks a lot, it helped indeed! Morten _

Re: [R] TukeyHSD

2003-01-31 Thread ripley
You implicitly (probably) used na.omit on your dataset, and it added an na.action attribute to the data frame. Try not doing that (e.g. bu using na.action = na.fail and search the archives for the same question about model.tables). Are there any missing values in your dataset? If so, you should w

[R] r-bugs web site temporarily down

2003-01-31 Thread Peter Dalgaard BSA
The machine that serves r-bugs.biostat.ku.dk has been taken off the net. A new machine is planned to replace it, but we need to reconfigure the DNS, which is not going to happen until Monday. The mail interface should still work. -- O__ Peter Dalgaard Blegdamsvej 3 c/ /'

[R] TukeyHSD

2003-01-31 Thread C-G Pettersson
Hello everybody! I´m working with a dataset from eleven field trails on barley fertilization. I use R 1.6.2 (Windows) It is quite easy to fit aov() objects to the dataset. The call: > (l1t4y.aov <- aov(Yield ~ Trial + Treatment, data=led1t4b)) Results in an object with this anova table:

[R] svm regression in R

2003-01-31 Thread Christoph Helma
Hallo, I have a question concerning SVM regression in R. I intend to use SVMs for feature selection (and knowledge discovery). For this purpose I will need to extract the weights that are associated with my features. I understand from a previous thread on SVM classification, that predictive mod

Re: [R] Decreasing my personal entropy ...

2003-01-31 Thread Douglas Bates
"Gunter, Bert" <[EMAIL PROTECTED]> writes: > R-Listers: > > A very minor -- and maybe silly -- question just for personal enlightenment. > > In S (either R or S-Plus, AFAIK) when one types or pastes a trellis graphics > command into the commands/console window, the graph is automatically > produ

Re: [R] Library for GUI?

2003-01-31 Thread Christian Schulz
library(tcltk) library(RGtk) christian - Original Message - From: "Francisco do Nascimento Junior" <[EMAIL PROTECTED]> To: "R-help" <[EMAIL PROTECTED]> Sent: Friday, January 31, 2003 3:02 PM Subject: [R] Library for GUI? Are there some library for to make graphic interface with buttons,

Re: [R] Library for GUI?

2003-01-31 Thread John Fox
Dear Francisco, See the tcltk package. John At 11:02 AM 1/31/2003 -0300, you wrote: Are there some library for to make graphic interface with buttons, fields, menus??? - John Fox Department of Sociology McMaster University Hamilton, Ontario,

Re: [R] Decreasing my personal entropy ...

2003-01-31 Thread John Fox
Dear Bert, At 08:45 AM 1/31/2003 -0500, Gunter, Bert wrote: R-Listers: A very minor -- and maybe silly -- question just for personal enlightenment. In S (either R or S-Plus, AFAIK) when one types or pastes a trellis graphics command into the commands/console window, the graph is automatically p

[R] Library for GUI?

2003-01-31 Thread Francisco do Nascimento Junior
Are there some library for to make graphic interface with buttons, fields, menus??? Tks, Francisco. ^^ Francisco Júnior, Computer Science - UFPE-Brazil "One life has more value that the world whole" ^^ ___

[R] RODBC & Sys.sleep

2003-01-31 Thread Christian Schulz
Hi, for a real-time tracking i fetch the data with Sys.sleep in a loop, but with the newest version of RODBC the "TDX MySQL Driver connect" pop every time intervall up. Is it possible to disable this, until now i have no success find anything ? MySQL 3.23.44 , RODBC 1-0-1,R.1.6.2,W2K many than

RE: [R] Varying texts in expression(paste())

2003-01-31 Thread Winfried Theis
Hi, I had a similar problem last year and Thomas Lumley helped me a lot on this. I needed to replace some variables by their values in a quite complicated text expression. Here is what Thomas suggested: In your example you wanted b and c to be elements of a vector. They actually have to be elem

[R] Decreasing my personal entropy ...

2003-01-31 Thread Gunter, Bert
R-Listers: A very minor -- and maybe silly -- question just for personal enlightenment. In S (either R or S-Plus, AFAIK) when one types or pastes a trellis graphics command into the commands/console window, the graph is automatically produced: e.g., trellis.device(...) xyplot(y~x) If one puts t

[R] Varying texts in expression(paste())

2003-01-31 Thread Morten Sickel
Hi, I am using R a lot to make plots relating to radioactivity, I am often using expression() to label the plots with nuclide names written with superscripts, e.g. expression(paste("Releases of ", { }^{99},Tc," (TBq/year)"))->ywtext But, is there any simple way to change the number and na

Re: [R] typo in read.table and count.fields?

2003-01-31 Thread Peter Dalgaard BSA
Claudio Agostinelli <[EMAIL PROTECTED]> writes: > It seems that in read.table and count.fields there is a typo in the quote > parameter. > > read.table(file, header = FALSE, sep = "", quote = "\"'", dec = ".", >row.names, col.names, as.is = FALSE, na.strings = "NA", >colCl

Re: [R] nearest neighbour interpolation

2003-01-31 Thread Barry Rowlingson
Roger Bivand wrote: Point-in-polygon functions are in the splancs package, but there you need the polygons. There's also nearest-neighbour functions in Splancs - the nn2d function returns the distances and indices of the nearest points in one set (the grid) to another set (the data points). '

[R] typo in read.table and count.fields?

2003-01-31 Thread Claudio Agostinelli
It seems that in read.table and count.fields there is a typo in the quote parameter. read.table(file, header = FALSE, sep = "", quote = "\"'", dec = ".", row.names, col.names, as.is = FALSE, na.strings = "NA", colClasses = NA, nrows = -1, skip = 0, check.names = T

Re: [R] Validation of clustering

2003-01-31 Thread Christian Hennig
Hi, I am not really sure what you want to know. If you simply want to know how to find the resulting clustering, consider help(clara.object). For validation and assessment of the clustering, you might use a silhouette plot, see help(plot.partition), help(partition.object). For more descriptions

Re: [R] nearest neighbour interpolation

2003-01-31 Thread ripley
On 31 Jan 2003, Jari Oksanen wrote: > On Fri, 2003-01-31 at 12:47, Roger Bivand wrote: > > On Thu, 30 Jan 2003, Matthew Oliver wrote: > > > > > Dear Help List, > > > > > > My name is Matt Oliver. I have been using R for about a year and find it very > > > helpful. However, I have a need for a f

Re: [R] nearest neighbour interpolation

2003-01-31 Thread Jari Oksanen
On Fri, 2003-01-31 at 12:47, Roger Bivand wrote: > On Thu, 30 Jan 2003, Matthew Oliver wrote: > > > Dear Help List, > > > > My name is Matt Oliver. I have been using R for about a year and find it very > > helpful. However, I have a need for a function that I cannot find. I am not very > > good

Re: [R] mgcv, gam

2003-01-31 Thread ripley
On Fri, 31 Jan 2003, Simon Wood wrote: > > I have some problems with gam in mgcv. Firts a detail: it would > > be nice igf gam would accept an na.action argument, but that not the > > main point. > - I find it hard to think of a sensible action except dropping the > associated data, but if you've

Re: [R] nearest neighbour interpolation

2003-01-31 Thread Roger Bivand
On Thu, 30 Jan 2003, Matthew Oliver wrote: > Dear Help List, > > My name is Matt Oliver. I have been using R for about a year and find it very > helpful. However, I have a need for a function that I cannot find. I am not very > good at programming so I thought I would ask the group. > Maybe th

Re: [R] mgcv, gam

2003-01-31 Thread Simon Wood
> I have some problems with gam in mgcv. Firts a detail: it would > be nice igf gam would accept an na.action argument, but that not the > main point. - I find it hard to think of a sensible action except dropping the associated data, but if you've a concrete suggestion I'm happy to add it to the