Sorry, I should have mentioned that
at the end of each script I have been using
dev.off().
-david paul
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
R1.8.0, Win2k:
When I run the code
lset( list( background = list(col = white)))
xyplot
(
y ~ TIME , data = foo.frame,
scales = list(alternating = FALSE),
ylim = c(.75,y.max),
panel = function(x, y, panel.number, ... )
{
panel.superpose(x =
: Sunday, October 19, 2003 12:48 PM
To: Paul, David A; '[EMAIL PROTECTED]'; 'Deepayan Sarkar'
Subject: Re: [R] problem with win.metafile( ): traceback()
On Sunday 19 October 2003 11:49, Paul, David A wrote:
For the first error message:
win.metafile(file = //.../plot1.wmf,
+ width = 8.5, height
R1.8.0, Win2k:
When I paste the code
win.metafile(file = //.../plot1.wmf,
width = 8.5, height = 6.25)
lset( list( background = list(col = white)))
xyplot( y ~ x | ID, data = Group1,
scales = list(alternating = FALSE),
ylim = c(.75,y.max),
panel = function(x, y,
-
From: Paul, David A
Sent: Sunday, October 19, 2003 12:28 PM
To: '[EMAIL PROTECTED]'
Subject: [R] problem with win.metafile( )
R1.8.0, Win2k:
When I paste the code
win.metafile(file = //.../plot1.wmf,
width = 8.5, height = 6.25)
lset( list( background = list(col = white
Suppose foo.frame is your dataframe. Use
foo.frame[j,]
to extract all of the elements in the jth row
of foo.frame. Also, you should get two excellent
books that will help your learning process:
The Basics of S and Splus, by Krause and Olsen,
Third Edition, Springer-Verlag
Modern Applied
paul
-Original Message-
From: Jesus Frias [mailto:[EMAIL PROTECTED]
Sent: Wednesday, September 17, 2003 10:35 AM
To: Paul, David A; R help
Subject: RE: [R] gnls( ) question
Hi Paul,
The message has to do with the gauss-newton algorithm and it
announces
library(tree)
?tree
should work. If you don't have the tree library,
you can download it off of CRAN at
http://cran.r-project.org/bin/windows/contrib
if you're using Windows, or go to
http://cran.r-project.org/src/contrib/PACKAGES.html
for the source code directly as gzipped .tar files.
Last week (Wed 9/10/2003, regression questions) I posted
a question regarding the use of gnls( ) and its dissimilarity
to the syntax that nls( ) will accept. No one replied, so
I partly answered my own question by constructing indicator
variables for use in gnls( ). The code I used to
I have been puzzling over how to fit some fixed effects models
to a set of data. My response function is
response - function(a, b, c, alpha1, alpha2, indicator, t, t2)
{
z = a +
b * (t) * exp(-alpha1 * t) +
indicator *c * (t2) * exp(-alpha2 * t2)
}
where
I am one of only 5 or 6 people in my organization making the
effort to include R/Splus as an analysis tool in everyday work -
the rest of my colleagues use SAS exclusively.
Today, one of them made the assertion that he believes the
numerical algorithms in SAS are superior to those in Splus
and R
My thanks to Drs. Armstrong, Bates, Harrell, Liaw, Lumley,
Prager, Schwartz, and Mr. Wang for their replies. I have
pasted my original message and their replies below.
After viewing http://www.itl.nist.gov/div898/strd/ as suggested
by Dr. Schwartz, it occurred to me that it might be educational
Use the nlme( ) package:
library(nlme)
The lme( ) and nlme( ) functions in this package are
designed to estimate mixed effects models. An
outstanding reference is Mixed Effects Models in S
and Splus, by Drs. Pinheiro and Bates.
Best,
david paul
-Original Message-
From: Jörg Peter
Win2k, R1.7.1:
I am currently working with some growth curve data from a
biotoxicology experiment. Each of 12 subjects had their blood
drawn at 0, 2, 4, 6, 8, and 10 weeks. For the purposes of the
project, it would be helpful if I were able to do the following:
a. Produce 12 panels,
Thank you VERY much! That does address my questions.
Respectfully,
david paul
-Original Message-
From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
Sent: Friday, August 29, 2003 1:54 PM
To: Paul, David A; '[EMAIL PROTECTED]'
Subject: Re: [R] Lattice plot questions
Let's consider
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
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
Undoubtedly a simple question:
I've looked at order() and sort() in the help pages for
R1.7.1. It doesn't appear that these functions are immediately
suited to doing the same thing as
PROC SORT DATA = BLAH;
BY X Y Z;
RUN;
in SAS. I have also checked Frank Harrell's Hmisc library.
You can use the well-known Taylor series approximation to the
variance of an arbitrary function:
Var( f(X) ) ~= Sum( s[i]^2*D2[i] ) + 2*Sum( Sum( s[i,j]*D[i]*D[j] ) )
where D2[i] is the second partial derivative of f(x) with respect
to the ith parameter and D[j] is the first partial derivative
R1.7.1/Win2k:
Apologies if this posts twice - the first message
was not in plain text.
I have looked in help.start() and tried typing
crossprod and %*% into the RGui to get an idea
for what R is using as internal algorithms for
its matrix computations/manipulations... to no
avail.
Could
The most commonly used dose-response functions for nonlinear calibration
curves are the four- and five-parameter logistic functions. The four-
parameter logistic is specified as
F(z) = delta + (alpha - delta)/(1 + (z/gamma)^beta)
so I'm not sure where you are getting your dose-response
R1.7, Win2k:
I have some Splus code that has allowed me, in the past,
to place multiple plots on the same graph:
plot(y1 ~ x1, data = foo1.frame, type = l, xlab=,ylab=,
xlim=c(...,...), ylim=c(...,...), axes=F)
par(new=T, xaxs = d, yaxs = d)
plot(y2 ~ x2, data = foo2.frame, col = 2)
Thanks to Dr. Thomas Hotz, Prof. Brian Ripley, Dr. Dennis
Murphy, and Dr. David Scott for their replies.
Trying an idea:
x - c(1,2,3,4)
y-c(2,4.2,5.9,9)
temp-data.frame(cbind(x,y))
attach(temp)
temp.lm - lm(y~x)
windows()
plot(temp.lm)
Hit Return to see next plot:
Hit Return to see next
Thank you both for pointing out that this is a lattice
plot (ie, R's version of Trellis graphics) and therefore
needs something other than par().
I was able to use IE6.0 to search for trellis and
find the relevant commands (after using help.start(),
of course). This brings up another question:
In Splus the code
test.lm - lm(y ~ x, data = test.data)
plot(test.lm)
generates a graphics window that contains
multiple graph sheets that one may choose
from via the page tabs at the bottom of
the window.
Is there a way to do this sort of thing in
R? As another example, I have some repeated
I have found Mixed Effects Models in S and Splus by
Drs. Pinheiro and Bates to be enormously helpful.
I highly recommend the book - it contains excellent
examples.
Best,
david paul
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Friday, June 13, 2003 2:54
In a contour plot having numeric values displayed
next to the contours, is it possible to modify these
numeric values so that they are bolded, enlarged,
and/or placed differently?
Much thanks in advance,
David Paul, Ph.D.
Battelle Memorial Institute
614.424.3176
I don't have the reference, but a biologist friend of mine
once showed me a refereed journal article that purported
to demonstrate numerical errors made by MSExcel. This
would have been Excel97 or Excel2000... In any case, the
journal's scope was biological in nature and the article
was of
to use and it avoids those ...you could try `help.search( blah
)'...
messages. Suggestions in this area would also be greatly appreciated.
Respectfully,
david paul
-Original Message-
From: Paul, David A
Sent: Monday, March 31, 2003 5:07 PM
To: '[EMAIL PROTECTED]'
Subject: [R] nonpos
I apologize, in advance, for cross-posting this to
the R listserv. I have submitted this query twice
to the S listserv (yesterday and this morning)and
neither post has made it, not sure why.
When I run the code
gls.1 -
gls(y ~ x, data = foo.frame,
weights = varPower(form = ~
Title: Exporting Splus Data to R
Using
data.dump(c(foo.frame.1,foo.frame.2, ...), file=DumpData, oldStyle=T)
in Splus generates the requisite file in my working Splus directory. However, when I use
library(foreign)
data.restore(file = C:/.../DumpData,verbose = TRUE, env =
31 matches
Mail list logo