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, D
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 par
I am desperately trying to impute missing data using 'imp.mix' but always
run into this yucky error message to which I cannot find the solution. It's
the first time I am using mix and I'm trying really hard to understand, but
there's just this one step I don't get...perhaps someone knows the answer
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
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
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
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"
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 heatm
> 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
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 indiv
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(x<5)
[1] 1 2 3 4
> which(x>2)
[1] 3 4 5 6 7 8 9 10
> which(x<5 && x>2)
numeric(0)
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
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 r
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(x<5)
> [1] 1 2 3 4
> > whic
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 to
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
http://biostat.
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 log
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 2
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 g
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:
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, an
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
https://stat.ethz.ch/ma
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 w
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 th
Ok, the time consumption could be a problem, I find it quite okay though.
My problem is to reach the actual messages.
If I use the option File/Attechments/Open, the window that opens just
contain a list of all headers from the last 24 hours (from which I just
choosed one - as I thought)
If I on
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 utf-
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 analys
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 "
> "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
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 latte
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 (
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?
> __
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 id
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! http://www.R-projec
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
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
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 ~40,0
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 cha
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 treatmen
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 o
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 e
[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 manipulat
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 best way to
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 somet
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,
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
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
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 Na/N
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
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
https://stat.e
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 charac
51 matches
Mail list logo
