We would like to announce the release of our software MatchIt, now
available on CRAN. MatchIt implements a variety of matching methods for
causal inference.
Abstract:
MatchIt implements the suggestions of Ho, Imai, King, and Stuart (2004)
for improving parametric statistical models by prepr
Sorry I didn't read the question properly. Please disregard, my mind was
elsewhere.
Tom
> -Original Message-
> From: Mulholland, Tom
> Sent: Thursday, 13 January 2005 10:52 AM
> To: Dr Carbon; r-help@stat.math.ethz.ch
> Subject: RE: [R] Finding seasonal peaks in a time series
>
>
You stand more chance if you do it yourself
https://stat.ethz.ch/mailman/listinfo/r-help
> -Original Message-
> From: Kevin Ita [mailto:[EMAIL PROTECTED]
> Sent: Thursday, 13 January 2005 2:25 AM
> To: R-help@stat.math.ethz.ch
> Subject: [R] Please unsubscribe me from you list
>
>
> Plea
You might find breakpoints in strucchange helpful
Tom
> -Original Message-
> From: Dr Carbon [mailto:[EMAIL PROTECTED]
> Sent: Thursday, 13 January 2005 6:19 AM
> To: r-help@stat.math.ethz.ch
> Subject: [R] Finding seasonal peaks in a time series
>
>
> I have a seasonal time series
I have often noted that "statistics can't prove a damn thing, but they can be
really useful in disproving something." Having spent most of 80s and half of
the 90s with the Australian Bureau of Statistics to find out how you collect
these numbers, I am disconcerted at the apparent disregard for m
Hi, there.
I have two questions about the diagnostics in multivarite statistics.
1. Is there any diagnostics tool to check if a multivariate sample is from
multivariate normal distribution? If there is one, is there any function
doing it in R?
2. Is there any function of testing if two multivari
There's something that you're not telling me. If you want something other than
your data, why use your data. If you have another set of data that you wish to
overlay on the image then you are going to have to scale one of the data
sources to match the other. I'm not sure where your problem is bu
> From: Weiwei Shi
>
> Hi, there:
> Thanks a lot for all people' prompt replies.
>
> In detail, I am facing a huge amount of data: over
> 10,000 and 400 vars. This project is very challenging
> and interesting to me. I tried rpart which gives me
> some promising results but not good enough. So I
On Wed, 12 Jan 2005, Berton Gunter wrote:
>R-Listers.
>
>The following is a rant originally sent privately to Frank Harrell in
>response to remarks he made on this list. The ideas are not new or original,
>but he suggested I share it with the list, as he felt that it might be of
>wider interest, n
Hi, there:
Thanks a lot for all people' prompt replies.
In detail, I am facing a huge amount of data: over
10,000 and 400 vars. This project is very challenging
and interesting to me. I tried rpart which gives me
some promising results but not good enough. So I am
trying randomForest and gbm now.
R-Listers.
The following is a rant originally sent privately to Frank Harrell in
response to remarks he made on this list. The ideas are not new or original,
but he suggested I share it with the list, as he felt that it might be of
wider interest, nonetheless. I have real doubts about this, and I
You can also check out:
http://www.i-pensieri.com/gregr/gbm.shtml
There are reference papers on there, too.
HTH,
Danny
On Wed, 12 Jan 2005, Peter Dalgaard wrote:
> Weiwei Shi <[EMAIL PROTECTED]> writes:
>
> > Hi, there:
> > I am wondering if I can find some detailed explanation
> > on gbm or
Weiwei Shi <[EMAIL PROTECTED]> writes:
> Hi, there:
> I am wondering if I can find some detailed explanation
> on gbm or explanation on examples of gbm.
What is gbm?
Green Belt Movement?
Georgie Boy Manufacturing?
I'm serious! Well, only sort of, but try Google on "gbm" and you'll
find those tw
I just got 25 hits from "www.r-project.org" -> search -> "R site
search". Might one or more of these help you? If they don't solve your
problem, I suggest you try "the posting guide!
http://www.R-project.org/posting-guide.html";. If that still doesn't
solve your problem, it should help
I have a seasonal time series. I want to calculate the annual mean
value of the time series at its peak
(say the mean of the three values before the peak, the peak, and the
three values after the peak).
The peak of the time series might change cycle slightly from year to year.
# E.g.,
nPts <-
Hi, there:
I am wondering if I can find some detailed explanation
on gbm or explanation on examples of gbm.
thanks,
Ed
__
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-
To all that replied, thanks... I have a clue where I can change the settings.
tnx,
Kurt Sys
>- Oorspronkelijk bericht -
>Van
: Gabor Grothendieck [mailto:[EMAIL PROTECTED]
>Verzonden
: woensdag
, januari
12, 2005 05:26 PM
>Aan
: r-help@stat.math.ethz.ch
>Onderwerp
: Re: [R] changing lan
Apparently the message below wasn't posted on R-help, so I'm sending it
again. Sorry if you received it twice.
--- bogdan romocea <[EMAIL PROTECTED]> wrote:
> Date: Tue, 11 Jan 2005 17:31:42 -0800 (PST)
> From: bogdan romocea <[EMAIL PROTECTED]>
> Subject: Re: [R] global objects not overwritten w
Dan Bolser wrote:
On Wed, 12 Jan 2005, Frank E Harrell Jr wrote:
Dan Bolser wrote:
Hello,
I am making some use of ROC curve analysis.
I find much help on the mailing list, and I have used the Area Under the
Curve (AUC) functions from the ROC function in the bioconductor project...
http://www.b
Please unsubscribe me from your list.
Thank you.
Kevin
-
The all-new My Yahoo! What will yours do?
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://sta
Bang wrote:
When I installed R 2.0.1 (replacing 1.9.0) for Windows, a code using
"model.response" began acting up. Here are the first several lines of a
code I had been tweaking for a spatial model (the code is mostly that of
Roger Bivand--I am adapting it to a slightly different data structure a
On Wed, 12 Jan 2005, Frank E Harrell Jr wrote:
>Dan Bolser wrote:
>> Hello,
>>
>> I am making some use of ROC curve analysis.
>>
>> I find much help on the mailing list, and I have used the Area Under the
>> Curve (AUC) functions from the ROC function in the bioconductor project...
>>
>> http
I don't know what SAS does, but transfer functions are essentially MA/AR
from an ARMA model, so you should be able to get what you want from the
various ARMA estimation tools in R.
Paul Gilbert
Samuel Kemp (Comp) wrote:
Hi,
Does anyone know of a function in R that can estimate the parameters of
Troels:
It would be best if you discussed this with a local statistician to make
sure that the data and analysis are properly addressing the scientific
issues. Perhaps that is why no one replied to your previous post. Also, this
is primarily a **statistical** issue, not really an ** R-issue **.
H
On Wed, 12 Jan 2005, Bang wrote:
> When I installed R 2.0.1 (replacing 1.9.0) for Windows, a code using
> "model.response" began acting up. Here are the first several lines of a
> code I had been tweaking for a spatial model (the code is mostly that of
> Roger Bivand--I am adapting it to a slight
You have some package in use which defines an S4 method for
model.response, most likely SparseM. You didn't say that, but therein
lies the problem. It's hard for us to guess why it is failing to find
its default method.
Incidentally,
mf <- lm(formula, data, na.action = na.action, method = "m
It uses the language set by LC_TIME: see ?Sys.setlocale and ?format.Date
which references it.
On Wed, 12 Jan 2005, Kurt Sys wrote:
Hi all,
I've got a small, practical question, which untill now I couldn't solve
(otherwhise I wouldn't mail it, right?) First of all, I'm talking about R
2.0.1 on a w
Kurt Sys pandora.be> writes:
:
: Hi all,
:
: I've got a small, practical question, which untill now I couldn't solve
: (otherwhise I wouldn't mail it, right?) First of all, I'm talking about
: R 2.0.1 on a winxp system (using the default graphical interface being
: 'Rgui').
: When I make plo
For the KS-test of normality with estimated parameters see
?lillie.test in package nortest.
Best,
Christian
On Wed, 12 Jan 2005, Christoph Buser wrote:
> Hi Kwabena
>
> I did once a simulation, generating normal distributed values
> (500 values) and calculating a KS test with estimated
> param
Kurt Sys <[EMAIL PROTECTED]> writes:
> Hi all,
>
> I've got a small, practical question, which untill now I couldn't
> solve (otherwhise I wouldn't mail it, right?) First of all, I'm
> talking about R 2.0.1 on a winxp system (using the default graphical
> interface being 'Rgui').
> When I make pl
Christoph Buser wrote:
Hi Kwabena
I did once a simulation, generating normal distributed values
(500 values) and calculating a KS test with estimated
parameters. For 1 times repeating this test I got about
1 significant tests (on a level alpha=0.05 I'm expecting about 500
significant tests by
When I installed R 2.0.1 (replacing 1.9.0) for Windows, a code using
"model.response" began acting up. Here are the first several lines of a
code I had been tweaking for a spatial model (the code is mostly that of
Roger Bivand--I am adapting it to a slightly different data structure and
the proble
Hi Kwabena
I did once a simulation, generating normal distributed values
(500 values) and calculating a KS test with estimated
parameters. For 1 times repeating this test I got about
1 significant tests (on a level alpha=0.05 I'm expecting about 500
significant tests by chance)
So I think if
I try again - perhaps it is analysis of covariance with treatment
(thio,ultiva) as two categories and antime as covariate. On the basis of
such a model, is then the probability of GCS <= 12 larger with thio treatment ?
Dear friends, forgive me a simple question, possibly related to quantreg
but
Hi all,
I've got a small, practical question, which untill now I couldn't solve
(otherwhise I wouldn't mail it, right?) First of all, I'm talking about
R 2.0.1 on a winxp system (using the default graphical interface being
'Rgui').
When I make plots, using dates on the x-axis, it puts the labels
I think you will find all the doc in the help files
> ?factor()
gets
The function 'factor' is used to encode a vector as a factor (the
terms 'category' and 'enumerated type' are also used for factors).
If 'ordered' is 'TRUE', the factor levels are assumed to be
ordered. For co
On Wed, 12 Jan 2005, Tobias Verbeke wrote:
On Tue, 11 Jan 2005 10:45:48 -0500
"Palos, Judit" <[EMAIL PROTECTED]> wrote:
Breslow-Day test
A statistical test for the homogeneity of odds ratios.
[..some definitions..]
Your message was not particularly clear, but if
you were looking for R code to do a
Dear help desk and R community,
I have a problem on how R2.0 handle the RAM, maybe a bug
In fact I used R1.9 since january 2004 with large data set using a
macosx G5 with 1G of ram without problem . Then I passed to 2.0 and
I found myself short in ram using virtual memory. I tried to use
the
On Wed, 2005-01-12 at 15:17 +0100, [EMAIL PROTECTED] wrote:
> hi,
>
> I am trying to grow a classification tree on some data, but I have a little
> problem. In order to do so I have to use a function like "tree" in R and on
> the
> internet help(tree) I get the following:
>
> "The left-hand-sid
On 12 Jan 2005 at 15:17, [EMAIL PROTECTED] wrote:
>
> hi,
>
> I am trying to grow a classification tree on some data, but I have a
> little problem. In order to do so I have to use a function like "tree"
> in R and on the internet help(tree) I get the following:
>
> "The left-hand-side (resp
hi,
I am trying to grow a classification tree on some data, but I have a little
problem. In order to do so I have to use a function like "tree" in R and on the
internet help(tree) I get the following:
"The left-hand-side (response) should be either a numerical vector when a
regression tree will
On Tue, 2005-01-11 at 17:39 +0100, Rau, Roland wrote:
> Dear R-Help-List,
>
> are there any plans to organize a "useR" conference in 2005?
>
> Best,
> Roland
As I understand it, no. The next one will be in 2006, so it will be
every other year, interleaved with the DSC meeting the odd years.
I
Your example is sequential, ignoring the tree-like structure of most
documents. Why not via a DOM or similar "XML-ish" structure?
While I'd never advocate general purpose XML as a user format, as you
note, that is what XSLT is for, and using XML as an electronic
internal document representation wo
On Tue, 11 Jan 2005 19:33:41 + Rick Ram wrote:
> Groundwork for the choice of break method in my specific application
> has already been done - otherwise I would need to rework the wheel
> (make a horribly detailed comparison of performance of break
> approaches in context of modelling post br
On Tue, Jan 11, 2005 at 04:24:11PM +0100, Lefebure Tristan wrote:
> example from a shell:
>
> echo -e "pdf(file=\"test.pdf\")\nplot(1:10,11:20)\ndev.off(dev.cur())\n">cmd.R
> R -s
> (write a file of command for R, and than feed R with it)
This may be on the verge of becoming offtopic, but let
Hi R UseRs,
I am interested in providing Reporting abilities to R and have initiated a
project called R Output MAnager (ROMA).
My starting point was my R2HTML package which provides (rough) HTML
exportations. I began with trying to mimic it for LaTeX but fastly did
realize that it was a bad idea
avneet singh wrote:
Hello:
I have been trying to impute missing values of a data
frame which has both numerical and categorical values
using the function transcan() with little luck.
Would you be able to give me a simple example where a
data frame is fed to transcan and it spits out a new
data fram
Dan Bolser wrote:
Hello,
I am making some use of ROC curve analysis.
I find much help on the mailing list, and I have used the Area Under the
Curve (AUC) functions from the ROC function in the bioconductor project...
http://www.bioconductor.org/repository/release1.5/package/Source/
ROC_1.0.13.t
On Tue, 11 Jan 2005 10:45:48 -0500
"Palos, Judit" <[EMAIL PROTECTED]> wrote:
> Breslow-Day test
> A statistical test for the homogeneity of odds ratios.
>
[..some definitions..]
Your message was not particularly clear, but if
you were looking for R code to do a Breslow-Day test,
Google found t
(1) I do read the posting guide (the fact that I missread o
missunderstood something does not imply not reading)
(2) I could change NAs to 0 (I know) but I have previously (older
versions of R and SQL*Plus) used the same select with the "right" output
(namely with 0s).
(3) AFAIK "strange" is not
PLEASE do read the help page, which says
nullstring: character string to be used when reading 'SQL_NULL_DATA'
character items from the database.
^^^
so this does not apply to numeric items.
You can of course easily change numeric NAs to 0s, if you want to.
On Wed, 12
[EMAIL PROTECTED] writes:
> Could somebody please tell me what does the r-help-bounce address do?
>
> When I try to respond to an r-help post, my mailer (Lotus Notes) generates
> the r-help bounce return address automatically in addition to the original
> sender address and the r-help address. I
There is something strange in R behaviour (perhaps).
I have run the same select in Oracle SQL*Plus (version 10.1.0.2.0) and
the output comes out with NULLs (which is what it ougth to be).
But in R I still get the same result with NAs (no matter I use
na.strings or nullstring arguments)
An output
Does SuSE have anything like Debian's "alien" for importing foreign packages?
If so, then you might use the Debian packages.
Otherwise, it isn't too hard to unpack in your home directory, and add
the appropriate " (load-file ...)" statement to your .emacs /
.xemacs/init.el
best,
-tony
On
Hi Robert,
It looks like there is no way to explicitly make an S4 object call a
function when it is garbage collected unless you resort to tricks with
reg.finalizer.
It turns out that Prof. Ripley's reply (thanks!!) had enough hints in it
that I was able to get the effect I wanted by using R's ex
Hi,
On Wed, 12 Jan 2005, Bigby wrote:
> Hello,
>
> numerical.summary(Data), which didnt work, he suggested we try library(s20x)
> first, which came up with an error on my console. I have version 2.0.1 of R
library(s20x) is a package written by the Department of Statistics at the
University of Au
maybe something like:
f <- y ~ x1 + x2
attr(terms(f), "term.labels")
but this wan't work if you have a more complex formula (e.g., f <- y ~
x1*x2 + I(x1^2)) and you want only c("x1", "x2").
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Publ
Hello:
I have been trying to impute missing values of a data
frame which has both numerical and categorical values
using the function transcan() with little luck.
Would you be able to give me a simple example where a
data frame is fed to transcan and it spits out a new
data frame with the NA valu
Daniel Almirall wrote:
R-list,
1. Given a formula (f) w variables referencing some data set (dat), is
there any easier/faster way than this to get the names (in character form)
of the variables on the RHS of '~' ?
dat <- data.frame(x1 = x1 <- rnorm(100,0,1), x2 = x2 <- rnorm(100,0,1), y = x1
+ x
Breslow-Day test
A statistical test for the homogeneity of odds ratios.
Homogeneity
In systematic reviews homogeneity refers to the
degree to which the results of studies included in a review are similar.
"Clinical homogeneity" means that, in studies included in a review, the
participants, i
On Tue, 11 Jan 2005, Francisca xuan wrote:
here is a function I wrote
cdfest=function(t,lambda,delta,x,y){
a1=mean(x< t)
a2=mean(x< t-delta)
a3=mean(y1< t)
s=((1-lambda)*a1+lambda*a2-a3)^2
s
}
when I try to integrate over t, I got this message:
integrate(cdfest,0,4,lambda=0.3,delta=1,x=x,y=
Hello,
I am making some use of ROC curve analysis.
I find much help on the mailing list, and I have used the Area Under the
Curve (AUC) functions from the ROC function in the bioconductor project...
http://www.bioconductor.org/repository/release1.5/package/Source/
ROC_1.0.13.tar.gz
However,
On Tue, 11 Jan 2005, Marc Schwartz wrote:
>On Tue, 2005-01-11 at 14:59 +, Dan Bolser wrote:
>> On 10 Jan 2005, Peter Dalgaard wrote:
>>
>> >Dan Bolser <[EMAIL PROTECTED]> writes:
>> >
>> >> Cheers. This is really me just being lazy (as usual). The latex
>> function
>> >> in Hmisc allows me to
R-list,
1. Given a formula (f) w variables referencing some data set (dat), is
there any easier/faster way than this to get the names (in character form)
of the variables on the RHS of '~' ?
dat <- data.frame(x1 = x1 <- rnorm(100,0,1), x2 = x2 <- rnorm(100,0,1), y = x1
+ x2 + rnorm(100,0,1))
Cserháti Mátyás wrote:
Dear all,
Thanks to those 3 people who sent me answers to my question. Got
the problem solved. Great!
Now, another question of mine is:
I would like to run an R script from the Linux prompt. Is there any way
possible to do this? The reason is, the calculation that I'm doin
Dear R-Help-List,
are there any plans to organize a "useR" conference in 2005?
Best,
Roland
+
This mail has been sent through the MPI for Demographic Rese...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listin
here is a function I wrote
cdfest=function(t,lambda,delta,x,y){
a1=mean(x< t)
a2=mean(x< t-delta)
a3=mean(y1< t)
s=((1-lambda)*a1+lambda*a2-a3)^2
s
}
when I try to integrate over t, I got this message:
integrate(cdfest,0,4,lambda=0.3,delta=1,x=x,y=y1)
Error in integrate(cdfest, 0, 4, l
On Tue, 11 Jan 2005, bogdan romocea wrote:
Dear useRs,
I have a function that creates several global objects with
assign("obj",obj,.GlobalEnv), and which I need to run iteratively in
another function. The code is similar to
f <- function(...) {
assign("obj",obj,.GlobalEnv)
}
fct <- function(...) {
On Tue, Jan 11, 2005 at 03:59:58PM +0100, Cserh?ti M?ty?s wrote:
> I would like to run an R script from the Linux prompt. Is there any way
> possible to do this? The reason is, the calculation that I'm doing takes a
> few hours, and I would like to automatize it.
>
> Or does it mean that I have
Hi,
Would anyone be aware of an R package implementing the functionality
found in Meeker's SPLIDA software written for S-Plus? I don't know if
anyone has tried to port the s/w to R directly, or if equivalent
functions are available within another package.
Thanks in Advance.
- Colin
Colin Cunni
R-help,
I'm using the RODBC package to retrieve data froma ODBC database which
contain NA's.
By using the argument nullstring = "0" in sqlQuery() I expect to
coerce them to numeric but still get NA's in my select.
I'm running on Windows XP
> version
_
platform i386-pc-mi
On Tue, 2005-01-11 at 14:59 +, Dan Bolser wrote:
> On 10 Jan 2005, Peter Dalgaard wrote:
>
> >Dan Bolser <[EMAIL PROTECTED]> writes:
> >
> >> Cheers. This is really me just being lazy (as usual). The latex
> function
> >> in Hmisc allows me to make a .ps file then grab a screen shot of
> that
I would suggest reading the posting guide,
(http://www.r-project.org/posting-guide.html) and give a reproducible
example, with the error message that you received. As is, I have no idea
what you are doing here, and certainly cannot run this code. You use
"..." as an argument to your functions (why
Hello,
My name is Graham, I am an engineering student and my lecturer wishes us to
get a numerical summary of some data. He said use the command
numerical.summary(Data), which didnt work, he suggested we try library(s20x)
first, which came up with an error on my console. I have version 2.0.1 of R
see man R
example from a shell:
echo -e "pdf(file=\"test.pdf\")\nplot(1:10,11:20)\ndev.off(dev.cur())\n">cmd.R
R -s Dear all,
>
> Thanks to those 3 people who sent me answers to my question. Got
> the problem solved. Great!
>
> Now, another question of mine is:
>
> I would like to run an R scrip
75 matches
Mail list logo
