We find it equally strange that you posted this!
The advice _is_ in the R-admin manual which the INSTALL file asks you to
read if you have any questions. It covers using enhanced BLAS libraries.
R builds out of the box on FC3, FC4 and Suse on AMD64. I use Goto's BLAS,
but ATLAS can be used
This is a consequence of how VC++ compiled gdal. R is protecting
itself against that. Since R did not crash, there is nothing to worry
about.
On Tue, 2 Aug 2005, Tony Gill wrote:
Hi all,
I just downloaded windows binaries of RGDAL (from
On Mon, 1 Aug 2005, Haibo Huang wrote:
Please refer to the following post.
Which is about Windows only, not Linux. (And on Windows, the answer given
is on the help page for memory.size. together with a better one.)
Ed
--- Mike Lawrence [EMAIL PROTECTED] wrote:
Date: Mon, 1 Aug 2005
On Tue, Aug 02, 2005 at 07:34:39AM +0100, Prof Brian Ripley wrote:
We find it equally strange that you posted this!
The advice _is_ in the R-admin manual which the INSTALL file asks you to
read if you have any questions. It covers using enhanced BLAS libraries.
R builds out of the box on
Hello!
I have two functions.
The first one prepares the arguments for the second one. What is the best
way to put all resoults of the first one into the second one? I tried
attach, however the object in the main enviorment have a priority over the
ones in list. An example is at the end.
On Tue, 2 Aug 2005, Göran Broström wrote:
On Tue, Aug 02, 2005 at 07:34:39AM +0100, Prof Brian Ripley wrote:
We find it equally strange that you posted this!
The advice _is_ in the R-admin manual which the INSTALL file asks you to
read if you have any questions. It covers using enhanced BLAS
?with will help you.
I would avoid using 'list' as a function name, as it will confuse people
and might confuse R too.
On Tue, 2 Aug 2005, [iso-8859-2] Alea }iberna wrote:
Hello!
I have two functions.
The first one prepares the arguments for the second one. What is the best
way to put
Thank you!
However, this does not do exacty what I want. I would like somehow to modify
only the function second.
BTW, I used list only to create a list, it is not one of my functions.
Thanks again,
Ales Ziberna
- Original Message -
From: Prof Brian Ripley [EMAIL PROTECTED]
To: Ales
On Tue, Aug 02, 2005 at 09:33:45AM +0100, Prof Brian Ripley wrote:
On Tue, 2 Aug 2005, Göran Broström wrote:
[...]
Incidentally, I have just tried building R on a Fujitsu Amilo amd64 with
debian-amd64 (unstable) and ATLAS. Both 'make' and 'make check' worked
without complaints with gcc-3.4.5,
On Tue, 2 Aug 2005, Ales Ziberna wrote:
Thank you!
However, this does not do exacty what I want. I would like somehow to modify
only the function second.
What happened when you tried my suggestion as
second - function(l, c) with(l, a + b + c)
? It does work for me.
BTW, I used list
On Tue, 2 Aug 2005, Göran Broström wrote:
On Tue, Aug 02, 2005 at 09:33:45AM +0100, Prof Brian Ripley wrote:
On Tue, 2 Aug 2005, Göran Broström wrote:
[...]
Incidentally, I have just tried building R on a Fujitsu Amilo amd64 with
debian-amd64 (unstable) and ATLAS. Both 'make' and 'make
I apologize, it seams did not interpret your first mail corectly!
Everything works now!
Thank you again!
Ales Ziberna
- Original Message -
From: Prof Brian Ripley [EMAIL PROTECTED]
To: Ales Ziberna [EMAIL PROTECTED]
Cc: R-help r-help@stat.math.ethz.ch
Sent: Tuesday, August 02, 2005 11:25
Chris Wallace [EMAIL PROTECTED] writes:
I am struggling with migrating some stata code to R
Thanks to all who replied. It was very helpful to see a combination
of more direct stata-R translations and more R-ish code. which.max()
solves my problem this time, but learning about split(),
hi all
i am busy teaching a regression analysis course to second year science
students. the course is fairly theoretical with all of the standard
theorems and proofs...
i would like to give the class a practical assignment as well. could you
suggest a good problem and the location of the data
On Tue, Aug 02, 2005 at 10:58:31AM +0100, Prof Brian Ripley wrote:
The problems we are seeing with gcc4 are mainly (but not
entirely) with gfortran: it seems not quite ready for production use.
-4.0.1 is already a considerable improvement over -4.0.0.
There is another project (www.g95.org)
Dear members,
assume given a function of type
test-function(formula, data , w){
..
glm1-glm(formula, family=poisson, data=data, weights=w)
..
}
and a simple example data frame as
test.frame-data.frame(x=1:10,y=(1:10)*2,a=(1:10)^3).
Let us now execute
test(y ~ x, test.frame,
Clark Allan wrote:
i would like to give the class a practical assignment as well. could you
suggest a good problem and the location of the data set/s?
it would be good if the data set has been analysed by a number of other
people so that students can see the different ways of tackling a
I don't think that is the best way to do what I guess you intended. Try
something like
test - function(formula, data , weights)
{
Call - match.call()
Call[[1]] - as.name(glm)
Call$family - quote(poisson)
glm1 - eval.parent(Call)
}
which is probably giving the
Prof Brian Ripley wrote:
This is a consequence of how VC++ compiled gdal. R is protecting
itself against that. Since R did not crash, there is nothing to worry
about.
To add to that: the RGDAL maintainer could work around this bug in
VC++. I haven't used that compiler, but with others
Question on graphs:
The default case for drawing a graph in R, is where a little space is left
on the x and y axis before the first tick i.e. even if I say xlim=c(0,1) --
there will be some space between the edge of the x-axis and where 0 is
placed. If I want 0 on the edge, how do I do it in
Thank you Dimitris Reid, you were very helpful.
Best wishes,
Arie.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
If you have a code that takes 2 weeks to run, than it might be a case of
inefficient algorithm design. I was able to go from overnight runs (SELDI
data analysis) to 20 minute runs by identifying single inefficient function
that took most of the time, and writing it in C.
Jarek
On Tue, 2005-08-02 at 18:20 +0530, Renuka Sane wrote:
Question on graphs:
The default case for drawing a graph in R, is where a little space is left
on the x and y axis before the first tick i.e. even if I say xlim=c(0,1) --
there will be some space between the edge of the x-axis and where
Le 02.08.2005 14:50, Renuka Sane a écrit :
Question on graphs:
The default case for drawing a graph in R, is where a little space is left
on the x and y axis before the first tick i.e. even if I say xlim=c(0,1) --
there will be some space between the edge of the x-axis and where 0 is
placed.
How about:
trapz = function(x, y)
{ # computes the integral of y with respect to x using trapezoidal
integration.
idx = 2:length(x)
return (as.double( (x[idx] - x[idx-1]) %*% (y[idx] + y[idx-1])) / 2)
}
Jarek
\===
Jarek Tuszynski,
Hi,
To find the area lying between the curve y = y(x) and 45 degree line (which,
assuming it goes through the origin, is y = x), you can use the following
function based on trapezoidal rule:
trap.rule - function(x,f) {sum(diff(x)*(f[-1]+f[-length(f)]))/2}
trap.rule(x,f=y-x)
This area will be
Any suggestions about drawing the ridge line on a response surface
represented as a contour plot?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer
Dear All,
is there a simple way to print a data.frame without its row.names?
example:
datum - as.Date(c(2004-01-01, 2004-01-06, 2004-04-12))
content - c('Neujahr', 'Hl 3 K.', 'Ostern')
df1 - data.frame(datum, content)
print(df1)
datum content
1 2004-01-01 Neujahr
2 2004-01-06 Hl 3 K.
3
hi all,
I wish to draw on the same graphic device 3 functions.
But i don't want them to be on different graph, i want to compare them on the
same
I don't need mfrow or mfcol, I need something else...
1 graph on 1 device inside this graph 3 ploted function.
I saw something unsing data.frame,
hi all,
I wish to draw on the same graphic device 3 functions.
But i don't want them to be on different graph, i want to compare them on the
same
I don't need mfrow or mfcol, I need something else...
1 graph on 1 device inside this graph 3 ploted function.
I saw something unsing data.frame,
Le 02.08.2005 15:45, Heinz Tuechler a écrit :
Dear All,
is there a simple way to print a data.frame without its row.names?
example:
datum - as.Date(c(2004-01-01, 2004-01-06, 2004-04-12))
content - c('Neujahr', 'Hl 3 K.', 'Ostern')
df1 - data.frame(datum, content)
print(df1)
datum content
Hi
Maybe matplot (see ?matplot) can help you.
Regards,
Christoph Buser
--
Christoph Buser [EMAIL PROTECTED]
Seminar fuer Statistik, LEO C13
ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND
phone: x-41-44-632-4673
What is it?
---
mod_R is a project dedicated to embedding the R interpreter inside
the Apache 2.0 (and beyond) web server for the purpose of writing web
applications in R. It's composed of three parts:
mod_R: the Apache 2.0 module that implements the glue to load the
Le 02.08.2005 15:56, [EMAIL PROTECTED] a écrit :
hi all,
I wish to draw on the same graphic device 3 functions.
But i don't want them to be on different graph, i want to compare them on the
same
I don't need mfrow or mfcol, I need something else...
1 graph on 1 device inside this graph 3
Romain Francois [EMAIL PROTECTED] writes:
Le 02.08.2005 15:56, [EMAIL PROTECTED] a écrit :
hi all,
I wish to draw on the same graphic device 3 functions.
But i don't want them to be on different graph, i want to compare them on
the same
I don't need mfrow or mfcol, I need something
Dear R-user,
I have a simple question, I just can't figure out a easy way to handle it.
My importing data x is like this:
COL1 COL2 id
1 12 49 1
2 70 120 1
3 58 124 1
51 14 13 2
52 88 100 2
53 90 134 2
I want to change the format of the data, i want to group data into
differenct part
At 16:05 02.08.2005 +0200, Romain Francois wrote:
Le 02.08.2005 15:45, Heinz Tuechler a écrit :
Dear All,
is there a simple way to print a data.frame without its row.names?
example:
datum - as.Date(c(2004-01-01, 2004-01-06, 2004-04-12))
content - c('Neujahr', 'Hl 3 K.', 'Ostern')
df1 -
And you can identify inefficient code fairly easily taking snapshots
from proc.time and computing elapsed time for sections of your code.
spencer graves
Tuszynski, Jaroslaw W. wrote:
If you have a code that takes 2 weeks to run, than it might be a case of
inefficient
use 'split'
x.1 - data.frame(COL1=1:50, COL2=50:1, id=sample(1:4,50,T))
x.2 - split(x.1, x.1$id)
str(x.2)
List of 4
$ 1:`data.frame': 10 obs. of 3 variables:
..$ COL1: int [1:10] 5 10 11 12 22 24 27 34 38 47
..$ COL2: int [1:10] 46 41 40 39 29 27 24 17 13 4
..$ id : int [1:10] 1
The problem with this is that if someone thinks it's useful, they
will likely slog through the obfuscation, then offer something much
simpler to do the same thing -- and then be shocked and confused if the
originator doesn't thank them for their efort! I've heard that S-Plus
allows
Hi all
i need to put on one graph 2 functions who's x axis is the same and y not.
I mean on horizontal the time, and on vertical left: pressure, on vertical
right: rpm of a motor, is R able to do that?
i've found this that i could adapt maybe (i don't need time series really?) :/ :
Heinz == Heinz Tuechler [EMAIL PROTECTED]
on Tue, 02 Aug 2005 17:46:07 +0200 writes:
...
Heinz I tried this, but then the column headers and column
Heinz contents are not aligned.
Use the tabulator if you need them aligned :
On 8/2/05, Spencer Graves [EMAIL PROTECTED] wrote:
And you can identify inefficient code fairly easily taking snapshots
from proc.time and computing elapsed time for sections of your code.
spencer graves
Using Rprof may be a better choice. See
?Rprof
Tuszynski,
Hi,
Search the list, you might be able to find the message but here it is.
Copied From Dr. Brian Ripley's post
x - 1:10
y - rnorm(10)
z - runif(10, 1000, 1)
par(mar=c(5,4,4,4) + 0.1) # Leave space for z axis
plot(x, y)
par(new=T)
plot(x, z, type=l, axes=F, bty=n, xlab=, ylab=)
Martin Maechler [EMAIL PROTECTED] writes:
Heinz == Heinz Tuechler [EMAIL PROTECTED]
on Tue, 02 Aug 2005 17:46:07 +0200 writes:
...
Heinz I tried this, but then the column headers and column
Heinz contents are not aligned.
##data
y1 - matrix(c(3,1,0,1,0,1,1,0,0,0,1,0,0,0,1,1,0,1,1,1),
nrow = 5, byrow = TRUE)
y2 - matrix(c
(3,0,10,3,3,0,0,1,1,0,0,0,0,0,1,0,1,0,0,2,1,0,1,1,0,2,1,1,4,1),
nrow = 5, byrow = TRUE)
y1 - as.data.frame(y1)
y2 - as.data.frame(y2)
rownames(y1) - rownames(y2) -
I have been using lmer in lme4() to analyse the effect of dropping a term from
the model as below
anova.name-anova(m1,m2)
Where m1 is an original model and m2 has one term removed. I can then create
my own type III tables for each variable in the model. However with many
variables this becomes
sp219 [EMAIL PROTECTED] writes:
I have been using lmer in lme4() to analyse the effect of dropping a term
from
the model as below
anova.name-anova(m1,m2)
Where m1 is an original model and m2 has one term removed. I can then create
my own type III tables for each variable in the model.
I am running R 2.0.1 on Windows trying to import a SAS dataset into R
using sas.get.
It is a small SAS dataset, but when I write the command, the SAS window
opens but never appears to complete. Am I doing something wrong?
Here is an example of my command:
temp - sas.get(J:/blah,name)
Jorge
On Tue, 2005-08-02 at 18:27 +0100, Gavin Simpson wrote:
##data
snip
coca(y1 ~ y2, method = symmetric, symmetric= TRUE)
sorry, this should have been:
coca.formula(y1 ~ y2, method = symmetric, symmetric= TRUE)
gives:
Error in model.frame(formula, rownames, variables, varnames, extras,
On Tue, 2 Aug 2005, Spencer Graves wrote:
And you can identify inefficient code fairly easily taking snapshots
from proc.time and computing elapsed time for sections of your code.
Or use the profiler, which makes it much easier. There was a Programmers'
Niche article about it in one
Hi all:
How does one read in an mps file and store it in a matrix as is. I have
pasted the format of the file. I have implemented an mps writer but would
like to have the capability to read in an mps file as well and store it in a
matrix with out using the 'linprog' routines in R. Is there
All,
I have been reading Dr. Harrell's excellent
Regression Modeling Strategies book and trying out
the exercises. I understand that contrast( ) is used
to obtain contrasts between two variables for given
levels of other nuisance variables; is there a way
to use contrast( ) to obtain, for
Dear everyone
I am a new user,would like to combine these code together by using a
loop,each function has three value as Tr = 1 - 3,how can i combine
together?
logitTr1 -logit[logit[,Study]logit[,Tr]==1,]
(number of row in each group (1-3) is difference but equal in colume)
fnTr1 -
Hi, is it possible to increase the memory limit to
infinite so that I don't need to worry about whether
it is enough or not? In S-plus, you can do this by
setting:
options( memory = as.integer( Inf ) )
is it possible to do this in R?
__
Hi, is there an excellent source of documentation for
installation, configuration of R packages
(Rmpi,snow,rsprng etc.) for parallel computing on
linus system?
Thanks
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
array chip wrote:
Hi, is it possible to increase the memory limit to
infinite so that I don't need to worry about whether
it is enough or not? In S-plus, you can do this by
setting:
options( memory = as.integer( Inf ) )
is it possible to do this in R?
You are on Windows, right?
See
Greg Tarpinian wrote:
All,
I have been reading Dr. Harrell's excellent
Regression Modeling Strategies book and trying out
the exercises. I understand that contrast( ) is used
to obtain contrasts between two variables for given
levels of other nuisance variables; is there a way
to use
Hello,
I've encountered the need to cut some chron objects of the form:
R mychron - chron(sort(runif(10, 0, 10)))
R mychron
[1] (01/01/70 16:36:20) (01/02/70 00:08:46) (01/03/70 16:54:49)
[4] (01/04/70 06:45:00) (01/07/70 06:21:24) (01/07/70 18:28:44)
[7] (01/08/70 00:47:05) (01/08/70
Jorge,
My understanding about sas.get is that it actually runs SAS in the
backend, output the SAS data, and then import it into R again.
It might be better to convert SAS to csv format and read it into R.
On 8/2/05, Jorge Sirgo [EMAIL PROTECTED] wrote:
I am running R 2.0.1 on Windows trying to
I'm using R for Mac Aqua version 2.1.1
If at the command line, I type:
ls(
and then a return, the command line prompt changes to
a plus sign indicating that the command has not been
completed. How do I break, or kill, the partial
command? Sometimes, I type a long command but miss a
bracket or
Wensui Liu wrote:
Jorge,
My understanding about sas.get is that it actually runs SAS in the
backend, output the SAS data, and then import it into R again.
He has SAS running locally.
It might be better to convert SAS to csv format and read it into R.
For that, the Hmisc package's
Hi,
Thank you all for the kind reply.
I recompiled R as the previous one turned profiling
off.
I am using package MAANOVA, running the matest
function which is a permutation test. The author did
warn that it takes a long time to run.
Here is one of the test results:
[ Rdata]# ./R CMD Rprof
Steve Adams steve_adams_sd at yahoo.com writes:
Hi, is there an excellent source of documentation for
installation, configuration of R packages
(Rmpi,snow,rsprng etc.) for parallel computing on
linus system?
In case you want it ready-to-run, Rmpi, SNOW, Rsprng, ... work out
of the box under
Thanks to all of you for your help.
As far as I see, the solution of Peter Dalgaard works exactly as I want.
All other solutions have limitations.
Heinz
At 19:19 02.08.2005 +0200, Peter Dalgaard wrote:
Martin Maechler [EMAIL PROTECTED] writes:
Heinz == Heinz Tuechler [EMAIL PROTECTED]
Assuming, as in your post:
set.seed(123)
mychron - chron(sort(runif(10, 0, 10)))
breaks - quantile(mychron)
# is one of these adequate?
cut(mychron, breaks)
cut(unclass(mychron), unclass(breaks), lab = FALSE)
On 8/2/05, Sebastian Luque [EMAIL PROTECTED] wrote:
Hello,
I've encountered
Hello everybody,
I am new to using Windows and R, and have been experimenting with various
packages.
I recently installed R 2.1.1 under Windows XP SP2, and tried installing the
latest versions of SciViews (0.7-6, and R package 0.8-7 found on the
sciviews.org website), and also Tinn-R 1.16.1.5
Gabor Grothendieck [EMAIL PROTECTED] wrote:
Assuming, as in your post:
set.seed(123)
mychron - chron(sort(runif(10, 0, 10)))
breaks - quantile(mychron)
# is one of these adequate?
cut(mychron, breaks)
cut(unclass(mychron), unclass(breaks), lab = FALSE)
Thank you Gabor, that showed me I
How about the following:
set.seed(123)
mychron - chron(sort(runif(10, 0, 10)))
(breaks - chron(pretty(quantile(mychron
[1] 01/01/70 01/03/70 01/05/70 01/07/70 01/09/70 01/11/70
spencer graves
Sebastian Luque wrote:
Gabor Grothendieck [EMAIL PROTECTED] wrote:
Assuming, as in your
Hello,
Can you get eigenvalues in addition to eigevectors using prcomp? If so how?
I am unable to use princomp due to small sample sizes.
Thank you in advance for your help!
Rebecca Young
--
Rebecca Young
Graduate Student
Ecology Evolutionary Biology, Badyaev Lab
University of Arizona
1041 E
I donn't think you can get it directly from prcomp ,but you can get it though
svd.In fact,the prcomp use the svd to do the principal components analysis .
=== 2005-08-03 10:06:56 您在来信中写道:===
Hello,
Can you get eigenvalues in addition to eigevectors using prcomp? If so how?
I
Rebecca Young wrote:
Hello,
Can you get eigenvalues in addition to eigevectors using prcomp? If so how?
I am unable to use princomp due to small sample sizes.
Thank you in advance for your help!
Rebecca Young
Hi, Rebecca,
From ?prcomp:
The calculation is done by a singular
Spencer Graves [EMAIL PROTECTED] wrote:
How about the following:
set.seed(123)
mychron - chron(sort(runif(10, 0, 10)))
(breaks - chron(pretty(quantile(mychron
[1] 01/01/70 01/03/70 01/05/70 01/07/70 01/09/70 01/11/70
I was looking for a way to cut it into specified units of time, rather
Can anyone suggest a simple way to calculate a 'hash code'
from an arbitrary R object?
hash(x) should return an integer or string
with the property that
if hash(x) != hash(y) then x and y are not identical
and the time to compute hash(x) should be quite short.
Any suggestions welcome
On Wed, 2005-08-03 at 12:05 +0800, Adrian Baddeley wrote:
Can anyone suggest a simple way to calculate a 'hash code'
from an arbitrary R object?
hash(x) should return an integer or string
with the property that
if hash(x) != hash(y) then x and y are not identical
and the time to
Dear all--
I am still forging my first arms with R and I am fighting with regexpr() as
well as portability between unix and windoz. I need to extract barcodes from
filenames (which are located between a double and single underscore) as well
as the directory where the filename is residing. Here is
Dear there,
I am wondering how to test whether a simple linear regression model
(e.g. y=1.05x) is significantly different from a 1 to 1 line (i.e. y=x).
Thanks.
Regards,
Jin
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch
Try this. The regular expression says to match
- anything
- followed by a double underscore
- followed by one or more digits
- followed by an underscore
- followed by anything.
The digits have been parenthesized so that they can be referred to in
the backreference \\1.Also use the R
Dear all,
I am trying to replicate some multilevel models with binary outcomes
using R's lmer and glmmPQL and Stata's gllmm, respectively.
The data can be found at http://www.uni-koeln.de/~ahf34/xerop.dta.
The relevant Stata output can be found at http://www.uni-
79 matches
Mail list logo