Hi everyone,
I'm looking for any examples of useful frailty models, in particular any
situation in which a cox proportional hazards model with frailty outperforms a
regular cox proportional hazards model with respect to prediction of the time
to event (or the X-year risk of an event). I have
Hi,
I am trying to creat a matrix with diagnal blocks, say, I have a matrix X
of any dimension (nxm) ,and would like to have:
X
Hello R Experts,
I want to aggregate parameters by week. But our production week ends Friday
night instead of Sunday Night which is the default value in R.
In order to solve the problem I want to substract two days from the current
data and than use the R function
Hi,
Can some one please explain me with an example what is
Symbolic Number Coding in R.
Thanks and Regards,
Pratap
-
[[alternative HTML version deleted]]
Short story: January 2006 I did some analysis in R-2.2.1 using nls. Repeating
the exercise in R-2.3.1 yesterday produced somewhat different results.
After some debugging I found that either nls is the problem or that mine
understanding of environments or scoping rules is lacking something.
On Thu, 21 Sep 2006, Thorsten Muehge wrote:
Hello R Experts,
I want to aggregate parameters by week. But our production week ends Friday
night instead of Sunday Night which is the default value in R.
The default in ISO8601, not just in R, but that is %W, not %U as used
below.
In order to
GS == Gavin Simpson [EMAIL PROTECTED]
on Thu, 21 Sep 2006 00:08:17 +0100 writes:
GS On Wed, 2006-09-20 at 18:56 -0400, Charles Annis, P.E. wrote:
Recommending a good book on statistics is like recommending a good book
on
sports: Which sports?
A good book for
If I remember well, there should be a package including the function
bdiag() making the job..but I am not able to remember its name..
However a quick search via RSiteSearch(bdiag) yields
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/40393.html
Hope this helps you,
vito
Tong Wang wrote:
Martin Maechler [EMAIL PROTECTED] writes:
GS I couldn't agree more. A number of my colleagues have bought Peter
GS Dalgaard's book to a) learn some R and b) learn some statistics. They
GS have found it very useful indeed.
Yes!
I'm pretty sure it has been the first book of its
Hello
You need adiag(), which
is now part of the magic package:
library(magic)
adiag(matrix(1:6,2,3),matrix(1:9,3,3))
[,1] [,2] [,3] [,4] [,5] [,6]
[1,]135000
[2,]246000
[3,]000147
[4,]00025
Prof Ripley,
You are absolutely correct, this code will not work at all - for
starters, M isn't correctly initialized, etc. I edited my code in
tinn-r and never ran it before posting... my apologies, I always seem
to be too quick off the mark to reply - despite the 4*runif(1)
suggestion in the
Hi
On 20 Sep 2006 at 14:19, Dave Evens wrote:
Date sent: Wed, 20 Sep 2006 14:19:04 -0700 (PDT)
From: Dave Evens [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject:[R] Spliting a huge vector
Dear R users,
I have a huge
Thanks for all of your information.
I will look into all of them
Rainer
Gavin Simpson wrote:
On Wed, 2006-09-20 at 17:03 +0200, Rainer M Krug wrote:
Hi
I heard so much about Emacs and ESS that I decided to try it out - but I
am stuck at the beginning.
Is there anywhere a beginners
On Wednesday 20 September 2006 17:16, Marc Schwartz (via MN) wrote:
On Wed, 2006-09-20 at 17:03 +0200, Rainer M Krug wrote:
Hi
I heard so much about Emacs and ESS that I decided to try it out - but I
am stuck at the beginning.
Is there anywhere a beginners manual for Emacs ESS to
On R 2.3.1 I was surprised to see the following results from
summary.default():
R x - scan()
1: 12148 10426 10912 9116 13226 11663 11781 10680 8457 10788 12605
10591 11040
14: 10815 12962 11644 10047 10478 10108 12353 11778 11092 11673 8758
11145 11495
27:
Read 26 items
R summary(x)
On Wed, 20 Sep 2006, Jeffrey Horner wrote:
Gamal Azim wrote:
Is it possible to remotely save an R session then
terminate R? Of course the destructive task after
'then' is rather straightforward by itself.
Yes. Sending the process a SIGUSR1 signal:
$ kill -USR1 pid
will save the global
On Wednesday 20 September 2006 22:21, Iuri Gavronski wrote:
I would like to buy a basic statistics book (experimental design,
sampling, ANOVA, regression, etc.) with examples in R. Or download it
in PDF or html format.
I went to the CRAN contributed documentation, but there were only R
Hi all,
I am very fond of Using R for Introductory Statistics
by John Verzani, 2005, Chapman Hall.
Regards
D. Trenkler
--
Dietrich Trenkler c/o Universitaet Osnabrueck
Rolandstr. 8; D-49069 Osnabrueck, Germany
email: [EMAIL PROTECTED]
__
Hi
I generate a new dataframe by doing:
npl.agg - aggregate(npl$DensPlants, list(year=npl$year, sim=npl$sim),
mean, na.rm=TRUE )
Now I want to plot it by using
coplot(npl.agg$x ~ npl.agg$year | npl.agg$sim, type=l)
but, as npl.agg$year is seen as a factor, the order of the points on the
That's an FAQ:
7.10 How do I convert factors to numeric?
It may happen that when reading numeric data into R (usually, when
reading in a file), they come in as factors. If f is such a factor
object, you can use
as.numeric(as.character(f))
to get the numbers back. More efficient, but
Dear useRs,
While exploring new R packages, I have found the Rattle.
This screenshot http://rattle.togaware.com/rattle-correlation.png
is very interesting
(others are in http://rattle.togaware.com/rattle-screenshots.html ).
Which function was used to produce this plot?
Is such plotting of the
Thanks a lot David.
Just a suggestion - it might help if this is added to the help of factor()
Rainer
David Barron wrote:
That's an FAQ:
7.10 How do I convert factors to numeric?
It may happen that when reading numeric data into R (usually, when
reading in a file), they come in as
Jeffrey Horner wrote:
$ kill -USR1 pid
will save the global environment in the file .RData, but you'll need to
remember the current working directory of the process to find it.
Remember? With a computer, you never need to remember!
$ ls -l /proc/$pid/cwd
is a symlink to the current
Hi
I am using R to fit a survival function to my data
(with a weibull distribution).
Data: Survival of individuals in relation to 4
treatments ('a','b','c','g')
syntax:
survreg(Surv(date2)~males2, dist='weibull')
But I have some problems interpreting the outcome and
getting the
Dear all
I am writing again with a question I posted a few weeks ago (to no avail). I
have a problem with SciViews for R. It's probably a slightly stupid question
but I cannot find a solution to a very elementary problem. I am using SciViews
0.8.9 on with R 2.3.1pat on a Windows XP Home
Dear Sir/Madam,
I am encountering one of those alien computer momements one finds every
so often in life. See the sequence below:
fish3.fis -read.csv(emperor2.csv, check.names = TRUE, strip.white =
TRUE)
colnames(fish3.fis)
[1] Month YearFishingArea
SumOfTotalCatch CPUE
Hi
is there any other editor for Linux, apart from emacs, which has the
completion feature (i.e. I enter libTAB and it offers me among other
options library)
Rainer
--
Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation
Biology (UCT)
Department of Conservation Ecology and Entomology
On Thu, 21 Sep 2006, Al-Kharusi, L. wrote:
Dear Sir/Madam,
I am encountering one of those alien computer momements one finds every
so often in life. See the sequence below:
fish3.fis -read.csv(emperor2.csv, check.names = TRUE, strip.white =
TRUE)
colnames(fish3.fis)
[1] Month
On Thu, 2006-09-21 at 13:07 +0100, Al-Kharusi, L. wrote:
Dear Sir/Madam,
I am encountering one of those alien computer momements one finds every
so often in life. See the sequence below:
fish3.fis -read.csv(emperor2.csv, check.names = TRUE, strip.white =
TRUE)
colnames(fish3.fis)
JGR support completion. It can runs under Linux,though JGR is more
than a editor.
On 9/21/06, Rainer M Krug [EMAIL PROTECTED] wrote:
Hi
is there any other editor for Linux, apart from emacs, which has the
completion feature (i.e. I enter libTAB and it offers me among other
options library)
On Thu, 2006-09-21 at 13:58 +0200, [EMAIL PROTECTED] wrote:
Dear all
I am writing again with a question I posted a few weeks ago (to no
avail).
STG
Perhaps the main reason you received no replies is that you failed to
read to the posting guide and act appropriately. Sciviews is *not* part
hi,
I am a system admin who has just set up R-2.3.1 in suse 9.1 (opteron)
for a previous windows user. Please forgive any obvious stupid questions :)
We are struggling with a number of things.
1. is there a gui interface ? How do we use it ? R --gui just gets us
to the command line.
2.
MJ Crawley Statistics: An Introduction Using R, received a good review in
Journal of the Royal Statistical Society, Vol 169.
Dan Bebber
Dr. Daniel P. Bebber
Department of Plant Sciences
University of Oxford
South Parks Road
Oxford OX1 3RB
UK
[[alternative HTML version deleted]]
Julie Harold [EMAIL PROTECTED] writes:
hi,
I am a system admin who has just set up R-2.3.1 in suse 9.1 (opteron)
for a previous windows user. Please forgive any obvious stupid questions :)
We are struggling with a number of things.
1. is there a gui interface ? How do we use it ? R
Hi
Thanks for the tip - sounds interesting.
I installed JGR, but I can't get it to work.
I installed it as mentioned on via install.packages() and if I try to
load it, I get the following:
library(JGR)
Loading required package: rJava
Loading required package: JavaGD
Loading required package:
Hi
I am not familiar with linux version but
On 21 Sep 2006 at 13:56, Julie Harold wrote:
Date sent: Thu, 21 Sep 2006 13:56:06 +0100
From: Julie Harold [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject:[R] newie help needed
First of all, since this is my first posting, I would like to thank
anybody who works/has worked on R, and made it such a beautiful open
source package!
My question is: how do I create a new Access database with R? I need a
channel before I can do anything, but if the mdb-file doesn't exist, I
Does anyone know how to get a meaningful estimate of the covariance
matrix for the first canonical variate coefficients for the right-hand
(X) side of a canonical regression? cancor returns coefficients but not
precisions. This would have to be subject to some constraint because of
Julie Harold j.m.harold at uea.ac.uk writes:
1. is there a gui interface ? How do we use it ? R --gui just gets us
to the command line.
Look at SciViews. Not a 100% menu driven GUI, but can be useful to new users.
3. can we define default directories - eg for data files ?
I think you
On 9/21/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
(i) the command area at the bottom cannot be found
(ii) the regular R window cannot be minimized/maximized anymore.
I cannot duplicate these problems with R 2.3.1, SciViews 0.8-9, Win XP(current).
jab
--
John Bollinger, CFA, CMT
Vladimir Eremeev [EMAIL PROTECTED] wrote in message
news:[EMAIL PROTECTED]
Dear useRs,
While exploring new R packages, I have found the Rattle.
This screenshot http://rattle.togaware.com/rattle-correlation.png
is very interesting
Which function was used to produce this plot?
Because of the replacement of the current NAS (Networked Disk Server),
there will be a mail server timeout
from 18:00 till about 20:00 Swiss time, i.e. UTC+2
This will also effect cran.CH.r-project.org,
our ftp server and http://ESS.r-project.org/
but *not* the subversion server
Suppose I have a square matrix P
P - matrix(c(.3,.7, .7, .3), ncol=2)
I know that
P * P
Returns the element by element product, whereas
P%*%P
Returns the matrix product.
Now, P^2 also returns the element by element product. But, is there a
slick way to write
P %*% P %*% P
Obviously,
Please stop posting your homework problems to R-help!
The posting guide quite explictly does not allow it.
M.Maechler, ETH Zurich
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On 9/21/2006 10:40 AM, Doran, Harold wrote:
Suppose I have a square matrix P
P - matrix(c(.3,.7, .7, .3), ncol=2)
I know that
P * P
Returns the element by element product, whereas
P%*%P
Returns the matrix product.
Now, P^2 also returns the element by element product. But,
On 9/21/06, Earl F. Glynn [EMAIL PROTECTED] wrote:
Which function was used to produce this plot?
library(ellipse)
?plotcorr
Look at the third example for the color plot.
That's pretty nice. I have been using symnum(). Are there any other
neat visualisations for correlation matrices?
Dear all, I have a vector like this:
z = rep(c(M,F),c(50,60))
How can I get the following frequency table:
Sex counts %
F60 54.5
M 50 45.5
I try:
DD- function(data,...)
{
n - nobs(data)
out - c(Frequency
On Thu, 21 Sep 2006, Anaid Diaz wrote:
Hi
I am using R to fit a survival function to my data
(with a weibull distribution).
Data: Survival of individuals in relation to 4
treatments ('a','b','c','g')
syntax:
survreg(Surv(date2)~males2, dist='weibull')
But I have some problems
BBands a écrit :
That's pretty nice. I have been using symnum(). Are there any other
neat visualisations for correlation matrices?
Although very simple, I find often image() sufficient.
Plus some ordering if needed.
__
R-help@stat.math.ethz.ch
Quoting Duncan Murdoch [EMAIL PROTECTED]:
On 9/21/2006 10:40 AM, Doran, Harold wrote:
Suppose I have a square matrix P
P - matrix(c(.3,.7, .7, .3), ncol=2)
I know that
P * P
Returns the element by element product, whereas
P%*%P
Returns the matrix product.
Now, P^2 also returns
At 12:03 20/09/2006, Thorsten Muehge wrote:
Hello R Experts,
how can I incorporate a variable in a data frame definition?
Thorsten, you have already had a reply about this, but I wonder
whether you _really_ want to do this. You do not say what the
scientific question is you are trying to
Dear R-friends
I fitted a regression using GLM and plotted fitted values agaist X axis.
Now I need plot upper and lower confidence intervals (95%) for this fit.
Acoording to ?lm help and can use something like :
[ predict.lm(mod, se.fit=T, interval = confidence, level=.95) ]
to
There is also mexp in the Matrix package and MatrixExp in the msm package.
Andreas
- Original Message -
From: Dimitrios Rizopoulos [EMAIL PROTECTED]
To: Duncan Murdoch [EMAIL PROTECTED]
Cc: Doran, Harold [EMAIL PROTECTED]; r-help@stat.math.ethz.ch
Sent: Thursday, September 21, 2006
On 9/21/06, Dimitrios Rizopoulos [EMAIL PROTECTED] wrote:
Quoting Duncan Murdoch [EMAIL PROTECTED]:
On 9/21/2006 10:40 AM, Doran, Harold wrote:
Suppose I have a square matrix P
P - matrix(c(.3,.7, .7, .3), ncol=2)
I know that
P * P
Returns the element by element product,
delete a entire vector of a dataframe
Hello,
i want to delete a vector and tried rm (t.d$V712). This did not work,
message was, could not find variable. I thought the $ defines the vectro
in a dataframe, when I just type t.d$V712 the content of this vector
is displayed.
Greetings, Thomas
And look at ?prop.table too.
On 21/09/06, lamack lamack [EMAIL PROTECTED] wrote:
Dear all, I have a vector like this:
z = rep(c(M,F),c(50,60))
How can I get the following frequency table:
Sex counts %
F60 54.5
M 50 45.5
Hello,
I`m a newbie to R so maybe this question is boring, but I have a large
table with several empty missing values, which come out as NA. How can
i ignore them or replace them by another number?
Greetings, Thomas
__
R-help@stat.math.ethz.ch
Dear all,
I'm working on MAC OS X 10.4.7 using R version 2.3.1 (June 2006)
and platform powerpc-apple-darwin8.6.0. Since I installed the new
version of R I cannot install packages any more (but it worked fine
before) and get the following message every time. I can, however,
install
You might want to look at the CrossTable function in the gmodels
package (in the gregmisc bundle).
On 21/09/06, lamack lamack [EMAIL PROTECTED] wrote:
Dear all, I have a vector like this:
z = rep(c(M,F),c(50,60))
How can I get the following frequency table:
Sex counts %
F
Dear everybody!
I am trying to make up a histogram of marks of a recent test I have
corrified.
Say I have a vector points and a vektor breaks to segregate the
marks. Finally there is a vector marks.
hist(points, breaks=breaks) leaves me with two problems:
At first the width of the bars of the
The cairoDevice package has recently been released on CRAN. cairoDevice
is an R graphics device based on the cairo
(http://www.cairographics.org) vector graphics library. It is
distinguished from many other devices in its ability to render high
quality anti-aliased graphics. It is available on
It depends a bit on what function you are using. For example,
set.seed(1)
xx - c(NA, rnorm(10))
mean(xx)
[1] NA
mean(xx, na.rm=TRUE)
[1] 0.1322028
Is how you would use this to compute a mean.
Harold
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf
Hola!
maybe someone can see what is wrong in the following graphics code.
The data needed for the example (printed with dump) is at the end of the
message.
I make 3 pages with multipanel plot, and want each page to have a subtitle
explaining the common x-axis, with the command:
Quoting Douglas Bates [EMAIL PROTECTED]:
On 9/21/06, Dimitrios Rizopoulos [EMAIL PROTECTED] wrote:
Quoting Duncan Murdoch [EMAIL PROTECTED]:
On 9/21/2006 10:40 AM, Doran, Harold wrote:
Suppose I have a square matrix P
P - matrix(c(.3,.7, .7, .3), ncol=2)
I know that
P * P
There is a function remove.vars in the gdata package, or you can do this:
dt-data.frame(a=(1:5),b=(6:10),c=(11:15))
dt
a b c
1 1 6 11
2 2 7 12
3 3 8 13
4 4 9 14
5 5 10 15
dt[[a]] - NULL
dt
dt
b c
1 6 11
2 7 12
3 8 13
4 9 14
5 10 15
If you know it is the first column that
Most R functions have arguments of the form na.action or na.rm that
allow you to specify how you treat NA's. In general, it's not a good idea to
replace NA's with numbers.
See also ?na.omit, ?na.action.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of
Hi Marc,
thanks for the answer, but that wasn't the problem.
I set my LD_LIBRARY_PATH in my .cshrc.
I solve my problem by compiling the library manually. It created me a
file (ROracle.so) of ~14mb (instead of the one of ~150K) that I replaced
with the other one.
I guess R compilation did
On Thu, 2006-09-21 at 19:58 +0200, Thomas Preuth wrote:
Hello,
I`m a newbie to R so maybe this question is boring, but I have a large
table with several empty missing values, which come out as NA. How can
i ignore them or replace them by another number?
Greetings, Thomas
On Thu, 2006-09-21 at 20:01 +0200, Thomas Preuth wrote:
delete a entire vector of a dataframe
Hello,
i want to delete a vector and tried rm (t.d$V712). This did not work,
message was, could not find variable. I thought the $ defines the vectro
in a dataframe, when I just type t.d$V712
Hi,
Is there some function, which generates this kind of n x n -matrix
from a vector?
rhset
[1] 1792 256 13312 512 1024 2048 8192 4096
m=matrix(nrow=length(rhset),ncol=length(rhset))
for(i in 1:length(rhset))
+ {
+ m[,i]=rhset
+
I'm not sure why it doesn't work, perhaps something to do with how you
have set the mar and oma parameters? But why not use
mtext(side=1,text=aleles(bp),outer=TRUE,col=snow3)
On 21/09/06, Kjetil Halvorsen [EMAIL PROTECTED] wrote:
Hola!
maybe someone can see what is wrong in the following
kone said the following on 9/21/2006 2:30 PM:
Hi,
Is there some function, which generates this kind of n x n -matrix
from a vector?
rhset
[1] 1792 256 13312 512 1024 2048 8192 4096
m=matrix(nrow=length(rhset),ncol=length(rhset))
for(i in 1:length(rhset))
+
Are you trying to install a package from source?
Try adding type='source' when you call install.packages().
The help page for install.packages() indicates the need for this on a Mac.
You might also check and see if ~/libs exists, and you have write
privileges, but take care of the source vs.
Hi Don,
type=source solved the problem.
Thanks very much for your help,
best,
Snaebjorn
On Sep 21, 2006, at 3:59 PM, Don MacQueen wrote:
Are you trying to install a package from source?
Try adding type='source' when you call install.packages().
The help page for install.packages()
Robert W. Baer, Ph.D.
Associate Professor
Department of Physiology
A. T. Still University of Health Science
800 W. Jefferson St.
Kirksville, MO 63501-1497 USA
- Original Message -
z = rep(c(M,F),c(50,60))
How can I get the following frequency table:
Dear R users,
Sorry, I made a mistake in my specification of
problem, should be:
# my vector is this
a.vector - seq(2, by=5, length=1000)
# so my cut values are
cut.values - c(30, 50, 100, 109, 300, 601, 803, 1000)
# and x should have an extra ) at the end
x - rep(1:length(cut.values),
Hi,
I'm updating the LMGene package from Bioconductor. Writing R Extensions
suggests
that all source files (the ones in the R directory) have a .R ending, so I
added it to the (one) source file.
The next time I installed and ran R, R didn't understand any of the
functions.
I tried various things
I am updating the Bioconductor package, LMGene. Thus I am modifying someone
else's package, editing or writing new documentation, etc.
When I modify the man files in LMGene and install the library, it doesn't
change the 'help' that R gives you. The 'help' you actually get in R is the
same as
Gentlemen,
(I am using R 2.2.1 in a Windows environment.)
I apologize but I did not fully comprehend all of your
answer. I have a dataframe called data1. I run
several liner regression using the lm function similar
to:
reg - ( lm(lm(data1[,2] ~., data1[,2:4])) )
I see from generous
?vcov
-- 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: r user [mailto:[EMAIL PROTECTED]
Sent: Thursday, September 21, 2006 2:11
On 9/21/06, John Tillinghast [EMAIL PROTECTED] wrote:
Hi,
I'm updating the LMGene package from Bioconductor. Writing R Extensions
suggests
that all source files (the ones in the R directory) have a .R ending, so I
added it to the (one) source file.
The next time I installed and ran R, R
On 9/21/2006 5:01 PM, John Tillinghast wrote:
Hi,
I'm updating the LMGene package from Bioconductor. Writing R Extensions
suggests
that all source files (the ones in the R directory) have a .R ending, so I
added it to the (one) source file.
The next time I installed and ran R, R didn't
Hmm... sounds like you're on Windows and have the Explorer setup such
that it is hiding file extensions. You can try to use list.files()
from R to see if the files actually have file extension.
If this is your problem, open My Computer - Tools - Folder
Options... and select tab View. Make sure
Yes, this was exactly the problem: I was using the unzipped package, not the
source.
Now it works!
-- Forwarded message --
From: Deepayan Sarkar [EMAIL PROTECTED]
Date: Sep 21, 2006 2:58 PM
Subject: Re: [R] Adding .R to source file keeps R from reading it?
To: John Tillinghast
Hi all,
How to store recursive resutls from a function for each step without using
global operators -? Thanks ahead.
Xiaohui Chen
Dept. of Statistics
UBC, Canada
_
Dont waste time standing in linetry shopping online. Visit
85 matches
Mail list logo