Hello,
Just look at the examples in ?persp.
There is a function called trans3d defined in it that will traduce your 3D
coordinates to 2D, and so you will be able to draw lines with lines function.
Romain.
Corey Bradshaw a écrit :
>I've created a perspective plot using 'persp' in the graphics p
Dear list!
When using the German special characters I didn't see all characters in the
correct way. Let's take the command
ff <- "äöüßÄÖÜ" # (for those ones who can't see this in the correct way \"a,
\"ö, \"ü
# \ss (?), \"A, \"O and \"U in LaTeX commands
'ff' will show the
Thanks, Romain. I've certainly used that to draw lines and points in the plots
produced by 'persp'; however, my problem is that I need to incorporate an
entirely new z function (not just a plane) onto the same plot (i.e., using the
same x and y values).
Corey
-Original Message-
From: [
Hi there, I need help on the following problem. Any suggestions will be
much appreciated.
Let x be an i.i.d. random sample of a continuous variable. For some fixed
numbers a1
1) the frequency of x falls in each interval (a1, a2], (a2, a3], (a3, a4],...
2) the conditional mean of x within each in
I tried lrm in library(Design) but there is always
some error message. Is this function really doing the
weighted logistic regression as maximizing the
following likelihood:
\sum w_i*(y_i*\beta*x_i-log(1+exp(\beta*x_i)))
Does anybody know a better way to fit this kind of
model in R?
FYI: one exa
While it is true that the large number of variables relative to
the number of observations restricts what can be inferred,
the situation is not as hopeless as Bert seems to suggest.
If it were, attempts at the analysis of expression array data
would be a waste to time. Methods developed to that
ge
I've created a perspective plot using 'persp' in the graphics package.
I'd like to add a second plane of z values to the existing plot, but I
cannot seem to do this using 'persp'. Is there an analogue to 'lines' or
'points' for perspectives?
Corey.
[EMAIL PROTECTED]
__
sorry, I had a typo there, it's datadist(b) for the
analysis of data frame "b".
--- Robert Balshaw <[EMAIL PROTECTED]> wrote:
> Not sure if this will help, but did you mean to use
> datadist(a) for
> the analysis of B?
>
> Rob
>
> > -Original Message-
> > From: [EMAIL PROTECTED]
> > [ma
Is there anything wrong with vectorized calculation? E.g.,
dd <- (M-10):M
r <- 0.5 * (1 - exp(-2 * (32 - dd) / 100))
round(r, 4)
Andy
> From: bogdan romocea
>
> Not sure if it's the best way, but you could do it this way:
> all.results <- vector(mode="numeric")
> for (i in 1:100)
> {
>
actually data frame "b" has 177 rows, the script ran
ok on the first 166 rows as a subset, but started to
break down if subset of the first 177 rows were used
as the dataset, or the first 166 rows plus 168th row,
the data in those rows are:
time status scores
165 172 0 -123.3
166
Hi, I encountered a weird problem when using the
Design and Hmisc problem. I have 2 data frame called
"a" and "b", both have 3 columns: "time", "status" and
"scores", a sample of the data frame is like:
data frame "a":
time status scores
1 21 1 99.61
2 38 0 101.11
3 51 0 1
"Liaw, Andy" <[EMAIL PROTECTED]> writes:
> Hi all,
>
> I've only now noticed that on our Opteron boxes running Suse Enterprise
> Server 8 that R-2.0.1 did not pass make chek-all. It bombed in
> tests/lapack.R at:
>
> > eigenok <- function(A, E, Eps = 1000 * .Machine$double.eps) {
> V <- E$v
Thanks for reminding me of that thread, Roger. I've forgotten all about
it...
Strangely, several people reported no problem with optimized BLAS, but I
tried linking against the latest version (0.96, threaded), and got the same
failure. I will try to see if I can upgrade GCC.
Best,
Andy
> From:
Again, I have to tell...
Take a look at Ubuntu
http://www.ubuntulinux.org/
When I get tired of MAC OSX, I installed the ppc version of Ubuntu on my
IBook. It's great and it's based on Debian.
On Mon, 2004-12-13 at 12:35 -0500, doktora v wrote:
> I'm using SUSE with success on intel laptop and
I am trying to find Multivariate Multiple Regression documentation for R. Could
somebody please indicate me some, if available at all.
Thanks
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/
Hi
Dr Carbon wrote:
How do I draw a rectangle across multiple plots on a device?
E.g.,
def.par <- par(no.readonly = TRUE)
par(mfrow = c(3, 1))
plot(1:10, rnorm(10), ylim = c(-4,4), type = "l")
plot(1:10, rnorm(10), ylim = c(-4,4), type = "l")
plot(1:10, rnorm(10), ylim = c(-4,4), type = "l")
rect(2
Claus Dethlefsen math.aau.dk> writes:
>
> Dear list
>
> I would like to be able to group terms in a formula using a function that I
> will call tvar(), eg. the formula
>
> Y ~ 1 + tvar(x:A) + tvar(z) + u + tvar(B) + tvar(poly(v,3))
>
> where x,u and v are numeric and A and B are fac
Hi all,
I am trying to set up r-mathlib so I can access R
functions from Ruby using ruby-rmathlib, and for the
life of me I can't find an rpm or source code for
r-mathlib (I can find debian packages ...). I'm using
Red Hat Enterprise Linux WS (V 3) (EM_64 bit
processor) and R version 1.9.1.
Than
Not sure if it's the best way, but you could do it this way:
all.results <- vector(mode="numeric")
for (i in 1:100)
{
...
this.run <- ...
all.results <- c(all.results,this.run)
}
At this point all.results contains the values of this.run from the
whole loop. I
?aov or ?manova
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
"The business of the statistician is to catalyze the scientific learning
process." - George E. P. Box
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of
>
I am looking for a simple implementation of a Chebyshev type I
low-pass filter in R. Searching through the R web site I only found
similar requests but never an answer pointing to a package providing
such standard filters (Chebyshev, Butterworth). Any idea where I
might find this?
Thanks, Chri
A little bit late, but I have to tell...
Take a look at Ubuntu
http://www.ubuntulinux.org
Is both comfortable to install and maintain; plus easy to expand.
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
Hiya,
I have been struggling to save the result from the FOR loop. What is the
best way to do it, as I need the result to merge with another dataset for
further analysis ?
for (dd in ((M-10):M)){
+ dist<-(32-dd)
+ r<-1/2*(1-exp(-2*dist/100))
+ map<-c(dd,round(r,4))
+ print(map)
+ next
+ }
Thanks
How do I draw a rectangle across multiple plots on a device?
E.g.,
def.par <- par(no.readonly = TRUE)
par(mfrow = c(3, 1))
plot(1:10, rnorm(10), ylim = c(-4,4), type = "l")
plot(1:10, rnorm(10), ylim = c(-4,4), type = "l")
plot(1:10, rnorm(10), ylim = c(-4,4), type = "l")
rect(2, -4, 3, 4)
par(de
It looks just like the classical F-test for lack-of-fit, using estimate of
`pure errors' from replicates, doesn't it? This should be in most applied
regression books. The power (i.e., probability of finding lack-of-fit when
it exists) of such tests will depend on the data.
Andy
> From: RenE J.V
Martin Maechler brought this up a while back (and added the regression
test). I believe it was a compiler problem and upgrading to gcc 3.4.1
fixed it.
See the thread starting here:
https://stat.ethz.ch/pipermail/r-devel/2004-July/030260.html
-roger
Liaw, Andy wrote:
Hi all,
I've only now notice
Hi all,
I've only now noticed that on our Opteron boxes running Suse Enterprise
Server 8 that R-2.0.1 did not pass make chek-all. It bombed in
tests/lapack.R at:
> eigenok <- function(A, E, Eps = 1000 * .Machine$double.eps) {
V <- E$vect
lam <- E$values
stopifnot(abs(A %*% V - V %*%
Halldor Björnsson vedur.is> writes:
:
: Hi,
: I am trying to understand the behaviour of the plot function.
: If I have
:
: novdate <- as.Date("2001/11/1") + (0:29)
: y <- 1:30
: b <- data.frame(novdate,y)
:
: then plot(b$novdate,b$y) will produce a plot where the x-ticmarks are
: giv
Sebastian:
I routinely deal with situations like that in one of two ways:
(1) Can you work with log(besselI)? If yes, that should solve the
problem.
(2) What do you do with the numbers returned from besselI? I
assume they are later used to compute numbers in a more sensible
On Tue, 14 Dec 2004, [ISO-8859-1] Halldor Björnsson wrote:
>
> Hi,
> I am trying to understand the behaviour of the plot function.
> If I have
>
> novdate <- as.Date("2001/11/1") + (0:29)
> y <- 1:30
> b <- data.frame(novdate,y)
>
> then plot(b$novdate,b$y) will produce a plot where the x-ti
Hi,
I am trying to understand the behaviour of the plot function.
If I have
novdate <- as.Date("2001/11/1") + (0:29)
y <- 1:30
b <- data.frame(novdate,y)
then plot(b$novdate,b$y) will produce a plot where the x-ticmarks are
given as dates (Nov 04, Nov 09 etc), but plot(b) will produce a plot
wh
Hu Chen <[EMAIL PROTECTED]> writes:
> Hi, I got a amd 64 + debian/pure64 system.
> R is installed by apt-get so it's no doubt that R is for 64bit.
> My question is, for those packages, for example, ipred, it's difficult
> to find in which official debian packages they are. How could I
> install th
I don't see how this is the `same' question at all. The object Martin has
doesn't look like a data frame (can you tell why?), so I believe is
unrelated to what you asked before.
Andy
ps: Why is this copied to [EMAIL PROTECTED]
> From: bogdan romocea
>
> I asked the same question a few weeks a
The problem, I think, boils down to the following:
> x <- matrix(c(3, 1, 2), ncol=3, dimnames=list(NULL, letters[1:3]))
> x
a b c
[1,] 3 1 2
> sort(x)
a b c
[1,] 1 2 3
> y <- rbind(x, c(1,1,2))
> y
a b c
[1,] 3 1 2
[2,] 1 1 2
> sort(y)
a b c
[1,] 1 1 2
[2,] 1 2 3
I.e., I'm gue
El Martes, 14 de Diciembre de 2004 12:30, jing tang escribió:
> Hi,
> I have several small question in R,
> 1) How to display all the variables in current workspace?
ls()
> 2) How to write a long command in two lines. Suppose one command line is
> long to be put within one line.
just push 'Enter
> "bogdan" == bogdan romocea <[EMAIL PROTECTED]>
> on Mon, 13 Dec 2004 12:26:46 -0800 (PST) writes:
bogdan> A simple for loop does the job. Why not write your own function?
movsd <- function(series,lag)
{
movingsd <- vector(mode="numeric")
for (i in lag:length(
Jari Oksanen wrote:
On 13 Dec 2004, at 19:53, doktora v wrote:
I'm looking to switch to Mac platform. Anyone had any experience
with that? I'm expecting on a power G4 laptop later this week hope
R behaves...
I have been a Linux user since 1999, and I got my first ever Mac
(iBook G4 laptop) la
> "Thomas" == Thomas Isenbarger <[EMAIL PROTECTED]>
> on Mon, 13 Dec 2004 14:27:56 -0600 writes:
Thomas> R-help folks:
Thomas> Thanks in the past for your help. I have another
Thomas> question that I hope has a simple answer. I have
Thomas> searched the R home pages
Hi,
Is this the mailing address of Chunlou Yung who used to work at Perform.com
-Suhail
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
39 matches
Mail list logo
