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
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
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
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, :
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)
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
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
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
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
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]
Hi,
I want to calculate the Lyapunov exponents of several time series (16 years
of weekly observations each with a strong seasonal component).
Ive got a reference that compares three methods: Feedforward Neural
Networks, Thin-Plate Splines and Response Surface (Ellner & Turchin 1995,
American
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
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
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
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
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
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.
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
> >
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
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
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
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
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
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
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
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
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
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
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.
>
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
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
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
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
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
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
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
_
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
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/ /'
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:
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
"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
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,
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,
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
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"
^^
___
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
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-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
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
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
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). '
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
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
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
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
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
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
> 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
58 matches
Mail list logo