My reaction when learning of a proposed patent on a new graph was: "oh
well, that's something I can forget about". Without a patent, code would
have been available in R in a very short period of time, the statistical
community would have been able to play around with it, see how it worked
on v
Dear Francisco,
1. Have a look at ?Devices which (under Windows at least)
lists a range of devices you can open before and
between plots.
2. Use par(ask=TRUE) before you start your plots
3. Christophe Declercq provided an excellent example on this
mailing list, under the topic "curves wit
Dear R helpers,
is there an equivalent of SAS's PROC LATTICE in R ?
I have a partially balanced square lattice design to study and I want to compare
lattice based analysis with the randomized complete block approach.
Thanks a lot for your help,
Arnaud
DISCLAIM
I have been trying to read a random sample of lines from a file into a
data frame using readLines(). The help indicates that readLines() will
start from the current line if the connection is open, but presented with
a closed connection it will open it, start from the beginning, and close
it when fi
On Wed, 27-Aug-2003 at 01:40PM -0400, Alvaro Muñoz wrote:
|> Drs. Harrell and O'Keefe,
|>
|>
|>
|> Thank you for your suggestions.
|>
|> Although it is at odds with your beliefs, University staff working
|> on licensing and technology transfer believe that a patent may be a
|> vehicle to achi
Hi everyone. R is new to me and I'm very impressed with its
capabilities but still cannot figure out how to do some basic things.
There seems to be no lack of documentation but finding what I need has
proven difficult. Perhaps you can help.
Here's what I'm after:
1. How do I create a new p
Alvaro Muñoz wrote:
Drs. Harrell and O'Keefe,
Although it is at odds with your beliefs, University staff working on
licensing and technology transfer believe that a patent may be a vehicle to
achieve a wide use. The audience of the proposed methods would be the end
users who are not sophisticated
Hello, everyone! I wonder if any R package can do the multivariate
Smirnov test. Specifically, let x_1,..,x_n and y_1,...,y_m be
multivariate vectors. I would like to test if the two samples are from
the same underlying multivariate distribution. Thanks in advance.
Jason
=
Jason G. Liao, Ph.
Thanks to all (and especially Prof. Tierney) for the response. The box we
are considering will spend probably over 90% of CPU time in R, so it's
comforting to know that R compiles and pass all the test (at least once) on
such platform.
(I switched my attention from Itanium to Opteron when I read
On 26 Aug 2003, Peter Dalgaard BSA wrote:
> Dirk Eddelbuettel <[EMAIL PROTECTED]> writes:
>
> > On Tue, Aug 26, 2003 at 03:17:19PM -0400, Liaw, Andy wrote:
> > > Has anyone tried using R on the the AMD Opteron in either 64- or 32-bit
> > > mode? If so, any good/bad experiences, comments, etc? W
In a recent talk ('Visions of the Past, Present & Future of Statistical
Graphics'),
I talked about, among other things, the lessons Minard's March on Moscow
graphic had
for modern statistical graphics, and illustrated aspects of power and
simplicity
in several programming languages where this gr
Hi,
There have been a number of reports of RMySQL crashing R when
attempting to connect to a MySQL server using dbConnect().
The problem appears to be in some binary versions of the MySQL
client library. Known instances include
(1) Red Hat MySQL binary RPM client library 3.23.32, but
upd
On Wed, 27 Aug 2003 13:40:59 -0400
Alvaro Muñoz <[EMAIL PROTECTED]> wrote:
> Drs. Harrell and O'Keefe,
>
>
>
> Thank you for your suggestions.
>
>
>
> Regarding your comments about the content of the paper, I respectfully
> disagree that "categorizing continuous variables is a fundamental vi
Drs. Harrell and O'Keefe,
Thank you for your suggestions.
Regarding your comments about the content of the paper, I respectfully
disagree that "categorizing continuous variables is a fundamental violation
of statistical graphics," nor are you to assume that all categorizations are
arbitrary.
Thank you all for the quick responses.
However, I'm not sure I unterstand the scaling matrix (denote S
henceforth) correcty. An observation x will be transformed by Sx into a
new vector space with the properties given by the description. What is
now the direction perpendicular to the seperating pla
RINNER Heinrich <[EMAIL PROTECTED]> writes:
> Dear R users!
>
> I am using R Version 1.7.1, Windows XP, package "foreign" (Version: 0.6-1),
> SPSS 11.5.1.
>
> There is one thing I noticed with "read.spss", and I'd like to ask if this
> is considered to be a feature, or possibly a bug:
> When rea
On Wed, 27 Aug 2003, RINNER Heinrich wrote:
> Dear R users!
>
> I am using R Version 1.7.1, Windows XP, package "foreign" (Version: 0.6-1),
> SPSS 11.5.1.
>
> There is one thing I noticed with "read.spss", and I'd like to ask if this
> is considered to be a feature, or possibly a bug:
> When rea
Hi, Ronaldo:
Have you talked with anyone in the Math Department in the Univ.Fed.de
Viçosa? They offer courses in Statistics there, and I would expect that
someone there could help you get copies of the articles of interest. I
wonder if such contacts might help you with other statistics-rel
On Wed, 27 Aug 2003, Martin Hoyle wrote:
> Dear R Users,
>
> I suppose this is a school boy question, but here it is anyway. I'm trying to
> re-create the residuals for a poisson GLM with simulated data;
>
> x<-rpois(1000,5)
> model<-glm(x~1,poisson)
> my.resids<-(log(x)- summary(model)$coeffic
As ?residuals.glm reveals, it's got an argument type:
type: the type of residuals which should be returned. The
alternatives are: `"deviance"' (default), `"pearson"',
`"working"', `"response"', and `"partial"'.
You calculated "response" residuals, R gives "deviance" residu
Please excuse my reference to Batelle: I confused who was asking and
who answering the question. Thanks, Peter, for the clarification and
for the alternative suggestions.
Spencer Graves
Peter Dalgaard BSA wrote:
Spencer Graves <[EMAIL PROTECTED]> writes:
The Battelle Institute surely should
Dear R users!
I am using R Version 1.7.1, Windows XP, package "foreign" (Version: 0.6-1),
SPSS 11.5.1.
There is one thing I noticed with "read.spss", and I'd like to ask if this
is considered to be a feature, or possibly a bug:
When reading character columns, character strings seem to get filled
Spencer Graves <[EMAIL PROTECTED]> writes:
> The Battelle Institute surely should have access to a library with
> such popular and prestigious journals as Biometrika and The American
> Statistians. If you don't have time for that, you surely should have
> money to purchase a copy from, e.g., "www
Sven Garbade <[EMAIL PROTECTED]> writes:
> Hi all,
>
> suppose I've got a vector y with some data (from a repeated measure
> design) observed given the conditions in f1 and f2. I've got a model
> with two unknown fix constants a and b which tries to predict y with
> respect to the values in f1 an
That's the linear model lm(y ~ I(1/f1) + f2), so yes, yes and
fuller answers can be found in most of the books and guides mentioned in
R's FAQ.
Note that how `good' the fit is will have to be relative, unless you
really can assume a uniform error with range 1, when you could do a
maximum-likelih
Dear R Users,
I suppose this is a school boy question, but here it is anyway. I'm trying to
re-create the residuals for a poisson GLM with simulated data;
x<-rpois(1000,5)
model<-glm(x~1,poisson)
my.resids<-(log(x)- summary(model)$coefficients[1])
plot(my.resids,residuals(model))
This shows tha
On Wed, 27 Aug 2003 11:04:21 -0300
"Ronaldo Reis Jr." <[EMAIL PROTECTED]> wrote:
> Can anybody send these articles for me?
>
> > NagelKerke, N. J. D. (1991) "A note on a general definition of the
> > coefficient of determination", Biometrika 78: 691-2.
The fitting functions lrm, psm, cph
Hi all,
suppose I've got a vector y with some data (from a repeated measure
design) observed given the conditions in f1 and f2. I've got a model
with two unknown fix constants a and b which tries to predict y with
respect to the values in f1 and f2. Here is an exsample
# "data"
y <- c(runif(10, -
On Wed, 27 Aug 2003, Petr Pikal wrote:
> Hallo
>
> On 27 Aug 2003 at 1:49, Eugene Salinas wrote:
>
> > Hi,
> >
> > I'm a recent R convert so I haven't quite figured out
> > the details yet...
> >
> > How do I select one variable by another one? Ie if I
> > want to draw the histogram of variable X
Marc Schwartz <[EMAIL PROTECTED]> writes:
> On Wed, 2003-08-27 at 09:21, [EMAIL PROTECTED] wrote:
> > Hi guys,
> >
> >Yesterday I posted my first couple of questions (see bottom of this
> > message) to this forum and I would like to thank you guys for all the
> > useful feedback I got. I just
The Battelle Institute surely should have access to a library with such
popular and prestigious journals as Biometrika and The American
Statistians. If you don't have time for that, you surely should have
money to purchase a copy from, e.g., "www.lindahall.org/docserv".
hope this helps. spenc
On Wed, 2003-08-27 at 09:21, [EMAIL PROTECTED] wrote:
> Hi guys,
>
>Yesterday I posted my first couple of questions (see bottom of this
> message) to this forum and I would like to thank you guys for all the
> useful feedback I got. I just would like to make some comments:
>
> 1. Exporting R
I think you've badly misinterpreted the purpose
of the R listserv with this request:
https://www.stat.math.ethz.ch/mailman/listinfo/r-help says
"The `main' R mailing list, for announcements about the
development of R and the availability of new code, questions
and answers about problems and s
On Tue, 26 Aug 2003, Ross Boylan wrote:
>
> 2) Modifications of existing packages.
> I modified the survival package (I'm not sure if that's properly called
> a "base" package, but it's close). I know in this particular case, if
> I'm serious, I probably should contact the package maintainer. Bu
On Wed, 27 Aug 2003, Prof Brian Ripley wrote:
> On Tue, 26 Aug 2003, Ross Boylan wrote:
>
> > So do people have any thoughts about whether introducing this is a step
> > forward or back?
>
> A useful analogue: we now distribute Fortran code not the original Ratfor.
>
As a footnote to Brian's comm
Can anybody send these articles for me?
> NagelKerke, N. J. D. (1991) "A note on a general definition of the
> coefficient of determination", Biometrika 78: 691-2.
>
> Cox, D. R. and Wermuth, N. (1992) "A comment on the coefficient of
> determination for binary responses", The Ame
Hi guys,
Yesterday I posted my first couple of questions (see bottom of this
message) to this forum and I would like to thank you guys for all the
useful feedback I got. I just would like to make some comments:
1. Exporting R graphs as vector graphics:
The best answer came from Thomas Lumley
Thank you so much for that fix (to my understanding).
I would be willing to add such an example to the help
page for future releases - though I'm sure others would
do it better - there are currently no examples where
INDICES is a name.
In fact in my real application it is more or less essentia
"Simon Fear" <[EMAIL PROTECTED]> writes:
> I tried to define a function like:
>
> fnx <- function(x, by.vars=Month)
> print(by(x, by.vars, summary))
>
> But this doesn't work (does not find x$Month; unlike other functions,
> such as
> subset(), the INDICES argument to "by" does not look for va
On Wed, 27 Aug 2003, Simon Fear wrote:
> I tried to define a function like:
>
> fnx <- function(x, by.vars=Month)
> print(by(x, by.vars, summary))
>
> But this doesn't work (does not find x$Month; unlike other functions,
> such as
> subset(), the INDICES argument to "by" does not look for vari
I tried to define a function like:
fnx <- function(x, by.vars=Month)
print(by(x, by.vars, summary))
But this doesn't work (does not find x$Month; unlike other functions,
such as
subset(), the INDICES argument to "by" does not look for variables in
dataset
x. Is fully documented, but I forget ev
On Wed, 27 Aug 2003 13:03:39 +1200 (NZST), you wrote:
>For real efficiency here, what's wanted is a variant of readLines
>where n is an index vector (a vector of non-negative integers,
>a vector of non-positive integers, or a vector of logicals) saying
>which lines should be kept.
I think that's
> "Jason" == Jason Fisher <[EMAIL PROTECTED]>
> on Tue, 26 Aug 2003 14:08:26 -0700 writes:
Jason> For anyone wishing to see R-Tcl/Tk-MySQL in action (Windows XP)...
Jason> http://moffett.isis.ucla.edu/gwplot/
Jason> Examples were by far the most useful learning tool
Ja
hist(X[Y>=A & Y<=B])
`An Introduction to R' explains such things, as do (in more detail) the
introductory texts (see the R FAQ).
On Wed, 27 Aug 2003, Eugene Salinas wrote:
> I'm a recent R convert so I haven't quite figured out
> the details yet...
>
> How do I select one variable by another
Hallo
On 27 Aug 2003 at 1:49, Eugene Salinas wrote:
> Hi,
>
> I'm a recent R convert so I haven't quite figured out
> the details yet...
>
> How do I select one variable by another one? Ie if I
> want to draw the histogram of variable X only for
> those individuals that also have a value Y in a
On Wed, 27 Aug 2003, Eugene Salinas wrote:
> I'm a recent R convert so I haven't quite figured out
> the details yet...
Usually it is good to read the manuals when you use a unfamiliar
software...
> How do I select one variable by another one? Ie if I
> want to draw the histogram of variable X
Eugene Salinas <[EMAIL PROTECTED]> wrote:
> How do I select one variable by another one? Ie if I
> want to draw the histogram of variable X only for
> those individuals that also have a value Y in a
> certain range?
>
> In STATA I would give something like:
>
> histogram X if ((Y>=A & Y<=B))
his
Hallo
On 26 Aug 2003 at 13:09, Jablonsky, Nikita wrote:
> Hi All,
>
> I am trying to match two character arrays (email lists) using either
> pmatch(), match() or charmatch() functions. However the function is
> "missing" some matches due to differences in the cases of some letters
try toupper or
Eugene,
R allows indexing with logical vectors, so your example would look
like
hist(X[(Y>=A) & (Y<=B)])
See the manual "An Introduction to R" for details.
HTH
Thomas
> -Original Message-
> From: Eugene Salinas [mailto:[EMAIL PROTECTED]
> Sent: 27 August 2003 09:49
> To: [EMAIL PROTEC
Hi,
I'm a recent R convert so I haven't quite figured out
the details yet...
How do I select one variable by another one? Ie if I
want to draw the histogram of variable X only for
those individuals that also have a value Y in a
certain range?
In STATA I would give something like:
histogram X if
On Tue, 26 Aug 2003, Ross Boylan wrote:
> I have two questions about packaging up code.
>
> 1) Weave/tangle advisable?
> In the course of extending some C code already in S, I had to work out
> the underlying math. It seems to me useful to keep this information
> with the code, using Knuth's tan
On Tue, 26 Aug 2003, Lily wrote:
> Seems crossval from library(bootstrap) can only be
> used for leave-one-out and k-fold cross validation?
> Here is a dumb question, suppose n=80, how to do
> exactly leave-50-out cross validation? K-fold cross
> validation is not eligible for this case since
> n/
I'm bored, but just to point out the obvious fact: to skip n lines in a
text file you have to read *all* the characters in between to find the
line separators.
I have known for 30 years that reading text files of numbers is slow and
inefficient. So do it only once and dump the results to a binar
If we are going to use unix tools to create a new dataset before calling
into R, why not simply use
cat my_big_bad_file | tail +1001 | head -100
to read lines 1000-1100 (assuming one header row).
Or if you have the shortlisted rownames in one file, you can use join
after sort. A working example
Seems crossval from library(bootstrap) can only be
used for leave-one-out and k-fold cross validation?
Here is a dumb question, suppose n=80, how to do
exactly leave-50-out cross validation? K-fold cross
validation is not eligible for this case since
n/ngroup is not an integer. Thanks!
___
> I do think that in an article you should also point out to people that
> there is a lot of numerical code available out there, written by people
> who know a lot more than we do about what they are doing. It's often
> easier than writing your own code and the results are better. One
> advantage o
Duncan Murdoch <[EMAIL PROTECTED]> wrote:
For example, if you want to read lines 1000 through 1100, you'd do it
like this:
lines <- readLines("foo.txt", 1100)[1000:1100]
I created a dataset thus:
# file foo.awk:
BEGIN {
s = "01"
for (i = 2; i <= 41; i++) s
57 matches
Mail list logo
