Hi
The October 2006 issue of R News is now available on CRAN under the
Documentation/Newsletter link.
Torsten
(on behalf of the R News Editorial Board)
___
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-announce
__
On Thu, 19 Apr 2007, Florian Koller-Meinfelder wrote:
Dear R-helpers,
I am looking for a segmentation package that gives some "tree identifier"
as output for every observation in the data set (my response variable is
binary). I have skimmed through "rpart", "ada" and "adabag": The output
"tree
On Fri, 2 Feb 2007, Henric Nilsson (Public) wrote:
> Torsten, consider the following:
>
>> ### ordinal regression
>> mammoct <- ctree(ME ~ ., data = mammoexp)
> Warning message:
> no admissible split found
>> ### estimated class probabilities
>> treeresponse(mammoct, newdata = mammoexp[1:5, ])
>
elated to the ordinal
> case (predicted probabilities > 1) and I've CC:ed the package's maintainer
> (Torsten Hothorn).
>
yep, you are right (as always)-- a bug introduced by a fix, grrr. Its a
little bit more complicated, but I'll make correct predictio
On Tue, 9 Jan 2007, Walter Durka wrote:
dear all,
Walter,
_please_ cc questions on contributed packages to the maintainer, since not
everybody follows r-help closely...
I want to perform a posthoc test for my ANCOVA:
a1<-aov(seeds~treatment*length)
With
summary(glht(a1, linfct = mcp(tr
On Tue, 9 Jan 2007, Bálint Czúcz wrote:
There is an improved version of the original random forest algorithm
available in the "party" package (you can find some additional
information on the details here:
http://www.stat.uni-muenchen.de/sfb386/papers/dsp/paper490.pdf ).
I do not know whether i
On Wed, 8 Nov 2006, Baoqiang Cao wrote:
> Dear All,
>
> I am using wilcox.test to test two samples, data_a and data_b, earch sample
> has 3 replicates, suppose data_a and data_b are 20*3 matrix. Then I used the
> following to test the null hypothesis (they are from same distribution.):
>
> wilc
On Fri, 3 Nov 2006, Hvidberg, Martin wrote:
> Dear list
>
> - Propoerly a beginner question, so bare with my frustration...
>
> I tried install the 'shapefiles' package into R 2.4.0 but it seems that the
> install had little effect...
>
>
> > install.packages(c("shapefiles"))
>
On Fri, 27 Oct 2006, Inman, Brant A. M.D. wrote:
> System: R 2.3.1 on a Windows XP computer.
>
> I am validating several cancer prognostic models that have been
> published with a large independent dataset. Some of the models report a
> probability of survival at a specified timepoint, usually
anything?
>
just wait for version 0.991-1 being available for windows (should happen
within the next hours) and have a look at ?glht
Best wishes,
Torsten
> Thanks again for your help.
>
> Regards,
> Ezhil
>
> --- Torsten Hothorn
> <[EMAIL PROTECTED]> wrote:
>
On Tue, 24 Oct 2006, A Ezhil wrote:
> Hi All,
>
> I have performed a 3-way ANOVA analysis for my
> experimental data using aov function. My simple R
> funtion for this is:
>
> 3aof <- function(x){
> m <- data.frame(R,S,T, x);
> anova(aov(x ~ R+S+T+R*S+R*T+S*T+R*S*T, m) )
> }
>
you can use
R> li
ous Inference for General Linear Hypotheses
Version: 0.991-1
Date: 2006-10-16
Author: Torsten Hothorn, Frank Bretz and Peter Westfall
Description: Simultaneous tests and confidence intervals
for general linear hypotheses in parametric models, including
linear, generalized linear, linear mixed eff
On Wed, 18 Oct 2006, Douglas Bates wrote:
On 10/18/06, Torsten Hothorn <[EMAIL PROTECTED]> wrote:
On Wed, 18 Oct 2006, SCHMERA Dénes wrote:
Dear All,
I would like to run a generalized linear mixed model with the software R (one
categorical predictor, one random factor, the distri
On Wed, 18 Oct 2006, SCHMERA D?nes wrote:
Dear All,
I would like to run a generalized linear mixed model with the software R (one
categorical predictor, one random factor, the distribution of the dependent
variable is binomial, and the link is logit). Thereafter, I would like to
perform mul
On Thu, 5 Oct 2006, Christos Hatzis wrote:
> Jue,
>
> On a second look, it appears that wilcox.test does report the
> offset-adjusted statistic U, as also mentioned in the help page.
>
> wilcox.test returns W=6 (instead of 12 as your example showed, unless
> "wilcox_test" is a different function)
On Tue, 19 Sep 2006, Andrew Zachary wrote:
> Not sure if anyone has posted on this problem ... I want to use rpart to
> build a binary tree on a relatively large dataset with ~1400 data points
> and 15 predictors. But I've noticed that rpart fails almost immediately
> in the call to C_s_to_rp, as
[...]
>
>> cochranq.test(K)
>
>Cochran's Q Test for Dependent Samples
>
> data: K
> Cochran's Q = 23.9298, df = 11, p-value = 0.01303
>
Cochran's Q fits into the `coin' framework and thus:
> K <- as.table(matrix(c(1,1,0,0, 1,1,0,1, 1,1,1,1, 1,1,1,1, 1,0,0,0,
+ 1,1,1,1, 1,1,1,1, 0,0,0
On Fri, 25 Aug 2006, Stefan Grosse wrote:
> Dear List,
>
Stefan,
> after updating the exactRanksumTests package I receive a warning that
> the package is not developed any further and that one should consider
> the coin package.
>
> I don't find the signed rank test in the coin package, only th
|___|
On Wed, 2 Aug 2006, Spencer Graves wrote:
> I'm not familiar with the correlation adjustment to Bonferroni you
> mention below, though it sounds interesting. However, I think there is
> something not right about
On Thu, 20 Jul 2006, Jacob van Wyk wrote:
> Hallo
>
> Is there an elegant way to do the following:
>
> Dataset consists of 2 variables: var1: some measurements, and var2: a
> grouping variable with two values, 1 and 2.
> There are (say) 10 measurements from group 1 and 15 measurements from group
On Fri, 14 Jul 2006, Daniele Medri wrote:
> Dear R.Users,
>
> using ctree() (from "party" library) on a data.frame, I want to append a
> column with the references for the groups/segments detected. While these
> nodes are easy readable in output, I need a vector for my obs.
>
Daniele,
do you me
On Mon, 19 Jun 2006, Bálint Czúcz wrote:
Dear R users!
Does someone know about any algorithms / packages in R, that perform
classification / regression / decision trees using multivariate
splits?
I have done some research, but I found nothing. Packages "tree" and
"rpart" seem only to be able
On Thu, 8 Jun 2006, Hans Gardfjell wrote:
Dear R-users,
Just like other users (as seen from previous posts on the list), I would
like to use female and male signs in plots. I found B. Ripley's post
about using Unicode characters. However, it doesn't works for me.
> text(locator(1),"\u2640")
On Sat, 27 May 2006, P Ehlers wrote:
> Looks like a bug (in mvt()?). Note that
>
> pmvnorm(lower = c(-Inf, -Inf), upper = c(Inf, 10))
>
> works as expected, as does replacing any of the 4 Infs with a finite
> value.
>
> Function sadmvn() in pkg:mnormt does give 1 with Inf lower/uppers.
>
yes, lo
Dear useRs,
we would like to inform you that we are approaching the regular
registration deadline (2006-04-15) of useR! 2006, the second R
user conference taking place in Vienna on June 15-17. So if you
plan to join the conference, this is the perfect time to register!
The pre-conference tutoria
Dear useRs,
Version 0.8-1 of the `party' package will appear on CRAN and its mirrors
in due course. This version implements two new tools:
o `mob', an object-oriented implementation of a recently suggested
algorithm for model-based recursive partitioning (Zeileis, Hothorn,
Horn
ract submission, registration and Vienna is
available at
http://www.R-project.org/useR-2006/
We look forward to receiving your abstract!
The organizing committee:
Torsten Hothorn, Achim Zeileis, David Meyer, Bettina Gruen,
Kurt Hornik and Friedrich Leisch
___
On Wed, 15 Feb 2006, [EMAIL PROTECTED] wrote:
> Hi all,
>
> I have being using wilcox.test to test for differences between 2 independent
> samples. I had understood the difference in location to be conventionally the
> difference in the sample medians however this is not the case when implemente
cal details are in
@article{party2006,
key = {566},
author = {Torsten Hothorn and Kurt Hornik and Achim Zeileis},
title = {Unbiased Recursive Partitioning: A Conditional Inference
Framework},
journal = {Journal of Computational and Graphical Statistics},
year = 2006,
tions of R. A web page offering more information on the `useR!'
conference, abstract submission, registration and Vienna is available at
http://www.R-project.org/useR-2006/
We will accept submissions until February 28, 2006.
Let the contributions roll in!
The organizing committee:
Tor
On Mon, 23 Jan 2006, Stephen wrote:
> Hi.,
>
> I would like to conduct a CHAID tree analysis - Chi-square Automatic
> Interaction Detector.
>
> >From what I can make out from MASS, tree, rpart, and a quick search it
> isn't available
> as a package.
>
> Is that correct or have I missed it?
corre
On Wed, 21 Dec 2005, Peter Dalgaard wrote:
> P Ehlers <[EMAIL PROTECTED]> writes:
>
> > > so a good guess at its definition is that it is obtained from W or one
> > > of the others by subtracting the mean and dividing with the SD.
> > >
> >
> > With the SD adjusted for ties, of course. (See, e.g.
On Thu, 28 Jul 2005, Eric Lecoutre wrote:
> Hi there,
>
> I often do receive some mails about this piece of code regarding
> Cochran-Armitage or Mantel Chi square.
This is a very late reply but maybe still interesting. The
conditional version of the Cochran-Armitage test for trend
for proportion
On Tue, 15 Nov 2005, Thomas Lumley wrote:
> On Tue, 15 Nov 2005, Claus Atzenbeck wrote:
>
> > On Tue, 15 Nov 2005, Peter Dalgaard wrote:
> >
> However, for Wilcoxon tests, the formula for effect sizes is:
> r = Z / sqrt(N)
>
> I wonder how I can calculate the Z-score in R for
On Wed, 9 Nov 2005, Julie Lejeune wrote:
> Hello,
> I would like to know what I do it to test correlation between ordered
> categorical variables.
> Tendancy Chi test?
a linear-by-linear association test is one possibilty, see `lbl_test' in
package `coin'.
Best,
Torsten
> Thank you,
> Melle
We are happy to inform you that the online abstract submission and
registration for `useR! 2006' is now available online from
http://www.R-project.org/useR-2006/
This second world meeting of the R user community will take place at
the Wirtschaftsuniversitaet Wien, Vienna, Austria, June 15 to 17
eRs. Please send submissions to Uwe Ligges <[EMAIL PROTECTED]>.
The deadline for submissions is
October, 8th, 2005
Keep the contributions rolling in!
The Editorial Board,
Doug Bates, Paul Murrell and Torsten Hothorn
__
R-help@stat.math.ethz
On Mon, 29 Aug 2005, Peter Dalgaard wrote:
> Torsten Hothorn <[EMAIL PROTECTED]> writes:
>
> > > Dear All,
> > >
> > > is there a stratified version of the Wilcoxon test (also known as van
> > > Elteren test) available in R?
> >
> > y
> Dear All,
>
> is there a stratified version of the Wilcoxon test (also known as van
> Elteren test) available in R?
you can plug it together using the `coin' infrastructure (see the
examples in the manual and vignette).
Torsten
> Thanks,
>
> Heinz Tüchler
>
> _
> Dear list,
>
> Please forgive me if this is a dumb question. I want to compare a set of
> paired data (pre and late). There are ties in pre and late. I searched
> on line, some document says use wilcox.exact in R instead of
> wilcox.test. Anyone has any experience using wilcoxon signed rank test
> I am trying to use a custom S4 object in my C code and I cannot get the
> access to its slots working.
>
>
>
> The following is a toy example, but even this crashes.
>
>
>
> In R I have:
>
>
>
> setClass("pd", representation(data="numeric"))
>
> x <- new("pd", data=1:5)
>
>
>
> test <- function(
On Thu, 4 Aug 2005, Peter Ho wrote:
> HI R-users,
>
> I am trying to repeat an example from Rayner and Best "A contingency
> table approach to nonparametric testing (Chapter 7, Ice cream example).
>
> In their book they calculate Durbin's statistic, D1, a dispersion
> statistics, D2, and a residu
f tutorial proposals: 2005-09-30
- early registration deadline: 2006-01-31
- submission deadline of abstracts: 2006-02-28
- registration deadline: 2006-05-31
We hope to meet you in Vienna!
The organizing committee:
Torsten Hothorn, Achim Zeileis, David Meyer, Betti
> Hi list,
>
> Does anybody know if R has functions to do the ANOVA permutation
> test? I googled and found R has the "vegan" package to do "ANOVA like
> permutation test for Constrained Correspondence Analysis". But does R
> have a function for general ANOVA-like permutation tests? Thanks in
regression models is available.
Best,
Torsten
__
Package: party
Title: A Laboratory for Recursive Part(y)itioning
Date: $Date: 2005/06/27 06:38:16 $
Version: 0.2-2
Author: Torsten Hothorn, Kurt Hornik and Achim Zeileis
On Thu, 2 Jun 2005, Liaw, Andy wrote:
> By post hoc I guess you meant pairwise comparisons. You might want to check
> out the npmc package on CRAN.
>
?oneway_test in `coin' has an example, essentially
### Length of YOY Gizzard Shad from Kokosing Lake, Ohio,
### sampled in Summer 1984, Hollande
__
Package: coin
Title: Conditional Inference Procedures in a Permutation Test Framework
Date: $Date: 2005/06/02 14:55:45 $
Version: 0.2-11
Author: Torsten Hothorn and Kurt Hornik, with contributions by
Mark van de Wiel and Achim Zeileis
Maintainer: Torsten Hothorn <[EM
> I am trying to do a Nonparametric Tukey-type multiple comparison
> post-hoc test to determine which groups are significantly different.
the manual page to `oneway_test' (in package `coin') has an example of
the Nemenyi-Damico-Wolfe-Dunn test taken from Hollander & Wolfe (1999).
Best,
Torsten
send submissions to Uwe Ligges <[EMAIL PROTECTED]>.
The deadline for submissions is
April, 10th, 2005
Keep the contributions rolling in!
The Editorial Board,
Doug Bates, Paul Murrell and Torsten Hothorn
__
R-help@stat.math.ethz
On Wed, 22 Dec 2004, Anne wrote:
> Sorry, sorry
> of course
> levels(testf)[c(2,1,3)]
>
> will do the job
>
Anne,
note that `red' is coded as `blue' now:
R> testf <- factor(c("red","red","red","blue","blue","white"))
R> testf
[1] red red red blue blue white
Levels: blue red white
R>
On Wed, 24 Nov 2004, Ramon Diaz-Uriarte wrote:
> Dear Cristoph,
>
> I guess you want to assess the error rate of a LDA that has been fitted to a
> set of currently existing training data, and that in the future you will get
> some new observation(s) for which you want to make a prediction.
> Then
ammeur
>
>
>
> -Message d'origine-
> De : Torsten Hothorn [mailto:[EMAIL PROTECTED]
> Envoyé : vendredi 13 août 2004 14:19
> À : Liaw, Andy; Laurent Houdusse
> Cc : '[EMAIL PROTECTED]'
> Objet : RE: [R] simtest for Dunnett's test
>
>
>
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
On Wed, 4 Aug 2004, Matt Nunes wrote:
> hello.
>
> I have a query about the Matrix package for R. I wrote some code a while
> ago using the Matrix package version 1.6.2 with an early version of R, to
> do some linear least squares for regression:
>
> xn
> [,1] [,2] [,3] [,4]
On Thu, 15 Jul 2004, [big5] house-ball wrote:
> Hi,
>
> I try to compute critical value for multivariate normal distribution, and I find the
> crit(fit, const = c(0, 1), d = 1, cov = 0.95, rdf = 0) which seems to compute
> critical value. However, I can't compute right critical value for multiva
On Fri, 25 Jun 2004, Douglas Grove wrote:
> I should have specified an additional constraint:
>
> I'm going to need to use this repeatedly on large
> vectors (length 10^6), so something efficient is
> needed.
>
give function `irank' in package `exactRankTests' a try.
Best,
Torsten
>
> On Fri
On Wed, 16 Jun 2004, Mahbub Latif wrote:
> Hi,
>
> I am looking for C/C++ codes for computing generalized
> inverse of a matrix. Can anyone help me in this
> regard?
>
experimental, undocumented and for square matrices only:
SEXP svd (SEXP x) {
/* experimental */
SEXP jobu, jobv, u, v,
On Wed, 26 May 2004, Chihiro Kuroki wrote:
>
>
> BTW, I have another strange example of simtest. I want to know
> why simtest returns these p-values.
>
> -- example 1 ---
> rm(list = ls())
> require(multcomp)
> y1 <- c(seq(3,7),seq(3,7))
> y2 <- c(rep(c(6,7,8,9),7))
> Hello,
>
> I am trying to get the P values from the output of a summary for lm.
>
> lm <- lm(y ~ age + sex)
> s <- summary(lm)
s$coefficients
gives you a matrix with the P-values in the fourth column
Torsten
>
> I thought that I might be able to get them using a combination of scan,
> grep a
>
> BTW, I have another strange example of simtest. I want to know
> why simtest returns these p-values.
>
> -- example 1 ---
> rm(list = ls())
> require(multcomp)
> y1 <- c(seq(3,7),seq(3,7))
> y2 <- c(rep(c(6,7,8,9),7))
> sort(runif(28),index=T) -> a
> y3 <- numeric(0)
On Thu, 20 May 2004, Chihiro Kuroki wrote:
> Hi, all:
>
> Two examples are shown below.
>
> I want to use the multiple comparison of Dunnett.
> It succeeded in upper case "example 1".
>
> However, the lower case "example 2" went wrong.
>
it was due to an error in the code underlying the `mvtnorm'
Yes, `pmvt' returns NaN without indicating this error. We need to check.
Thanks for the report (and *please* cc emails reporting problems with
packages to the maintainer!),
Torsten
On Thu, 20 May 2004, Chihiro Kuroki wrote:
> Hi, all:
>
> Two examples are shown below.
>
> I want to use the mult
On Tue, 13 Apr 2004, Hui Han wrote:
> Hi,
>
> I am doing feature selection for my dataset. The following is
> the extreme case where only one feature is left. But I got
> the error below. So my question is that do I have to use
> more than one features?
>
> sample.subset
> udomain.edu hpclass
>
On Sun, 4 Apr 2004, David L. Van Brunt, Ph.D. wrote:
> Playing with randomForest, samples run fine. But on real data, no go.
>
> Here's the setup: OS X, same behavior whether I'm using R-Aqua 1.8.1 or the
> Fink compile-of-my-own with X-11, R version 1.8.1.
>
> This is on OS X 10.3 (aka "Panther")
On Wed, 31 Mar 2004, Hui Han wrote:
> Dear all,
>
> Can anybody give me some hint on the following error msg I got with using
> randomForest?
>
> I have two-class classification problem. The data file "sample" is:
> --
> udomain.edu udomain.
On Sat, 6 Mar 2004, Padmanabhan, Sudharsha wrote:
>
> Hello,
>
> I need help in performing a Van_der_Waerden normal scores test in R. I
> have two arrays of scores(final on therapy scores from drug and placebo) and
> want to use the normal scores procdeure to test for significance.
> (observation
Hi,
I could not find any information on how `read.spss' deals with date
information. As an example, I created a file containing two variables,
one numeric (values = (1, 2)) and one of type "Datum" in SPSS (german
version with values "11.02.2003" and "03.04.1999" and I get in R:
SPSSfile = url("h
On Tue, 17 Feb 2004, Mahmoud K. Okasha wrote:
> Hello R-users,
>
> I would like to generate two-way contingency tables with zero in one cell. I tried
> to use the function r2dtable but I could not force one cell to have zero value.
>
r2dtable samples from the conditional distribution of the tab
> Stuart V Jordan <[EMAIL PROTECTED]> writes:
>
> > > mvrnorm(n = 1000,B,V)
> > Error in mu + eS$vectors %*% diag(sqrt(pmax(ev, 0)), p) %*% t(X) :
> > non-conformable arrays
> > > mvrnorm(n = 1000,t(B),V)
> > Error in mu + eS$vectors %*% diag(sqrt(pmax(ev, 0)), p) %*% t(X) :
> > n
> Dear R-listers,
>
> I am trying to compute simultaneous confidence intervals with simint from the
> package multcomp. 230 measures (abundance) have been taken in 23 sites (factor) of a
> data.frame (donnees: a file can be sent on request, saved with
> save(donnees,file="donnees")). I would lik
We are happy to inform you that the topics of the keynote lectures to be
presented at the first R user conference `useR! 2004' in Vienna (May
20-22th) are now available from the conference web page at
http://www.ci.tuwien.ac.at/Conferences/useR-2004/
The submission process for oral and p
On Thu, 11 Dec 2003, Jeffrey Chang wrote:
> Hello everybody,
>
> I'm running R 1.8.1 on both Linux and OS X compiled with gcc 3.2.2 and
> 3.3, respectively. The following call seems to freeze the interpreter
> on both systems:
> > chisq.test(matrix(c(233, 580104, 3776, 5786104), 2, 2),
> simulat
>
> > Jason,
> >
> > For many bibliography styles, the command \citeyear{key} will work. If this
> > doesn't work for the style you are using, you can investigate style-specific
> > methods or consider other styles. I find that natbib is good for
> > author-year formats.
>
> Yup. \usepackage[roun
> Hello everybody,
>
> I'm seeing some strange behavior on R 1.8.1 on Intel/Linux compiled
> with gcc 3.2.2. The p-value calculated from the chisq.test function is
> incorrect for some input values:
>
>
> > chisq.test(matrix(c(0, 1, 1, 12555), 2, 2), simulate.p.value=TRUE)
>
> Pearson's
> On 28-Nov-03 Torsten Hothorn wrote:
> > yes, thats my understanding too. The "enumerative techniques" as
> > you call it condition on the data actually observed and determine
> > the null distribution of the associated test statistic from the data.
> >
On Thu, 27 Nov 2003, Spencer Graves wrote:
> Hi, Torsten:
>
> Thanks for the reference to library(exactRankTests). That seems
> like a reasonable alternative to "prop.test" with small samples.
>
> However, aren't "exact tests" and the related bootstrap
> methodology what Deming calle
> Hello,
>
> I'm looking for some guidance with the following problem:
>
> I've 2 samples A (111 items) and B (10 items) drawn from the same unknown
> population. Witihn A I find 9 "positives" and in B 0 positives. I'd like to
> know if the 2 samples A and B are different, ie is there a way to fin
> Hi,
>
> is it correct that i need ~ 2GB RAM that it's
> possible to work with the default setting
> ntree=500 and a data.frame with 100.000 rows
> and max. 10 columns for training and testing?
>
no. You may parallelize the computations: perform 5 runs of RF with `ntree
= 100' (or less) and sav
> Hello,
> Is there like the kruskal wallis test in relation to ANOVA (no
> restrictions on normallity and variance homogenity) something (in R)
> for MANOVA?
You may know Fortunato Pesarin's book
@book{multivaria:2001,
key = {386},
author= {Fortunato Pesarin},
title = {Mu
ference, abstract submission, registration and Vienna is available at
http://www.ci.tuwien.ac.at/Conferences/useR-2004/
We hope to meet you in Vienna!
For the organizing committee:
Torsten Hothorn, Achim Zeileis and David Meyer.
___
[EMAIL P
On Mon, 10 Nov 2003, Christoph Lehmann wrote:
> Hi
> Since a stepwise procedure for variable selection (as e.g. in SPSS) for
> a LDA is not implemented in R and anyway I cannot be sure, that all the
> required assumptions for e.g. a procedure using a statistic based on
> wilks' lambda, hold (such
> I've never seen anything written about multiple comparisons,
> as in the multcomp package or with TukeyHSD, but using a glm.
> Do such procedures exist? Are they sensible?
> Are there any packages in R that implement such comparisons?
since version 0.4-0 in `multcomp':
0.4-0 (13.08.2003)
> 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
On Mon, 20 Oct 2003, michael watson (IAH-C) wrote:
> Hi
>
> I wondered what the standard procedure was for upgrading from one version of R to
> the next. I currently have R 1.7.1 and want the latest release, R 1.8.0. I am
> running SUSE linux 8.2. The main thing is that I want to keep all of
are given at
http://www.bioinf.uni-hannover.de/R-Kurse/
Best,
Torsten
___
| |
| Dr. rer. nat. Torsten Hothorn
On Wed, 10 Sep 2003 [EMAIL PROTECTED] wrote:
> Hello, there:
> I got data matix with missing values. I want to calculate any possible
> pairwise Spearman correlation rho for each column. Is there a function just
> like cor(x, y, use="complete.obs") for Pearson correlation?
> Thanks in advance!
>
On 26 Aug 2003, Stefan [ISO-8859-1] Böhringer wrote:
> How can I extract the linear discriminant functions resulting from a LDA
> analysis?
>
> The coefficients are listed as a result from the analysis but I have not
> found a way to extract these programmatically. No refrences in the
> archives
On Thu, 21 Aug 2003, Gavrilov, Pavel M wrote:
> Hello. I have been working with GeoDA, and have created a spatial weights
> file for my data. I am now looking to use R to run regressions on this
> data. However, I don't know and can't figure out how to get my data into R
> to run these regressi
> Have I been sleeping in class?
>
> rw1071 from CRAN, windows XP
>
> incidencia is made by a call to tapply
>
> > class(incidencia)
> [1] "array"
> > incidencia <- unclass(incidencia)
> > class(incidencia)
> [1] "array"
>
>
`unclass' only removes the `class' attribute. However, arrays do no
Version 0.4-0 of the `multcomp' package has just been uploaded to CRAN.
Known bugs have been fixed and, inspired by some questions posted on
r-help, methods for `lm' and `glm' objects for the generics `simint' and
`simtest' have been implemented. Thanks to all bug-reporters!
Torsten
>From the C
On Tue, 12 Aug 2003, Ciril Rozman wrote:
> Sirs,
> I have recently been interested in your Maxstat.
For questions related to add-on packages please at least cc to the
maintainer since some might not read r-help regulary.
> I have computed
> with my own programme the ranks (by using the Kaplan-
On Thu, 14 Aug 2003, Spencer Graves wrote:
> This seems to identify a possible bug in R 1.7.1 under Windows 2000:
>
> > tstDf <- data.frame(y = 1:11, x=1:11)
> > fit <- nls(y~a/x, data=tstDf, start=list(a=1))
> > predict(fit, se.fit=TRUE)
> [1] 7.0601879 3.5300939 2.3533960 1.7650470 1.4120
> Why do
> x<-b%*%ginv(A)
> and
> x<-solve(A,b)
> give different results?.
they do (in cases the solution to A x = b is unique):
R> A <- matrix(c(0,-4,4,0),nrow=2,ncol=2)
R> A
[,1] [,2]
[1,]04
[2,] -40
R> b <- c(-16,0)
R> x1 <- ginv(A) %*% b <- NOT b %*% ginv(A)
R> x1
ula
X5 ~ .
R> rpart(thisformula, data = df.treino)
n= 5
node), split, n, deviance, yval
* denotes terminal node
1) root 5 3.032904 -0.3392065 *
Torsten
> -Original Message-
> From: Torsten Hothorn [mailto:[EMAIL PROTECTED]
> Sent: Tue
On Tue, 22 Jul 2003, Luis Miguel Almeida da Silva wrote:
> Dear helpers
>
> I want to use rpart several times in a loop to build a classification tree. My
> problem is that rpart needs a formula as argument and for that the variables need to
> have names and this doesn't happen in my case. Ev
On Thu, 17 Jul 2003, Angel wrote:
>
> Hi,
> How can I remove elements from a vector and them put them back in place??
> An example (very simple, my vector/operations are much larger/complicated):
> Got a vector, lets say:
> a<-c(1,2,40,10,3,20,6);
you need to save the index of the elements o
>
> Hi,
>
> I am trying to get an error estimation
> for a classification done using lda.
>
> The examples work fine, however I don't get
> my own code to work.
>
> The data is in object d
>
> > d
> class hydrophobicitycharge geometry
> 1 2 6490.0400 1434.97
> This question concerns rpart's facility for user-defined functions that
> accomplish splitting.
>
> I was interested in modifying the code so that in each terminal node,
> a linear regression is fit to the data.
>
if you just want to have a linear model in each terminal node instead of
the mea
> Dear collegues,
>
> Using maxstat I am getting the following:
>
> > blood <- maxstat.test(Surv(SUPER, FV)~ZAP,data=zap70, smethod="LogRank")
> Error in maxstat(y = structure(c(24.4301369863014, 26.4164383561644,
> 18.7835616438356, : couldn't find function "cscores"
>
`cscores' is provided by
> I have the following problem. It is not of earthshaking importance,
> but still I have spent a considerable amount of time thinking about
> it.
>
> PROBLEM: Is there any way I can have a single textfile that contains
> both
>
> a) data
>
> b) programcode
>
> The program should act on the dat
1 - 100 of 118 matches
Mail list logo