Dear Ted,
Yes, that makes sense, and I hadn't thought of it -- I was thinking in
terms of a nonparametric estimate of the hazard function. Spencer Graves
makes a similar point. Andy Liaw was kind enough to point out to me that
the muhaz function is in the muhaz package. As it turns out, muhaz
Sorry, but there was an error in the link that I posted below.
It should have been
http://www.vsn-intl.com/genstat/downloads/datald.htm
---
Date: Tue, 4 Nov 2003 19:54:06 -0500
From: Gabor Grothendieck <[EMAIL PROTECTED]>
To: <[EMAIL PROTECTED]>, <[EMAIL PROTECTED]>
Cc: <[EMAIL PROTECTED]>
On Tue, 2003-11-04 at 17:15, Paul Sorenson wrote:
> Thanks to those who pointed me at the solutions to the legend
> overprinting the bars. I took the "easy" way of rescaling the y axis,
> picking the scaling factor for stacked bars is somewhat problematic
> but sufficient for my application.
>
>
On 05-Nov-03 John Fox wrote:
> Dear Monica,
>
> I'm not sure what the muhaz function is (it's not in the survival
> package), but regardless, unless I'm seriously mistaken, there's no
> information to estimate the hazard function if you haven't observed
> any events.
>
> I hope that this helps,
>
If you are willing to assume an exponential distribution, then you can
get a one-sided confidence bound on the exponential parameter. The
likelihood is the probability of what you observe. Write that down and
figure it out. I'm sorry I don't have time to say more right now, but
I've done it
Dear Monica,
I'm not sure what the muhaz function is (it's not in the survival package),
but regardless, unless I'm seriously mistaken, there's no information to
estimate the hazard function if you haven't observed any events.
I hope that this helps,
John
At 08:10 PM 11/4/2003 -0500, Monica L.
On Tue, 2003-11-04 at 13:12, Prof Brian Ripley wrote:
> Look into external pointers. That is how we have tackled this, e.g. in
> the ts package.
I got the following to work. Any comments? I indicated some areas of
uncertainty in the comments. I'm also unsure about the differences, if
any, betw
On Wed, 05 Nov 2003 09:22:34 +1000
Ross Darnell <[EMAIL PROTECTED]> wrote:
> I have fitted a logistic regression model
>
> > failed.lr2$call
> lrm(formula = failed ~ Age + task2 + Age:task2, data = time.long,
> na.action = na.omit)
Use Age*task2 and omit na.action as na.omit is the default
On Tue, 4 Nov 2003 16:19:44 -0800 , you wrote:
>Hello,
>
>I am running a program in r that calls a function, which calls another
>function, which calls another etc. These functions are of the form:
>
>example<- function(x,y,z)
>
>{x, y, and z are defined within curly braces like this}
>
>Here'
It would seem that R 1.8 doesn't like the version of libiconv that
Panther has installed. Can we compile this for Aqua ourselves or will
there be a binary for we who upgrade sooner rather than later?
[ludwig:~] gall% R
dyld: mkdir version mismatch for library:
/usr/local/lib/libiconv.2.dylib (c
Thank you for your responses.
Here's the workaround. In SPSS, I saved the file EXvideo.sav as
Excel. Still in SPSS, I read in this Excel file, and saved it
under a new name, EXvideo-SPSStoExceltoSPSS.sav, as an SPSS file.
This file I was able to read into R with no problem, using read.spss.
Th
see ?fit.contrast in library gregmisc.
Cheers,
Simon.
Simon Blomberg, PhD
Depression & Anxiety Consumer Research Unit
Centre for Mental Health Research
Australian National University
http://www.anu.edu.au/cmhr/
[EMAIL PROTECTED] +61 (2) 6125 3379
> -Original Message-
> From: Igor Royt
On Tue, 4 Nov 2003, Jacob Wegelin wrote:
>
> Is there any documentation on what kind of SPSS file can and cannot be
> read by read.spss? Alternatively, how can one modify or "clean" an SPSS
> file to make it readable by read.spss? What properties must a *.sav file
> before read.spss can read it?
Dear All,
I would like to ask if it is possible to estimate a hazard function
using the muhaz command when all the data is right-censored. My data
has information of the number of weeks people has been unemployed but
all of them are unemployed at the date of the survey, that is, I cannot
obser
Could anyone please explain how to set up contrasts between means in R. I want to know
if "before I conduct an experiment and believe the mean for 1 and 2 will be different
from means 3 and 4, Is this true?" That is what I have to prove or disprove, I
thought that contrasts would be the way to
David Baird has created a free data conversion utility
called dataload. See
http://www.vsn-intl.com/genstat/downloads/dataload.htm
It can translate SPSS files into other formats
including several that are readable by R (including csv
and rda). I don't use SPSS myself but I have used
data
R help list subscribers,
I am a new user of R. I am attempting to use R to explore a set of
equations specifying the dynamics of a three trophic level food chain. I
have put together this code for the function that is to be evaluted by
LSODA. My equations Rprime, Cprime, and Pprime are meant
Jake -
Very puzzling that a freshly-written copy of the file gives the
same problem.
As to your second question, my recollection is that read.spss()
will be in package 'foreign'. The absolute reference is the
source code itself, and this is often quite readable, whether
it's commented or not.
Hello,
I am running a program in r that calls a function, which calls another
function, which calls another etc. These functions are of the form:
example<- function(x,y,z)
{x, y, and z are defined within curly braces like this}
Here's my question. To start the main function, I input as an i
Here is an answer to a 1999 post. I didn't find a direct answer anywhere
on the Web, perhaps because it is "obvious" once one sees it.
Suppose you have data from a longitudinal study, where each subject was
measured *up to* four times, with missing measurements, so that the data
look like this:
Jake -
The error message and warnign message shown below say something
is wrong with this file's SPSS system-file header. If you are
really able to open this one in SPSS, do so, change maybe a
column name or row name or two, and save it again under a
different file name. See if read.spss() cho
Is there any documentation on what kind of SPSS file can and cannot be
read by read.spss? Alternatively, how can one modify or "clean" an SPSS
file to make it readable by read.spss? What properties must a *.sav file
before read.spss can read it?
The file in this example is 270KB, with 5 rows an
I have fitted a logistic regression model
> failed.lr2$call
lrm(formula = failed ~ Age + task2 + Age:task2, data = time.long,
na.action = na.omit)
using the Design package functions and would like to generate a
nomogram from this model.
the datadist information is generated and stored in
>
In tseries, look for ?adf.test & ?pp.test.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Erin Hodgess
Sent: Wednesday, November 05, 2003 8:20 AM
To: [EMAIL PROTECTED]
Subject: [R] Cointegration
Do any packages exist for cointegration, please?
Do we ne
Is this a homework assignment ?
- tom blackwell - u michigan medical school - ann arbor -
On Tue, 4 Nov 2003, Feng Zhang wrote:
> Hey, R-listers
>
> When we say a function f(t) is smooth, does this mean that
> f has infinite differentials with respect to t?
>
> Or any other formal definit
Thanks to those who pointed me at the solutions to the legend overprinting the bars.
I took the "easy" way of rescaling the y axis, picking the scaling factor for stacked
bars is somewhat problematic but sufficient for my application.
I have another couple of barplot questions:
- Can I
Do any packages exist for cointegration, please?
Do we need them, if the answer to the previous is no, please?
Thanks,
Erin
mailto: [EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Hey, R-listers
When we say a function f(t) is smooth, does this mean that
f has infinite differentials with respect to t?
Or any other formal definition on this?
Thanks for your points.
Fred
[[alternative HTML version deleted]]
__
[EMAIL PR
> From [EMAIL PROTECTED] Wed Nov 5 11:31:52 2003
>
> but I cannot find any function to find the avalable region names in
> a map database
>
They are in the .N file for the appropriate map.
Try:
library(maps)
data(worldMapEnv)
file.show(paste(Sys.getenv(worldMapEnv), "world.N", sep="/"))
These
On 5 Nov 2003 at 8:21, Ray Brownrigg wrote:
> Except the Bolivia map must be a local one (from [EMAIL PROTECTED]);
> it is not in library(mapdata).
>
> Ray Brownrigg
Sorry, I should have typed
> library(maps)
> library(mapdata)
> map("worldHires","bolivia")
which indeed shows a correctly-lo
On 4 Nov 2003 at 13:14, Stephane DRAY wrote:
That the matrix is symmetric is sufficient to guarantee real
eigenvalues (proof is easy). I don't know about any general
conditions for asymmetric matrices, and doubt there are.
But in many structured situation you could be able to show similarity
"Bill Shipley" <[EMAIL PROTECTED]> writes:
> Hello. I am trying to determine whether I should be using ML or REML
> methods to estimate a linear mixed model. In the book by Pinheiro &
> Bates (Mixed-effects models in S and S-PLUS, page 76) they state that
> one difference between REML and ML is
Look into external pointers. That is how we have tackled this, e.g. in
the ts package.
On Tue, 4 Nov 2003, Ross Boylan wrote:
> I have a likelihood I would like to compute using C++ and then
> optimize. It has data that need to persist across individual calls to
> the likelihood. I'd apprecia
I have a likelihood I would like to compute using C++ and then
optimize. It has data that need to persist across individual calls to
the likelihood. I'd appreciate any recommendations about the best way
to do this. There are several, related issues.
1. Use the R optimizer or a C optimizer?
Beca
At 14:40 04/11/2003, [EMAIL PROTECTED] wrote:
>
>I'm not familiar with the definition of an indefinite matrix.
I see it in Harville, D.A., Matrix algebra from a statistician's
perspective. 1997, New-York: Springer-Verlag. 630.
it is a matrix A that is neither positive definite nor positive defi
I'm not familiar with the definition of an indefinite matrix.
(I assume you mean the eigenvalues are real...)
If A is a symmetric real valued matrix then the eignevalues of A are real.
There is a short proof of this due to Styan that can be found, among other
places, in Searle's 'Matrix Algeb
Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
>
> On Tue, 4 Nov 2003 [EMAIL PROTECTED] wrote:
>
> > On 3 Nov 2003 at 15:04, [EMAIL PROTECTED] wrote:
> >
> > >
> > > I believe in R you need
> > >
> > > library(maps)
> > >
> > > before using usa(). Of course, the library has to be installed fir
Hello list,
Sorry, these questions are not directly linked to R.
If I consider an indefinte real matrix, I would like to know if the
symmetry of the matrix is sufficient to say that their eigenvectors are real ?
And what is the conditions to ensure that eigenvectors are real in the case
of an a
> >
> For independent samples there is a simple technique described by Larsen,
> Am. Stat, May 1992
> for turning summary statistics back into a data set with an equivalent
> ANOVA summary.
>
> This is implemented and described in my STAT2DAT macro,
> http://www.math.yorku.ca/SCS/sasmac/stat2d
> Hi,
>
> I would like to interface a C code into R. Is it possible to use in the C code,
> functions from a R package (for instance, to use pmvnorm within loops in the C code
> and to call the result in a R function)?
Probably the most easiest way is copying the C-sources from the original
pa
I am still waitting to hear from someone are you still in buissness when i access the
web site it says account susspended. looking for payless4smokes.com i am a faithful
costomer i would like to know what is going on . thanks dawn corvin- Original
Message -
From: DAWN CORVIN
To:
Sent:
Neil -
Maybe also the "Function and variable index", pages 94-96, and
the "Concept index", pages 97-98 in "An Introduction to R",
cran.r-project.org/doc/manuals/R-intro.pdf.
- tom blackwell - u michigan medical school - ann arbor -
On Tue, 4 Nov 2003, RINNER Heinrich wrote:
> The
Dear Bill,
I am not a lme-expert, but I believe the Pinheiro&Bates' book is rather
clear here.
However you know that a lme model is, for instance
fixed= y~x1+x2 and random=y~x1|group
and you can fit it by ML or REML.
If you are interested in testing for x2 by means the LRT (namely by
comparing
There is a very nice Search Engine, with Keywords arranged by topic, that
comes with R.
If you are under Windows (and have Html help installed), you can find this
in the Help menu, under "Html help".
(The relevant file can also be found as
...\doc\html\search\SearchEngine.html.)
Maybe this helps,
Hello. I am trying to determine whether I should be using ML or REML
methods to estimate a linear mixed model. In the book by Pinheiro &
Bates (Mixed-effects models in S and S-PLUS, page 76) they state that
one difference between REML and ML is that « LME models with different
fixed-effects struc
Timur Elzhov <[EMAIL PROTECTED]> writes:
> On Tue, Nov 03, 2003 at 17:06:37 +0100, you wrote:
>
> >> REAL(x)[0] = 10;
> >> REAL(y)[0] = 20;
> >> SETCADR(R_fcall, x);
> >> SETCADR(R_fcall, y);
> >> ...
> >> .Call("f", c, new.env()) returns
> >> `.Primitive("c")(20)'
> >
That's usually what help.search() is for. For example,
help.search("file", package = "base")
-roger
Neil Osborne wrote:
Hi,
I'm new to R and I'm finding it quite a chore trawling through the R
documentation to find a function to carry out simple atomic tasks. Is
any one aware of R help docu
Dear R-Users:
Which is the best way to fee memory inside a function to
make room for some other process that uses most of the
RAM available and latter retrive this data again to
finish the function?
Thank you very much for your help
Kenneth Cabrera
--
___
Hi,
I'm new to R and I'm finding it quite a chore trawling through the R
documentation to find a function to carry out simple atomic tasks. Is any
one aware of R help documentation that is aranged in functional categories
for e.g.:
String manipulation
File I/O
Dataframe, List manipulation
etc,
On Tue, Nov 03, 2003 at 17:06:37 +0100, you wrote:
>> REAL(x)[0] = 10;
>> REAL(y)[0] = 20;
>> SETCADR(R_fcall, x);
>> SETCADR(R_fcall, y);
>> ...
>> .Call("f", c, new.env()) returns
>> `.Primitive("c")(20)'
> Hmmm. This is definitely not right:
>
> SETCADR(R_fcall
Dear Bernd,
As Brian Ripley points out, the problem is that you can't break a quoted
string over input lines.
It is nevertheless awkward that you have to provide the recode directives
in one long line. I'll think about alternatives for the next version of the
car package. One simple possibilit
I strongly encourage you tho check this book:
The basics of S and S-Plus / Andreas Krause, Melvin Olson
It has a very detailed and easy to follow info on how to link C code with S
(which is nearly identical to R). The R documentation is certainly very
complete, but it is written for people that a
On Tue, 4 Nov 2003 [EMAIL PROTECTED] wrote:
> On 3 Nov 2003 at 15:04, [EMAIL PROTECTED] wrote:
>
> >
> > I believe in R you need
> >
> > library(maps)
> >
> > before using usa(). Of course, the library has to be installed first.
> >
>
> Yes, I new that. When I did help.search("usa") I HAD m
Nothing to do with recode: you cannot have line breaks inside quoted
strings.
> "c(101,25,167,45,75)=25;
Error: syntax error
> c(104,51)=51"
Error: syntax error
On Tue, 4 Nov 2003, Bernd Weiss wrote:
> Dear all,
>
> it seems to be that 'recode' can't handle any line breaks in its
On Tue, 4 Nov 2003, Jeffrey Todd Lins wrote:
>
> I recently attempted to read a .txt file using both read.table("
> ",header=TRUE) and read.delim(" ",header=TRUE) and received the
> following message
Not from those commands, though.
> Error in edit.data.frame(get(subx, envir = parent), ...) :
On 3 Nov 2003 at 15:04, [EMAIL PROTECTED] wrote:
>
> I believe in R you need
>
> library(maps)
>
> before using usa(). Of course, the library has to be installed first.
>
Yes, I new that. When I did help.search("usa") I HAD maps
attached. It dosen't have any function usa(). R is more interna
Dear all,
it seems to be that 'recode' can't handle any line breaks in its code.
The following command causes no problem:
datameta$smpid.r <-
recode(datameta$smpid,"c(101,25,167,45,75)=25;c(104,51)=51")
But if I type ...
datameta$smpid.r <-recode(datameta$smpid,
"c(101,25,167,45,75)=2
Hello,
I'm not sure, but if I understand your question well, we have had such
questions already plenty of times. I guess, Google/STFW/RTFM/writing R
extensions. It's easy, it's well-explained, ...
Kurt.
Peyrard Nathalie wrote:
Hi,
I would like to interface a C code into R. Is it possible to
I recently attempted to read a .txt file using both read.table(" ",header=TRUE)
andread.delim(" ",header=TRUE)and received the following message
Error in edit.data.frame(get(subx, envir = parent), ...) :
symbol print-name too long
I am able to also create a variable x<
I think "Writing R Extensions" may help you with that...
On Tue, 4 Nov 2003, Peyrard Nathalie wrote:
> Date: Tue, 04 Nov 2003 10:19:00 +0100
> From: Peyrard Nathalie <[EMAIL PROTECTED]>
> To: [EMAIL PROTECTED]
> Subject: [R] interfacing C into R and R packages
>
> Hi,
>
> I would like to interf
thanks to Paulo Justiniano, I saved an error but still can't build R
(make -k repeat the same error on each package) :
**
-- Making package base
adding build stamp to DESCRIPTION
installing demos
Hi,
I would like to interface a C code into R. Is it possible to use in the C code, functions from a R package (for instance, to use pmvnorm within loops in the C code and to call the result in a R function)?
Nathalie
__
[EMAIL PROTECTED] mailing list
Hi
Claus Gwiggner wrote:
> > there is no Oracle Client on my machine (Redhat Linux 9.0/686i,
> > R1.8.0). I cannot thus compile the ROracle package.
You can in principle: Copy libs and the tnsconfig from the Oracle server
manually (no fun, but at least for Ora8i I managed this in some 30 min).
Thanks to Andy and Thomas,
Reading help(hclust) more carefully would have done it but
sometimes you do not
see the wood for the trees...
So hc$merge does exactly what I want.
I have never been aware of the command str to get the structure of
an R-object. It
seems pretty useful to me.
Thanks,
Sorry for re-posting this message; the first one was lost in another
thread.
>Hi,
> there is no Oracle Client on my machine (Redhat Linux 9.0/686i,
> R1.8.0). I cannot thus compile the ROracle package.
> I read the readme.client file and did:
> I unziped the package ROracle and copied it into t
Thanks to Andy and Thomas,
Reading help(hclust) more carefully would have done it but sometimes you do not
see the wood for the trees...
So hc$merge does exactly what I want.
I have never been aware of the command str to get the structure of an R-object. It
seems pretty useful to me.
Thanks,
A
You could compile it from the sources: there is a repository of the
sources in the same CRAN area as the current sources.
Or you could correct your code.
On Mon, 3 Nov 2003, Ernesto Calvo wrote:
> I am looking for the binary rw1051. I have some code that will only run in version
> 1.5.1 and I
67 matches
Mail list logo