,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
'' is your data) do what you want? You can then save the
result
to a txt or csv file in whatever manner suits you; sink(), write(),
write.table(),
write.csv(),
cheers,
Rolf Turner
##
Attention
or non-parametric bootstrapping)
would seem to be
quite feasible, but I thought I'd look for an ``analytic'' approach
first.
Thanks for any help that anyone can give.
cheers,
Rolf Turner
##
Attention
:
if(interactive()) {
# Do the non-batch thing
} else {
# Do the batch thing.
}
Haven't tested this, but.
cheers,
Rolf Turner
##
Attention:\ This e
On 20/03/2009, at 9:17 AM, Lyman, Mark wrote:
Thanks for the suggestion Rolf. Unfortunately, this uses the BATCH
profile when I use for example R CMD INSTALL. I would like the BATCH
profile to be used only when R CMD BATCH is used. Sorry, about the
vagueness of my question.
I have found a way
?
Is it possible to twiddle the control parameters, for either or both
functions,
so as to obtain identical results?
Thanks.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged
On 20/03/2009, at 12:25 PM, Peter Palenchar wrote:
I have a list:
x-c(abcg, bcah, defi)
I want the third and fourth thing in each item of the list (bc,
ca,
ef).
?substr
##
Attention:\ This e-mail message is
going on in the actual hypothesis testing mechanism.
HTH.
cheers,
Rolf Turner
##
Attention:
This e-mail message is privileged and confidential. If you
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
that is sufficiently Unix-like for the following roll-your-own to work:
cfsl - function(file) {
# cfsl -- check for symbolic link
xxx - system(paste(ls -l,file),intern=TRUE)
if(substr(xxx,1,1)==l) TRUE else FALSE
}
HTH
cheers,
Rolf Turner
? *Are* there ``R binary'' files lurking about that
I am somehow
not seeing? Why won't read.xls() work on this data set?
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid
On 25/03/2009, at 12:09 PM, Frank E Harrell Jr wrote:
snip
(2) Scrolling down to ``Byar and Green prostate cancer data''
appeared
to get
me to the right place. But I couldn't see any signs of any ``R
binary
files''.
Please look again. It's under the heading R. Unfortunately
with a digit.
And it can't have white space in it.
There are probably other rules, but essentially anything *sensible*
as a variable name can be used as a variable name.
cheers,
Rolf Turner
##
Attention:\ This e
On 27/03/2009, at 2:52 PM, Marc Schwartz wrote:
On Mar 26, 2009, at 8:40 PM, Rolf Turner wrote:
On 27/03/2009, at 2:04 PM, Mike Miller wrote:
Importing data with a header row using read.delim, one variable
should be
named @5HTT but it is automatically renamed to X.5HTT, presumably
interface --- Unsubscribe or edit options
Doesn't seem too difficult to me.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
, the spatstat package provides
quite
a few facilities for analyzing quadrats/quadrat counts. This package
might
prove to be of some use to you.
cheers,
Rolf Turner
##
Attention:\ This e-mail message
Hi.
I am trying to do histograms in lattice, and I want to get both
counts and percents in the same plot. To try to be clearer ---
there are 3 levels to my factor; I'd like to get a 2 x 3 array
of plots where the top row consist of histograms by counts and
the bottom consists of (the
Thanks very much for your input.
On 30/03/2009, at 12:29 PM, Felix Andrews wrote:
Modifying your example...
library(lattice)
set.seed(42)
XX - data.frame(y=runif(300,0,10),a=factor(sample(letters[1:3],300,
TRUE,c(0.5,0.3,0.2
XX - rbind(XX,XX)
0.098 2.040
so unless there's something that I am misunderstanding (always a serious
consideration) Wacek's apply method looks to be about 1.4 times
*faster* than
the do.call/pmax method.
cheers,
Rolf Turner
On 30/03/2009, at 3:55 PM, Bert Gunter wrote
that means.
Bottom line --- I can't watch the video.
But that's the story of my life. ***Nothing*** ever works for me! :-)
Except R, which *mostly* works.
cheers,
Rolf Turner
On 31/03/2009, at 12:49 PM, Ted Harding wrote:
On 30-Mar-09 23:04:40, Jim Porzak wrote
This whole thing is an April Fool's joke. Isn't it?
***Please***!!! (Let it be an April Fool's joke.)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
On 2/04/2009, at 8:37 AM, Jason Rupert wrote:
Is there any syntax in R that allows a switch-type condition to
be used?
switch(variable){
case CONSTANT_VALUE;
break;
default:
break;
}
?switch
##
Attention:\
infinite values.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Uhhh, Bill, he wanted E to be recoded as ``Treat 3''.
And he didn't say ***what*** he wanted to happen to D.
Fancy. A chance to ``correct'' Bill Venables for a second time
in two days! :-)
cheers,
Rolf
On 2/04/2009, at 12:39 PM, bill.venab...@csiro.au wrote:
Here
On 2/04/2009, at 11:27 PM, Uwe Ligges wrote:
snip
Not many of us are clairvoyants
snip
Hey! Wish I'd said that! :-)
cheers,
Rolf
P. S. Great minds think alike.
R.
,].
And that's what you get. Look carefully --- it is not correct to say
that when
you looked at d2 *no* data were removed; d2 is equal to d1 with its
first row
removed. I.e. your operation didn't remove what you wanted, but it
removed
*something*.
cheers,
Rolf
to plain number:
-10, 0 , 10, 20, ...etc
Look at the ``scipen'' argument of options().
A perfunctory experiment indicates that options(scipen=1) should work
in your case.
cheers,
Rolf Turner
, 10, 30, ...etc
^
in the x-axis tick marks.
axis(side=1,at=0) might give what you want.
But I get a 0 tick mark plotted automatically when I try a simple
example.
This could be due to OS difference I guess.
cheers,
Rolf Turner
On 4/04/2009, at 10:37 AM, Emmanuel Charpentier wrote:
Le vendredi 03 avril 2009 à 14:17 -0400, stephen sefick a écrit :
I am starting to use R for almost any sort of calculation that I
need.
I am a biologist that works in the states, and there is often a need
to convert from standard
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
to have become an obsession with some
of the
R community, particularly many of the senior members --- bugs the
living
Drambuie out of me.
cheers,
Rolf Turner
##
Attention:\ This e
,
Rolf Turner
set.seed(1)
alpha - numeric(100) # Vector to store coefficient estimates
X1 - numeric(100)
X2 - numeric(100)
a - 0 # True parameters
B1 - .5
B2 - .85
for (i in 1:100){
x1 - rnorm(1000)
x2 - rnorm(1000)
epsilon - rnorm(1000)
LP - a + B1*x1 + B2*x2 + epsilon # Linear
On 8/04/2009, at 2:31 PM, Roslina Zakaria wrote:
Hi Rolf,
I would like to extend the problem that I asked you before
regarding the newton method using 4 functions with 4 parameters.
My functions involve the modified bessel function of the first kind
which I can type them without any
There is no bug.
Both of your examples work fine for me. You may have a corrupt
version of ``iris'' somewhere in your search path before ``datasets''.
What does find(iris) tell you? What does names(iris) tell you?
cheers,
Rolf Turner
On 9/04/2009, at 8:16 AM, Jose
called
in the right circumstances.
Sounds like a job for namespaces (about which I was *complaining*
recently;
:-) ) although I don't understand these well enough to be sure that they
could help here. Peter will know, I have every confidence!
cheers,
Rolf Turner
and the discrete Fourier transform, the latter
being what
is calculated by fft(). (NOTE: ***NOT*** ``FFT'' --- R is case
sensitive.)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged
On 17/04/2009, at 10:21 PM, Duncan Murdoch wrote:
Benjamin Tyner wrote:
Many thanks Duncan. Perhaps this merits a more explicit note in the
documentation?
The quote I gave is from the documentation. How could it be more
explicit?
This is unfortunately typical of the attitude of R-core
On 19/04/2009, at 9:45 PM, Duncan Murdoch wrote:
On 18/04/2009 8:47 PM, Rolf Turner wrote:
On 17/04/2009, at 10:21 PM, Duncan Murdoch wrote:
Benjamin Tyner wrote:
Many thanks Duncan. Perhaps this merits a more explicit note in the
documentation?
The quote I gave is from the documentation
On 19/04/2009, at 8:59 PM, Patrick Burns wrote:
Rolf Turner wrote:
On 17/04/2009, at 10:21 PM, Duncan Murdoch wrote:
Benjamin Tyner wrote:
Many thanks Duncan. Perhaps this merits a more explicit note in the
documentation?
The quote I gave is from the documentation. How could
end up with:
f
[,1] [,2] [,3]
[1,]016
[3,]361
Note: 1 is an arbitrary cut-off value.
d[apply(d,1,function(x){any(x=1)}),apply(d,2,function(x){any(x=1)})]
cheers,
Rolf Turner
W_t according to phi(B)W_t = theta(B)a_t and then set
Z_t = W_t + mu
where mu = phi_0/phi(1).
HTH
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
to
Sigma[1,1] - (Sigma[1,2]^2)/Sigma[2,2]
Knowing that, you can use pnorm to calculate Pr(Xx | Y=y).
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
On 2/10/2008, at 11:02 AM, Ravi Varadhan wrote:
I think it is meaningful to ask for a non-trivial Pr (X x, Y=y)
when you
are writing down the likelihood for parameter estimation. This is
commonly
the case in likelihood estimation in bivariate failure time
models. If one
interprets
.
(c) xlim is not a function; it is an *argument* (to the plot3d()
function.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
/Rhelp02a/archive/137154.html
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
() and loess().
cheers,
Rolf Turner
##
Attention:
This e-mail message is privileged and confidential. If you are not the
intended recipient please delete the message and notify the sender.
Any views
.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
.'' :-) )
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch
On 8/10/2008, at 5:16 PM, Deepayan Sarkar wrote:
On 10/7/08, Rolf Turner [EMAIL PROTECTED] wrote:
I am trying to do a histogram lattice plot and I would like the
histogram to be filled with a different colour in each panel.
Note: I want every bar in each histogram to be the same colour
On 8/10/2008, at 6:48 PM, Erin Hodgess wrote:
Hi R People:
I am looking at the Braun/Murdoch book, A First Course in
Statistical Programming in R, and I have a question about a function
there. It's on page 52, Example 4.5; the sieve of Erastosthenes.
There is a line:
primes - c()
Is there
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
discussion and provides a number of
references.
(Feller volume 1 seems to be a good place to start, as usual.)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid
what you want? If not, please
get back to us with more detail.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
there be an FAQ discussing the difference between ifelse()
and if{ } else { }? I couldn't find one, and this issue seems to
arise over and OVER again on this list.
Unless their noses are rubbed in it, many people seem not to grok the
distinction.
cheers,
Rolf Turner
a
line linking the points (m.residuos).
For crying out loud don't write ``wanna''. In spoken English it may
be acceptable (to some people) to *pronounce* ``want to'' as ``wanna''.
It is *not* acceptable to write it that way.
Rolf Turner
On 14/10/2008, at 9:02 AM, Michael Just wrote:
Hello,
How can I add notes (i.e. text) to a sink output?
sink(test.txt)
#This text will describe the test
summary(x)
sink()
How can I add that text above to the sink output?
?cat
?
On the contrary. This is *exactly* what it means.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org
I can't help with your inquiry generally, but I can address the
issue of ``optional braces''.
On 16/10/2008, at 8:34 AM, Ted Byers wrote:
snip
I take it FUN = function(x, grp) quantile(x$DATA, names=FALSE) is the
function definition for a function called FUN. I would guess,
On 16/10/2008, at 10:03 AM, jim holtman wrote:
try putting as.character in the call:
x = read.csv(as.character(V4[[i]]), header = FALSE
No. This won't help. V4 is a column of the data frame optdata,
and hence is a vector. Not a list! Use single brackets --- V4[i] ---
and all will be
,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
On 16/10/2008, at 10:47 AM, Wacek Kusnierczyk wrote:
Rolf Turner wrote:
On 16/10/2008, at 12:28 AM, Henrique Dallazuanna wrote:
Try this:
lapply(1:n, rnorm)
That has nothing to do with what the inquirer *asked*.
probably because it's not what the inquirer knew he *should* have
asked
On 16/10/2008, at 10:44 AM, Erin Hodgess wrote:
Dear R people:
Is there a way to perform simple polynomial multiplication; that is,
something like
(x - 3) * (x + 3) = x^2 - 9, please?
I looked in poly and polyroot and expression. There used to be a
package that had this, maybe?
The
, 16 Oct 2008, Rolf Turner wrote:
On 16/10/2008, at 10:03 AM, jim holtman wrote:
try putting as.character in the call:
x = read.csv(as.character(V4[[i]]), header = FALSE
No. This won't help. V4 is a column of the data frame optdata,
and hence is a vector. Not a list! Use single brackets
package ``Oarray'' from CRAN.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
On 16/10/2008, at 2:27 PM, Moshe Olshansky wrote:
Hi Rolf,
Thank you for making me aware of the existence of PolynomF package.
By the way, your solution needs a small modification:
a - polynom(c(-3,1)) and not polynom(-3,1) and similar for b.
Indeed; I typed the code into the email
- as.function(l[[4]])
foo(3)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https
On 17/10/2008, at 8:22 AM, useR wrote:
Let's try to bring this discussion back again after Frank made
very funny remark!
Frank's remark was *serious*. Take it seriously.
cheers,
Rolf Turner
is that the system is exactly singular!!!
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
to make parity errors in respect of looking at upper and lower tails
when
dealing the cumulative distribution of a discrete random variable.
cheers,
Rolf Turner
##
Attention:\ This e-mail message
On 22/10/2008, at 8:18 AM, Ted Byers wrote:
snip
... even with all the power and utility of Excel ...
snip
Is this some kind of joke?
cheers,
Rolf Turner
##
Attention:\ This e
:30])
num [1:11] -0.315 -0.693 -0.771 0.448 0.204 ...
Try x[20:30,1,drop=FALSE]
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
. You will need to clarify it
considerably.
If you do so, you may be able to pose a meaningful question, and
perhaps answer
it yourself.
cheers,
Rolf Turner
On 23/10/2008, at 11:59 AM, Tom.O wrote:
This might not be the correct forum for this question
contributions.
HTH
cheers,
Rolf Turner
##
Attention:
This e-mail message is privileged and confidential. If you are not the
intended recipient please delete the message and notify the sender.
Any views or opinions
a clever
mate; he couldn't; he asked some clever mates, who finally came up
with the proof
that I have attached in pdf form.
(Anyone who's interested: If the attachment gets removed by the list,
send me email
and I'll send the pdf file to you directly.)
cheers,
Rolf
to the literature,
including Peter's own substantial contributions.
HTH
cheers,
Rolf Turner
##
Attention:
This e-mail message is privileged and confidential. If you are not the
intended recipient please delete
to the literature,
including Peter's own substantial contributions.
HTH
cheers,
Rolf Turner
##
Attention:
This e-mail message is privileged and confidential. If you are not the
intended recipient please delete
file to you directly.)
cheers,
Rolf Turner
P. S. ``WLU'' (in the pdf file) means Wilfrid Laurier University
(Waterloo, Ontario, Canada).
##
Attention:
This e-mail message is privileged and confidential
as it is, to make up for
this lack
of knowledge. Learn some statistics first, before you try to do a
statistical
analysis.
cheers,
Rolf Turner
P. S. If just want to test for a difference between two distributions/
populations
on the basis of independent samples from
19
101 1 10100 110
by(samples$Effort,samples$Region.Label,length)
samples$Region.Label: 1
[1] 10
I.e. I get 10 as you expected.
cheers,
Rolf Turner
##
Attention
- data.frame(hester.=h,value=z)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing
46 whereas the response y has been
trimmed
down to length 44 by eliminating any rows of the data where any of
the variables
involved are missing. Hence a problem.
The solution of the problem requires some code re-writing by the
maintainer of geepack.
cheers,
Rolf
On 30/10/2008, at 10:46 AM, [EMAIL PROTECTED] wrote:
I usually just run my R programs at the R command prompt but for my
latest one I want to save any output that gets written to the
screen so
I am
trying to use R CMD BATCH and send the output to an output file. I
realize I could use sink
condemned to eternal perdition for doing so however. They
would say that you extract things from S4 objects *only* using the
extractor functions provided. How one knows what extractor functions
actually *are* provided remains a mystery to me, but there it is! :-)
cheers,
Rolf
hanging around in your search path?
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org
, or you
could say that R-squared is -100%. (!!!)
Bottom line --- the R-squared concept makes no sense in this context.
The R-squared concept is at best dubious, and should be used, if at all,
only in the completely orthodox setting.
cheers,
Rolf Turner
analysis. Doing tests
mindlessly without a clear understanding of exactly what it is that
you are testing is a recipe for nonsense if not disaster.
cheers,
Rolf Turner
##
Attention:\ This e-mail message
to optim().
Remember that the ... arguments must be given in the name=value form.
Remember that optim() usually doesn't like it if fn() returns Inf or
NA or NaN.
That's about all I can come up with in the way of ideas. Good luck.
cheers,
Rolf Turner
P. S
with `counts' here
Is there a more efficient and/or idiomatic way to do this?
No.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
would be very much appreciated!
Look at the last line of your data frame.
cheers,
Rolf Turner
##
Attention:
This e-mail message is privileged and confidential. If you are not the
intended recipient please
.
cheers,
Rolf Turner
P.S. The foregoing URL discusses the issues in terms of a *single*
observation. Note that if X_1, ..., X_n are i.i.d. Poisson(lambda)
then Y = X_1 + ... + X_n is Poisson(n*lambda).
Note that the ``exact'' confidence limits are very conservative
Boy are you confused. This has nothing at all to do with substitution.
Instead do
test - with(fc,ave.fc[match(diff_mirs_list,Probe)])
cheers,
Rolf Turner
On 7/11/2008, at 11:46 AM, Iain Gallagher wrote:
Hello list.
I have a vector of values:
eg
head
of tests
into a
single ``suitcase''. (``Portmanteau'' is an old-fashioned word for
suitcase.)
HTH.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
. This is the approach that is (almost?) universally adopted in
constructing prediction
intervals in (ARMA) time series analysis.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged
problem?
Thanks.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https
On 12/11/2008, at 4:28 AM, Terry Therneau wrote:
I've read the back and forth this morning, and I have to side with
Vince.
Well, I've read back and forth this morning and I have to side
with Berwin Turlach --- whose postings were, I thought, extremely
well
On 12/11/2008, at 11:29 AM, Peter Dalgaard wrote:
Ben Bolker wrote:
Sometime soon when I have the time and energy I will start
campaigning for an additional drop argument to subset that
does what one expects (!!??) with subsetted factor variables ...
Not that one again! For at least one
On 13/11/2008, at 7:50 AM, sarahkm wrote:
Hi-
I'm a biologist trying to figure out the growth rate of salamanders in
different ponds. I collected individuals from various populations at
different dates, and using the size and date collected, I want to
figure out
the growth curve of each
Please read and understand the response that I already sent you
on this issue.
Rolf Turner
On 14/11/2008, at 10:44 AM, Edna Bell wrote:
Dear R gurus:
Here is the following from Montgomery's Design and Analysis of
Experiments, 5th edition.
str(rout1.df)
'data.frame': 16 obs
. Don't use F. It causes trouble.
(Note: You *cannot* have spurious objects name FALSE (and not equal
to FALSE) hanging around; R won't let you. That's why you use FALSE and
not F.)
cheers,
Rolf Turner
401 - 500 of 2038 matches
Mail list logo