Thomas Lumley <[EMAIL PROTECTED]> writes:
> If you are going to be doing a lot of this sort of thing it might be worth
> looking at the parallel processing facilities in the `snow' package.
> There's a description of their use in another gene expression problem in
> the new R Newsletter.
Actuall
The function 'wilcox.test' in R and S gives (almost) identical results (see
below). 'qwilcox' however, does not:
> qwilcox(p,5,5)
p: 0.025 0.975
R> 3 22
S> 18 37
I originally wanted to ask a questions, but then I found the answer. Given
the confu
On 06/11/03 03:28, Jonck van der Kogel wrote:
>Hi all,
>I have (yet another) question about a function in R. What I would like
>to do is test for the presence of a certain value in a vector, and have
>the positions that this value is at returned to me.
>For example, let's say I have a vector:
>x
On 06/10/03 18:30, Morgan Hough wrote:
>Sorry for the probable repeat post but I can only search the list up to
>2002 (is there a better way?).
Yes, see my search page below.
>I am using the RH9 RPM from CRAN but
>packages like AnalyzeFMRI say that tcltk is not found. Do I need to do
>more to ge
?which
HTH
Edith Hodgen
Statistician and Data Manager
New Zealand Council for Educational Research
Phone: +64-4-384 7939 x 812
Fax: +64-4-384 7933
[EMAIL PROTECTED]
Web site http://www.nzcer.org.nz
>>> Jonck van der Kogel <[EMAIL PROTECTED]> 11/06/2003 13:28:38 >>>
Hi all,
I hav
Sorry for the probable repeat post but I can only search the list up to
2002 (is there a better way?). I am using the RH9 RPM from CRAN but
packages like AnalyzeFMRI say that tcltk is not found. Do I need to do
more to get Tk GUIs working on RH9 or does the RPM not have tcltk support
built in (shou
Hi all,
I have (yet another) question about a function in R. What I would like
to do is test for the presence of a certain value in a vector, and have
the positions that this value is at returned to me.
For example, let's say I have a vector:
x <- c(1,1,2,2,3,3,4,4)
Now I would like a function t
Hi, A, this is G -- please let me know how to reach you!
[[alternate HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
I did the following, and had it plot as needed:
>legend(3.25, -2.5, leg.txt, pch=15, col= levels(as.factor(codes(cv.wshed.grp
Problem was in my (lack of) understanding of codes()
But my second question is this: is there a way to color code the points
on the first two PCA axes by group id wit
On Tue, 10 Jun 2003, Xiao-Jun Ma wrote:
> I have a gene expression matrix n (genes) X p (cases), where n = 8000 and p
> = 100. I want to fit each gene as univariate in a coxph model, i.e., fitting
> 8000 models. I do something like this:
>
> res <- apply(data, 1, coxph.func)
>
> which takes about
I have a gene expression matrix n (genes) X p (cases), where n = 8000 and p
= 100. I want to fit each gene as univariate in a coxph model, i.e., fitting
8000 models. I do something like this:
res <- apply(data, 1, coxph.func)
which takes about 4 min, not bad. But I need to do large numbers of
per
It means what it says! The residuals from the manova fit have a
degenerate distribution: that's a problem with bootstrapping.
I don't think you've done this correctly: assuming you are intending to
bootstrap residuals you seem to have resampled the independent variable
and not added back the m
Hi,
Does anyone know what the error message mean?
> Boot2.Pillai <- function(x, ind) {
+ x <- as.matrix(x[,2:ncol(x)])
+ boot.x <- as.factor(x[ind, 1])
+ boot.man <- manova(x ~ boot.x)
+ summary(manova(boot.man))[[4]][[3]]
+ }
>
> man.res <- manova(as.matrix(pl.nosite) ~
+
Hi all,
I thought that I would follow up on this thread with some "refined"
information.
As I mentioned in a prior post in this thread, I posted a query to the
OOo folks regarding the inability to replicate the IEEE representation
issue in OOo Calc.
Recall the following results:
OOo Calc 1.0.2
I'm using R 1.6.2 on a Windows 2000 machine.
I've plotted the results of an MDS run labeled by a numerical ID, and
color coded by a group code:
plot(cv.mds.spr$points, type="n", main="Non-Metric Multidimensional
Scaling of SprRun CV Watersheds")
text(cv.mds.spr$points, labels = as.character(cv.ws
On Tue, 10 Jun 2003, Rafael Bertola wrote:
> I´ve a data set and i want fit a kernel density estimate to the data.
> but using the k-nearest neighbour method.
> How i do this with R.
Well, you define the exact algorithm you want to use, and then you program
it, R being a fully-featured programmi
I´ve a data set and i want fit a kernel density estimate to the data.
but using the k-nearest neighbour method.
How i do this with R.
thanks
--
[EMAIL PROTECTED]
--
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/
> "MM" == Martin Maechler <[EMAIL PROTECTED]>
> on Tue, 10 Jun 2003 18:12:36 +0200 writes:
MM> Ping, you found another bug in silhouette.default() --
MM> which can happen when there's one cluster with exactly
MM> one observation.
MM> I'll let you know more, once I have
Ping,
you found another bug in silhouette.default() -- which can
happen when there's one cluster with exactly one observation.
I'll let you know more, once I have a complete fix.
Martin Maechler <[EMAIL PROTECTED]> http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum LEO C16
The shapefiles and EMME/2 packages are available at:
http://www.odot.state.or.us/tddtpau/R.html
Rmap looks like a great package but it doesn't work on the OS that I use at
work (Windows NT).
Ben Stabler
>-Original Message-
>From: Hisaji Ono [mailto:[EMAIL PROTECTED]
>Sent: Tuesday, June
An exponential distribution is a gamma distribution, and as far as fitting
the MLE of the coefficients all gammas give the same MLEs. (You can
specify the dispersion and hence that the gamma is exponential when asking
for summaries, anova, etc.)
On Wed, 11 Jun 2003, Masayoshi Hayashi wrote:
> I
On Tuesday 10 June 2003 17:31, Masayoshi Hayashi wrote:
(B> I am learning glm function, but how do you fit data using exponential
(B> distribution with glm?
(B
(BThe Gamma family is parametrised in glm() by two parameters:
(Bmean and dispersion; the "dispersion" regulates the shape.
(B
(BS
I am learning glm function, but how do you fit data using exponential
(Bdistribution with glm?
(BIn the help file, under "Family Objects for Models", no ready made option
(Bseems available for the distribution as well as for other distributions
(Bsatisfying GLM requirements not listed there.
(
> Hello to all-
> 1. When I run a regression which implements the
> augmented Dickey-Fuller
> test, I am confused about the names given to the regressors
> in the output.
> I understand what "xGE" stands for in a standard "lm" test
> involving an
> independent variable GE for instance
First, *which* SOM routine, from which package?
Credit where credit is due and all that
If you mean the one in class(VR), it uses the standard R random number
generation. You don't need to `pass the random-seed as an argument
somehow', just set the seed as you would anyway (e.g. via set.seed
Hi all,
I have a question about the SOM routine. You can either supply the
initial representatives for the lattice yourself or else they are
chosen randomly from the dataset. Is it possible to pass the
random-seed as an argument somehow, when choosing the random
initialisation of the lattice?
A
Hello to all-
1. When I run a regression which implements the augmented Dickey-Fuller
test, I am confused about the names given to the regressors in the output.
I understand what "xGE" stands for in a standard "lm" test involving an
independent variable GE for instance, but if I lags a
Hi.
Currently I couldn't your uploaded packages in CRAN at Austria.
Where could I get them currently?
And do you know Rmap(http://www.maths.lancs.ac.uk/Software/Rmap/) which
supports shape files and other geo-spatial data formats with map projection?
Regards.
- Original Message -
I have been writing some S4 classes and have a problem about how I might
pass a signature to "c()".
Take the following example:
setClass("collection", representation("list", date="POSIXt"))
x <- new("collection", list(1,2,3), date=Sys.time())
y <- new("collection", list(4,5,6), date=Sys.time())
Hello
Im a beginner on R and I would like how to develop a ROC statistic to
evaluate a GAM model with a binomial distribution (Im using mgcv package)
Thanks in advance
--
David Nogués Bravo
Functional Ecology and Biodiversity Department
Pyrenean Institute of Ecology
Spanish Research Council
Av
Eigenvectors are defined only up to a scalar constant (assuming distinct
eigenvalues). However, your `by hand' answer does not pass the simple
test Av = lambda v for some lambda. So you cannot reproduce incorrect
answers in R!
Your example is unusual: A is of rank 1.
On 9 Jun 2003, Christoph
31 matches
Mail list logo