Last year there was a request seeking functions to plot Venn diagrams in
R. Seeing no reply and for other reasons needing this, I wrote a quick
routine. The general problem of creating a Venn diagram with overlap
areas proportional to the number of counts in each overlap is over
determined. An
Many thanks to Patrick Connolly and Deepayan Sarkar for their quick
responses to my query. In essence,
grid graphics cannot be used in conjunction with base R graphic functions
such as mtext(). Deepayan suggests using bold() in the ylab in place of font
= 2; the bold() does produce the desired f
On Monday 03 February 2003 03:47 pm, Yang, Richard wrote:
> Dear all;
>
> I wish to create a graphic object combing an xyplot() and an mtext().
You cannot. xyplot() uses grid for all its graphics, and grid graphics cannot
be used in conjunction with base R graphics functions.
> My
> code looks
Dear all;
I wish to create a graphic object combing an xyplot() and an mtext(). My
code looks like following,
gmv <- {
trellis.device("windows", bg="white", width = 7, height = 7)
xyplot(Mvol ~ Age | Nl * Th , data = Hft1,
prepanel = function(x, y) prepanel.loess(x, y
Patricia Maforte dos Santos wrote:
>
> Hi,
>
> I'm having difficulties in the manipulation of the function table.
> Necessary to have access its elements and I am not obtaining. For example,
> if I have a data.frame with name "df" and make:
> t < - data.frame(table(df))
Dou you really want to ge
Dear Michael,
You could use plot(x, as.numeric(y) - 1).
John
At 01:10 PM 2/3/2003 -0700, Michael Lynn Fugate wrote:
Hi,
When I use the function plot() to plot a categorical variable with two
levels versus a continuous variable, the two levels of the categorical
variable are plotted at the valu
Michael Lynn Fugate wrote:
>
> Hi,
>
> When I use the function plot() to plot a categorical variable with two
> levels versus a continuous variable, the two levels of the categorical
> variable are plotted at the values of 1.0 and 2.0. I would like them to
> be plotted at the values of 0.0 and
> "Michael" == Michael Lynn Fugate <[EMAIL PROTECTED]> writes:
> Hi, When I use the function plot() to plot a categorical
> variable with two levels versus a continuous variable, the
> two levels of the categorical variable are plotted at the
> values of 1.0 and 2.0. I would l
Hi,
When I use the function plot() to plot a categorical variable with two
levels versus a continuous variable, the two levels of the categorical
variable are plotted at the values of 1.0 and 2.0. I would like them to
be plotted at the values of 0.0 and 1.0. How can I do this?
Example:
x <- 1:
Thanks to Brian D. Ripley and Uwe Ligges. It works now.
Zhu Wang
Statistical Science Department
Southern Methodist University
3225 Daniel Avenue, Heroy 123
Dallas, TX 75275
Tel: (214)768-2453
Fax: (214)768-4035
>
>
> > -Original Message-
> > From: Uwe Ligges [mailto:[EMAIL PR
How can I extend the signature of a generic function, because the following
example is unclear for me.
setClass("A",representation(data="numeric"));
setGeneric("func",function(obj,...){
res <- standardGeneric("func");
res@cal <- match.arg();
return(res);
});
func1.A <- function(obj){
pr
This has nothing to do with R as such; I'm simply trying to exploit
the vast resource of knowledge and expertise that resides in the R
community.
For my sins, I am teaching a course on Game Theory this term. Game
Theory I know from ***nothing***. I am trying to learn the subject
as I go along,
There are jackknife functions about, but this is not jackknifing. Unless
popstat is itself vectorized (meaning I think that it can take list of
datasets, or perhaps a matrix) I doubt if anything is better than (a).
Remember Jackson's rules of programming (which are quoted in `S
Programming'). D
Hi all,
I'm implementing a population genetic statistic that requires repeated
re-estimation of population parameters after a single observation has
been left out. It seems to me that one could:
a) use looping in R,
b) use a vectorized approach in R,
c) loop in a dynamically loaded c-function,
d
On Mon, 3 Feb 2003, John Bjornar Bremnes wrote:
> The solution turned out to be simple:
But wrong. That's not the right class. Thomas Lumley's suggestion is
better, and mine is correct but less elegant.
> > x <- 0:100
> > class(x) <- "POSIXct"
> > x
--
Brian D. Ripley, [EMA
On Mon, 3 Feb 2003, John Bjornar Bremnes wrote:
> A colleague of mine would like to know how seconds since 1970
> (represented as integers) can be converted to date-time objects?
structure(as.double(x), class=c("POSIXt", "POSIXlt"))
will do it in the current version.
--
Brian D. Ripley,
Thomas Lumley <[EMAIL PROTECTED]> writes:
> On Mon, 3 Feb 2003, John Bjornar Bremnes wrote:
>
> > A colleague of mine would like to know how seconds since 1970
> > (represented as integers) can be converted to date-time objects?
>
> ISOdate(1970,1,1)+seconds
I think you might want
ISOdate(197
Hello, The other day I wrote about my SD() function did not work with NA
in the input and it did not recognise the na.rm=T parameter.
BDR ask me if I had another SD function loaded. I did not wrote any sd
function but with ls() there it was sd <- function(x)
sqrt(var(as.vector(x)), I found out t
Steve Howie <[EMAIL PROTECTED]> writes:
> Hi
>
> I've installed the "R" system under Solaris 8 on a SunFire V880. I
> configured it using the "--with-x" option. Installation seemed to go
> fine, but when I start R with the X-Windows GUI option,
>
> R --gui=X11
>
>
On Mon, 3 Feb 2003, John Bjornar Bremnes wrote:
> The solution turned out to be simple:
>
> > x <- 0:100
> > class(x) <- "POSIXct"
> > x
>
Yes, but that relies on the internal representation of POSIXct objects,
which is probably unwise.
-thomas
__
On Mon, 3 Feb 2003, John Bjornar Bremnes wrote:
> A colleague of mine would like to know how seconds since 1970
> (represented as integers) can be converted to date-time objects?
ISOdate(1970,1,1)+seconds
-thomas
__
[EMAIL PROTECTED] mailing l
The solution turned out to be simple:
> x <- 0:100
> class(x) <- "POSIXct"
> x
John Bjornar
John Bjornar Bremnes wrote:
>
> A colleague of mine would like to know how seconds since 1970
> (represented as integers) can be converted to date-time objects?
>
> thanks
>
> --
> John Bjornar Bremne
Hi
I've installed the "R" system under Solaris 8 on a SunFire V880. I
configured it using the "--with-x" option. Installation seemed to go
fine, but when I start R with the X-Windows GUI option,
R --gui=X11
I just get put into the command-line mode.
Any ideas what I'm missing?
Thanks
S
Peter Dalgaard BSA <[EMAIL PROTECTED]> writes:
> 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.
We got swamp
It makes sense and it works!
Thanks and best regards,
José A. S. Alegria
-Original Message-
From: Marc Schwartz [mailto:[EMAIL PROTECTED]]
Sent: segunda-feira, 3 de Fevereiro de 2003 4:18 PM
To: José Santos Alegria; [EMAIL PROTECTED]
Subject: RE: [R] Overlaying a moving average cur
Hi,
I'm having difficulties in the manipulation of the function table.
Necessary to have access its elements and I am not obtaining. For example,
if I have a data.frame with name "df" and make:
t < - data.frame(table(df))
and later trying to make a comparison of the type t[1,1 ] < 0.1, generates
a
>-Original Message-
>From: [EMAIL PROTECTED]
>[mailto:[EMAIL PROTECTED]] On Behalf Of José
>Santos Alegria
>Sent: Monday, February 03, 2003 5:30 AM
>To: [EMAIL PROTECTED]
>Subject: [R] Overlaying a moving average curve on top of a barplot
>
>
>I'm using standard barplot (Windows version 1
On Mon, Feb 03, 2003 at 10:00:46AM -0500, Liaw, Andy wrote:
> > From: Martin Maechler [mailto:[EMAIL PROTECTED]]
> >
> > > "Giovanni" == Giovanni Petris <[EMAIL PROTECTED]>
> > > on Mon, 3 Feb 2003 08:36:06 -0600 (CST) writes:
> >
> > Giovanni> Has anybody else experienced somethi
A colleague of mine would like to know how seconds since 1970
(represented as integers) can be converted to date-time objects?
thanks
--
John Bjornar Bremnes
Norwegian Meteorological Institute (met.no)
Research and Development Department
P.O.Box 43 Blindern, N-0313 Oslo, Norway
Phone: (+47) 229
On Mon, 3 Feb 2003 [EMAIL PROTECTED] wrote:
> All,
>
> I was wondering if someone had come across the problem of producing partial
> regression plots for glm objects in R?
>
> When using Splus in the past I have passed glm objects to the plot.gam
> function.
> To my knowledge this functionality is
Sorry for this off-topic subject.
I am fighting for running SJava under Windows.
SJava_0.64 (compiled by Simon Urbanek, thanks), R 1.6.2, Java JDK 1.4.0_02,
Windows XP pro:
> library(SJava)
> .JavaInit()
Error in .JavaInit() : Couldn't start Java Virtual Machine: Cannot find the
Omegahat interfac
> From: Martin Maechler [mailto:[EMAIL PROTECTED]]
>
> > "Giovanni" == Giovanni Petris <[EMAIL PROTECTED]>
> > on Mon, 3 Feb 2003 08:36:06 -0600 (CST) writes:
>
> Giovanni> Has anybody else experienced something like the
> example below?
> not recently.
> Could it be that your v
> "Giovanni" == Giovanni Petris <[EMAIL PROTECTED]>
> on Mon, 3 Feb 2003 08:36:06 -0600 (CST) writes:
Giovanni> Has anybody else experienced something like the example below?
not recently.
Could it be that your version of the `grid' package (which is
loaded by `lattice') or `latti
Has anybody else experienced something like the example below?
Any clues about where I could start looking?
Thank you in advance,
Giovanni
> version
_
platform sparc-sun-solaris2.7
arch sparc
os solaris2.7
system sparc, solaris2
Dear Jonathan,
The all.effects function in the car package will identify the high-order
terms in a linear or generalized-linear model and absorb terms marginal to
them (e.g., for a two-way interaction, the main effects marginal to the
interaction and the constant). The plot method for the objec
Hello Rabaie,
how about the contributed package: "RQuantLib"?
HTH,
Bernhard
-Original Message-
From: Rabaie Remon [mailto:[EMAIL PROTECTED]]
Sent: 02 February 2003 19:24
To: [EMAIL PROTECTED]
Subject: [R] Option pricing
Dear sir:
I want just ask if there is in R any package to use l
Dear Toby,
I think that you mean "partial-residual" (i.e., component+residual) plots
rather than "partial-regression" (i.e., added-variable) plots. In either
event, see the cr.plots and av.plots functions in the car package; both
have methods for GLMs.
I hope that this does what you need,
Joh
This discussed almost daily.
You *print* a trellis/lattice plot, and no printing is done in your
example.
On Mon, 3 Feb 2003 [EMAIL PROTECTED] wrote:
> Something I don't understand ... (!)
It's in all good books on S/R
> With the lattice library loaded, I have a loop
>
> for( Z in ...
Something I don't understand ... (!)
With the lattice library loaded, I have a loop
for( Z in ... ) {
...
xyplot(y~x | t, xlab=..., ylab=... )
}
and no plot appears on the R graphics device.
However, when I run the commands within {...}
separately for each instance of Z, I get the
plot d
My problem is not with "plot()" but with "barplot()"! I guess it may have to do with
the fact that the barplot's bars have a non-negligible width and the moving average
line not! Is it?
José A. S. Alegria
-Original Message-
From: Morten Sickel [mailto:[EMAIL PROTECTED]]
Sent: segu
Hi
On 1 Feb 2003 at 7:17, John Fox wrote:
> Dear Rex,
>
> You can use if.else, which is vectorized:
>
> > v.new <- ifelse(v < 0, -v/2, v)
> > v.new
> [1] 5.0 5.0 1.5 1.5 7.0 8.0 10.0
or simply
abs(v/((v<0)+1))
it should be a little bit quicker
>
> (By the
Jose Santos Alegria wrote:
>I'm using standard barplot (Windows version 1.6.2 of R) to represent a
certain weekly
>metric "v" and I would like to properly overlay on top of it its moving
average "mean.8"
>(window of 8 weeks). I must be doing something wrong since the moving
average (using >"line
I'm using standard barplot (Windows version 1.6.2 of R) to represent a certain weekly
metric "v" and I would like to properly overlay on top of it its moving average
"mean.8" (window of 8 weeks). I must be doing something wrong since the moving average
(using "lines") doesn't overlay properly, i
library(akima)
plan <- interp.old(x,y.z)
image(plan)
is very close.
This is in MASS (the book) with examples.
On Mon, 3 Feb 2003, Davies, Jonathan wrote:
> I've been told that you may be able to help.
>
> I have a complex linear model with multiple two-way interaction terms. Is
> there a way t
I've been told that you may be able to help.
I have a complex linear model with multiple two-way interaction terms. Is
there a way to view the interaction terms in a 3d plot equivalent to the S
functions:
> plan<-interp(x,y,z)
> image(plan)
where x and y are my explanatory variables and z my re
Hi
1) I am using the function ca() of the package multiv to make a
correspondance analysis of a matrix of categorical datas which are
numericaly coded. I would like to be sure that it is considered as a matrix
of categorical datas and not numerical. I cannot find any explicit mention
of that
- Original Message -
From: "ptremb17" <[EMAIL PROTECTED]>
length(varnames) = 9
length(c(F,T,F,F,F,F,F,F,F,F))) = 10
To: "R-HELP" <[EMAIL PROTECTED]>
Sent: Saturday, February 01, 2003 3:24 AM
Subject: [R] Read.table problem
> Hi !
>
> I am new to R, and using the MAC version onto Mac OS 9
47 matches
Mail list logo