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 achieve a
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
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
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
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
OK, I am try to produce an array of ranks. I have a set of data (s)
that looks like this
rt subj
312 dave
467 dave
411 dave
383 kim
398 kim
...
Now I want to make a column that is an array of the ranks of rt by
subject. The closest of gotten is using the following.
r
You need to save the connection object returned by file() and then use that
object in other functions.
You need to change the appropriate lines to the following (at least):
con - file(c:/data/perry/data.csv,open=r)
cline - readLines(con,n=1)
close(con)
(I don't know if more changes are needed
You are using the connection the wrong way. You need to do something like:
fcon - file(c:/data/perry/data.csv, open=r)
for (iline in 1:slines) {
isel - isel + 1
cline - readLines(fcon, n=1)
...
}
close(fcon)
BTW, here's how I'd do it (not tested!):
strvec - rep(,slines)
selected -
Need some help with R?
An R workshop, led by John Maindonald, is planned for the 1st October
at the University of Southern Queensland, Toowoomba, Australia. The
workshop is being organised by the Queensland branch of the
Statistical Society of Australia as part of their conference.
The
Need some help with R?
An R workshop, led by John Maindonald, is planned for the 1st October
at the University of Southern Queensland, Toowoomba, Australia. The
workshop is being organised by the Queensland branch of the
Statistical Society of Australia as part of their conference.
The workshop
Duncan Murdoch [EMAIL PROTECTED] wrote:
One complication with reading a block at a time is what to do when you
read too far.
It's called buffering.
Not all connections can use seek() to reposition to the
beginning, so you'd need to read them one character at a
I have been reading this discussion (or debate, depends on your point of
view) with great interest in the last few days.
On Thu, 28 Aug 2003, David Scott wrote:
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,
I need to install some geostatistics packages that are not in the binaries,
for example gstat. I downloaded rm171 and installed, and I even added
precomplied packages successfully. Unfortunately, I do not know how to add
packages available as gstat. If you can spare a minute, could you give me a
Thanks Hedderik, Andrew, David and of course, you Mark. You have all
been very helpful. I got good answers two 1. and 2. and arsenal of
ideas for 3. Mark's solution for 3. is exactly what I was thinking
about. I'm indeed very grateful.
Best,
-Francisco
On Wednesday, August 27, 2003, at
I am calculating LOD scores in genetic epidemiology where the number of
non recombinants is 11 and recombinants is 1. I created a function to
calculate the z scores for the following thetas.
theta - c(0,0.01,0.05,0.10,0.20,0.45,0.5)
The z scores for the respective thetas in this cases are:
Could you clarify that you want to do this on the Carbon MacOS port, and
not the Darwin MacOS X port?
Since gstat uses Unix-like configure tools, I think this would need a
_lot_ of expertise, something not worth doing for the last release of the
Carbon port.
If you are running MacOS X, why
Am Mittwoch, 27. August 2003 22:22 schrieb Jason Liao:
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
Talking about Excel, you can produce excellent graphs in Excel. Yes you
have to work at it, but you can get there. The problem is that they are
not the default. My gut feeling is that R will make more of an impact in
the presentation of graphics than any implementation in Excel.
So even if a
Hi
On 27 Aug 2003 at 20:32, John Christie wrote:
OK, I am try to produce an array of ranks. I have a set of data (s)
that looks like this
rtsubj
312 dave
467 dave
411 dave
383 kim
398 kim
...
Now I want to make a column that is an array of the ranks of rt by
David Scott wrote:
but if any of us are to do
likewise I guess we will either pay up or secretly write code and play
around with diamond graphs while hidden in the basement.
Okay, lets stand up and be counted: who has been writing diamond graph
code? Mine's 60 lines.
Like others I welcome the
RINNER Heinrich wrote:
In R:
library(foreign)
test - read.spss(test.sav, to.data.frame=T)
test
XCHAR
1 a
2 ab
3 abcde
levels(test$XCHAR)
[1] a ababcde
Shouldn't it rather be a ab abcde (no blanks)?
I think, that should be no problem since the blanks in XCHAR may be
easily
Marc Schwartz wrote:
In general, if you want to leave the existing device open and have a new
device open for a new plot, you simply call the device name that you
want to open (ie. under Linux you would use X11() ) to open a new
plotting device on the display. See ?Devices for more details.
X11()
Dear All
I am trying to replicate a numerical application (not computed on R) from an
article. Using, ks.test() I computed the exact D value shown in the article
but the p-values I obtain are quite different from the one shown in the
article.
The tests are performed on a sample of 37 values
Thomas Petzoldt [EMAIL PROTECTED] writes:
Marc Schwartz wrote:
In general, if you want to leave the existing device open and have a new
device open for a new plot, you simply call the device name that you
want to open (ie. under Linux you would use X11() ) to open a new
plotting device
You appear to be applying the KS test after estimating parameters. The
distribution theory is for an iid sample from a known continuous
distribution (and does not otherwise depend on the distribution). Since
your H_0 is not pre-specified, that distribution theory is not correct.
(Some
Thanks to Brian Ripley, Douglas Bates and Thomas Petzoldt for their
comments.
I agree that it is not really a problem, as you can easily use sub() after
read.spss() to get rid of the blanks (I had already done that).
On the other hand, it might be important to _know_ about the fact that
On Thu, 28 Aug 2003, Prof Brian Ripley wrote:
You appear to be applying the KS test after estimating parameters. The
distribution theory is for an iid sample from a known continuous
distribution (and does not otherwise depend on the distribution). Since
your H_0 is not pre-specified, that
With the Shifted Exponential test, H_0 is data is a sample coming from a
Shifted Exponential distribution with shift=30 and lambda = 0.001566907
You appear to be applying the KS test after estimating parameters.
I do in order to define H_0 as explained above
The distribution theory is for an
On Thu, 28 Aug 2003, franck allaire wrote:
With the Shifted Exponential test, H_0 is data is a sample coming from a
Shifted Exponential distribution with shift=30 and lambda = 0.001566907
You got that after looking at the data. H_0 has to be specified before
looking at the data.
You appear
Hi,
I am a beginner user of R. I have a trivial question I am almost
ashamed I cannot figure it out does not matter how many times I
am reading the help.
I have a table in .txt format, tab delimited. I can read it with
read.delim() with no problems.
Afterwards I would like to use boxplot
Hi all,
I would like to call a C program from R function. I tried to use the
.C() function
without success. I need a very simple example (such as the hello
program) to understand
how it works. Could you give a such example?
(I use the 1.7.1 Version of R with linux)
Does it work in the case
Monica Palaseanu-Lovejoy [EMAIL PROTECTED] writes:
I am a beginner user of R. I have a trivial question I am almost
ashamed I cannot figure it out does not matter how many times I
am reading the help.
I have a table in .txt format, tab delimited. I can read it with
read.delim() with
Martin Olivier wrote:
Hi all,
I would like to call a C program from R function. I tried to use the
.C() function
without success. I need a very simple example (such as the hello
program) to understand
Do you want to call a C program (executable) or a shared library? If you
want to run a
Hi,
Thank you so much for all your rapid answers. I am impressed.
What i didn't know was that i have to assign my data to an object
to work further on. It was not clear from the help (at least for me)
that 'data()' itself is calling data already in R packages. All of you
make that clear.
Monica Palaseanu-Lovejoy [EMAIL PROTECTED] writes:
Hi,
I am a beginner user of R. I have a trivial question I am almost
ashamed I cannot figure it out does not matter how many times I
am reading the help.
I have a table in .txt format, tab delimited. I can read it with
If you read the data into a data frame, you should be able to simply pass the name of
the data frame in a call to boxplot.
my.data - read.delim(mytext.txt)
boxplot(my.data)
If you only want a boxplot of column 5
boxplot(my.data[,5])
See ?boxplot for other options to make the boxplot look the
On Thu, 28 Aug 2003, Martin Olivier wrote:
Hi all,
I would like to call a C program from R function. I tried to use the
.C() function
without success. I need a very simple example (such as the hello
program) to understand
how it works. Could you give a such example?
There is one in
Dear All,
A couple of questions about the nls package.
1. I'm trying to run a nonlinear least squares
regression but the routine gives me the following
error message:
step factor 0.000488281 reduced below `minFactor' of
0.000976563
even though I previously wrote the following command:
Have you considered normal probability plots (qqnorm) to identify
outliers? These will identify much more, of course, including the need
for transformations, mixtures of distributions, etc.
hope this helps. spencer graves
Monica Palaseanu-Lovejoy wrote:
Hi,
Thank you so much for all your
giovanni caggiano [EMAIL PROTECTED] writes:
Dear All,
A couple of questions about the nls package.
1. I'm trying to run a nonlinear least squares
regression but the routine gives me the following
error message:
step factor 0.000488281 reduced below `minFactor' of
0.000976563
even
Please use an informative subject line. The r-help archives at
www.r-project.org - search - R site search indexes that, and I find
answers to today's problems in the r-help discussions of yesterday's
questions.
Only yesterday, I was got essentially that error message. I solved
it by
Hi,
Figure 13 of Emmanuel Paradis's R for Beginners was produced by termplot
working on an aov object. The lower right-hand plot is labelled Cook's
distance plot, and I'd really like to produce a similar type of figure,
but in a totally different context. (I'm not even sure what this kind of
giovanni caggiano [EMAIL PROTECTED] writes:
A couple of questions about the nls package.
1. I'm trying to run a nonlinear least squares
regression but the routine gives me the following
error message:
step factor 0.000488281 reduced below `minFactor' of
0.000976563
even though I
On Thu, 28 Aug 2003, Frank Gibbons wrote:
Hi,
Figure 13 of Emmanuel Paradis's R for Beginners was produced by termplot
working on an aov object. The lower right-hand plot is labelled Cook's
distance plot, and I'd really like to produce a similar type of figure,
but in a totally
On Thu, 28 Aug 2003, Frank Gibbons wrote:
Hi,
Figure 13 of Emmanuel Paradis's R for Beginners was produced by termplot
working on an aov object.
No, it was produced by plot() working on a aov object, as its caption
indicates. The termplot() is Figure 14.
The lower
I agree with what you said about using trace = TRUE when you are
having trouble getting nls to converge. This allows you to see what
is happening to the parameters during the iterations and that it is
often quite instructive; as is plotting your data and thinking about
whether you should expect
A Guide for the Unwilling S User would help orient you
to how R works. It is meant to do that as quickly and painlessly
as possible.
Patrick Burns
Burns Statistics
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling S User)
Monica
Thanks to all who responded, and so promptly too: it works exactly as you
describe.
Figure 13 of Emmanuel Paradis's R for Beginners was produced by termplot
working on an aov object.
No, it was produced by plot() working on a aov object, as its caption
indicates. The termplot() is Figure 14.
The MacOSX (carbon port) command turned out to be, surprisingly enough,
macintosh()
Thanks!
On Thursday, August 28, 2003, at 07:19 AM, Peter Dalgaard BSA wrote:
Thomas Petzoldt [EMAIL PROTECTED] writes:
Marc Schwartz wrote:
In general, if you want to leave the existing device open and have a
I am working with a RedHat 9 system with gcc 3.2.2, makeinfo 4.3, and
perl 5.8.
I have just downloaded the latest r-patched with rsync and
successfully compiled the code and passed the tests conducted with
'make check'.
After that, 'make dvi' gives the following result:
make dvi
make[1]:
Have you looked at Modern Applied Statistics with S by Venables and
Ripley?
hope this helps. spencer graves
Frank Gibbons wrote:
Thanks to all who responded, and so promptly too: it works exactly as
you describe.
Figure 13 of Emmanuel Paradis's R for Beginners was produced by
termplot
On Thu, 28 Aug 2003, Frank Gibbons wrote:
Thanks to all who responded, and so promptly too: it works exactly as you
describe.
Figure 13 of Emmanuel Paradis's R for Beginners was produced by termplot
working on an aov object.
No, it was produced by plot() working on a aov object,
I'm (relatively) new to R myself, and recently found
the documentation by Dr. William Cleveland et al at
http://cm.bell-labs.com/cm/ms/departments/sia/project/trellis/
http://cm.bell-labs.com/cm/ms/departments/sia/project/trellis/software.writi
ng.html
It has been helpful, and covers S/Splus
Dear R users:
Is there a function or easier way to randomize the data between rows
(considering a row as a whole unit)? For example,
original data:
A1 A2 A3 A4
B1 1 2 3 4
B2 5 6 7 8
B3 9 10 11 12
B4 13 14 15 16
randomized data:
A1 A2 A3 A4
B4 13 14 15 16
B2 5 6
A - array(1:16, dim=c(4,4))
dimnames(A) - list(LETTERS[1:4], letters[1:4])
A[sample(1:4), ]
A - array(1:16, dim=c(4,4))
dimnames(A) - list(LETTERS[1:4], letters[1:4])
A[sample(1:4), ]
a b c d
C 3 7 11 15
B 2 6 10 14
D 4 8 12 16
A 1 5 9 13
hope this helps. spencer graves
[EMAIL PROTECTED]
mtext does not appear to be rendering a 'bold' expression.
Is there another parameter to set? Thx.
example (does not create bold (font=2) on plot)
mtext( font=2, expression(paste(y, = , x + z), side=3 )
__
[EMAIL PROTECTED] mailing list
In my experience, transformations of the type Doug just described has
often made sums of squares (or log(likelihood)) contours more parabolic,
thereby increasing the accuracy of the simple normal approximations to
the distributions of parameter estimates. It is wise to check these
things, as
57 matches
Mail list logo