.
Linux. It's best for ***everything***.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped
On 7/09/2007, at 10:17 AM, (Ted Harding) wrote:
On 06-Sep-07 18:42:32, Philip James Smith wrote:
Hi R-ers:
I need to compute the distance between 2 street addresses in
either km or miles. I do not care if the distance is a shortest
driving route or if it is as the crow flies.
Does
On 5/09/2007, at 9:37 AM, Greg Snow wrote:
The best option is to use a bar chart or dot chart instead of a pie
chart.
Right on, Red Freak!!! :-)
cheers,
Rolf Turner
On 5/09/2007, at 11:16 AM, Lisa Hu wrote:
To make it specific, I need to simulation Y with inverse beta
distribution,
that is, Y~inverseF(X), where F is the CDF of beta distribution.
THANKS
On 9/4/07, Lisa Hu [EMAIL PROTECTED] wrote:
Dear All,
I need to use the inverse of some
?
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
cogently expressed!
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
__
R-help
and
indicates a certain
lack of insight which *could* do harm in more complex situations.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti
()
for(i in 1:count) rslt[[i]] - read.csv(paste(etc.))
do.call(rbind,rslt)
}
cheers,
Rolf Turner
##
Attention:\ This e-mail message
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch
is a vector of candidate
values for
the ``join point'', e.g. b - seq(lo,hi,length=100), where lo and hi are
suitably chosen lower and upper bounds on the location of the join
point.
cheers,
Rolf Turner
bsr - function(x,y
on Excel, available at the URL:
http://www.stat.uiowa.edu/~jcryer/JSMTalk2001.pdf
Executive summary: Friends don't let friends use Excel for statistics.
cheers,
Rolf Turner
On 24/08/2007, at 12:51 PM, Moshe Olshansky wrote:
This won't work since it produces a matrix (try this).
On the contrary, Patrick's solution is correct. I tried it. It
works just fine.
cheers,
Rolf Turner
What should work
On 22/08/2007, at 8:44 AM, Douglas Bates wrote:
In some ways it might be better to remove the dependence on coda and
write the classes and methods in the lme4 package. That way I can use
S4 classes and irritate all those people who rail against S4 classes
and methods (and you know who you
On 22/08/2007, at 8:49 AM, Emmanuel Levy wrote:
Hi,
I was wondering the what would be the (most efficient) way to generate
a sequence
of sequences, i mean:
if I have 1,2 and 3.
I'd like to generate a sequence of length N*3 (N ~ 1,000,000 or more)
Where random permutations of the
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
!
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
__
R-help@stat.math.ethz.ch mailing
'' --- as
Mathew *initially* (and quite reasonably) expected it to do.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti
ed., Wiley, 2000) for a good introduction to complex demodulation.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped
that anyone may be able to give.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped
on
unobserved).
Read the documentation for the various function *very* carefully!
Drop me a line if you want more info.
Best wishes,
Ted.
Also, consider mlest() in the mvnmle package.
On 15-Aug-07 21:16:32, Rolf Turner wrote:
I have a data matrix X (n x k, say) each row of which
On 15/08/2007, at 9:53 AM, Jake Verschuyl wrote:
I have some bird flight height data that follows a gamma
distribution. The
data (x) goes from 0 to 700 meters (n=1055). The calculated
parameters
calculated from the fitdistr(x) are (shape = 5.1379, rate =
0.017541), and
therefore
cleverer than I will come to your rescue.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped
On 13/08/2007, at 11:27 AM, Wolfgang Huber wrote:
Dear Rolf and Binabina
perhaps this is of use to some:
Colour for Presentation Graphics. Ross Ihaka.
www.stat.auckland.ac.nz/~ihaka/colour/color.pdf
Choosing Color Palettes for Statistical Graphics
Achim Zeileis and Kurt Hornik.
to Nair's original question; it is
clear
that he is interested in a two-independent-sample t-test.
cheers,
Rolf Turner
##
Attention:\ This e-mail
?append
cheers,
Rolf Turner
On 28/07/2007, at 2:12 AM, Henrique Dallazuanna wrote:
Hi, I don't know if is the more elegant way, but:
X-c(1,2,3,4,5)
X - c(X[1], 0, X[2:5])
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O
On 27/07/07
family form --- and hence is not estimated
by ppm(), at least not
directly.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
https
and p2) in
order that the values of q1 and q2 lie between 0 and 1.
I.e. for some values of rho your goal will be impossible to achieve.
cheers,
Rolf Turner
Typo in my previous email:
sqrt((1-p1)*(1-p2)) should have read sqrt(p1*(1-p1)*p2*(1-p2))
I think!
cheers,
Rolf
##
Attention:\ This e-mail message is privileged
?tapply
##
Attention:\ This e-mail message is privileged and confidenti...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
I was recently asked by one of new colleagues how to combine, or catentate,
or concatentate two factors. If x and y are factors, doing c(x,y) appears
to coerce x and y to numeric mode before catenating them. So what does one
do if one wants the result to be a factor whose levels are the union
much hand-cuffed by the officious ITS policies here as to what
I can install on my Mac. (Effectively, nothing.)
``cheers'' (Hah!),
Rolf Turner
[EMAIL PROTECTED
On 22/5/07 6:48 AM, hadley wickham [EMAIL PROTECTED] wrote:
Hi Rolf,
Two possible solutions:
* DISPLAY=0.0.0.0:0 R - and then X11() should work without having to use
xterm
* install.packages(CarbonEL); library(CarbonEL); quartz()
It is clear that life is determined to frustrate
is that it would be a
hard problem.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo
As I have previously asked, in response to a similar
question: Is this a homework problem?
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch
I must say I agree with Richard O'Keefe who wrote:
I wrote:
R documentation comments really belong
in [.Rd] files where help() can find them.
Barry Rowlingson [EMAIL PROTECTED] replied:
R documentation comments belong in .Rd files at the moment, but how
joyous would
method that I wrote
got dispatched, but the function method didn't?
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help
/library
export R_LIBS
HTH.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
This sounds so simple it must be a homework problem, no?
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
Ricardo Rios wrote:
Hi Rolf Turner, I have a statistical model, it model need this
numbers for calculate the probability. This numbers must be random.
For example I need that
magicfunction(3)
[1] 0.3152460 0.5231614 0.1615926
magicfunction(3)
[1] 0.6147933 0.3122999 0.0729068
?sink
__
R-help@stat.math.ethz.ch 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, minimal, self-contained, reproducible code.
,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting
,
Rolf Turner
__
R-help@stat.math.ethz.ch 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, minimal, self-contained
, what's
the point of having this property for folders?)
Can anyone enlighten me?
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch
by the *relative* path name,
with no problem.
Thanks to all.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
)
cheers,
Rolf Turner
[EMAIL PROTECTED]
Original message:
Dear R helpers,
To illustrate my problem, here is a simplified example. I want to write
a table to a file similar to:
x a
1
.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
! '' as it should.
I would appreciate enlightenment. Ta.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https
?legend
__
R-help@stat.math.ethz.ch 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, minimal, self-contained, reproducible code.
,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org
,
apparently even to the cognoscenti.
cheers,
Rolf Turner
[EMAIL PROTECTED]
P. S. A library, as someone --- Martin Maechler??? --- has repeatedly
tried to explain to the list
Try ``model = FALSE'' rather than ``model = F'' and see if it makes a
difference. You make have an unwanted variable named ``F'' lurking
somewhere.
(In general it is a *bad* idea to use ``F'' when ``FALSE'' is
intended.)
cheers,
Rolf
interpretation of the dollar sign you have
to use *single* quotes.
grep '^dog$' /usr/share/dict/words
*does* work. (Try it!)
cheers,
Rolf Turner
[EMAIL PROTECTED
,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch 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
be a reasonable place to start.
Bon chance.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https
the calculations for you; all
you would be doing is reinventing the wheel.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https
?system
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
are called.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
readable and avoids unintended consequences.
Another infelicity in your code: ``plot=F''. Use ``plot=FALSE''.
(Note that the symbols ``F'' and ``T'' are assignable, *unlike*
``TRUE'' and ``FALSE''.)
cheers,
Rolf Turner
(spiffy,0,3)
0.2353400 with absolute error 2e-09
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
summary(object)$cov.unscaled
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Peter Dalgaard wrote:
Rolf Turner [EMAIL PROTECTED] writes:
summary(object)$cov.unscaled
You need to multiply that with sigma. However, vcov(object) is easier.
Well, I thought unscaled meant unscaled --- the plain
unvarnished covariance matrix! I figure that multiplying
'' is to allow for 0-origin.
data2 - ts(data2,start=0)
???
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing
Well, I like do.call():
puddy.tat - do.call(c,lapply(paste(a,1:10,sep=),get))
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
Peter Dalgaard wrote:
How about
dim(M) - c(4,248,74)
mn - apply(M, c(2,3), mean)
Hey! That's sexy! Much better than my kludgy
suggestion!
cheers,
Rolf
__
be done is via
m.new - t(apply(array(t(m.old),dim=c(74,4,992/4)),c(1,3),mean))
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
, then
the p-value is
uniformly distributed on [0,1].
This should be kept in mind when assessing the ``instability''
of p-values.
cheers,
Rolf Turner
.
cheers,
Rolf Turner
[EMAIL PROTECTED]
Original message:
I have a question regarding normality testing with the nortest
package. I have to admit, my question is so general
plausible. But
still not very plausible. Or maybe major American cities
are even worse than we had been led to believe.
cheers,
Rolf Turner
__
R-help@stat.math.ethz.ch
Jarek wrote:
At one company I was working for, I had to run all the licenses of
all the software I had on my machine, through the legal department.
When they read GNU Public License (GPL) their only comment was: We
have no idea what that license means. Do not touch any software using
it.
,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R
Duncan Murdoch wrote (amongst other things):
Statistical computing is not easy, so how could R be? Who has ever
claimed it is? Any package that makes statistical computing appear to
be easy is probably giving you wrong answers half the time, or is
extremely limited in scope.
for any hints.
cheers,
Rolf Turner
[EMAIL PROTECTED]
P. S. Version info:
version
_
platform sparc-sun-solaris2.9
arch sparc
os
See argument ``fixed'' in help(arima).
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
charts are terrible.
Don't mess with them. Period.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman
a
helpful, humble, ***practical*** response, like that from
Detlef Steuer, from the Microsoft(dung) Monolith.
cheers,
Rolf Turner
[EMAIL PROTECTED
will also (almost surely!)
tell you that the data are not Poisson.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch
be missing
a point or three. As I said, I don't really grok environments.
Given that what I want to do makes any kind of sense at all, can
someone start me off in the right direction?
cheers,
Rolf Turner
Brian Ripley wrote:
Can I also suggest local()? This does a similar thing in a perhaps
more natural way.
Sorry, I'm not with you. I'm slow, and as I said, I don't
really grok environments.
Let's look at a toy example. Suppose I want to create a function foo:
function(x){x^n}
and
Thomas Lumley wrote:
foo - local({
n - 4
function(x) {x^n}
})
Bewdy! Thanks a lot. Slowly but slowly ( :-) ) I learn things!
cheers,
Rolf
.
Is there a sexier way?
cheers,
Rolf Turner
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R
, it
might be nice to have some sort of warning in the help file to keep
young players from getting confused.
cheers,
Rolf Turner
[EMAIL PROTECTED
?
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R
.
cheers,
Rolf Turner
[EMAIL PROTECTED]
P. S. The help on optim() says the Nelder-Mead method is
``relatively slow''. However for this problem optim() seems to come
in about 10 seconds when
,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting
, but I would like to deepen my
understanding. There's probably something about the ``expression''
concept that I'm not grokking here
Thanks for any insight.
cheers,
Rolf Turner
[EMAIL PROTECTED
,
Rolf Turner
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Peter Dalgaard wrote:
Well, come the S4 revolution and dots will cause trouble no more...
Yes but an infinite number of other things will cause infinitely more
trouble. God save us from S4.
cheers,
Rolf Turner
,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
into a
Word document. If so, would you be willing to share it with me and
my student?
If so, please be gentle in your explanation. I am not myself (repeat
***NOT***) a user of Word!
Thanks.
cheers,
Rolf Turner
.
Can anyone help me out?
cheers,
Rolf Turner
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R
apologies for confusing everyone with my initial mis-statement of
the problem.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https
me out?
cheers,
Rolf Turner
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting
Please ignore this message; I apologise for the annoyance.
cheers,
Rolf
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
I think you want to use the function try(); see ?try.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https
to send a separate email
about this.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman
.
It seems to be the complicated strings which form the entries of x that
are messing things up;
y - c(b,c,NA,a,NA)
y[order(y,na.last=FALSE)]
gives
[1] NA NA a b c
is it ought.
cheers,
Rolf Turner
,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
deriv(expression(sum(x)), x)
does not make any sense.
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing
Martin C. Martin wrote:
Rolf Turner wrote:
deriv(expression(sum(x)), x)
does not make any sense.
Good point. But this does:
deriv(expression(sum(log(a*x))), a)
where a is a scalar.
Okay --- I see what you're getting at now.
But I think
1 - 100 of 245 matches
Mail list logo