Dear Roland,
Have you looked at Zelig (http://gking.harvard.edu/zelig)?
Several professors in my department are going to use it to teach
R to political science undergraduates and graduate students this
fall. We just presented it at the Political Methodology
meeting, and will present it again at t
Part of that decision may depend upon how big the dataset is and what is
intended to be done with the ID's:
> object.size(1011001001001)
[1] 36
> object.size("1011001001001")
[1] 52
> object.size(factor("1011001001001"))
[1] 244
They will by default, as Andy indicates, be read and stored as do
I think I agree with Prof. Lumley: For me, R would likely be the
platform of choice for most new statistical algorithm development. It
has a steep learning curve, and I would therefore not recommend it for
anyone who would not likely be involved in this kind of activity (unless
my targe
Hi Roland and folks:
Roland's and Susanna's comments were very interesting. I think we need to keep in mind
that unlike SPSS, R is more than a statistical package. It's a tool. SPSS does indeed
make life easy, particularly for data preprocessing (or data "cleaning"), something
that is perhaps n
On Sat, 14 Aug 2004, Matthew Walker wrote:
> I have read (in a very recent email to r-help) that R's methods package
> does not support private methods. However, I also looked at the source
> for the function "is":
>
> > is
> function (object, class2)
> {
> cl <- .class1(object)
> if (mis
If I'm not mistaken, numerics are read in as doubles, so that shouldn't be a
problem. However, I'd try using factor or character.
Andy
> From: Dan Bolser
>
> I store an id as a big number, could this be a problem?
>
> Should I convert to at string when I use read.table(...
>
> example id's
>
.class1 is an 'internal' method in the namespace of the package 'methods'.
It can be accessed with the ::: operator. A good start for documentation
may be Vol 3/1 of R-News.
>methods:::.class1
function (x)
class(x)[[1]]
--Matt
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAI
I store an id as a big number, could this be a problem?
Should I convert to at string when I use read.table(...
example id's
1001001001001
1001001001002
...
1002001002005
Bigest is probably
1011001001001
Ta,
Dan.
__
[EMAIL PROTECTED] mailing li
Hi,
I'm trying to understand R's object oriented abilities.
I have read (in a very recent email to r-help) that R's methods package
does not support private methods. However, I also looked at the source
for the function "is":
> is
function (object, class2)
{
cl <- .class1(object)
if (miss
There are quite a few helpful free contributions--under contributed go
figure-- on the R site.
"Using R for Data Analysis and Graphics" by John Maindonald (PDF [702kB],
data sets and scripts are available at JM's homepage).
I found this a big help( Thanks to John Maindonald) as well as trying to
s
Hello to all,
There are several functions to perform scatter diagrams with labels in the
library 'ade4'.
#
library(ade4)
bob <- as.data.frame((array(1:5, c(4,2
row.names(bob)
row.names(bob) <- paste("row",1:4,sep=".")
row.names(bob)
par(mfrow=c(1,2))
s.label(bob, su
"An Introduction to R" by Venables, Smith, and the R Development Core Team is a
nicely written little paperback that is even easier on the pocketbook ($14 on
Amazon.com) than Dalgaard's . It's a nicely written basic introduction the the
language that makes no attempt to teach statistics.
Cheers,
B
There are a number of resources but I believe the best for starters is
_Introductory Statistics with R_, by Peter Dalgaard. It's inexpensive ($35)
and very well written, easy to read and understand even if you don't know a
lot of statistics. In fact the book's intention is to teach you stats by
w
An Introduction to R is pretty good
(http://cran.r-project.org/doc/manuals/R-intro.pdf). And there are
many others available at
http://cran.r-project.org/other-docs.html
-roger
T. Murlidharan Nair wrote:
Hi !!
Is there a good tutorial for the R language ? I really find it hard to
find methods
Hi !!
Is there a good tutorial for the R language ? I really find it hard to
find methods in R.
Thanks and Cheers ../Murli
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-proj
You might find it easier to do this type of thing with an open-source image
analysis program such as ImageJ: http://rsb.info.nih.gov/ij/
Rob
- Original Message -
From: "Doug Bourne" <[EMAIL PROTECTED]>
To: <[EMAIL PROTECTED]>
Sent: Friday, August 13, 2004 1:10 PM
Subject: [R] cell count
I have a bunch of jpg files that I need to count cells on. Some cells are green and
some cells are red. Is there a function like locator that I could use to manually
count the cells? I would like to be able to click on a cell, have it marked so I
don't count it again and have a counter keep t
Marc and Partha,
Thank you both for your suggestions of using par("mfrow"), that's exactly
what I needed to do. In fact I already had it in the code (which I didn't
write originally) and simply had to put my postscript statement in front of
it. :)
It works great now! Thanks for the quick respon
Have you looked at ?par("mfrow") ?
"Johnson, Heather" <[EMAIL PROTECTED]>
Sent by: [EMAIL PROTECTED]
08/13/2004 12:28 PM
To: "'[EMAIL PROTECTED]'" <[EMAIL PROTECTED]>
cc:
Subject:[R] Question from Newbie on PostScript and Multiple Plots
Hi,
As I'm pr
On Fri, 2004-08-13 at 11:28, Johnson, Heather wrote:
> Hi,
>
> As I'm pretty new to R I hope this question isn't too basic.
>
> I am currently looping through my dataset and for each iteration am
> producing three separate plots. When I output these plots to the screen
> they are nicely grou
Hi,
As I'm pretty new to R I hope this question isn't too basic.
I am currently looping through my dataset and for each iteration am
producing three separate plots. When I output these plots to the screen
they are nicely grouped as three plots per page, however, when I try to send
it to a Po
Chris,
I understood the following: you want to try every single covariate in a
cox model with only one covariate, then take the best ones according to
p-value.
Assume your columns look like:
stop status event x1 x2 x3 etc
You want to add column 3 (x1), then 4, etc.
I suggest a for() loop:
z<-NU
"by" is perhaps what you looking for
try to write "?by" for help
>I have a panel data set (25 years) and I would like to break
>it up by year, and run a few tests. How is this done in R?
>In STATA it is quite simple (keep if year =1970). If someone
>could give me an example of how to do it f
You probably want something like
subset(dataset, year == 1970)
and loop through the years.
Or, perhaps even classier, would be
datalist <- split(dataset, as.factor(year))
results <- lapply(datalist, myAnalysisFunction)
-roger
Justin G. Gardner wrote:
Hello R list.
I am sure this is an easy question
>First, I do not know how to specify such a link function in R.
>Second, if I can specify such alink, I could use (in place of d*j), the
>smooth baseline estimated after doing a Cox regression. But I don't know
>how to fit (for instance) a piecewise constant baseline hazard with a
>Poison glm
I fo
On Fri, 2004-08-13 at 09:22, Luis Rideau Cruz wrote:
> R-help
>
> Is there any option to get closer the x-axis and names.arg from barplot?
>
> Thank you
Using mtext() you can do something like the following:
data(VADeaths)
# Now place labels closer to the x axis
# set 'axisnames' to FALSE so t
My deep appreciation to everyone who responded to my question. I am
digesting your proposals (a little distracted by my governor's dramatic
resignation in New Jersey). It seems that I am able to implement the
algorithm in R (with the framework kindly provided by Thomas). I will
post the code when i
"Justin G. Gardner" <[EMAIL PROTECTED]> writes:
> Hello R list.
>
> I am sure this is an easy question, please forgive my
> ignorance of R.
>
>
> I have a panel data set (25 years) and I would like to break
> it up by year, and run a few tests. How is this done in R?
> In STATA it is quite
Thanks Mayeul,
I actually would like to test each variable individually and use those have low
p-value to build a classifier (not in cox model). Therefore, I need to write a
function to subset those low p-value variables, instead of putting them as covariates.
Any ideas?
Chris
-Original M
Hello R list.
I am sure this is an easy question, please forgive my
ignorance of R.
I have a panel data set (25 years) and I would like to break
it up by year, and run a few tests. How is this done in R?
In STATA it is quite simple (keep if year =1970). If someone
could give me an example
R-help
Is there any option to get closer the x-axis and names.arg from barplot?
Thank you
Luis Ridao Cruz
Fiskirannsóknarstovan
Nóatún 1
P.O. Box 3051
FR-110 Tórshavn
Faroe Islands
Phone: +298 353900
Phone(direct): +298 353912
Mobile: +298 580800
Fax: +298
On Fri, 13 Aug 2004, Mayeul KAUFFMANN wrote:
>
> > you can always use parametric models and a
> > full likelihood (but you may have to program them yourself).
> > Prof Brian Ripley
>
> I started trying this but I could not make the counting process notation
> work on this.
> (Andersen, P.K. and Gi
Clécio da Silva Ferreira wrote:
Hi
I'm a Brazilian student and I'd like to work with R.
But, I had some problems to install R for windows with shortcut of the
statisticals packages in tool bar.
What do you mean with "shortcut of the statisticals packages in tool
bar." What are the "statistical
> > coxph does not use any information that are in the dataset between
event
> > times (or "death times") , since computation only occurs at event
times.
> This is the consequence of the use of partial likelihood in the Cox
model.
>[...]
>You need to make more assumptions, such as a
>smooth baselin
Hi
I'm a Brazilian student and I'd like to work with R.
But, I had some problems to install R for windows with shortcut of the
statisticals packages in tool bar.
I'll be glad if you could help me to install it without miss the additionals
packages I had install (more than 50!!)
Thanks!
Cléc
> Yes, but my control group is not always in the first colonne.
> So, I can have to a matrix like this:
> 1 0 0 -1 0
> 0 1 0 -1 0
> 0 0 1 -1 0
> 0 0 0 -1 0
>
see the `base' argument to `contrMat'.
Torsten
> Laurent Houdusse
> Analyste Programmeur
>
>
>
> -Message d'origine--
Ok,
I use simtest because i retrieve the same result
But i don't know how to do the comparaison of P
Laurent Houdusse
Analyste Programmeur
-Message d'origine-
De : Liaw, Andy [mailto:[EMAIL PROTECTED]
Envoyé : vendredi 13 août 2004 14:24
À : Laurent Houdusse; '[EMAIL PROTECTED]'
Objet
> From: Laurent Houdusse
>
> However, on several tests, my results are correct:
> The t-value seems always to correspond to the q'value calculated by
> SigmaStat...
>
> So,What do I have to use to find the same results as
> SigmaStat (value and
> comparaison)?
Seems like you don't have much
Yes, but my control group is not always in the first colonne.
So, I can have to a matrix like this:
1 0 0 -1 0
0 1 0 -1 0
0 0 1 -1 0
0 0 0 -1 0
Laurent Houdusse
Analyste Programmeur
-Message d'origine-
De : Torsten Hothorn [mailto:[EMAIL PROTECTED]
Envoyé : vendredi 13
On Fri, 13 Aug 2004, Liaw, Andy wrote:
> Before you do any of that, you should realize the fact that simtest does
> _not_ do Dunnett's test: It can use `Dunnett' _contrasts_ to do comparisons
> with a control, but the actual procedure is different. You are unlikely to
> get the same result as o
However, on several tests, my results are correct:
The t-value seems always to correspond to the q'value calculated by
SigmaStat...
So,What do I have to use to find the same results as SigmaStat (value and
comparaison)?
Laurent Houdusse
Analyste Programmeur
-Message d'origine-
De :
Before you do any of that, you should realize the fact that simtest does
_not_ do Dunnett's test: It can use `Dunnett' _contrasts_ to do comparisons
with a control, but the actual procedure is different. You are unlikely to
get the same result as other packages that perform standard Dunnett's tes
Hi!
I use simtest fonction of multcomp package to compile a Dunnett's test.
I have 10 treatments and one control group, so i create a matrix with:
m<-matrix(0,10,11)
m[1,1]<--1
m[1,2]<-1
m[2,1]<--1
m[2,3]<-1
m[3,1]<--1
m[3,4]<-1
m[4,1]<--1
m[4,5]<-1
m[5,1]<--1
m[5,6]<-1
m[6,1]<--1
m[6,7]<-1
m[7,1
"Luis Rideau Cruz" <[EMAIL PROTECTED]> writes:
> R-help,
>
> I wish to replace NULL elements(or missing) in the following list :
...
> If I use
> lapply(z2,function(x) ifelse(is.null(x),0,as.numeric(x)))
...
> i.e,only the first element in z2
ifelse() gives a result that is the same length
Hi Luis,
try to use,
lapply(z2, function(x) if(is.null(x)) c(0,0) else x)
I hope this helps.
Best,
Dimitris
Dimitris Rizopoulos
Doctoral Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/396887
Fa
R-help,
I wish to replace NULL elements(or missing) in the following list :
> z2
$cod
mean sd
62.56190 12.65452
$haddock
mean sd
36.61490 11.50365
$ling
mean sd
86.17949 20.43587
$saithe
meansd
50.275847 5.453606
$whiting
NULL
$"norway pout"
>How can i plot an array and instead of having on the x labels the
>indexes of the array I want to display an other String array of the
>same length
Do this:
> plot(myarray,xaxt="n",xlab="")
> axis(1,at=1:length(myarray),lab=my.vector.of.names)
__
[EM
On Fri, 13 Aug 2004 [EMAIL PROTECTED] wrote:
> Do you know if there is any function in R available which calculates an R^2
> for a robust MM-regression?
It's very easy to do, *but* why would you want a non-robust measure for a
robust fit?
[Just 1 - sum(residuals(fit)^2)/sum((y-mean(y))^2), for
> you can always use parametric models and a
> full likelihood (but you may have to program them yourself).
> Prof Brian Ripley
I started trying this but I could not make the counting process notation
work on this.
(Andersen, P.K. and Gill, R.D. (1982). Cox's regression model for counting
process
I've tried doing this by specifying x=TRUE, which provides me with a single
imputation, that has been useful. However, the help file possibly suggests that I
should get a flat-file matrix of n.impute imputations, presumably with indexing. I'm
a bit stuck using alternatives to aregImpute, as ne
How can i plot an array and instead of having on the x labels the
indexes of the array I want to display an other String array of the
same length
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Dear Jason,
Of course you can do almost eveything in R, but I have rarely seen someone
prototyping something in Java to then implement it in R. KD-trees are for
performance, and interpreted R is not really fast. So interfacing some
foreign code might really be the better choice here.
Opportunity
Joern Kamradt wrote:
Hi,
I am using the rgl package for 3D display. Unfortunately, I am not able
to get the snapshot command running.
I tried the following:
> example(rgl.surface)
rgl.sr> data(volcano)
rgl.sr> y <- 2 * volcano
rgl.sr> x <- 10 * (1:nrow(y))
rgl.sr> z <- 10 * (1:ncol(y))
rgl.sr> y
Dear R-users,
Do you know if there is any function in R available which calculates an R^2
for a robust MM-regression?
With best regards,
Carsten Colombier
Dr. Carsten Colombier
Economist
Group of Economic Advisers
Swiss Federal Finance Administration
Bundesgasse 3
CH-3003 Bern
phone +41 31 322
54 matches
Mail list logo