Hi: I have a large data set that I'm testing and I'm finding that
preplot.gam is taking a very long amount of time to compute (like, more
than 20 minutes). My machine is 32-bit, Debian unstable, 4GB memory,
dual Xeon 3GHz. While the data set is very large, the gam() procedure is
able to
Dear useRs,
Release 1.0.0 of the new R package 'randomSurvivalForest' is now available
on CRAN and its mirrors. The package implements Ishwaran and Kogalur's
Random Survival Forests algorithm for right censored survival data. The
algorithm is closely patterned after Breiman's random forests,
Timothy Rye wrote:
I'm interested in clustering my data using the Gower Similarity Coefficient,
and I was wondering if R is capable of using that metric
Timothy Rye
RSiteSearch(Gower) points to a number of relevant messages in the
archives. It also reveals gdist() in the mvpart package,
Rashmi Mathur wrote:
Hello,
How do I split a y-axis to plot data on different scales?
Eg:
x - 1:10
y - c(-0.01,0.79,0.74,0.55,-0.67,0.32,-0.47,-0.05,723,759)
plot(x,y)
I'd like to show these data on the same plot, but the way it's written, all
contrast in the first 8 data points is
[EMAIL PROTECTED]
__
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
and provide commented, minimal, self-contained, reproducible code.
I have been trying for a LONG time to figure out how to do what I
imagine is a fairly simple graph rendering issue. Can anyone help me
figure out how to take two LARGE graphs G1 and G2 (not random) and set
their node/edge attributes separately and uniquely such that when they
are joined, G1
AJ,
hmmm, which package are you using? This might help answering the question,
which is quite hard to decypher in this form. At least for me.
Gabor
On Mon, Aug 14, 2006 at 02:25:17AM -0800, adlai burman wrote:
I have been trying for a LONG time to figure out how to do what I
imagine is a
Dear mailing list
I have a data in two columns and how can i convert it to one row . thank you
in advance
inpute
1 2
3 4
5 6
7 8
9 1
out put
1 2 3 4 5 6 7 8 9 1
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
On 8/14/06, yohannes alazar [EMAIL PROTECTED] wrote:
I have a data in two columns and how can i convert it to one row . thank you
in advance
inpute
1 2
3 4
5 6
7 8
9 1
out put
1 2 3 4 5 6 7 8 9 1
An example follows:
input - matrix(1:10,5,2)
input
[,1] [,2]
[1,]16
Dear all,
I need to calculate tr(A B), tr(A B A B) and similar quantities **fast** where
the matrices A, B are symmetrical. I've searched for built-in functions for
that purpose, but without luck. Can anyone help?
Thanks in advance
Søren
[[alternative HTML version deleted]]
probably something like:
mat - matrix(c(1:9, 1), 5, byrow = TRUE)
c(t(mat))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
input - matrix(1:10,5,2)
input
[,1] [,2]
[1,]16
[2,]27
[3,]38
[4,]49
[5,]5 10
as.vector(input)
[1] 1 2 3 4 5 6 7 8 9 10
Neuro
From: yohannes alazar [EMAIL PROTECTED]
To: r-help r-help@stat.math.ethz.ch
Subject: [R] column to row
Date: Mon, 14
I would suspect that something simple like
sum(diag(crossprod(A,B)))
would be quite competitive...
url:www.econ.uiuc.edu/~rogerRoger Koenker
email[EMAIL PROTECTED]Department of Economics
vox: 217-333-4558University of Illinois
fax:
check the following:
A - matrix(rnorm(100*100), 100, 100); A - A + t(A)
B - matrix(rnorm(100*100), 100, 100); B - B + t(B)
sum(diag(A %*% B))
sum(A * B)
system.time(for(i in 1:1) out - sum(diag(A %*% B)))
system.time(for(i in 1:1) out - sum(A * B))
I hope it helps.
Best,
Dimitris
Chuck == Chuck Cleland [EMAIL PROTECTED]
on Mon, 14 Aug 2006 04:53:55 -0400 writes:
Chuck Timothy Rye wrote:
I'm interested in clustering my data using the Gower
Similarity Coefficient, and I was wondering if R is
capable of using that metric
Timothy Rye
Martin == Martin Maechler [EMAIL PROTECTED]
on Mon, 14 Aug 2006 14:30:39 +0200 writes:
Chuck == Chuck Cleland [EMAIL PROTECTED]
on Mon, 14 Aug 2006 04:53:55 -0400 writes:
Chuck Timothy Rye wrote:
I'm interested in clustering my data using the Gower
Similarity
Hello
One sometimes (quite often really ) marvel at the choice of defaults in
R's graphic engine.
For some obscure reason, mtext uses the typographical descender (bottom
of letters) to align text. That is: gG will end up slightly higher
that GG. Depending on the font, Q might end up higher
Dear list,
I am new to lattice plots.
I want to make a barchart with 10 and more levels.
I need to use a grey scale for printing purposes.
The problem is that with 10 or more levels in factors it is very
difficult to distinguish each level in the plot and legend, since the
greys are very
For importance it's probably best to stick with absolute values of
coefficients, instead of value of the penalty parameter for which the
coefficients changed to non-zero.
Friedman skipped a lot of details on his rule ensemble in that talk, due to
time constraint. In his implementation he was
You could try something like:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/0334.html
replacing the time-stamp with save.image() or save(). Not sure how it works
in BATCH though...
HTH,
Andy
From: John Morrow
Hello fellow R'ers, I have a simple calculation with a very
large data set
That's what I've tried before, on three dual-Xeon boxes, so I know it worked
(as documented a that time).
Andy
From: Paul Y. Peng
Luke Tierney just reminded me that makeCluster() can take a
number greater than the number of machines in a cluster. It
seems to be a solution to this
I need to save data in fixed-width format without headers and reading the
help archive leads me to believe that sprintf is pretty much the only way to
do this. My question is, is there anyway to change the output so the text
in each column is left justified instead of right justified? My code
I have used lme() on data from a between-within subjects experiment. The correct
ANOVA table is known because this is a textbook example (Experimental Design by
Roger Kirk Chapter 12: Split-Plot Factorial Design). The lme() F-values differ
from
the known results. Please help me understand why.
Andreas Svensson [EMAIL PROTECTED] writes:
Hello
One sometimes (quite often really ) marvel at the choice of defaults in
R's graphic engine.
For some obscure reason, mtext uses the typographical descender (bottom
of letters) to align text. That is: gG will end up slightly higher
On 8/14/06, Rafael Duarte [EMAIL PROTECTED] wrote:
Dear list,
I am new to lattice plots.
I want to make a barchart with 10 and more levels.
I need to use a grey scale for printing purposes.
The problem is that with 10 or more levels in factors it is very
difficult to distinguish each level
On Mon, 14 Aug 2006, roger bos wrote:
I need to save data in fixed-width format without headers and reading the
help archive leads me to believe that sprintf is pretty much the only way to
do this. My question is, is there anyway to change the output so the text
in each column is left
Your lme statement is OK. To get the usual split-plot anova, your aov
statement should be
fit2 - aov(y ~ a*b*c + Error(s), data = d)
This gives the same F-values as lme.
On Mon, Aug 14, 2006 at 10:27:19AM -0400, William Simpson wrote:
I have used lme() on data from a between-within subjects
Your lme statement is OK. To get the usual split-plot anova, your aov
statement should be
fit2 - aov(y ~ a*b*c + Error(s), data = d)
No, this gives wrong F-values. By wrong I mean it does not agree with the
published table.
Table 12.10-2, page 559:
Number of obs =
--- Dimitrios Rizopoulos
[EMAIL PROTECTED] wrote:
try the following:
mdf - data.frame(us.state, count, year, month)
mdf[order(mdf$year, mdf$month), ]
Thansk to Dimitris and Dieter. This has helped since
seems to have shown me a way around the problem. It
just means that I have to sort
Hello Folks-
Is there a way to create confidence bands with 'glmmPQL' ???
I am trying to decide a way to best plot the model which we created with the
glmmPQL function in R. I would like
to plot my actual averaged data points within 95 % confidence intervals from
the model. Plotting
the
Didn't get any useful response to the following question. Trying again.
I can't seem to get computationally stable estimates for the following
system:
Y=a+bX+cX^2+dX^3, where X~N(0,1). (Y is expressed as a linear
William Simpson [EMAIL PROTECTED] writes:
Your lme statement is OK. To get the usual split-plot anova, your aov
statement should be
fit2 - aov(y ~ a*b*c + Error(s), data = d)
No, this gives wrong F-values. By wrong I mean it does not agree with the
published table.
Well, it's the
Thanks very much Peter!
Your lme statement is OK. To get the usual split-plot anova, your aov
statement should be
fit2 - aov(y ~ a*b*c + Error(s), data = d)
No, this gives wrong F-values. By wrong I mean it does not agree with the
published table.
Well, it's the model that is
Hi
The reason I use expression is to get italics. Perhaps this is possible
some other way?
/A
Peter Dalgaard wrote:
Andreas Svensson [EMAIL PROTECTED] writes:
Hello
One sometimes (quite often really ) marvel at the choice of defaults in
R's graphic engine.
For some obscure reason,
Hi,
I have a tiny question concerning .Options$max.print
I have to set up this value to a greater value than 1 because I
want to concatenate my output of a function to one single string (for
connivence).
I did this via .Options$max.print - 64000 or options(max.print=64000)
Then I
Andreas Svensson wrote:
Hi
The reason I use expression is to get italics. Perhaps this is possible
some other way?
1. Yes, see ?mtext:
plot(1:10)
mtext(Hello World, 3, font=3)
2. expression() has to use bounding boxes for several reasons, in
particular fractions spring to mind at
The pro's and con's of using scale breaks were discussed by
Cleveland (1985) The Elements of Graphing Data (Wadsworth, pp. 85-91,
149). I don't know what Cleveland said about this is the second edition
of this book, but I believe there are times when scale breaks are
appropriate,
HansJB == Hans-Joerg Bibiko [EMAIL PROTECTED]
on Mon, 14 Aug 2006 17:56:30 +0200 writes:
HansJB Hi, I have a tiny question concerning
HansJB .Options$max.print
HansJB I have to set up this value to a greater value than
HansJB 1 because I want to concatenate my output
Hi,
I am usiing R 2.3.1 under Linux kernel 2.6.11 with libreadline/libhistory 5.1.
Bothe R and readline were compiled without quircks with gcc 3.3.3 and g77.
Every time I try to edit command line by using del key the following error
happens:
8--
ipacks -
At 14:21 13/08/2006, Dirk Enzmann wrote:
I receive R-help messages in digest form that makes it difficult to
answer a post. I noticed that my answer is not added to the thread
(instead, a new thread is started) although I use the same subject
line (starting with Re: ) as the original post. Is
At 13:23 13/08/2006, Michael Zatorsky wrote:
Hi,
I'm working on producing a simple cumulative frequency
distribution.
Thanks to the help of the good people on this list I
now have four vectors that I'd like to join/relate
into a table. e.g.
v1 - myHistogram$breaks # classes
v2 -
Are you saying you did not receive my reply stamped 8/10/2006 2:12
AM in my Sent folder, or that my reply was not useful? In case the
former is correct, my comments were as follows:
Have you tried writing a function to compute SS = sum of squares
deviations between the the left and
Hello,
does anyone know something about missing data in the MCMC pack?
thank you,
Mariagiulia
Mariagiulia Matteucci
Dipartimento di Scienze Statistiche “Paolo Fortunati”
Università di Bologna
Via Belle Arti, 41
40126 Bologna (ITALY)
e-mail: [EMAIL PROTECTED]
TEL: +39 051 264182
FAX: +39 051
Hi,
is possible to make a R script to run under a console without open the R
environment?
Something like this example.R
#!/usr/bin/R
function(name=Put here your name) {
print(name)
}
In a console I make
./example.R name=Ronaldo Reis Júnior
then program print my name.
It is possible?
Thanks
Spencer Graves [EMAIL PROTECTED] writes:
p.s. I don't see how it's obvious that 'a=-c'.
Y=a+bX+cX^2+dX^3, where X~N(0,1). (Y is expressed as a linear combination
of the first three powers of a standard normal variable.) Assuming that
E(Y)=0 and Var(Y)=1, one can obtain the following
Quoting John Kane [EMAIL PROTECTED]:
--- Dimitrios Rizopoulos
[EMAIL PROTECTED] wrote:
try the following:
mdf - data.frame(us.state, count, year, month)
mdf[order(mdf$year, mdf$month), ]
Thansk to Dimitris and Dieter. This has helped since
seems to have shown me a way around the
Dear List,
I am new to TS-Modeling in R. I would like to fit an ARMA(1,1) model
for a balanced panel, running Y on a full set of unit and year dummies
using an arma(1,1) for the disturbance:
y_it=unit.dummies+yeardummies+e_it
where: e_it=d*e_it-1+u_it+q*u_it-1
How can I fit this model in R?
Hello,
I am quite new to R and although I have read some material
(apparently not enough) I'm getting very frustrated. I'm trying to do
some simple circular statistics in R (CircStats, primarily circ.mean)
and I'm confused about the results. What is a negative radian or
degree
Isabelle Deguise ideguise at interchange.ubc.ca writes:
Hello,
x - c(5.009684,37.814266, 295.722970, 0.00,
326.463366,242.678840)
radx - rad(x)
circ.mean(radx)
[1] -0.4283351
deg(-0.4283351)
[1] -24.54179
I would just like this number to be converted to the
On Mon, 14 Aug 2006, Alexandre Aguiar wrote:
Hi,
I am usiing R 2.3.1 under Linux kernel 2.6.11 with
libreadline/libhistory 5.1.
That is known to be buggy. Either patch it or downdate to 5.0.
Bothe R and readline were compiled without quircks with gcc 3.3.3 and g77.
Every time I try
Hi,
Instead of having to program a loop to load several workspaces in a directory,
it would be nice to store the filenames in a list filelist and then to apply
load to this list
lapply( filelist, load)
Unfortunately, although it seems that R is loading the files, the contained
objects are not
--- Dimitrios Rizopoulos
[EMAIL PROTECTED] wrote:
Quoting John Kane [EMAIL PROTECTED]:
--- Dimitrios Rizopoulos
[EMAIL PROTECTED] wrote:
try the following:
mdf - data.frame(us.state, count, year, month)
mdf[order(mdf$year, mdf$month), ]
Thansk to Dimitris and Dieter.
A reproducible example here would help (please see posting guide). A guess:
is your filelist a list of (quoted) character strings? Correct pathnames to
the files with correct separators for your OS?
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the
On Mon, 2006-08-14 at 14:12 -0300, Ronaldo Reis-Jr. wrote:
Hi,
is possible to make a R script to run under a console without open the R
environment?
Something like this example.R
#!/usr/bin/R
function(name=Put here your name) {
print(name)
}
In a console I make
./example.R
On Mon, 14 Aug 2006, Berton Gunter wrote:
A reproducible example here would help (please see posting guide). A guess:
is your filelist a list of (quoted) character strings? Correct pathnames to
the files with correct separators for your OS?
I think the issue is (from the help page)
Usage:
Le 14.08.2006 19:12, Ronaldo Reis-Jr. a écrit :
Hi,
is possible to make a R script to run under a console without open the R
environment?
Something like this example.R
#!/usr/bin/R
function(name=Put here your name) {
print(name)
}
In a console I make
./example.R name=Ronaldo Reis
At 12:21 PM 8/8/2006, [EMAIL PROTECTED] wrote:
We look forward to hearing your comments and inputs on R+ ... please
feel free to suggest a final name for our commercially supported R.
I don´t understanding this mail.
Is possible exist a commercial version of R?
If R source is licensed for GPL
Em Seg 14 Ago 2006 17:58, Marc Schwartz (via MN) escreveu:
http://kavaro.fi/mediawiki/index.php/Using_R_from_the_shell
I made a small change to the wrapper example by implementing dynamic
allocation of memory and sizing the command line buffer to 32768 that is the
real maximum size (at least
# This works fine:
a - 1
b - 2
c - 3
E - expression(a * exp(b*X) + c)
X - c(0.5, 1.0, 2.0)
eval(E)
[1] 5.718282 10.389056 57.598150
D(E, b)
a * (exp(b * X) * X)
eval(D(E, b))
[1] 1.359141 7.389056 109.196300
# But if (a,b,c) are replaced with (A[1], A[2], A[3]),
Je serai absent(e) du 14/08/2006 au 21/08/2006.
Je répondrai à votre message dès mon retour.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Note that it's commercially supported, not just commercial. Also,
commercial does not mean closed source. Open source (and even GPL) and
commercial need not be mutually exclusive. Remember that the free in
free software is as in freedom rather than as in beer (as far as GPL is
concern, at
I think this is the sort of problem which is most elegantly handled by
computing on the language. Here is an INelegant solution:
A - c(1, 2, 3)
for(i in 1:3)assign(paste('A',i,sep=''),A[i])
E - expression(A1 * exp(A2*X) + A3) ## could also use substitute() here,
I think
## instead of
Dear Lister,
I understand merge() can be used to join 2 data frames based on 1 variable.
But how about merge based on more than 2 variables?
Thank you so much!
--
WenSui Liu
(http://spaces.msn.com/statcompute/blog)
Senior Decision Support Analyst
Health Policy and Clinical Effectiveness
Wensui Liu wrote:
Dear Lister,
I understand merge() can be used to join 2 data frames based on 1 variable.
But how about merge based on more than 2 variables?
Thank you so much!
Just specify the 2 (or more) variable names in a column vector for by)
merge(dat1, dat2, by= c(VarA, VarB))
what if the names are different in 2 data frames?
On 8/14/06, Simon Blomberg [EMAIL PROTECTED] wrote:
Wensui Liu wrote:
Dear Lister,
I understand merge() can be used to join 2 data frames based on 1 variable.
But how about merge based on more than 2 variables?
Thank you so much!
Hi all,
I was trying a probit regression using polr() and got this message,
Error in model.matrix.default(Terms, m, contrasts) :
cannot allocate vector of length 828310236
The data is about 20M (a few days ago I asked a question about large file,
thank you for responses, then I use MS
Then instead of by, use by.x and by.y to specifiy the variable names
separately for both data frames. See ?merge, especially the examples.
Wensui Liu wrote:
what if the names are different in 2 data frames?
On 8/14/06, Simon Blomberg [EMAIL PROTECTED] wrote:
Wensui Liu wrote:
67 matches
Mail list logo