At least in R 1.8.0 it is a list matrix:
class(jank)
[1] matrix
typeof(jank)
[1] list
The help page is not quite correct, as it does not mention what happens
when you create a data frame with a *list* as a column.
How did you create this data frame? Columns cortrange and logcortrange
are
On Sat, 18 Oct 2003, Ernie Adorio wrote:
Dear R experts,
1. How can I turn off the display of rulers in image() command?
Set axes=FALSE in image(). In this case - as in many others - running
example() - here example(image) - on the function does point in the right
direction. The function
PaulSch == Schwarz, Paul [EMAIL PROTECTED]
on Wed, 15 Oct 2003 12:09:11 -0700 writes:
PaulSch I am converting some S-PLUS scripts that I use for
PaulSch creating manuscript figures to R so that I can take
PaulSch advantage of the plotmath capabilities. In my
PaulSch S-PLUS
Richard A. O'Keefe [EMAIL PROTECTED] writes:
Scott Norton [EMAIL PROTECTED] wrote:
I have a list of character vectors. I'm trying to see if
there is a way (in a single line, without a loop) to pull out
the first element of all the vectors contained in the list.
You
Hi,
Following plot is displaying fine, (starting arount the 10. september),
except that the xaxp parameter has no effect.
I'd like to have a tickmark every 7 days...
plot(timeline,
subset(myd, TYPE==A)$list1,
ylim=c(100*floor(min(subset(myd, TYPE==A)$list1)/100-1),
Dear R-users,
I found a strange problem
working with products of two matrices, say:
a - A[i, ] ; crossprod(a)
where i is a set of integers selecting rows. When i is empty
the result is in a sense random.
After some trials the right answer
(a matrix of zeros) appears.
---
I agree completely.
In fact, I have about 5000 observations, which should
be enough.
I was using 200 samples because of RAM limitations and
I'm afraid to think about what amount of RAM I'll
need to fit an aov() for such data.
--- John Fox [EMAIL PROTECTED] wrote:
Dear Alexander,
If I
--- Deepayan Sarkar [EMAIL PROTECTED] wrote:
On Thursday 16 October 2003 19:03, Alexander
Sirotkin \[at Yahoo\] wrote:
Thanks for all the help on my previous
questions.
One more (hopefully last one) : I've been very
surprised when I tried to fit a model (using
aov())
Hello.
I have a problem using sort() in windows 2000 and windows NT 4.0, running R
1.8.0 on both. I want to sort a vector of characters names, where I have
used Scandinavian letters, like 'Æ', 'Ø', and 'Å' (for those who cannot
display these letters this question seems rather meaningless, i
Hi!
How can I fill with colors a portion of a graph (e.g.: I want fill in red the
area within two confidence intervals)?
Thank you very much
Cristian
~~
Cristian Pattaro
~~
Unit of Epidemiology Medical Statistics
I still have R 1.7.1, and the problem appears there as well:
a - A[i, ] ; crossprod(a)
[,1] [,2] [,3] [,4] [,5]
[1,] 1.195616e-301 7.042305e-302 9.563047e-302 2.281448e-302 2.198017e-302
[2,] 6.905419e-302 1.204915e-301 3.382433e-302
Hi all,
By default, the heatmap function gives an image with a dendrogram added
to the
left side and to the top. Is it possible to only add the dendrogram to
the left side
and let the order of the columns unchanged ?
I tried
heatmap(mat, col=rbg,Rowv=res.hclust$order,Colv=1:dim(mat)[[2]]).
In
[EMAIL PROTECTED] wrote:
Hi!
How can I fill with colors a portion of a graph (e.g.: I want fill in red the
area within two confidence intervals)?
You can construct the coordinates of the polygon that fills this
region, then use 'polygon' to fill it. Here:
# first set up the plot - we want the
Ivar Herfindal wrote:
Hello.
I have a problem using sort() in windows 2000 and windows NT 4.0,
running R 1.8.0 on both. I want to sort a vector of characters names,
where I have used Scandinavian letters, like 'Æ', 'Ø', and 'Å' (for
those who cannot display these letters this question seems
Or
do.call(cbind,x)[1,]
which of course makes a whole new copy of x and
gives you a nasty warning as well, but does not
use a conceptual `for` loop. Which I think was the
original question, to which AFAIK the answer is no, there is no
easy subscripting construct such as x[[1:3]][1] that will do
Hi All,
I've the following data frame with 54 rows and 4 colums:
x
Ratio Dose Time Batch
R.010mM.04h.NEW0.02 010mM 04h NEW
R.010mM.04h.NEW.1 0.07 010mM 04h NEW
...
R.010mM.24h.NEW.2 0.06 010mM 24h NEW
R.010mM.04h.OLD0.19 010mM 04h OLD
Sorry, I just figured it out: x[x$Batch == 'OLD',] instead of x[x$Batch ==
'OLD']. I didn't know this has to be in the same format then x[1:20,] where I
already used the comma.
sorry for posting the previous message ...
Arne
-Original Message-
From: [EMAIL PROTECTED]
On Fri, 17 Oct 2003 [EMAIL PROTECTED] wrote:
I've the following data frame with 54 rows and 4 colums:
x
Ratio Dose Time Batch
R.010mM.04h.NEW0.02 010mM 04h NEW
R.010mM.04h.NEW.1 0.07 010mM 04h NEW
...
R.010mM.24h.NEW.2 0.06 010mM 24h
On Wed, 15 Oct 2003 15:47:59 +0200
Uwe Ligges [EMAIL PROTECTED] wrote:
Shawn Way wrote:
I'm looking for design and hmisc version 2.0 for R 1.8 for windows. I've
found design 2.0 in the downloads for R1.7 but not hmisc.
I've also checked Dr. Harrell's site and it only goes to 1.6
Hi,
thanks for your replies regarding the problem to select a sub data frame by
expression. I start getting an understanding on how indexing works in R.
thanks for your replies,
Arne
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: 17
A couple of comments:
o Methods such as decision trees do not need to expand factors into columns
of 1df contrasts, so the memory requirement is vastly different. The models
produced is also very, very different.
o Why would you want all possible interactions of 10 variables, 6 of
which are
Dear Alexander,
At 01:29 AM 10/17/2003 -0700, Alexander Sirotkin \[at Yahoo\] wrote:
I agree completely.
In fact, I have about 5000 observations, which should
be enough.
I was using 200 samples because of RAM limitations and
I'm afraid to think about what amount of RAM I'll
need to fit an aov()
From: Gamal Abdel-Azim [mailto:[EMAIL PROTECTED]
While trying to expand the memory/object size in R, I noticed
that R might be using
only heap memory. Is this true? Are all objects in R created
in the heap not
allocated? It's not logical that this is the case!! Otherwise
the whole R
On 17 Oct 2003 at 1:33, Alexander Sirotkin [at Yahoo] wrote:
You mentioned in an earlier post that at least one of your factors
have 40 levels. If you use the default contrast, contrast.traetment,
the design matrix for this factor will be dominated by zeros. Maybe
you shoukd look at tha CRAN
If I take the following simple data:
YEAR MONTH DAY WEIGHT.KG
2003 10 6 1.2
2003 10 12 1.2
2003 10 16 1.3
and format the date data and plot it:
dates - strptime(paste(DAY,MONTH,YEAR),%d%m%Y)
plot(c(min(dates),max(dates)),c(0,max(WEIGHT.KG)),
xlab=Date,ylab=Weight (kg),type=n)
Thank you very much...
Shawn Way
-Original Message-
From: Frank E Harrell Jr [mailto:[EMAIL PROTECTED]
Sent: Friday, October 17, 2003 6:38 AM
To: Uwe Ligges
Cc: [EMAIL PROTECTED]; [EMAIL PROTECTED]
Subject: Re: [R] Design and Hmisc
On Wed, 15 Oct 2003 15:47:59 +0200
Uwe Ligges [EMAIL
Somehow R creates `a' as a matrix with 0 rows and 5 columns. I don't know
how crossprod() or other linear algebra functions deals with such a
degenerate matrix.
I'd suggest R Core to add checks for strictly positive dimensions in such
functions.
(Also, I find it strange that A[1,] is a vector,
On Fri, 17 Oct 2003 11:58:49 +0200, Uwe Ligges [EMAIL PROTECTED]
dortmund.de wrote:
Ivar Herfindal wrote:
Hello.
I have a problem using sort() in windows 2000 and windows NT 4.0,
running R 1.8.0 on both. I want to sort a vector of characters names,
where I have used Scandinavian letters,
On Fri, 17 Oct 2003, Alexander Sirotkin [at Yahoo] wrote:
I did study ANOVA and I do have enough observations.
200 was only a random sample of more then 5000 which I
think should be enough. However, I'm afraid to even
think about amount of RAM I will need with R to fit a
model for this
On Fri, 17 Oct 2003, Liaw, Andy wrote:
From: Gamal Abdel-Azim [mailto:[EMAIL PROTECTED]
While trying to expand the memory/object size in R, I noticed
that R might be using
only heap memory. Is this true? Are all objects in R created
in the heap not
allocated?
To the extent that
On Fri, 17 Oct 2003, Liaw, Andy wrote:
Somehow R creates `a' as a matrix with 0 rows and 5 columns. I don't know
how crossprod() or other linear algebra functions deals with such a
degenerate matrix.
I'd suggest R Core to add checks for strictly positive dimensions in such
functions.
Yes.
You can set any locale you like, and I suspect your machines are in
different locales (I believe older versions of Windows, including NT4, had
limited support for locales).
On Fri, 17 Oct 2003, Ivar Herfindal wrote:
On Fri, 17 Oct 2003 11:58:49 +0200, Uwe Ligges [EMAIL PROTECTED]
On Fri, 17 Oct 2003, Liaw, Andy wrote:
Somehow R creates `a' as a matrix with 0 rows and 5 columns. I don't know
how crossprod() or other linear algebra functions deals with such a
degenerate matrix.
I'd suggest R Core to add checks for strictly positive dimensions in such
functions.
On Friday 17 October 2003 02:20, Martin Maechler wrote:
PaulSch == Schwarz, Paul [EMAIL PROTECTED]
on Wed, 15 Oct 2003 12:09:11 -0700 writes:
PaulSch I am converting some S-PLUS scripts that I use for
PaulSch creating manuscript figures to R so that I can take
PaulSch
Dear all,
I am trying to compute a matrix of Pearson's `r' or Spearman's `rho'
rank correlation coefficients using rcorr (Hmisc) the following way:
mx-rcorr(x, type=spearman)[1]
but then ...
is.matrix(mx)
[1] FALSE
Even if I use as.matrix the result is not better.
What can I do?
Thank you
On Fri, 17 Oct 2003 16:36:47 +0200
Luca De Benedictis [EMAIL PROTECTED] wrote:
Dear all,
I am trying to compute a matrix of Pearson's `r' or Spearman's `rho'
rank correlation coefficients using rcorr (Hmisc) the following way:
mx-rcorr(x, type=spearman)[1]
Instead of [1] use $r. Or
On Friday 17 October 2003 03:33, Alexander Sirotkin \[at Yahoo\] wrote:
One more (hopefully last one) : I've been very
surprised when I tried to fit a model (using
aov())
for a sample of size 200 and 10 variables and
their interactions.
That doesn't really say
Regarding the very few tests I did (RMySQL versus RODBC using a MySQL ODBC
driver, but I do not remember details here), RMySQL is faster. It should be
great, if you need to access a MySQL database from R, to try both and decide
by yourself. If you do that, I am very interested by the results.
Hi,
I have a string representing an environment:
bob
And an environment
bob
environment: 0x3901234ac
How do write a function that takes the string and returns the
environment?
Crispin
This email is confidential and intended solely
Philippe Grosjean wrote:
Regarding the very few tests I did (RMySQL versus RODBC using a MySQL ODBC
driver, but I do not remember details here), RMySQL is faster. It should be
great, if you need to access a MySQL database from R, to try both and decide
by yourself. If you do that, I am very
I am not seeing this on Linux. The x axis marks are at midnight GMT,
hence 1am BST on my system.
On Fri, 17 Oct 2003, Jacob Etches wrote:
If I take the following simple data:
YEAR MONTH DAY WEIGHT.KG
2003 10 6 1.2
2003 10 12 1.2
2003 10 16 1.3
and format the date data and plot it:
Currently, I am intrested in parameter fitting of Generalized Lambda
Distribution.And I have found two packages in R related to Gld,Davies and gld. What's
a pity that no method in Davies deals with fitting of gld,and starship used in
package:gld is quite time-consuming when sample size is
Crispin Miller wrote:
Hi,
I have a string representing an environment:
bob
And an environment
bob
environment: 0x3901234ac
How do write a function that takes the string and returns the
environment?
get(bob)
Uwe Ligges
Crispin
This
Is get(bob) what you are looking for?
It is the usual way to go from the name of an R object (as a character
string) to the actual object.
On Fri, 17 Oct 2003, Crispin Miller wrote:
Hi,
I have a string representing an environment:
bob
And an environment
bob
environment: 0x3901234ac
AndyL == Liaw, Andy [EMAIL PROTECTED]
on Fri, 17 Oct 2003 09:10:16 -0400 writes:
AndyL One of the good thing about R (and S in general, I
AndyL guess) is that if a function does mostly what you
AndyL want, except for some small things, you can just make
AndyL another copy of
On Thursday 16 October 2003 10:03, Stefán Hrafn Jónsson wrote:
Dear R community.
I have two problems with figures. First deals with short vector on the
x-axis and the second with two-panel barchart.
1) For demonstration I create the following pseudo data for three years,
2001:2003. The
WordofMouthConnection.com Search Awareness System
This email is a website-generated message, but it is not spam.
An acquaintance of yours recently conducted a search on your email address in our
online community, WordofMouthConnection.com. It could be a friend, a family member,
co-worker,
I am creating lists of vectors withing a loop. I also would like to
change the dim attribute to the vectors in order to make them
matrices.
I have tried the following, but it doesn't work...
sim - c('simMeans','simVars','simWeights')
indexTable - table(modelIndex)
for (i in sim) {
+
We are currently looking for a Director Biostatistics and Data Management in
our Cambridge MA facility.
We are looking for 7+ years of Statistical Analysis in a Biotech /
Pharmaceutical environment with Phase 2 and Phase 3 clinical trial data.
Experience with clinical protocol design,
The problem is related to time zones. The easiest way to
handle this is to avoid using POSIXt and use chron instead
so you don't have to worry about them.
require(chron)
day - 6:16
dts - dates(paste(10, day, 03, sep=/))
plot(dts,day)
abline(v=dts)
---
From: Jacob Etches [EMAIL PROTECTED]
WordofMouthConnection.com Search Awareness System
This email is a website-generated message, but it is not spam.
An acquaintance of yours recently conducted a search on your email address in our
online community, WordofMouthConnection.com. It could be a friend, a family member,
co-worker,
Ernie Adorio wrote:
Dear R experts,
1. How can I turn off the display of rulers in image() command?
Are rulers the same as axes? If so, try
image(...,axes = FALSE)
See example(image) for detials.
2. Rather than typing my commands at the command line, how can I input a file
contents aside
On Fri, Oct 17, 2003 at 01:16:00PM -0400, Gabor Grothendieck wrote:
The problem is related to time zones. The easiest way to
handle this is to avoid using POSIXt and use chron instead
so you don't have to worry about them.
require(chron)
day - 6:16
dts - dates(paste(10, day, 03,
Hello everybody,
I would like to perform a test for normality (without specifying the
mean a variance) on the sample data (80 observations). I found that
Lilliefors test is appropriate. Does anybody have it programmed already,
or is there a function for this test in R?
Thank you very much,
On 13 Oct 2003 at 8:45, Martin Maechler wrote:
kjetil == kjetil halvorsen [EMAIL PROTECTED]
on Sun, 12 Oct 2003 09:55:00 -0400 writes:
kjetil Hola! I have the following in a .Rd file:
kjetil \eqn{\mbox{coef} = c(\mbox{coef}[1],\ldots, \mbox{coef}[n]) }
kjetil
From: Dirk Eddelbuettel [EMAIL PROTECTED]
On Fri, Oct 17, 2003 at 01:16:00PM -0400, Gabor Grothendieck wrote:
The problem is related to time zones. The easiest way to
handle this is to avoid using POSIXt and use chron instead
so you don't have to worry about them.
I do see the described behavior, on three systems, linux R 1.8.0, Mac
OS X R 1.8.0, and Solaris R 1.7.1.
Plot 1 is different than plot 2; in plot 1 the points are offset to
the left of the axis tick marks.
datet - as.POSIXct(dates)
## 1
plot(datet,WEIGHT.KG)
## 2
plot(datet,WEIGHT.KG,xaxt='n')
Hi,
I wish to compile the R source on SuSE,
but am unable to find the gcc for it.
Can anyone send me a pointer of where
they got it from.
Thanks,
Dipti
--
/\ Dipti Kamdar
\\ \ Solution
On 12 Oct 2003 at 21:45, Uwe Ligges wrote:
[EMAIL PROTECTED] wrote:
I am running
Rcmd check (Windows XP, rw1080 from cran)
on a new package.
This reports undocumented code objects for 14 functions,
which all have their .Rd files!
What might be happening?
1) You
On 17 Oct 2003 at 13:59, Martina Pavlicova wrote:
There is shapiro.test in package ctest, which have much better power
properties than Lillefors test. So there is no need to have
Lilliefors test in R, except for archeological interest.
Kjetil Halvorsen
Hello everybody,
I would like to
So someone forgot to specify the timezone, if the current one was not
wanted.
However, I don't see how timezones can account for a 24hour difference
as originally reported.
On Fri, 17 Oct 2003, Don MacQueen wrote:
I do see the described behavior, on three systems, linux R 1.8.0, Mac
OS X R
On Fri, 17 Oct 2003 [EMAIL PROTECTED] wrote:
On 12 Oct 2003 at 21:45, Uwe Ligges wrote:
[EMAIL PROTECTED] wrote:
I am running
Rcmd check (Windows XP, rw1080 from cran)
on a new package.
This reports undocumented code objects for 14 functions,
which all have
At Friday 02:20 PM 10/17/2003 -0400, Gabor Grothendieck wrote:
[material deleted]
Time zones are not part of the problem yet POSIXt forces this
extraneous complication on you. chron has no time zones in the
first place and therefore allows you to work in the natural frame
of the problem, avoiding
From: Tony Plate [EMAIL PROTECTED]
I also see the usefulness of a time-zone-free time/date class,
but why does chron need to be moved to the base to be useful here?
Because other software makes use of times in the base. Package
writers figure that what is in the base is the most available
Dipti Kamdar wrote:
Hi,
I wish to compile the R source on SuSE,
but am unable to find the gcc for it.
Can anyone send me a pointer of where
they got it from.
Thanks,
Dipti
Why do you think installing gcc is a topic related to R-help?
You can install gcc from rpms that are on the SuSE
Thanks for all the responses.
After re-examining my data I came to realize that
second order interactions would be enough in my
particular case. With second order instructions I
managed to fit a model with less then 512MB RAM.
Thanks to everybody.
--- John Fox [EMAIL PROTECTED] wrote:
Dear
I am confused by the following difference in the behavior of R and S. Any
clarification would be greatly appreciated.
Jan Erik Backlund
Dow AgroSciences, LLC.
[EMAIL PROTECTED]
R : Copyright 2003, The R Development Core Team
Version 1.8.0 (2003-10-08)
sw
Fertility Agriculture
On Fri, 17 Oct 2003, Backlund, Jan Erik (JE) wrote:
I am confused by the following difference in the behavior of R and S. Any
clarification would be greatly appreciated.
sw[[2,1]] in R is short for sw[[2]][[1]], which in the case of a data
frame is sw[1,2], as your example shows.
See ?[[
[[ operates recursively, so
sw[[2,1]] is the same as sw[[2]][[1]]
Giovanni
Date: Fri, 17 Oct 2003 12:30:07 -0400
From: Backlund, Jan Erik (JE) [EMAIL PROTECTED]
Sender: [EMAIL PROTECTED]
Cc:
Precedence: list
I am confused by the following difference in the behavior of R and
Version 1.4 of package g.data is available on CRAN. This upgrade is
necessary for it to work under R-1.8.0, and is fully backward compatible.
Description: Create and maintain delayed-data packages (DDP's). Data stored in
a DDP are available on demand, but do not take up memory until
Yes. The timezone is not the whole problem.
What one would really like is that plot understands that it is
being given daily data and acts accordingly, in the same
way that plot already understands that its being given a factor
object or a dendogram, etc. and produces the right plot.
The OO
I've been working on a new package and I have a few questions regarding the
behaviour of the nlm function. I've been (for better or worse) using the nlm
function to fit a linear model without suppling the hessian or gradient
attributes in the objective function. I'm curious as to why the nlm
Dear R users,
Does anyone know why the following two ways to calculate
correlation variance give different answers? I also obtain
different answers when I use, say, spearman method in
cor(). The problem does not happen in R 1.7.1 (pearson
correlation only, of course in R 1.7.1).
set.seed(1234)
73 matches
Mail list logo