Sam McClatchie noaa.gov> writes:
> Can anyone tell me how to produce vector graphics in R that can be
> ungrouped for editing in CorelDraw?
Postscript should do. In theory, EMF will do too, but it quirky (at least with
my version of CorelDraw).
Dieter
Hi Amna,
I believe you are looking for these functions
?hclust
[with method = "ward"]
?kmeans
Best regards,
Stephen
--- amna khan <[EMAIL PROTECTED]> wrote:
> Hi Sir
>
> How to perform cluster analysis using Ward's method and K- means
> clustering?
>
> Regards
>
> --
> AMINA SHAHZADI
>
Hi,
I am trying to create a plot with pairs()
using a gplot.hexbin() for each pair.
For pairs I can provide a custom upperPanel function:
pairs(iris[1:4], panel=mypanel)
and mypanel() calls plot.hexbin():
library(hexbin) # Bioconductor
mypanel <- function(x, ...){
hb <- hex
Weiwei Shi wrote:
> Hi,
>
> I am looking for an algorithm (better if it is implemented in R) which
> can do the following:
>
> from the following list:
> a,b,c,d
> a,b,c
> b,c
> a,b,c,d,e
>
> to calculate
> a,b,c,d: 2
> a,b,c: 3
> a,b: 3
> a,c: 3
> b,c: 4
> b,c,d:2
>
> here, the order is not im
On 10/18/07, Ravi Varadhan <[EMAIL PROTECTED]> wrote:
> Dear Prof. Ripley,
> You are right. I failed to pay attention to the following remark in the
> help page:
> "Note that the storage used by DQRDC and DGEQP3 differs."
> Although this does not matter for solving a linear system (i.e. qr.solv
Dear all,
I have a question with respect to counting process formulation of the
coxph(survival) model.
I have two groups of observations for which I have partitioned each
observation into two distinct time intervals, namely, entry day till day 13,
and day 13 till death or censorship day (of co
Hi Paul,
Why do you want to do that? Maybe you could publish table along side
your graph if you want people to be able to read the raw numbers.
Hadley
On 10/16/07, H. Paul Benton <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> Just a quick one, hopefully. I have a histogram made from the method
> 'h
I'm a little confused about what I'm doing wrong using the setSqlTypeInfo
function to modify character size output:
> setSqlTypeInfo("MySQL",
+ list(double="double",
+ integer="integer",
+ character="varchar(30)",
+ logical=
Jim Lemon wrote:
> Hi Paul,
> It looks like you have a 1548 bar histogram, which means that you
> can't even resolve the bars on the standard display, much less
> annotate them.
umm yes but I know that I only want to see everything from 0 : 150 so
it's ends up only being 300bar's which I agree
Here is an approach that will work as long as there are less than 32
unique characters in the strings. It uses logical operations to make
the comparison
>library(bitops)
> x <- scan(textConnection('a,b,c,d
+ a,b,c
+ b,c
+ a,b,c,d,e'), sep='\n', what='')
Read 4 items
> x <- strsplit(x, ',')
> x
[[
Dear R People:
I am trying to put Rmpi with R-2.6.0 on an Ubuntu operating system.
When I do the installation, all is well.
But when I go into R and fire up the library, this is what I get:
> library(Rmpi)
Error in dyn.load(file, ...) :
unable to load shared library
'/home/cms/Desktop/R-2.6
Hi, there:
I think one possible solution from me is
for each of b2, I calculate the combinations using combn; then collapse to
calculate by using unlist and table.
For example,
b2[[1]] can give ab, ac, bc, abc,...
b2[[[2]] can give ab, ac, bc, abc
...
then unlist
then table to calculate the f
The rule is search on any combinations whose size is more than 1, like
ab, ac, abc, ...
The storage can be like b2
> b2
[[1]]
[1] "a" "b" "c" "d"
[[2]]
[1] "a" "b" "c"
[[3]]
[1] "b" "c"
[[4]]
[1] "a" "b" "c" "d" "e"
On 10/18/07, Nordlund, Dan (DSHS/RDA) <[EMAIL PROTECTED]> wrote:
> > -Ori
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Weiwei Shi
> Sent: Thursday, October 18, 2007 3:32 PM
> To: jim holtman
> Cc: [EMAIL PROTECTED]
> Subject: Re: [R] alignment algorithm or pattern frequency calculation
>
> for example,
>
> a,b: 3 whic
Writing up a comparision of stat's packages for
work. With the amount that's been written on R
I figured there'd be some sort of guidelines or
benchmarks published but Google and CRAN have
gotten me nowhere.
Any links to quantity-of-data writeups or
data/hardware benchmarks with R would be apprci
the naive way is to exhausive search, but too bad for large and real situation.
On 10/18/07, Weiwei Shi <[EMAIL PROTECTED]> wrote:
> for example,
>
> a,b: 3 which means a and b appear together 3 times in my input list.
>
> On 10/18/07, jim holtman <[EMAIL PROTECTED]> wrote:
> > It would be helpful
for example,
a,b: 3 which means a and b appear together 3 times in my input list.
On 10/18/07, jim holtman <[EMAIL PROTECTED]> wrote:
> It would be helpful if you explained what the numbers mean.
>
> On 10/18/07, Weiwei Shi <[EMAIL PROTECTED]> wrote:
> > Hi,
> >
> > I am looking for an algorithm
System:
Linux kernel 2.6.15 Ubuntu dapper
R version 2.5.1
ESS 5.2.11 on Emacs 21.4.1
-
Colleagues
Having read the posts on producing perfect graphs in R, or using
inkscape to edit R graphics output, I have a related question.
Lately I am publishing in a journal that is ve
It would be helpful if you explained what the numbers mean.
On 10/18/07, Weiwei Shi <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I am looking for an algorithm (better if it is implemented in R) which
> can do the following:
>
> from the following list:
> a,b,c,d
> a,b,c
> b,c
> a,b,c,d,e
>
> to calculate
What is 'foo'? It appears that foo does not take 2 arguments:
> foo <- function(x) x
> foo(1,2)
Error in foo(1, 2) : unused argument(s) (2)
>
On 10/18/07, Aydemir, Zava (FID) <[EMAIL PROTECTED]> wrote:
> Hi
>
> what does this mean when i get above error? When I call function foo(),
> specifying
Hi,
I am looking for an algorithm (better if it is implemented in R) which
can do the following:
from the following list:
a,b,c,d
a,b,c
b,c
a,b,c,d,e
to calculate
a,b,c,d: 2
a,b,c: 3
a,b: 3
a,c: 3
b,c: 4
b,c,d:2
here, the order is not important.
Thanks.
--
Weiwei Shi, Ph.D
Research Scientist
I am not sure I understand your question. I get the following output
from your code and it does not look like act.surv.time is repeated:
> m
treat strata time censorTime act.surv.time censoring
[1,] 1 1 0.8331923 0.654723930.65472393 1
[2,] 1 1 0.33
hie
i'm tryimg to generate two survival data using the following code (I know
its ugly ) but it seems to repeat two of the variables can any one tell me
whats the porblem.
n=20
n1=n/2
n2=n/4
a11=1 ;a12=1.4 ;a21=16 ;a22=a12 * a21
t1<-array(1,c(n1))
Hi
what does this mean when i get above error? When I call function foo(),
specifying parameters a and b as below
a<-5
b<-6
> y <- foo(a,b)
> Error in foo(a, b) : unused argument(s) (6)
Thanks
Zava
This is not an offer (or solici
On 10/18/07, Sebastian Weber <[EMAIL PROTECTED]> wrote:
> Hello all lattice experts!
>
> I'm currently trying to plot inset figures in my lattice graphics. The
> plot is composed out of two panels and in each panel I want to show an
> inset figure which is a closeup of the data in each panel. Thus
Dear Prof. Ripley,
You are right. I failed to pay attention to the following remark in the
help page:
"Note that the storage used by DQRDC and DGEQP3 differs."
Although this does not matter for solving a linear system (i.e. qr.solve()),
the Q and R representations are indeed different.
Thank
On 10/18/07, Jim Lemon <[EMAIL PROTECTED]> wrote:
> Sergey Goriatchev wrote:
> > Hello,
> >
> > Let me describe what I have and what I want to do:
> >
> > I have two (7x6) matrices (call them A and B) that have same row and
> > column names (rows=species, columns=variables) but contain numerical
>
On 10/18/07, amna khan <[EMAIL PROTECTED]> wrote:
> Hi Sir
>
> How to perform cluster analysis using Ward's method and K- means clustering?
For beginning, try to perform it using the GUI Rcmdr.
Regards,
Liviu
__
R-help@r-project.org mailing list
https:
On Thu, 18 Oct 2007, Ravi Varadhan wrote:
> Hi,
>
> This is in response to Simon's observation about QR decomp being way too
> slow for the "badx" matrix posted by Art Owen. This is due to the use of
> LINPACK routine DQRDC. QR decomp is much faster when LAPACK routine is
> used.
>
>> system.tim
I am trying to compute bootstrap confidence intervals for a sample
using R 2.5.1 for Windows.
I can get "Normal", "Basic", "Percentile", "BCa" and "ABC" from
boot.ci() and boot() in the Davison & Hinkley "boot" package.
But I can't figure out how to use tilt.boot() to get the "tilting"
confide
[Resend with proper subject line]
I am trying to compute bootstrap confidence intervals for a sample
using R 2.5.1 for Windows.
I can get "Normal", "Basic", "Percentile", "BCa" and "ABC" from
boot.ci() and boot() in the Davison & Hinkley "boot" package.
But I can't figure out how to use tilt.b
Hi Sir
How to perform cluster analysis using Ward's method and K- means clustering?
Regards
--
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
You want a symbol, so
for (i in 1:10) text(7, order[i], as.name(greeks[i]))
(or as.symbol) gets you from character string to symbol.
BTW1: please check the code you post: yours doesn't work.
BTW2: the "o" is "omicron": you can see a list of the names at
http://www.stat.auckland.ac.nz/~paul/R/C
On Thu, 2007-10-18 at 10:18 -0500, Bret Collier wrote:
> All,
>
> I have been digging around in the help files and found bsamsize in
> Hmisc, but I am wondering if i am using it right.
>
> So, here is the question: given a binomial response (success/failure)
> for 2 groups (treatment/control)
Hi,
I am wondering which classification algorithms in R have implemented
class weight for imbalanced data prediction?
Thanks.
--
Weiwei Shi, Ph.D
Research Scientist
GeneGO, Inc.
"Did you always know?"
"No, I did not. But I believed..."
---Matrix III
___
One of the examples here shows two Y axes each with a label:
library(zoo)
example(plot.zoo)
On 10/18/07, John Theal <[EMAIL PROTECTED]> wrote:
> I have the following R code to create a plot with two y axes. I am
> essentially trying to
> plot a price series with a volume series on the same graph
On Thu, 18 Oct 2007, Ralf Goertz wrote:
> S Ellison, Donnerstag, 18. Oktober 2007:
>>> I think there is reason to be surprised, I am, too. ...
>>> What am I missing?
>>
>> Read the formula and ?summary.lm more closely. The denominator,
>>
>> Sum((y[i]- y*)^2)
>>
>> is very large if the mean value
Hi all users,
Can anyone explain how to do monte carlo error propagation simulation using R.
Thanks in advance.
With regards,
Abu
_
100’s of Music vouchers to be won with MSN Music
__
R-
All,
I have been digging around in the help files and found bsamsize in
Hmisc, but I am wondering if i am using it right.
So, here is the question: given a binomial response (success/failure)
for 2 groups (treatment/control) and I want to estimate the necessary
sample size (n) to determine if
Two follow up questions from post regarding bootstrapping contrasts
for a RM ANOVA:
1. Because my data are repeated measures, isn't it true that I want to
shuffle my column heading individually for each subject (row) by
resampling WITHOUT replacement?
2. For my DV (latency) I have 3 columns of da
Yes, a scatterplot would be better. But for the sake of interest, is
it possible to
get the label for the secondary y-axis to appear? I have been
playing around with
it all day in an attempt to get it to work and have been
unsuccessful. I would like
to know if it is possible and/or what I
A zoo variable is a vector or matrix with an index so you can't
mix types (factors and numeric) in a single zoo variable; however,
you can represented the factor numerically:
library(zoo)
set.seed(1)
a <- data.frame(
nn = letters[1:4],
dd = as.Date("2007-08-01") + 1:4,
x =
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/98227.html
Your Rprofile has the setting
options(stringsAsFactors = TRUE)
If you override it globally by using
options(stringsAsFactors = FALSE)
it will give you what you want but observe Gabor's
caveat. I don't know of the other solution re s
I haven't followed this thread very closely, but it sounds like it may be
related to the somewhat arcane idea of "Desirability functions" for multiple
responses (mostly in th e experimental design context, as I recall). There
were some papers on this in TECHNOMETRICS a couple of decades ago. The
JM
***
XLSolutions Corporation R/Splus Training: November-December, 2007
***
NOVEMBER-DECEMBER: we added 7 new courses to our R/Splus training
webpages
See full s
Thomas Hoffmann wrote:
> Dear R-helpers
>
> I would like to do a fit of the form: y = a (x+c)**b, where a, b and c
> are unknown.
>
> Does anybody know how to do it?
If they did but didn't answer, would it matter?
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.ht
On 10/18/07, John Theal <[EMAIL PROTECTED]> wrote:
> I have the following R code to create a plot with two y axes. I am
> essentially trying to
> plot a price series with a volume series on the same graph. (i.e. to
> compare price with volume). I can label the first
Why not use a scatterplot? Y
Hi all,
I was trying to convert a data frame to a zoo object so I can use some
time series functions like lag(). But it seems then everything became a
factor, so I have to convert it back to numeric to run the correct
regressions. Is there a way to avoid it? Here is an example:
##
> Any suggestions?
When generating a bitmap (png is then the best solution) you can
control the resolution with "width" and "height" parameters (in
pixels):
png("myimage.png",width=1800,height=1800)
Then you just need to raster yor graph in an adequate resolution (more
pixels means bigger file,
Hi,
This is in response to Simon's observation about QR decomp being way too
slow for the "badx" matrix posted by Art Owen. This is due to the use of
LINPACK routine DQRDC. QR decomp is much faster when LAPACK routine is
used.
> system.time(qr(badx, LAPACK=T))
[1] 1.11 0.03 1.14 NA NA
>
S Ellison, Donnerstag, 18. Oktober 2007:
> >I think there is reason to be surprised, I am, too. ...
> >What am I missing?
>
> Read the formula and ?summary.lm more closely. The denominator,
>
> Sum((y[i]- y*)^2)
>
> is very large if the mean value of y is substantially nonzero and y*
> set to 0
Prof Brian Ripley stats.ox.ac.uk> writes:
>
> Please do check the archives: this has been mentioned several times,
> including yesterday, with a workaround (and it is fixed in R-patched).
>
> And I presume you meant the *graphics device* menu, not the console menu.
>
Thank you and yes, I mea
Please do check the archives: this has been mentioned several times,
including yesterday, with a workaround (and it is fixed in R-patched).
And I presume you meant the *graphics device* menu, not the console menu.
On Thu, 18 Oct 2007, Juan Lewinger wrote:
> Juan Lewinger usc.edu> writes:
>
>>
On 10/18/2007 9:42 AM, Gang Chen wrote:
> Hi Liviu,
>
> Thank you very much for the response! I knew that would work within
> R, but was just wondering why
>
> > R CMD INSTALL nlme
>
> does not work on the shell terminal. Any clue?
Just what the message said: on your system, nlme is not a v
I have tried to print a table of greek alphabet names and symbols
without success. I can print one character at a time but can't seem to
find a way to automate an entire list of the symbols. Some of the code I
have tried is below. I have searched on help, worked the examples in the
December 2002 R
Hi Liviu,
Thank you very much for the response! I knew that would work within
R, but was just wondering why
> R CMD INSTALL nlme
does not work on the shell terminal. Any clue?
Thanks,
Gang
On Oct 17, 2007, at 6:27 PM, Liviu Andronic wrote:
> On 10/17/07, Gang Chen <[EMAIL PROTECTED]> wrot
Dear R stat experts,
I am looking for an R function/code that can perform the EM algorithm to
estimate the parameters of a mixture of two negative binomial distributions.
I know that there are functions written that use the EM to estimate a
mixture of normal distributions and also within the pack
I have the following R code to create a plot with two y axes. I am
essentially trying to
plot a price series with a volume series on the same graph. (i.e. to
compare price with volume). I can label the first
y axis successfully, but the problem is in labeling the 2nd y-axis.
Essentially, t
Juan Lewinger usc.edu> writes:
> Also, 1) for the emf generated image look at the little circles at the
bottom of each boxplot: the fill color
> has bled out. Why?
> 2) From the console menu I cannot save a plot as PDF (no file is saved when
I try)
In re to 2) I didn't mention (sorry, didn't
Just saw this in the Windows Secret blog:
Get official and unofficial fixes for Excel
By Brian Livingston
Despite the hotfix that Microsoft recently released for Excel 2007, as
I described on Oct. 11, some math errors that you should know about
still lurk in both Excel 2007 and Excel 2003.
I'
>I think there is reason to be surprised, I am, too. ...
>What am I missing?
Read the formula and ?summary.lm more closely. The denominator,
Sum((y[i]- y*)^2)
is very large if the mean value of y is substantially nonzero and y*
set to 0 as the calculation implies for a forced zero intercept. In
R-Help,
I have a 3 dimensional array with dimensions:
> dim(test)
[1] 241 104 12
I copy-paste an example of the object
> test[1:5,1:5,1:2]
, , ar = 1996
slagnr
puntar12 34 5
I1 60.30303 36.50610 0 0.0 0
I10 58.19512 35.48209 0 75.0 0
I11 0.
On 18-Oct-07 11:56:48, erkan yanar wrote:
>
> Ist there a possibility to write a R-Script using something like
>
>#!/usr/bin/R
> and then alle the requestet commands?
>
> Of course "R CMD BATCH" exists, but there was (for me) no
> possibility to write something like a HERE-script.
> regards
> er
Dear R-helpers,
I need to insert an R (2.6.0) generated plot containing semi-transparent colors
into PowerPoint (2002). When I directly paste it from the clipboard or insert
it as (enhanced) Metafile (I'm on Windows XP) the semi-transparent colors don't
show. When I insert it at as a Bmp, Png o
I found a great site with many Perl/tk examples (it's in German):
http://gd.tuwien.ac.at/languages/perl/Hajji-Perlkurs/part4/tkperl.html
It's quite simple to convert these examples to R's tcltk. One
that gave me some trouble was:
# (please open the file to get the context)
# (...)
# a checkbutto
Thx
erkan
--
über den grenzen muß die freiheit wohl wolkenlos sein
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, mi
Achim Zeileis, Donnerstag, 18. Oktober 2007:
> On Thu, 18 Oct 2007, Toffin Etienne wrote:
>
> > Hi,
> > A have small technical question about the calculation of R-squared
> > using lm().
> > In a study case with experimental values, it seems more logical to
> > force the regression line to pass th
[erkan yanar]
>Ist there a possibility to write a R-Script using something like
>#!/usr/bin/R
>and then alle the requestet commands?
You could use:
#!/usr/bin/env Rscript
(or if you happen to be sure of the location of the program):
#!/usr/local/bin/Rscript
See ?Rscript for more inform
Ist there a possibility to write a R-Script using something like
#!/usr/bin/R
and then alle the requestet commands?
Of course "R CMD BATCH" exists, but there was (for me) no
possibility to write something like a HERE-script.
regards
erkan
--
über den grenzen muß die freiheit wohl wolkenlos
On Thu, 18 Oct 2007, Toffin Etienne wrote:
> Hi,
> A have small technical question about the calculation of R-squared
> using lm().
> In a study case with experimental values, it seems more logical to
> force the regression line to pass through origin with lm(y ~ x +0).
> However, R-squared value
On 18/10/2007 7:02 AM, Etienne Toffin wrote:
> Hi,
> A have small technical question about the calculation of R-squared
> using lm().
> In a study case with experimental values, it seems more logical to
> force the regression line to pass through origin with lm(y ~ x +0).
> However, R-squared
Dear David and Gilles,
Currently, Anova() in the car package doesn't handle lmer objects. I'm
working on a default method for Anova that will work with objects that
respond to coef() and vcov() and produce Wald-like tests. This should
almost work for lmer objects, replacing coef() with fixef(), bu
Hi,
A have small technical question about the calculation of R-squared
using lm().
In a study case with experimental values, it seems more logical to
force the regression line to pass through origin with lm(y ~ x +0).
However, R-squared values are higher in this case than when I
compute th
Hi,
A have small technical question about the calculation of R-squared
using lm().
In a study case with experimental values, it seems more logical to
force the regression line to pass through origin with lm(y ~ x +0).
However, R-squared values are higher in this case than when I
compute th
I looked at this a while ago - what are you trying to do? ie what are you
correlating?
From: [EMAIL PROTECTED] on behalf of Jorge Ivan Velez Valbuena
Sent: Thu 18/10/2007 12:20 AM
To: r-help
Subject: [R] Correlation matrices comparison
Hi all,
Somebody knows
Sergey Goriatchev wrote:
> Hello,
>
> Let me describe what I have and what I want to do:
>
> I have two (7x6) matrices (call them A and B) that have same row and
> column names (rows=species, columns=variables) but contain numerical
> values that differ (same experiment done twice and values calc
> > On 10/17/07, Ravi Varadhan <[EMAIL PROTECTED]> wrote:
> >
> >> What if simultaneously maximizing f(x,y) and g(x,y) is an incompatible
> >> objective?
> >>
> >> Modifying Duncan's example slightly, What if:
> >>
> >> f(x,y) = -(x-y)^2 and
> >> g(x,y) = -(x-2)^2-(y-x-1)^2?
> >>
> >> Here:
> >> (1
Dear useRs,
a new version of the 'relations' package has appeared on CRAN. New
features include:
o support for fuzzy relations added
o support for sets moved to separate 'sets' package
o new SD fitters for the S ("symmetric") and M ("matches") families
o fitters for Cook-Seiford
Dear R-helpers
I would like to do a fit of the form: y = a (x+c)**b, where a, b and c
are unknown.
Does anybody know how to do it?
Thanks
Thomas
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the post
In mgcv::gam smooths are represented using basis functions. The `k' argument
of `s' specifies the number of basis functions to use. Its default value is
10. If you have fewer than k unique covariate values then the basis can not
be set up. So you need to reduce k to at most the number of unique
Hello,
Let me describe what I have and what I want to do:
I have two (7x6) matrices (call them A and B) that have same row and
column names (rows=species, columns=variables) but contain numerical
values that differ (same experiment done twice and values calculated
and put in different matrices).
Hello all lattice experts!
I'm currently trying to plot inset figures in my lattice graphics. The
plot is composed out of two panels and in each panel I want to show an
inset figure which is a closeup of the data in each panel. Thus in each
panel I call xyplot again to plot the inset which I place
On Wed, 2007-10-17 at 19:57 -0600, James wrote:
> On Oct 17, 2007, at 4:36 PM, James wrote:
>
> > On Oct 17, 2007, at 10:18 AM, Waterman, DG ((David)) wrote:
> >
> >> I agree. Avoid the lines like:
> >> iv = c( iv, min(i, j) )
> >>
> >> I had code that was sped up by 70 times after fixing the
83 matches
Mail list logo
