On Wed, Jun 08, 2005 at 08:39:21PM -0400, [EMAIL PROTECTED] wrote:
I was wondering if someone familiar
with survival analysis can help me with
the following.
I would like to fit a Weibull curve,
that may be dependent on a covariate,
my dataframe labdata that has the
fields cov, time,
He,
I used the coxph function, with four covariates.
Let's say something like that
model.1 - coxph(Surv(Time,Event)~X1+X2+X3+X4,data=DATA)
So I obtain the 4 coefficients B1,B2,B3,B4 such that
h(t) = h0(t) exp(B1*X1+ B2*X2 + B3*X3 + B4*X4).
When I use the function on the same data
Hello,
On my barplot, I have on the axis y, the names of the rows of my matrix r_mat!
is it possible to change the color of these names on my barplot?
barplot2(t(r_mat),beside=TRUE,horiz=TRUE,plot.grid=TRUE,xlab=R,font.lab=4,las=2,xlim=c(0,1))
thanks,
Sabine
Le 09.06.2005 10:48, Navarre Sabine a écrit :
Hello,
On my barplot, I have on the axis y, the names of the rows of my matrix r_mat!
is it possible to change the color of these names on my barplot?
barplot2(t(r_mat),beside=TRUE,horiz=TRUE,plot.grid=TRUE,xlab=R,font.lab=4,las=2,xlim=c(0,1))
It is a very nice touch that optim() offers SANN (simulated annealing)
as a random search algorithm.
The R community already has genoud - an implementation of a genetic
algorithm for search.
Wouldn't it be neat if optim() would additionally offer method=GA
where it internally uses code from
On Thu, 9 Jun 2005 [EMAIL PROTECTED] wrote:
He,
I used the coxph function, with four covariates.
Let's say something like that
model.1 - coxph(Surv(Time,Event)~X1+X2+X3+X4,data=DATA)
So I obtain the 4 coefficients B1,B2,B3,B4 such that
h(t) = h0(t) exp(B1*X1+ B2*X2 + B3*X3 + B4*X4).
When
Dear @ll Friends,
After trying to install some libraries such as fSeries, I receive this error:
Fehler: unable to create temp directory
'C:/Programme/R/rw2010pat/library\file12763'
Also after trying to install other packages ,I receive similar errors, for
example, installation of
hello,
to use the function plot3d, i should use the package R.basic!
plot3d {R.basic}
If people know exactly a site to load this package, please give me the URL!
Thanks
Sabine
-
[[alternative HTML version deleted]]
http://www.maths.lth.se/help/R/R.classes/
in section 2.2 Manual download:
Download the Windows binary distribution (*.zip) or the non-Windows source
distribution (*.tar.gz)
Just Click on It
Grüße
Andreas
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of
Le 09.06.2005 11:48, Navarre Sabine a écrit :
hello,
to use the function plot3d, i should use the package R.basic!
plot3d {R.basic}
If people know exactly a site to load this package, please give me the URL!
Thanks
Sabine
Hello Sabine,
Do you plan to ask that same question everyday.
On Thu, 9 Jun 2005 11:48:32 +0200 (CEST) Navarre Sabine wrote:
hello,
to use the function plot3d, i should use the package R.basic!
plot3d {R.basic}
Better question than your previous two...
If people know exactly a site to load this package, please give me the
URL!
...but not
Quoting Prof Brian Ripley [EMAIL PROTECTED]:
On Thu, 9 Jun 2005 [EMAIL PROTECTED] wrote:
He,
I used the coxph function, with four covariates.
Let's say something like that
model.1 - coxph(Surv(Time,Event)~X1+X2+X3+X4,data=DATA)
So I obtain the 4 coefficients B1,B2,B3,B4 such that
h(t) =
Navarre Sabine wrote:
hello,
to use the function plot3d, i should use the package R.basic!
plot3d {R.basic}
I have just googled for
plot3d R.basic
and the first hit pointed me to the correct URL
Uwe Ligges
If people know exactly a site to load this package, please give me the
I do not see R.basic on CRAN, but digging a bit turns up Henrik's R.classes
bundle, which you can get by
install.packages(R.classes, contriburl=http://www.maths.lth.se/help/R;)
After that, just do
library(R.basic)
example(plot3d)
Andy
From: Navarre Sabine
hello,
to use the function
It is a file sustem problem on your computer (Windows, unstated).
Either you do not have permission to write in that directory or you have
some rogue files there. If you do have permission, go in and delete all
folders starting with `file', then reboot Windows and try again.
If you do not
On Thu, 9 Jun 2005 [EMAIL PROTECTED] wrote:
Quoting Prof Brian Ripley [EMAIL PROTECTED]:
On Thu, 9 Jun 2005 [EMAIL PROTECTED] wrote:
He,
I used the coxph function, with four covariates.
Let's say something like that
model.1 - coxph(Surv(Time,Event)~X1+X2+X3+X4,data=DATA)
So I obtain
Sabine,
It helps us to help you if you tell us what you want to do, preferably
with a code example, and what system/version of R you have: type version.
On my PC I get:
version
_
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major2
Now you made me curious!
Installed the package myself as well!
Maybe you can send an example graph to Romain, to be included in the
graph gallery. Then we can all see what plot3d does!
Thanks,
Sander.
Sander Oom wrote:
Sabine,
It helps us to help you if you tell us what you want to do,
Hi
I have recently been experimenting with the use of kriging, primarily
through Goldensoftware's Surfer. Having had a look around at the
packages available for R I note a couple that are available so have
initially added fields.
What I am interested in doing is batch processing ~300 datasets
Dear Ronaldo,
As you discovered, there are no predict methods in the sem package, either
for tsls or sem objects. You might have a look at the systemfit package,
which does support predict().
More generally, you should be able to substitute out identities and base
predictions on the reduced form
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Mike J Smith
Sent: 09 June 2005 09:58
To: r-help@stat.math.ethz.ch
Subject: [R] krig.image help
Hi
I have recently been experimenting with the use of kriging, primarily
through Goldensoftware's
hi,
My code:
data-matrix(data=c(0.425,0.5,0.75,0.125,0.25,0.475,0.375,0.25,0.625,0.5,0.1,0.125,0,0.25,0.25),nrow=3,ncol=5,byrow=TRUE,
dimnames=list(c(Good,Medium,Bad),c(Content,Logistic,Trainer,Supply,User
contribution)))
plot(corresp(data,nf=2),xlim=c(-1,1),ylim=c(-1,1));
The plot is
On 6/9/05, Navarre Sabine [EMAIL PROTECTED] wrote:
hi,
My code:
data-matrix(data=c(0.425,0.5,0.75,0.125,0.25,0.475,0.375,0.25,0.625,0.5,0.1,0.125,0,0.25,0.25),nrow=3,ncol=5,byrow=TRUE,
dimnames=list(c(Good,Medium,Bad),c(Content,Logistic,Trainer,Supply,User
contribution)))
Hi there,
I am trying to export a regression output to Latex. I am using the xtable
function in the xtable library. Doing
myfit - lm(myformula, mydata)
print.xtable(xtable(myfit), file=myfile)
only returns the estimated coefficients and the correspondent standard erros,
t-statiscs and
On Thu, 2005-06-09 at 14:40 +0200, Navarre Sabine wrote:
hi,
My code:
data-matrix(data=c(0.425,0.5,0.75,0.125,0.25,0.475,0.375,0.25,0.625,0.5,0.1,0.125,0,0.25,0.25),nrow=3,ncol=5,byrow=TRUE,
dimnames=list(c(Good,Medium,Bad),c(Content,Logistic,Trainer,Supply,User
contribution)))
On Thu, 2005-06-09 at 14:40 +0200, Navarre Sabine wrote:
My code:
data-matrix(data=c(0.425,0.5,0.75,0.125,0.25,0.475,0.375,0.25,0.625,0.5,0.1,0.125,0,0.25,0.25),nrow=3,ncol=5,byrow=TRUE,
dimnames=list(c(Good,Medium,Bad),c(Content,Logistic,Trainer,Supply,User
contribution)))
Hi Philippe and R community,
I am trying to use some functions from the svDialogs package but I get some
werid errors I do not understand:
library(svDialogs)
m_list - as.list(1:10)
guiDlgDoubleList(m_list, m_list)
Error in guiDlgDoubleList(m_list, m_list) :
Not yet implemented!
On Thu, 9 Jun 2005, Watalu, Y. (aka Wataru) wrote:
Hi,
I'm also wondering which expression the survreg() uses
for Weibull regression. Referring to help(survreg) and
help(survreg.distributions), I guess survreg() fits the
following model.
survreg() uses a different parametrization, say
F(x,
On Thu, 9 Jun 2005 [EMAIL PROTECTED] wrote:
I only have an other question:
since h(t) = h0(t) exp(B1*X1+ B2*X2 + B3*X3 + B4*X4)
represent the hazard at time t.
In a linear prediction,
what Value = B1*(X1-mean(X1)) + B2*(X2-mean(X2)) +
represent?
coxph() parametrizes the model so
Hello there:
Dear Dr. Murdoch,
I'm a statistician at Princess Margaret Hospital. Could you please help
me with the bootstrapping?
Can boot function (along with boot.ci) handle multivariate statistics?
My intention is that I want to group the outcome data and evaluate the
proportion of
Here is a link to see what's possible.
http://www.combinatorics.org/Surveys/ds5/VennSymmEJC.html
Venn did 4 sets with ellipses (see wikipedia.)
(There was an article in the last year (I think) in one of my math
journals that presented someone who made these complex Venn diagrams as
artwork and
Hi,
I'm trying to use transform my mercator locations into utm but I'm doing
something wrong because only x is transformed ... see
xy.sp
SpatialPoints:
loni lati
[1,] -8.85 38.16
[2,] -9.19 37.99
[3,] -9.11 37.97
[4,] -9.06 38.15
[5,] -9.03 37.87
[6,] -9.14 37.81
[7,] -9.09
data from multilevel units,first sample the class ,and then the student in
calss.following is the 2-level model. and the level-1 model deals with the
student,and the level-2 model deals with the class level the students belong to.
Level-1 Model
Y = B0 + B1*(ZLEAD) + B2*(ZBUL) +
Dear R users:
I am trying to add a label on axis(4) but I don't find
any reference in R archives. Could you please help me
out?
Thanks in advance.
The followings are the data set and what I have done:
Sample data:
time y1 y2
2 9 24.396
4 7.667 19.082
10 4.667 73.984
See ?mtext.
--sundar
Hai Lin wrote:
Dear R users:
I am trying to add a label on axis(4) but I don't find
any reference in R archives. Could you please help me
out?
Thanks in advance.
The followings are the data set and what I have done:
Sample data:
time y1 y2
2 9 24.396
4
Dear R-Users,
I have two questions:
a)
in a directory there are 3 files:
[1] Data.~csvKopie von Data.~csv VorlageTradefile.csv
The command dir( fold, pattern = \.csv ) gives back *all* the 3 files
With dir( fold, pattern = \\.csv ) I get back only VorlageTradefile.csv.
I don't
On Thu, 2005-06-09 at 09:26 -0700, Hai Lin wrote:
Dear R users:
I am trying to add a label on axis(4) but I don't find
any reference in R archives. Could you please help me
out?
Thanks in advance.
The followings are the data set and what I have done:
Sample data:
time y1 y2
2
Hello,
I've written a function that plots a few functions in a diagram.
The xlim and or ylim is not always the same, and set automatically by R.
A legend is part of this object.
Now the problem is: where to put the legend? Me would help a function that
returns the limits and scaling of the axis.
Hi,
The command dir( fold, pattern = \.csv ) gives back *all* the 3 files
With dir( fold, pattern = \\.csv ) I get back only VorlageTradefile.csv.
I don't understand this behaviour, IMHO the regex expression \.csv
becomes the string .csv and \\.csv becomes \.csv. So the first
string should
The documentation for lme includes the following:
lme(fixed, data, random, correlation, weights, subset, method,
na.action, control, contrasts = NULL)
...
Arguments:
...
random: optionally, any of the following: (i) a one-sided formula
Carsten Steinhoff wrote:
Hello,
I've written a function that plots a few functions in a diagram.
The xlim and or ylim is not always the same, and set automatically by R.
A legend is part of this object.
Now the problem is: where to put the legend? Me would help a function that
returns the
On Thu, 2005-06-09 at 19:05 +0200, Carsten Steinhoff wrote:
Hello,
I've written a function that plots a few functions in a diagram.
The xlim and or ylim is not always the same, and set automatically by R.
A legend is part of this object.
Now the problem is: where to put the legend? Me
Dimitri Joe wrote:
Hi there,
I am trying to export a regression output to Latex. I am using the xtable
function in the xtable library. Doing
You mean the xtable *package*.
myfit - lm(myformula, mydata)
print.xtable(xtable(myfit), file=myfile)
only returns the estimated coefficients and
On 6/9/05, Hans-Peter [EMAIL PROTECTED] wrote:
Dear R-Users,
I have two questions:
a)
in a directory there are 3 files:
[1] Data.~csvKopie von Data.~csv VorlageTradefile.csv
The command dir( fold, pattern = \.csv ) gives back *all* the 3 files
With dir( fold, pattern =
(2) Have you seen Pinheiro and Bates (2000) Mixed-Effects Models in S
and S-PLUS (Springer)? This is the definitive documentation on
library(nlme). I've learned a lot from this book, and it might help you
even if you use lme4.
i know this book and also want to read it,but the problem
Dear All,
I am trying to specify the following fixed effects model for lme:
y ~ constant1 - beta1*(x - beta2)
where y is the response, x is the independent variable, and the
operators above are real arithmetic operations of addition, subtraction,
and multiplication. I realize that this
Ravi Varadhan wrote:
You can use locator(1) as an argument in the legend function to manually
position the legend, i.e. after your plot is plotted, the system will wait
for you to indicate (by clicking the mouse) where the legend should be
placed.
legend(locator(1), ...)
Hope this helps,
Ravi.
On 6/9/05, Eric Hack [EMAIL PROTECTED] wrote:
Dear All,
I am trying to specify the following fixed effects model for lme:
If you have a linear fixed-effects model you should use lm, not lme.
y ~ constant1 - beta1*(x - beta2)
where y is the response, x is the independent variable,
I'm seeing some inconsistent behavior when re-assigning values in a data
frame. The first assignment turns all of the 0s in my data frame to 2s,
the second fails to do so.
df1 - data.frame(a = c(NA, 0, 3, 4))
df2 - data.frame(a = c(NA, 0, 0, 4))
df1[df1 == 0] - 2 ## Works
df2[df2 == 0] - 2
On Thu, 9 Jun 2005, McGehee, Robert wrote:
I'm seeing some inconsistent behavior when re-assigning values in a data
frame. The first assignment turns all of the 0s in my data frame to 2s,
the second fails to do so.
df1 - data.frame(a = c(NA, 0, 3, 4))
df2 - data.frame(a = c(NA, 0, 0, 4))
On 6/9/2005 1:05 PM, Carsten Steinhoff wrote:
Hello,
I've written a function that plots a few functions in a diagram.
The xlim and or ylim is not always the same, and set automatically by R.
A legend is part of this object.
Now the problem is: where to put the legend? Me would help a function
Thanks for the response. It is actually a repeated measures study, I
just mention the fixed effects specification because I think I know the
random effect specification, i.e.:
Random = ~ 1|subject
And thanks for the tip about the nonlinear model and the S-plus list. I
will check out nlme and
Hi all,
I originally posted this to the bioconductor group, but maybe it's better
suited to the r-help...
I'm using som() to partition samples of gene expression data into clusters.
The point is to classify control vs. experimental cases (sample clustering).
The original matrix was 22283 x
On Wed, 2005-06-08 at 14:15 +0200, Thomas Petzoldt wrote:
...
Hello Ross,
I see that your question was related to S4, but I just noticed a
solution based on the R.oo package so I thought I would add a solution
based on the proto package too. We had similar problems several times
ago and (to
On Wed, 2005-06-08 at 09:13 +0200, Uwe Ligges wrote:
[extensive deletions. Discussion concerned my desire to have a function change
the value
of an object in a way that had effects outside of the function, without
returning the object.]
You have to think about scoping rules and it
will be
To add on Mark's comments, another neat option to intreractivelly place the
label in the desired location is to use locator(), i.e:
text( locator(1),This is Y Axis 2, srt = 270, xpd = TRUE)
Once you execute this command, left click on the desired position for the
label in the plot and enjoy
Ross Boylan wrote:
I believe your example assumes that foo is updating the outer a by
cheating and directly modifying enclosing environments. (I figure it
also needs to get the name of its actual argument to do this, which
would also involve slightly dirty tricks.) This does seem to be the
On 6/9/2005 3:04 PM, Ross Boylan wrote:
On Wed, 2005-06-08 at 09:13 +0200, Uwe Ligges wrote:
[extensive deletions. Discussion concerned my desire to have a function change
the value
of an object in a way that had effects outside of the function, without
returning the object.]
You have to
Ross Boylan wrote:
The paper also notes that call-by-value vs call-by-reference, which is
the root of the slot update problem, is really orthogonal to the
FOOP/COOP distinction. It's easy to imagine FOOP with
call-by-reference. R.oo uses references.
Yes, it can be found on
On Thu, 9 Jun 2005, Thomas Lumley wrote:
On Thu, 9 Jun 2005, McGehee, Robert wrote:
I'm seeing some inconsistent behavior when re-assigning values in a data
frame. The first assignment turns all of the 0s in my data frame to 2s,
the second fails to do so.
But they differ in several ways, so
Dear Gorjanc,
It is a good idea to make an CTV for genetics. I made an informal summary
to Hum Genomics a short while ago and quoted below.
... some packages for genetic data analysis, which fall into several
categories: basic data manipulation (genetics) phylogenetic analysis
(PHYLOGR, ape),
On Thu, 9 Jun 2005, Ernesto Jardim wrote:
Hi,
I'm trying to use transform my mercator locations into utm but I'm doing
something wrong because only x is transformed ... see
Two points: I think you are getting the transform you specify, from
+proj=merc to +proj=utm with a different x
On 6/8/05, Mike R [EMAIL PROTECTED] wrote:
snip
From the command line, I'd like to be able to start an
interactive R session and at the same time, be able
to specify on the command line a project-specific or
task-specific script (R code) that is to be executed
at the beginning of the session.
test of thread behaviour
this is the initial post
__
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
this is the first reply to the initial thread
subject line and header unchanged
__
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
this is the second reply to the initial thread
subject line has been changed
__
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
this is the third reply to the initial post
subject line has been changed
__
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
this is the fourth reply to the initial post
subject line edited ( Re was removed )
__
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
this is fifth and final (?) reply to the initial post
if this post gets added to the original thread,
then i've found no way to initiate a new thread by
clicking the reply button
had no idea that this behaviour occurred. wonder
how many archives (other email lists) are corrupted
because of
okay, one last post.
gmail categorizes email by threads
gmail's behaviour is different than this R-archive:
https://www.stat.math.ethz.ch/pipermail/r-help/2005-June
in fact, gmail has categorized the posts as i had intended
but then gmail, i think, calls them conversations, not threads
now
Thanks for your help, that worked.
David
On 6/8/05, Martin Maechler [EMAIL PROTECTED] wrote:
MM == Martin Maechler [EMAIL PROTECTED]
on Wed, 8 Jun 2005 18:57:55 +0200 writes:
David == David Finlayson [EMAIL PROTECTED]
on Wed, 8 Jun 2005 09:24:54 -0700 writes:
David
I have some disconnected boundary data from a finite element ocean model
and I'd like to make a plot.
Maptools looks promising, but since my data is not in a shapefile or a
map, I'm unclear on what the best way to approach the problem.
geom[1:10,]
lon lat depth
1 -75.42481
this is the seventh reply to the initial post
subject line is unaltered
I added a ficticious (i think) email to the To: field
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hi,
Try allt.som$visual
Edgar Acuna
UPRM
On Thu, 9 Jun 2005, Ken Termiso wrote:
Hi all,
I originally posted this to the bioconductor group, but maybe it's better
suited to the r-help...
I'm using som() to partition samples of gene expression data into clusters.
The point is to classify
I am doing an analysis that requires me to calculate correlations for a matrix
of 15,000 rows x 50 columns. For each row I want to calculate the correlation
to all other rows and then for each row, find the n (say 10) most correlated
rows. If read in the 15,000 x 50 data from file and pass it
Hi
I am using R for logistic regression and finding it very useful.
However, I wondered if anyone could point me to any course or notes on
this subject using R.
All help most welcome.
Stephen
--
Internal Virus Database is out-of-date.
Checked by AVG Anti-Virus.
[[alternative
Hi,
I am thinking to use one of them but not sure which one is better. I
think Rpy cannot call python from R while the PRPython can in
two-directional calling. Am I right?
thanks,
--
Weiwei Shi, Ph.D
Did you always know?
No, I did not. But I believed...
---Matrix III
On Tue, 7 Jun 2005, David Forrest wrote:
...
http://fawn.unibw-hamburg.de/cgi-bin/Rwiki.pl?UsePictures
This page demonstrates the use of pictures and the use of a public image
server for hosting graphics files.
Also, there's a page:
Hi All,
Is it possible to control the placement R's device windows,
such as:
x()
Or more generally, to passing X-resources?
x( geometry=500x200+0+0 )
Has anyone by chance worked out the possibilities for controlling
R windows in .Xresources or .Xdefaults ?
Thanks in advance,
Mike
Hi,
Is it possible to redirect the staff that normally goes to the R console
window into a file. sink() does this for sdterr and stdout. But I need
something that redirects everything that appear in the console window
(including the top-level commands).
I want to achieve the same effect as
Will work, but as you mentioned there ought to be an easier way :-)
-Original Message-
From: Mike R [mailto:[EMAIL PROTECTED]
Sent: Thursday, June 09, 2005 10:29 PM
To: Vadim Ogranovich
Subject: Re: [R] Redirect console to file
On 6/9/05, Mike R [EMAIL PROTECTED] wrote:
On
Hello Marco,
For the first error, the message is clear: not implemented yet!.
Several dialog boxes are not done yet, but the functions already exist,
mainly as placeholders for future development.
Regarding the second, there was a bug in the function (corrected in
SciViews 0.8-6 that I have
Dear all
I am tryiing to fit the following lme with an ARMA correlation structure:
test - lme(fixed=fev1f~year, random=~1|id2, data=pheno2,
correlation=corARMA(value=0.2, form=~year|id2), na.action=na.omit)
But I get the following error message:
Error in getGroupsFormula.default(correlation,
83 matches
Mail list logo
