Thanks for the response. I ***think*** I'm making a bit of
progress
On 29/07/2008, at 10:14 AM, Douglas Bates wrote:
On Sun, Jul 27, 2008 at 9:06 PM, Rolf Turner
<[EMAIL PROTECTED]> wrote:
What I *don't* understand is the correlation structure of
On 30/07/2008, at 6:12 AM, dxc13 wrote:
useR's,
I am trying trying to find out if there is a faster way to do a
certain
computation. I have successfully used FOR loops and the apply
function to
do this, but it can take some time to fully compute, but I was
wondering if
anyone may know
27;' searches ``backwards'' --- with the convention
that the most recent command is the ``top'' of the file.
And then when you've found the command that you want, you can edit it
(with vi syntax) before pressing to re-issue the command.
No h
a box!!!
But be that as it may, you can eliminate the box
by putting ``bty="n"'' in the call to legend.
I have no idea what you mean by ``box-shadow''.
cheers,
Rolf Turner
###
r Bloomfield's
book ``Fourier Analysis of Time Series --- An Introduction'' (2nd
ed.), Wiley Series in
Probability and Statistics, 2000.
cheers,
Rolf Turner
P.S. BTW doing ``library(stats)'' is silly;
ion of what ``hooks'' are)
so I won't go there.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
e n is the series length.
I.e. the fft in R (and in S/Splus) does not apply any normalizing
factor,
so that the inverse transform only ``inverts'' up to a constant
multiple.
cheers,
Rolf Turner
On 5/08/2008, at 1:31 AM, Frank E Harrell Jr wrote:
I wonder why we don't just use the exact nonparametric confidence
interval for the median, which is just as easy to compute. Also,
it will be asymmetric if the data are skewed, as it should be.
The boxplot.stats() fun
create a pdf file ``.pdf''.
(It'll say ``Hmm ... looks like a package'' before it starts
processing.)
With a table of contents and an index and ***everything***!!! :-)
cheers,
Rolf Turner
###
On 8/08/2008, at 8:39 AM, carol white wrote:
Hi,
How to save an R object for example a matrix or vector and not all
objects created in a session (which is usually performed by
save.image or q("yes"))?
Two approaches:
(1) Clean up your workspace so that it contains only the objects you
x27;t using
a computer,
but get R to do the arithmetic for you.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
ended to be
used, rather than
trying to force your own inappropriate paradigm upon R. You'll find
it more efficient
to do things R's way rather than trying to force R to do things your
way. R gets
it right. [Almost] always!
cheers,
ily be ``generalized'' and bundled up in a function if
so desired.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and conf
It is safer and takes 6 fewer key strokes to use ``if(X)''. So why not
do that, eh?
cheers,
Rolf Turner
P. S. The example that triggered this post actually used ``if(X==T)''.
This is doubly dangerous since ``T'' can be overwritten. If it was
There was a thread on this recently. Solutions were posted to allow you
to join interpolated points to ``real'' ones using a different line
type.
See
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/127669.html
cheers,
Rolf Turner
On 18/08/2008, at 2:
m to their neighbours by dotted lines (where real
points
are joined to each other by solid lines).
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{
On 19/08/2008, at 10:15 AM, hadley wickham wrote:
On Mon, Aug 18, 2008 at 3:54 PM, Rolf Turner
<[EMAIL PROTECTED]> wrote:
On 19/08/2008, at 1:04 AM, stephen sefick wrote:
The real data are counts of aquatic insects at distinct locations
on a
river continuum on the way down a rive
Read the help for load(); note the ``Value'' section.
cheers,
Rolf Turner
On 19/08/2008, at 11:51 AM, Juan Manuel Barreneche wrote:
Hello list! i write because i'm having a strange problem with the
"save"
function, here an example:
x <
They are *NOT* equal to x, x^2, and x^3.
Try:
> m <- cbind(1,poly(1:10,3))
> round(t(m)%*%m,digits=8)
If you don't want orthogonal polynomials, use e.g. poly
(x,degree=3,raw=TRUE).
cheers,
ganappi
foregoing code by doing
u <- matrix(c(1201,688),2,2,byrow=TRUE)
v <- c(44308,4371)
m <- u/v
But as I indicated, that's probably not what you really want to do.
You need to get
yo
On 20/08/2008, at 4:19 PM, suman Duvvuru wrote:
I would like to know how to convert a string with characters to all
uppercase or all lowercase? If anyone could let me know if there
exists a
function in R for the conversion, that will be very helpful.
?tolower
#
e of .Random.seed be
made equal to that read in from the file?
***PLEASE*** don't tell me to use set.seed() instead, or something
like that.
I ***know*** about set.seed() --- I wasn't born yesterday, y'know!
Please just
accept that I want to do what I
On 21/08/2008, at 10:21 AM, Charles C. Berry wrote:
On Thu, 21 Aug 2008, Rolf Turner wrote:
For reasons that are best known to myself [ ;-) ] I have a value
of .Random.seed
saved (via dput()) in a file ``.Random.seed.save''.
In my .Rprofile I have the lines:
.Random.se
On 21/08/2008, at 3:24 PM, Charles C. Berry wrote:
On Thu, 21 Aug 2008, Rolf Turner wrote:
But why does this happen on *my* system, and not on Chuck's???
Because I goofed by rerunning ALL the lines in .Rprofile rather
than just the last two.
When I run just the
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz
to R
documentation?
Yes.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.or
?RSiteSearch
On 26/08/2008, at 1:43 PM, michal33 wrote:
Hi there,
I am trying to run LSD.test(model)
I used the following commands:
attach(model)
m1<- glm(ttl.m ~ site+year, family=quasipoisson, data= model)
df<-df.residual(m1)
MSerror<-deviance(m1)/df
The following command did not work:
co
control)) # You do not need to specify the
method explicitly.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
___
See FAQ 7.22
On 27/08/2008, at 11:07 AM, John Sanabria wrote:
Hi,
I have the following script:
t.R ---
grafica <- function() {
v <- read.csv('preprocessed/komolongma.ece.uprm.edu.active',sep=',')
x <- as.ts(v$active)
bitmap(file="output.png")
densityplot(~x,col='blue',main='Densi
On 27/08/2008, at 12:32 PM, michal33 wrote:
Thanks for the answers, this one worked:
library(Agricolae)
No, it couldn't have. There is no such package as ``Agricolae''.
There ***is*** a package called ``agricolae''. :-)
cheers,
You need to do:
system("perl presentPerformance.pl",intern=TRUE)
It does pay to read the help, you know.
cheers,
Rolf Turner
On 28/08/2008, at 8:36 AM, kevinchang wrote:
Dear R users,
I am trying to call a Perl subrout
at have you done
to offend the gods? :-)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.or
NA
FALSE 6 NA
Consequently the result of ifelse(a>1,6:7,NA) is c(NA,7,6,7,NA).
What's the problem?
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and co
On 28/08/2008, at 3:00 PM, Charilaos Skiadas wrote:
On Aug 27, 2008, at 10:40 PM, Rolf Turner wrote:
On 28/08/2008, at 2:02 PM, James Milks wrote:
The title says it all. Does anyone know of a way to save your
packages when you upgrade to a new version of R? This may seem
petty, but
he reshape package.
cheers,
Rolf Turner
On 18/08/2009, at 4:00 PM, rajclinasia wrote:
Hi Everyone,
I have a dataset like this
mean sd0%
25% 50%
75% 100% n
BODY TEMPERATURE 36.41099 0.4015699 35.1
vector(unname(a)))
}
E.g.:
library(spatstat)
set.seed(42)
x <- runif(20)
y <- runif(20)
i <- sample(1:20,100,TRUE)
X <- ppp(x=x[i],y=y[i],marks=rnorm(100))
Y <- foo(X)
plot(unmark(X))
plot(Y,add=TRUE,cols="red")
HTH
cheers,
Rolf Turner
P.S. Note that
e lines of the help file(s) would be useful.
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
If the poster asks a reasonable question, answer it, rather than
engaging in
bloody-minded obscurantism.
cheers,
Rolf Turner
-- Bert
Bert Gunter
Genentech Nonclinical Biostatisics
-Original Message-
From: r-help-boun...@r-project.org [mail
o bxp() that anyone
might be forgiven for thinking ``There must be a way to do what
I want; I just haven't twigged to the correct way of putting
these arguments together.'' Deliberately steering a new user
into such a misapprehension is unforgivable.
ory almost surely will
not be
saved, so it is highly probable that you are SOL. Unless you closed
the R session (and save the workspace) before the power cut. Bummer.
Good luck.
cheers,
Rolf Turner
##
A
would very likely be pretty
daunting to the OP. If he could do it, he probably wouldn't have
asked the question in the first place.
The response to which I objected was SERIOUSLY MISLEADING. And
therefore objectionable. Full stop.
cheers,
Rolf Turner
On 20/08/2009,
misleading and would have had the
OP tearing his hair out reading and re-reading the help
files and wondering what he was missing.
cheers,
Rolf Turner
One should not aim at being possible to understand, but at
being impossi
simply,
provided that one is clever enough.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
rix would be the standard errors of the parameter
estimates.
Is the pop quiz over?
You might also have asked him ``Would you like coffee with that?'' :-)
cheers,
Rolf Turner
##
's very bad graphical practice
and
will mislead the viewer with probability 1.
(2) If you bloody-mindedly ***insist*** on doing it, see:
http://wiki.r-project.org/rwiki/doku.php?id=tips:graphics-base:2yaxes
c
,
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
v), lty=1)
dev.off()
}
}
Clearly that doesn't work. I'm not sure how to make R see the iv
and dv
strings as variables. Advice?
?get
cheers,
Rolf Turner
##
Attention:\ This e-mail messag
On 27/08/2009, at 8:00 AM, Duncan Murdoch wrote:
On 26/08/2009 3:20 PM, Noah Silverman wrote:
Hi,
Is there a way to build up a vector, item by item. In perl, we can
"push" an item onto an array. How can we can do this in R?
I have a loop that generates values as it goes. I want to end up
this:
x <- list()
i <- 0
for(item in list) {
if(threshold met) {
i <- i+1
x[[i]] <-
}
}
x <- unlist(x)
cheers,
Rolf Turner
?reshape
On 28/08/2009, at 11:37 AM, Richardson, Patrick wrote:
I have a dataset that I'm trying to rearrange for a repeated
measures analysis:
It looks like:
patient basefev1 fev11h fev12h fev13h fev14h fev15h fev16h fev17h
fev18h drug
201 2.46 2.68 2.76 2.50 2.30 2.14 2
elp to you.
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/
Fortune candidate?
cheers,
Rolf Turner
On 30/08/2009 6:08 PM, Noah Silverman wrote:
Normally, I would just write a huge "for each" loop, but have read
that
is hugely inefficient with R.
On 31/08/2009, at 10:47 AM, Duncan Murdoch wrot
Many thanks to Gabor Grothendieck and William Dunlap who both
solved my problem for me, right rapidly!
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9
and x.pre are counts, rather than the raw
data, this
does not appear to make any sense. You are getting the ecdf-s of the
counts
rather than of the raw data which is a toadally different story. And
toadally
misleading. Re-think.
cheers,
Rolf Turner
On 7/09/200
symbols()-created circle to pass through points which
really
do lie on a circle you need to make the aspect ratio of the plot
equal to 1:
plot(c(-5,0,0,5), c(0,5,-5,0),asp=1)
symbols(0,0, circles=c(sqrt(25)), inches=FALSE, add=TRUE)
OM!
)
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/listinfo
ould someone please give some hints on the way to find the
derivative on the curve at some points ?
See
?smooth.spline
and
?predict.smooth.spline
cheers,
Rolf Turner
##
Attention:\
E.g.:
set.seed(42)
x <- rnorm(100)
y <- runif(120)
qqnorm(x)
points(qqnorm(y,plot=FALSE),col="red")
If you ***really*** (???) want a line graph, then do
qqnorm(sort(x),type="l")
lines(qqnorm(sort(y),plot=FALSE),col="red")
e word biology which is what I dont want *
(1) RTFM
(2) paste(a,b,sep="")
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
___
e that the components of props
may be
different from call to call (sometimes there's only x data,
sometimes there's no
ylab, etc.).
Thanks in advance for your time and help.
do.call(plot,props)
See ?do.call .
kidneys. As I may have said on this list previously, I think
fortune("machine learning") is relevant.)
cheers,
Rolf Turner
--
N
On 9/14/09 10:12 AM, Ravi Varadhan wrote:
Noah,
It may be just me - but how does "any"
On 15/09/2009, at 7:42 AM, Ted Harding wrote:
Hi Folks,
I don't often grumble,
Well Ted, that's just a flaw in your character! I guess
you can't help it! :-)
cheers,
Rolf
##
Attentio
^2 is meaningless for such models. They are fitted via maximum
likelihood, not least squares.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{drop
#x27; contain the string ``n.'' (literal ``.'')
somewhere before the string ``mytype''.
cheers,
Rolf Turner
On Tue, Sep 15, 2009 at 5:52 PM, Sebastien Bihorel
wrote:
Dear R-users,
I am trying to use the grep function to test whether a particu
27;s *test statistic* using the ``='' (which is the
``most extreme'' point of the null hypothesis, the point ``closest'' to
the alternative, i.e. the point least likely to lead one to reject
the null.
Thus confusion amongst the young and naive is minimize
b*u2 and choose a and b so that a^2 + b^2 = 1
and
(y,x2) = (y,x1). Note that ``(v1,v2)'' means the inner (dot) product
of v1 and v2.
``Choosing'' a and b involves solving a quadratic equation.
To get things in orthocomplements of things, use the Gramm-Sch
t;)
How can I suppress that output is printed?
If you don't want it printed, then why the
are you (explicitly!) using print??? Words fail me!!!
cheers,
Rolf Turner
##
Attention:\ This e-mail mes
On 22/09/2009, at 10:13 AM, Martin Batholdy wrote:
Am 21.09.2009 um 23:59 schrieb Rolf Turner:
On 22/09/2009, at 9:52 AM, Martin Batholdy wrote:
hi,
I use xtable to convert data.frames to html tables.
But when I use the print-command I always get the whole output
printed
even if I
this a homework question?
(b) Have you figured it out yet?
(c) Hint: You have spikes at +/- 40 in a range from -50 to 50. You
*want* spikes
at 10 and 90 Hz. Could it be that you haven't set your frequency
vector ``f''
quite right? :-)
cheers,
Rolf Tur
On 23/09/2009, at 12:14 PM, Henrique Dallazuanna wrote:
Try this:
treeName <- "PINE"
write.table(d, file = sprintf("C:/%sheight.txt", treeName))
Why not use the (less esoteric) paste() function?
cheers,
knowing R
> backwards and forwards, I did read ?read.table and ?raw. But your suggestion
> is not at all helpful Peter:
Etc. etc.
This sort of reaction is completely inappropriate for this list.
Peter's answer was perfectly reasonable.
If you don't like the
On 12/02/2010, at 9:39 AM, Peter Ehlers wrote:
> R-people,
>
> Duncan Murdoch's response in
>
> https://stat.ethz.ch/pipermail/r-help/2010-February/227869.html
>
> reminded me of something I had been meaning to ask.
>
> A while ago I started using
>
> for(i in seq_len(v)) {}
>
> in pr
insemius's later protestation
that his advice was *not* a version of RTFM.)
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
ht ball park. I think!
This is all posited on the distribution of the estimate of log(theta) being
normal (or ``Gaussian''). Whether this is a justifiable assumption in your
setting is questionable.
Some simulation experiments might be illuminating.
cheers,
R
On 18/02/2010, at 3:25 AM, Christine SINOQUET wrote:
> I want to unsuscribe from the R list.
> Thanks.
Go right ahead!!! Don't let us stop you!
cheers,
Rolf Turner
##
Attention:\ This e-m
having an effective
environment in
which to conduct statistical analyses, there can be no question that the
R/S/S-Plus
group win hands down.
cheers,
Rolf Turner
##
Attention:
This e-mail message is privile
eld entries are surrounded
by quotes, shouldn't they be left as character? Even if they are
all F's and T's?
Furthermore using F's and T's to represent TRUE's and FALSE's is
bad practice anyway. Since FALSE and TRUE are reserved words it
would make sense f
ued responses you shouldn't be using
the Poisson
family; the quasi family *might* be appropriate --- if you know what
you're doing.
> Can anyone offer some enlightenment?
I hope you feel enlightened.
cheers,
Rolf Turner
ces are ***two*** dimensional. It's an
array.
(b) Use abind() from the abind package.
(c) RSiteSearch("bind arrays") would have led you straight to this
but just RSiteSearch("bind") makes it the 57-th item so you probably
wouldn't have found it that way.
On 3/03/2010, at 12:14 PM, Peter Ehlers wrote:
> On 2010-03-02 15:32, Rolf Turner wrote:
>>
>> On 3/03/2010, at 11:22 AM, Luis Felipe Parra wrote:
>>
>>> Hello I have a 2x10x200 matrix and I would like to bind to it another 2x10
>>> matrix in order to
ot objects ``rel'' and ``len''
in your workspace (and these are different from the columns of ``con'')
then the attach()---detach() procedure will use the objects in your
workspace. They ``mask'' the columns of ``con''. The with() procedure
is not beset with this p
Do your own homework!
It looks pretty trivial; what's your problem?
cheers,
Rolf Turner
On 4/03/2010, at 2:39 PM, Pitmaster wrote:
>
> Hi guys... I have problem with this excersise...
>
> Consider the pressure data frame again.
>
> (a
and
(b) Between discrete and continuous distributions.
Given that you understand those differences you will see that all three
answers are correct.
cheers,
Rolf Turner
>
> for a discrete uniform distribution
scientific progress, I'd suggest that you
> ignore that sentence (and any conclusion drawn subsequent to it).
Surely a fortune candidate.
cheers,
Rolf Turner
##
Attention:
This e-mail message is pr
On 9/03/2010, at 11:17 AM, Mike Prager wrote:
> Rolf Turner wrote:
>>
>> I solved the problem by putting in a colClasses argument in my
>> call to read.csv(). But I really think that the read functions
>> are being too clever by half here. If field entries
Define clearly what you mean by ``mean slope''. Write down a mathematical
expression for this. The answer is then obvious.
No R function required.
cheers,
Rolf Turner
##
Attention:
This e-m
When you do a[!is.na(a)] you get a ***vector*** --- not a matrix.
``Obviously''!!! The non-missing values of a cannot be arranged in
a 5 x 6 matrix; there are only 26 of them. So (as my late Uncle
Stanley would have said) ``What the hell do you expect?''.
The ``trick''
On 15/03/2010, at 3:07 PM, S Ellison wrote
(in response to a very confused question about superimposing density
curves on barplots):
> So the bad news is that not a lot of what you're doing is right.
Sounds like a fortune to me! :-)
cheers,
peared.
> Error information: The instruction at "0x6abf12cf" referenced memory at
> "0x0286fff8". The memory could not be "read".
> Anybody met the similar problem? How to avoid this?
> I appreciate your suggestions.
I suggest that
On 15/03/2010, at 3:37 PM, S Ellison wrote:
> I begin to think that R needs a _mis_fortunes package...
I am tempted to nominate *that* as a fortune! :-)
cheers,
Rolf Turner
##
Attention:\ This e-m
u could use
crossdist() --- and then average these distances by using apply() on the
result.
If this doesn't answer your question, can you please clarify your question?
cheers,
Rolf Turner
##
Attentio
o that it looks like this...
> [,1] [,2] [,3]
> [,1] 43 1
> [,2] 52 2
> [,3] 36 1
>
> Does anyone know a quick and straightforward way to do this?
Does
t(mat[nrow(mat):1,])
float
On 1/07/2009, at 12:34 AM, Craig P. Pyrame wrote:
... it's probably not a good idea to submit bug
reports just because I misunderstand what R does.
Gotta be a fortune!!!
cheers,
Rolf
##
Attention:\
atrix(~x)
(where x is your factor) is pretty straightforward. How much simpler
can it get? The power of the formula concept is working for you here;
it would be hard to do without it.
cheers,
Rol
hat a more enlightening
error message gets issued. Probably not worth the effort since few
users
would ever entertain the idea of using "[<<-".
cheers,
Rolf Turner
##
Attention:\ Thi
eading and
confusing.
But if you insist --- and be it on your own head --- see
http://wiki.r-project.org/rwiki/doku.php?id=tips:graphics-base:2yaxes
Don't say I didn't warn you!
cheers,
Rolf Turner
nse
at all.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.org mailing list
https://stat.ethz.
On 8/07/2009, at 6:28 AM, Sunny wrote:
I am new with R and want do some analysis with a point vector data
file. Any
help is appreciate. Sunny
See the vignette ``Handling shapefiles in the spatstat package'' from
the
spatstat package on CRAN.
cheers,
ently ``replied'' to an R-help digest, and included
megabytes
of totally irrelevant material in your post. It took me several
minutes to delete it.
STEP ONE: LEARN HOW TO USE EMAIL!!!
(b) Do the following:
myList <- list();
for(i in n-1)
{
myList[[i]] <-c(as.nume
27;, i.e. of reading the Posting Guide.
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid...{{dropped:9}}
__
R-help@r-project.or
501 - 600 of 2104 matches
Mail list logo