Dear Prof. Ripley,
thanks for your help. Everything is working fine!
With best regards
Michael Wolf
-Ursprüngliche Nachricht-
Von: Brian D Ripley [mailto:[EMAIL PROTECTED]
Gesendet: Mittwoch, 15. Dezember 2004 09:01
An: Wolf, Michael
Betreff: Re: [R] Problem with German special
the method is described in the article:goodman leo A.,a modified multiple
regression approch to analysis of dischotomous variables,american sociological
review 33(hebruary):28-46
thank you in advance:)
__
[EMAIL PROTECTED] mailing list
Hi Corey
Did you try to use par(new=TRUE) before ploting the second
persp graph?
Cheers
Petr
On 15 Dec 2004 at 10:44, Corey Bradshaw wrote:
I've created a perspective plot using 'persp' in the graphics package.
I'd like to add a second plane of z values to the existing plot, but I
cannot
Hi all,
I need to fit a trend surface to a 3D-dataset with the surf.gls() function
but my data contain missing values. There is just 2 NA in a total of 360
data points but I can't suppresse the variable associated with these NA
values to the analysis. The problem is surf.gls can't manage with
Hello,
Does someone has experience working with time-series where the ocurrencies do
not happen at regular time spaces? In one day I can have 20 ocurrencies and in
the following day I can have the double. ts datatype must have regular time
spaces. What is the best way to put the data if I want
Hello all!
I recently switched mail program to Mozilla Thunderbird, running on W2k.
Everything works fine, except the MIME-digests from this list. The
decoding doesn´t work properly.
I had contact with Martin Maechler some time ago, and he suggested a try
on this list for ideas on how to do the
I think this will do what you want:
# Need this function to remove spaces from term labels later on
removeSpace - function(string) gsub([[:space:]], , string)
# Specify which terms are in a tvar group
# (could remove spaces separately)
tvar - unname(sapply(c(x:A, z, B, poly(v,3)),
Hi,
I have ~40,000 rows in a database, each of which contains an id column and
20 additional columns of count data.
I want to cluster the rows based on these count vectors.
Their are ~1.6 billion possible 'distances' between pairs of vectors
(cells in my distance matrix), so I need to do
jmoreira at fe.up.pt writes:
: Does someone has experience working with time-series where the ocurrencies
do
: not happen at regular time spaces? In one day I can have 20 ocurrencies and
in
: the following day I can have the double. ts datatype must have regular time
: spaces. What is the
Dear Dan,
I would think about transforming your columns in such a way (square
root, log?) that methods operating on n*p matrices and assuming
roughly elliptical within-clusters distributions such as kmeans or
clara, or, after dimension reduction, EMclust or fixmahal can be applied.
Maybe you can
On Wed, 15 Dec 2004, Diethelm Wuertz wrote:
[EMAIL PROTECTED] wrote:
Hello,
Does someone has experience working with time-series where the ocurrencies do
not happen at regular time spaces?
package fBasics from Rmetrics (www.rmetrics.org) has S4 timeDate and
timeSeries objects similar
I did not find this in the archive (hope it isn't there...):
the current release of R (2.0.1) for MacOS (10.3.6) seems not to handle
german special characters like '' correctly:
f - ''
can be entered at the prompt, but echoing the variable yields
[1] \303\274 (I think the unicode of the
Kerry Bush wrote:
I tried lrm in library(Design) but there is always
some error message. Is this function really doing the
weighted logistic regression as maximizing the
following likelihood:
\sum w_i*(y_i*\beta*x_i-log(1+exp(\beta*x_i)))
Does anybody know a better way to fit this kind of
model in
CG Pettersson wrote:
Hello all!
I recently switched mail program to Mozilla Thunderbird, running on W2k.
Everything works fine, except the MIME-digests from this list. The
decoding doesn´t work properly.
I had contact with Martin Maechler some time ago, and he suggested a try
on this list for
There's also nlm for non-linear newton type minimization,
best, ingmar
On 12/15/04 3:48 PM, Clark Allan [EMAIL PROTECTED] wrote:
hi all
other than optim, optimise, and some other related optimisation
functions are there any optimisation packages in R?
Dear ronggui,
If memory serves me right, this is a relatively early presentation for
sociologists of loglinear models for contingency tables (in which all the
variables are dichotomous). You can fit such a model in R as a generalised
linear model using glm() or via the loglin() function. The
On Wed, 15 Dec 2004, BUHARD O wrote:
Hi all,
I need to fit a trend surface to a 3D-dataset with the surf.gls()
function but my data contain missing values. There is just 2 NA in a
total of 360 data points but I can't suppresse the variable associated
with these NA values to the analysis.
Brian D Ripley wrote:
You wrote your mail in UTF-8. R does not support UTF-8, and that is both
documented and announced on startup in such a locale (at least on OSes
with standard-conforming implementations):
thanks for clarifying this point.
nevertheless:
1. the mail was (on purpose) sent in
When I am running R, the backspace key deletes the whole word instead of
one character.
Before I start R, and after I exit R, the backspace key works as it
should in my xterm terminal window.
This gets really annoying since I make a lot of mistakes while typing
and don't always want to retype the
Turns out that this problem is a propagation of a bug from a
Depends package - gdata doesn't import stats in its NAMESPACE
file, although it makes use of objects defined in stats like
reorder and na.omit.
I'm curious as to the reason why 'R CMD check' did not detect
this at the source package,
My guess is `area under the ROC curve'. There's the roc package in
BioConductor that I believe can compute this.
Andy
From: Spencer Graves
What's AUC? If you mean AIC (Akaike Information
Criterion), and
if you fit logistic regression using glm, the help file
says that glm
At 17:07 15/12/2004, Spencer Graves wrote:
Dear R-helper,
I would like to compare the AUC of two logistic regression models (same
population). Is it possible with R ?
Thank you
Roman Rouzier
Roman
If I understand your question You have 2 ROC curve from same dataset. In
this case you can use a
Joe Nocera wrote:
I believe that Roman is referring to AUC as the Area Under Curve from a Receiver
Operating Characteristic.
If this indeed your quantity of interest - it can be calculated in R. You can
download
code at:
http://www.bioconductor.org/repository/release1.5/package/Win32/
and/or
Why the last which in the following example doesn't
work? Is there a simple way to identify the indices of
array elements that meet multiple criteria?
Thanks.
YC Tao
x-1:10
which(x5)
[1] 1 2 3 4
which(x2)
[1] 3 4 5 6 7 8 9 10
which(x5 x2)
numeric(0)
which(x 5 x 2)
[1] 3 4
From the help page for Logical Operators:
?!
'' and '' indicate logical AND and '|' and '||' indicate
logical OR. The shorter form performs elementwise comparisons in
much the same way as arithmetic operators. The longer form
evaluates left to right
Dear R gurus,
I have the following model:
appcov.aov - aov(yield ~ prevyield + trt + block)
where prevyield is a continuous numeric covariate and trt and block are
factors (yes, I did factor()!)
Now, when I do a TukeyHSD, my diff's are all screwed up!
For instance:
treatment mean for treatmen E
Here's what I get:
summary(lm(yield ~ prevyield + trt + block))
Call:
lm(formula = yield ~ prevyield + trt + block)
Residuals:
Min 1Q Median 3Q Max
-22.616 -9.254 2.051 10.687 19.421
Coefficients:
Estimate Std. Error t value Pr(|t|)
(Intercept) 38.5797
Hi.All and R developers:
When I look into the R source code, I have a question.Since R has its
own data structure(i.e. SEXP),how does it convert the result to the normal
output after it has computed? For example,when I input,
abs(-3)
I learned that in R's execution, the expression is
[EMAIL PROTECTED] wrote:
Hello,
Does someone has experience working with time-series where the ocurrencies do
not happen at regular time spaces?
package fBasics from Rmetrics (www.rmetrics.org) has S4 timeDate and
timeSeries objects similar
to those in SPlus for irregular time series
I'm pretty new to R and I have a stats problem that's got me baffled. I'm
trying to carry out a repeated measures test on Poission distributed data,
where half the subjects in each block (4 blocks) were treated, and half were
controls. Measurements were carried out before and after the
It sounds like clara in package cluster might help.
Regards,
Matt Wiener
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Dan Bolser
Sent: Wednesday, December 15, 2004 6:37 AM
To: R mailing list
Subject: [R] Massive clustering job?
Hi,
I have
On Wednesday 15 December 2004 01:10, Corey Bradshaw wrote:
Thanks, Romain. I've certainly used that to draw lines and points in
the plots produced by 'persp'; however, my problem is that I need to
incorporate an entirely new z function (not just a plane) onto the
same plot (i.e., using the
hi all
other than optim, optimise, and some other related optimisation
functions are there any optimisation packages in R?__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
You wrote your mail in UTF-8. R does not support UTF-8, and that is both
documented and announced on startup in such a locale (at least on OSes
with standard-conforming implementations):
gannet% env LANG=en_GB.utf8 R
R : Copyright 2004, The R Foundation for Statistical Computing
Version 2.0.1
CG == CG Pettersson [EMAIL PROTECTED]
on Wed, 15 Dec 2004 16:43:38 +0100 writes:
CG Ok, the time consumption could be a problem, I find it quite okay
though.
CG My problem is to reach the actual messages.
CG If I use the option File/Attechments/Open, the window that opens just
You don't say what sort of optimisation you have in mind. If you
are looking for something that will handle non-standard problems,
you can have a look at 'genopt' from S Poetry.
Patrick Burns
Burns Statistics
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A
I was going to say ``Why not just use glm()?'', but when I tried the
example given in the original message I got a different but similarly
nervous-making warning:
Warning in eval(expr, envir, enclos) : non-integer #successes
in a binomial glm!
Looking into the code I found that
Dear R-helper,
I would like to compare the AUC of two logistic regression models (same
population). Is it possible with R ?
Thank you
Roman Rouzier
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
What's AUC? If you mean AIC (Akaike Information Criterion), and
if you fit logistic regression using glm, the help file says that glm
returns an object of class glm, which is a list containing among other
things an attribute aic. For example, suppose you fit a model as follows:
Hi all,
I encountered a weird problem when using the
Design and Hmisc libraries in S-Plus (it worked well
in R). I have a data frame called b, which
has 3 columns: time, status and
scores, a sample of the data frame is like:
data frame b:
time status scores
1 27 0 -126.7
2 24 0
AUC means Area Under the Curve and is a common summary statistic for
repeated measures experiments (e.g. repeated measurements of serum
concentration of a drug in an individual)in PK/PD studies
(pharmacokinetic/pharmacodynamic).
I think the poster may actually mean nonlinear regression for
I believe that Roman is referring to AUC as the Area Under Curve from a
Receiver
Operating Characteristic.
If this indeed your quantity of interest - it can be calculated in R. You can
download
code at:
http://www.bioconductor.org/repository/release1.5/package/Win32/
and/or
libR seems to include a main() function. Should it?
I'm hoping this simpler question will get more of a response than my
previous post about R a la carte :)
The evidence that libR includes main() is two-fold. First, when I use
nm, I see the _main defined. Second, when I run a program linked
Because you have one too many `'.
Andy
From: Y. C. Tao
Why the last which in the following example doesn't
work? Is there a simple way to identify the indices of
array elements that meet multiple criteria?
Thanks.
YC Tao
x-1:10
which(x5)
[1] 1 2 3 4
which(x2)
[1] 3 4 5
Good afternoon,
I ran heatmap and hclust on the same matrix x (strictly, I ran
heatmap(x), and hclust(dist(t(x))), and realized that the two
dendrograms were slightly different, in that the left-right
arrangement of one pair of subclusters (columns) was reversed in the
two functions (but all
Hierarchical clustering does NOT give an ordering for the clusters. It only
gives the clustering, so order is not invariant.
Sean
- Original Message -
From: Min-Han Tan [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: Wednesday, December 15, 2004 3:55 PM
Subject: [R] hclust and heatmap -
Damián Cirelli [EMAIL PROTECTED] writes:
Dear R gurus,
I have the following model:
appcov.aov - aov(yield ~ prevyield + trt + block)
where prevyield is a continuous numeric covariate and trt and block
are factors (yes, I did factor()!)
Now, when I do a TukeyHSD, my diff's are all
I have interactions where I shouldn't so nevermind, I'm a dumb ass.
Thanks again
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Hi there,
Turns out it was a simply 'par' problem
I just needed to establish par(new = TRUE) after the first 'persp' plot
and prior to coding the 2nd 'persp' plot.
Thanks for your suggestions.
Corey
-Original Message-
From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
Sent: Thursday,
49 matches
Mail list logo