Hi
In your case it would be wiser to do everything in one step.
plot(c(1.5,2.5), c(media_cob_veloc1.5, media_cob_veloc2.5), ...)
then add your arrows calls
With par(new = T) you will get completely new plot which superpose your
old one. It is like if you take two pictures on the same frame of
Duncan Murdoch wrote:
> On 20/06/2010 6:36 PM, Muenchen, Robert A (Bob) wrote:
>>
>>> -Original Message-
>>> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
>>> On Behalf Of Ivan Calandra
>>> Sent: Sunday, June 20, 2010 3:47 PM
>>> To: r-help@r-project.org
>>> Su
On Mon, 21 Jun 2010, cahyo kristiono wrote:
Dear all..
I have a respon variable y. Predictor variable are x1, x2, x3, x4, x5
(1)?What is?the syntax to get paramater estimation of?ZINB Model by
Newton Raphson (not BFGS)
I guess you are talking about the function zeroinfl() from the package
"p
Dear R users:
In stat, there is a "stpower" function for power analysis and sample-size
determination in survival
models. Is there a counterpart in R?
Thanks
Yao Zhu
Department of Urology
Fudan University Shanghai Cancer Center
Shanghai, China
[[alternative HTML version deleted]]
Hi R-users,
I just started learning R. I have a project on lot quality assurance sampling
(LQAS). In this project I have to develop LQAS plans to make decision on
stopping / continuing a programme. The LQAS plans is based on cluster
sampling: selection of k clusters (villages) of m children ea
Instead of nb2mat try
as.spam.listw(nb2listw(cell2nb(...)))
this will coerce the adjacency matrix into a sparse matrix representation
saving lot of memory.
Nikhil
On Sun, Jun 20, 2010 at 10:27 PM, Daniel Malter wrote:
>
> Hi, thanks much. This works in principle. The corrected code is below:
I have the following script to convert Longitude and Latitude into
lambert projection ( adapted from Geneland manual). How do I get the
proejction with units in metres rather than whatever they are
currently? alos is it possible to add a fualse XY so that all values
are positive ( all my latt
Hi Jeff,
I am not aware of any problems with POSIXct and ggplot2, so I'm not
sure what you're referring to. I don't ask for reproducible examples
to be malicious - they make it easier for me to figure out exactly
what's wrong and recommend a fix. Without a simple reproducible
example, all I can
That's a bit disingenuous, Hadley. There have been reports of problems with
POSIXct and ggplot2 for at least six months, particularly with values that have
no time portion. You have promised a fix before, but l haven't seen it, so I
convert to Date to work around the bug.
"Hadley Wickham" wro
Here is my solution:
tapply(zzz$DEATH1,zzz$GENDER,sum, na.rm=TRUE)
Dear R People:
Here is a little section of a data frame:
> zzz[1:10,]
DATE GENDER AGE Co DEATH1
3945 2009-04-16 M 24Botulinic 23
3851 2009-04-16 M 35 Constitutional 23
8
Hi, thanks much. This works in principle. The corrected code is below:
a <- nb2mat(cell2nb(nrow(x),ncol(x),torus=T), style="B")
g <- delete.vertices(graph.adjacency(a), which(x!=1)-1)
plot(g)
clusters(g)
the $no argument of the clusters(g) function is the sought after number.
However, the funct
Dear R People:
Here is a little section of a data frame:
> zzz[1:10,]
DATE GENDER AGE Co DEATH1
3945 2009-04-16 M 24Botulinic 23
3851 2009-04-16 M 35 Constitutional 23
8495 2009-04-16 F 49 Constitutional 27
10967 2009-04-16
You would probably need to talk to the people maintaining the mirrors as you
would need access to the download statistics.
I presume that the IP addresses of the people who download are also stored
somewhere, so you could possibly georeference the download statistics.
Flame mode on
Then you
At 09:01 PM 6/20/2010, Ted Harding wrote:
On 20-Jun-10 19:49:43, Hadley Wickham wrote:
Whales are a different kettle of fish! They are much more directly
observable, in principle, than are R-users. For one thing, a whale
has to come to the surface to breathe every so often, and if you
are in a sh
On 21-Jun-10 01:30:26, David Winsemius wrote:
> On Jun 20, 2010, at 9:14 PM, David Winsemius wrote:
>> On Jun 20, 2010, at 8:17 PM, (Ted Harding) wrote:
>>
>>> On 20-Jun-10 19:54:02, David Winsemius wrote:
On Jun 20, 2010, at 1:38 PM, Ekaterina Pek wrote:
> Hi, Ted.
> Thanks for your r
Dear all..
I have a respon variable y. Predictor variable are x1, x2, x3, x4, x5
(1) What is the syntax to get paramater estimation of ZINB Model by Newton
Raphson (not BFGS)
(2) What syntax to plot probability of observed & predicted of ZINB
Thx.
Regards
Krist.
[[alternative HTML versio
On Jun 20, 2010, at 9:14 PM, David Winsemius wrote:
On Jun 20, 2010, at 8:17 PM, (Ted Harding) wrote:
On 20-Jun-10 19:54:02, David Winsemius wrote:
On Jun 20, 2010, at 1:38 PM, Ekaterina Pek wrote:
Hi, Ted.
Thanks for your reply. It helped. I have further a bit of
questions.
It may be
On Jun 20, 2010, at 8:17 PM, (Ted Harding) wrote:
On 20-Jun-10 19:54:02, David Winsemius wrote:
On Jun 20, 2010, at 1:38 PM, Ekaterina Pek wrote:
Hi, Ted.
Thanks for your reply. It helped. I have further a bit of questions.
It may be that lm(log(b) ~ log(a)) is, from a substantive point of
On Jun 20, 2010, at 5:52 PM, Schmidt Martin wrote:
Dear R users
A quite simple question: how can I get the mean duration
(persistence) of daily data, when the data set looks as follows:
[1] 7 7 7 9 9 5 7 5 5 5 9 5 5 6 6 6 1 1 1 2 2 4 4 4 4 7 7 3 3 2 4 4
7 7 7 7 7
[38] 7 7 5 5 9 9 5 5 5 1
Hi Michael
Thank you very much for the intel regarding eta^2. It is pretty much the sort
of thing that I am wanting.
Found a good paper regarding all this:
Estimating an Effect Size in One-Way Multivariate Analysis of Variance (MANOVA)
H. S. Steyn Jr; S. M. Ellisa
Multivariate Behavioral Re
On 20-Jun-10 19:49:43, Hadley Wickham wrote:
>> I've given thought in the past to the question of estimating the R
>> user base, and came to the conclusion that it is impossible to get
>> an estimate of the number of users that one could trust (or even
>> put anything like a margin of error to).
>
At 08:44 PM 6/20/2010, RICHARD M. HEIBERGER wrote:
You have only one Age value in each Lot. Hence the Age is aliased
with one of the
5 degrees of freedom between lots. You can see it in the picture
xyplot(Conc ~ Age|Lot, data=blots, pch=16)
There is no opportunity for a slope within each of
Are you sure your data was the Date class?
> x <- read.table('clipboard', header=TRUE)
> x
DATE DEATH
1207 2009-04-16 2009-05-06
1514 2009-04-16 2009-05-06
2548 2009-04-16 2009-05-08
3430 2009-04-16 2009-05-09
3851 2009-04-16 2009-05-09
3945 2009-04-16 2009-05-09
7274 2009-04-16 20
You have only one Age value in each Lot. Hence the Age is aliased with one
of the
5 degrees of freedom between lots. You can see it in the picture
xyplot(Conc ~ Age|Lot, data=blots, pch=16)
There is no opportunity for a slope within each of the groups.
Rich
[[alternative HTML version
I am sure this can be written in a much more elegant and faster code.
One way I can think of, is to modify cell2nb code and create a new
function that creates the neighbour lists of only cells that are not
0. These are best directed to R-sig-Geo list.
However, the following might work.
li
Thanks for the help!
> a.df[1:10,c(1,6,9)]
DATE DEATH DEATH1
1207 2009-04-16 2009-05-06 20 days
1514 2009-04-16 2009-05-06 20 days
2548 2009-04-16 2009-05-08 22 days
3430 2009-04-16 2009-05-09 23 days
3851 2009-04-16 2009-05-09 23 days
3945 2009-04-16 2009-05-09 23 days
7274 2
Hi,
Are you sure these are date objects and not strings? For example
> d1<-runif(10,0,100)
> d2<-runif(10,0,200)
> df<-data.frame(d1=d1+as.Date("2010-6-20")+d1,d2=as.Date("2009-6-20")+d2)
> df
d1 d2
1 2010-09-23 2009-06-30
2 2010-06-25 2009-08-21
3 2010-10-10 2009-08-04
4
I'm using R 2.10.1 with the latest version of all packages (updated today).
I'm confused as to why I'm getting a hard singularity in a simple set
of experimental data:
> blots
ID Lot Age Conc
1 1 A 3 4.44
2 2 A 3 4.56
3 3 B 41 4.03
4 4 B 41 4.57
5 5 C 229 4.49
6
On 20-Jun-10 19:54:02, David Winsemius wrote:
> On Jun 20, 2010, at 1:38 PM, Ekaterina Pek wrote:
>> Hi, Ted.
>> Thanks for your reply. It helped. I have further a bit of questions.
>>
>>> It may be that lm(log(b) ~ log(a)) is, from a substantive point of
>>> view, a more appropriate model for wh
Dear R People:
I have a data frame with the two following date columns:
> a.df[1:10,c(1,6)]
DATE DEATH
1207 2009-04-16 2009-05-06
1514 2009-04-16 2009-05-06
2548 2009-04-16 2009-05-08
3430 2009-04-16 2009-05-09
3851 2009-04-16 2009-05-09
3945 2009-04-16 2009-05-09
7274 2009-04-16
Use a variable of class "factor"
> s <- c("ABC", "ACCDEDF", "ACCGEDF", "ACCGEGF", "ACCDEDF", "ACCGEGF")
> fs <- factor(s)
> levels(fs)
[1] "ABC" "ACCDEDF" "ACCGEDF" "ACCGEGF"
> unclass(fs)
[1] 1 2 3 4 2 4
attr(,"levels")
[1] "ABC" "ACCDEDF" "ACCGEDF" "ACCGEGF"
On Sun, Jun 20, 2010 at
Hi,
I have been a matlab user and is learning R.
I want to convert a large list of strings to a list of unique numeric
ids to reduce storage space.
For example,
there is a string list (there are duplicates)
ABC
ACCDEDF
ACCGEDF
ACCGEGF
.
ACCDEDF
ACCGEGF
and I want to have a correspondi
I'm not sure what exactly you want, but maybe ?rle can help you.
-Peter Ehlers
On 2010-06-20 15:52, Schmidt Martin wrote:
Dear R users
A quite simple question: how can I get the mean duration (persistence)
of daily data, when the data set looks as follows:
[1] 7 7 7 9 9 5 7 5 5 5 9 5 5 6 6
On 20/06/2010 6:36 PM, Muenchen, Robert A (Bob) wrote:
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Ivan Calandra
Sent: Sunday, June 20, 2010 3:47 PM
To: r-help@r-project.org
Subject: Re: [R] Popularity of R, SAS, SPSS, Stata.
Hi all, I am sorry if this is a very basic quesion, but I have no experience
with analyzing spatial data and could not find the right function/package
quickly. Any hints would be much appreciated. I have a matrix of spatial
point patterns like the one below and want to find the number of independe
>-Original Message-
>From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org]
>On Behalf Of Muenchen, Robert A (Bob)
>Sent: Sunday, June 20, 2010 6:43 PM
>To: Hadley Wickham; ted.hard...@manchester.ac.uk
>Cc: r-help@r-project.org
>Subject: Re: [R] Popularity of R, SAS, SPS
>-Original Message-
>From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org]
>On Behalf Of Hadley Wickham
... What about snowball
>sampling with R-help as an initial frame?
That's an interesting idea! I could put together a Two-item web survey:
1. What stat package do
>-Original Message-
>From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
>On Behalf Of Ivan Calandra
>Sent: Sunday, June 20, 2010 3:47 PM
>To: r-help@r-project.org
>Subject: Re: [R] Popularity of R, SAS, SPSS, Stata...
>
>Bob,
>
>I have no idea whether it is realistic
Thanks all for the advice,
I'm using the time with parts of a second (to 6 decimal places) hashed, and so
I ran a loop : (my attempt at repeatable code!)
> library(digest)
> op <- options(digits.secs=6)
> a <- NA
> for (i in 1:5) { a[i] <- digest(Sys.time(), algo='crc32') }
> options(op)
an
Hi Ottar,
It's impossible to tell what the problem is without a reproducible
example (http://gist.github.com/270442)
Hadley
On Sun, Jun 20, 2010 at 4:38 PM, Ottar Kvindesland
wrote:
> I have a problem that puzzles me a bit today. When loading off data from a
> database and plotting using ggplot
On Jun 20, 2010, at 6:12 PM, David Winsemius wrote:
On Jun 20, 2010, at 5:40 PM, Nick Matzke wrote:
Wow, did I stump everyone with this? Mostly I'd just like to build
a class without getting the error below.
Hardly.
a) LHS= left hand side
b) It is not an error message in the first place
>-Original Message-
>From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org]
>On Behalf Of Ted Harding
>Sent: Sunday, June 20, 2010 3:42 PM
>To: r-help@r-project.org
>Subject: Re: [R] Popularity of R, SAS, SPSS, Stata...
>
>
>I've given thought in the past to the question
How about getting statistics of downloads of the R-base from the different CRAN
mirrors ?
This should (in principle) allow one to estimate the total # of people who
intended to use R at some point in their life.
It may even be possible to analyze those numbers for temporal trends since the
da
Thanks Rainer,
I actually solved this by just adding a few fonts - from the x11-fonts ports.
I'm not sure about choosing the best ones, but this seems to have fixed the
"Cairo" device, although my capabilities has not changed at all.
We thought we could avoid the need to install X11 at all, but
On Jun 20, 2010, at 5:40 PM, Nick Matzke wrote:
Wow, did I stump everyone with this? Mostly I'd just like to build
a class without getting the error below.
Hardly.
a) LHS= left hand side
b) It is not an error message in the first place, and for another
thing, why are you creating instanc
Dear,
I want to compute coefficient of determination (R-squared) to complement AIC
for model selection of
multivariable GLM.
However, I found this is not a built-in function in glm. neither is it
available through reviewing the question in the R-help archive.
Please kindly help and thanks a lot.
Dear R users
A quite simple question: how can I get the mean duration (persistence)
of daily data, when the data set looks as follows:
[1] 7 7 7 9 9 5 7 5 5 5 9 5 5 6 6 6 1 1 1 2 2 4 4 4 4 7 7 3 3 2 4 4 7 7 7 7 7
[38] 7 7 5 5 9 9 5 5 5 1 7 9 9 9 9 9 5 5 5 5 5 3 4 4 4 4 4 4 4 4 2 2 2 2 2 2 8
Wow, did I stump everyone with this? Mostly I'd just like
to build a class without getting the error below.
Cheers!
Nick
Nick Matzke wrote:
Hi all! I'm trying to build a simple class, but I'm getting a weird error.
E.g. if I do:
==
gregion = setClass("gregion", contains = "data.fra
I have a problem that puzzles me a bit today. When loading off data from a
database and plotting using ggplot2, I wish to present data as a time series
with time of day. The code is text-book like as shown below:
# Fetch data
con <- dbConnect(dbDriver("MySQL"), user="user", password = "pwd",
dbnam
On Jun 20, 2010, at 1:38 PM, Ekaterina Pek wrote:
Hi, Ted.
Thanks for your reply. It helped. I have further a bit of questions.
It may be that lm(log(b) ~ log(a)) is, from a substantive point of
view,
a more appropriate model for whetever it is than lm(b ~ a). Or it may
not be. This is a s
> I've given thought in the past to the question of estimating the R
> user base, and came to the conclusion that it is impossible to get
> an estimate of the number of users that one could trust (or even
> put anything like a margin of error to).
I find it hard to believe that it should be harder
Bob,
I have no idea whether it is realistic, but if you look for the papers that
used R or SAS (or anything), you might get better results by searching for the
way R and SAS are cited.
It looks to me that what I'm saying is not clear, so here an example.
To cite R in a paper you have to write
On 20-Jun-10 19:07:21, Muenchen, Robert A (Bob) wrote:
>>I wonder if there are any capture-recapture type methodologies for
>>estimating open-source software usage? Another idea would be to
>>combine with some other known numbers, e.g. book sales, conference
>>attendance etc. You'd need personal i
>I wonder if there are any capture-recapture type methodologies for
>estimating open-source software usage? Another idea would be to
>combine with some other known numbers, e.g. book sales, conference
>attendance etc. You'd need personal information to link the data sets
>together.
>
>Hadley
This
Hi all, good morning,
My question is not really R related rather a practical problem and wondering if
stat-gurus here can show some light how that can be solved with some
statistical/mathematical tool.
I have some 10 items on which 10,000 viewers put their views based on some 12
attributes say
>-Original Message-
>From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org]
>On Behalf Of David Winsemius
>Sent: Sunday, June 20, 2010 1:05 PM
>To: Stefan Grosse
>Cc: r-help@r-project.org
>Subject: Re: [R] Popularity of R, SAS, SPSS, Stata...
>
>
>On Jun 20, 2010, at 10:
Hello everybody
I have some questions:
1)Â Â Â Â Â For generating data from Generalized pareto distribution, I used
POT package and I have example that used shape=2, but it canât be appropriate.
2)Â Â Â Â Â Exact density of histogram of theta in this program doesnât
display.
i<-1
X<-rep(0,
> I agree with all your points. What I have so far is nowhere near the big
> picture, but it's a start. When you install some software it asks if you
> mind it reporting usage stats back to its home site. I know that sort of
> thing has been discussed before on R-help. I'd love to see that added so
Hi, Ted.
Thanks for your reply. It helped. I have further a bit of questions.
> It may be that lm(log(b) ~ log(a)) is, from a substantive point of view,
> a more appropriate model for whetever it is than lm(b ~ a). Or it may
> not be. This is a separate question. Again, Spearman's rho is not
> de
>-Original Message-
>From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org]
>On Behalf Of Stefan Grosse
>Sent: Sunday, June 20, 2010 10:25 AM
>To: r-help@r-project.org
>Subject: Re: [R] Popularity of R, SAS, SPSS, Stata...
>
>Am 20.06.2010 15:31, schrieb Muenchen, Robert
On Jun 20, 2010, at 10:24 AM, Stefan Grosse wrote:
Am 20.06.2010 15:31, schrieb Muenchen, Robert A (Bob):
I've been fiddling around with various ways to estimate the
popularity
of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-
PLUS,
R-PLUS and Revolution R. It's not an eas
Am 20.06.2010 15:31, schrieb Muenchen, Robert A (Bob):
> I've been fiddling around with various ways to estimate the popularity
> of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-PLUS,
> R-PLUS and Revolution R. It's not an easy task. You can see what I've
> come up with so far at
On Sun, Jun 20, 2010 at 09:09:58AM -0500, li li wrote:
> Dear all,
> Does anyone know where to download the package 'Biobase' in R.
> Thank you!
>Hannah
http://lmgtfy.com/?q=R+package+biobase
--
Three out of two people have difficulties with fractions.
On Jun 20, 2010, at 8:29 AM, Hadassa Brunschwig wrote:
Hi,
I am aware that this question might be a basic one.
I did browse the help and archives but I still haven't understood
how to do the following:
I run R locally and would like to read in data from a server
for which there is a username
On Jun 20, 2010, at 10:09 AM, li li wrote:
Dear all,
Does anyone know where to download the package 'Biobase' in R.
Thank you!
It's a BioConductor bundle. Read the BioC documentation.
--
David Winsemius, MD
West Hartford, CT
__
R-help@r-pr
Dear all,
Does anyone know where to download the package 'Biobase' in R.
Thank you!
Hannah
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Hi All,
I've been fiddling around with various ways to estimate the popularity
of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-PLUS,
R-PLUS and Revolution R. It's not an easy task. You can see what I've
come up with so far at http://r4stats.com/popularity . I'm sure people
will h
I just checked, and it doesn't seem to have failed:
> model.lrm24$fail
[1] FALSE
Any other possibilities?
From: Christos Argyropoulos
Sent: Sun, June 20, 2010 7:23:34 AM
Subject: RE: [R] "Unable to fit" error message from the lrm function in the rms
library
On 20-Jun-10 12:02:58, Ekaterina Pek wrote:
> Hi,
> It's rather statistical question than R-question.
>
> There is some linear correlation but the picture "plot(a,b) +
> abline(lm(b~a))" is quite crowded in the left lower corner. The
> picture "plot(log(a), log(b)) + abline(lm(log(b)~log(a))" is m
Hi,
I am aware that this question might be a basic one.
I did browse the help and archives but I still haven't understood
how to do the following:
I run R locally and would like to read in data from a server
for which there is a username and password. That is,
how do I open a connection to a serv
Gabor Grothendieck [Sat, Jun 19, 2010 at 11:19:35AM
CEST]:
> On Sat, Jun 19, 2010 at 4:56 AM, Johannes Huesing
> wrote:
> > A rule would be a data structure containing the ID of the rule, the
> > rule in human readable language, an expression evaluating variables
> > within the environment of t
Hi,
It's rather statistical question than R-question.
There is some linear correlation but the picture "plot(a,b) +
abline(lm(b~a))" is quite crowded in the left lower corner. The
picture "plot(log(a), log(b)) + abline(lm(log(b)~log(a))" is much
nicer ("Milky Way"). Is it correct to use the secon
Hi,
The error message you are getting (probably) means that the algorithm did not
converge. Did you check for convergence? (Look at the "fail" element of the
returned lrm object)
Christos
__
73 matches
Mail list logo
