= number of columns and rows.
Probably the best approach to this in general is to do everything in lattice
and use grid viewports to position the plots. But I admit this requires
climbing a learning curve.
HTH.
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help
Thanks Deepayan. I knew it was simple ... (but I doubt I would ever have
figured it out).
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
Sent: Thursday, November 15, 2007 8:02 PM
To: Bert Gunter
Cc: R Help
Subject: Re
it sloppily.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of mysimbaa
Sent: Tuesday, November 20, 2007 2:30 PM
To: r-help@r-project.org
Subject: [R] display basename
Hello everybody,
I'm sorry for asking
Of course. I was just putting it inline to be as minimalist as possible...
I'm in a minimal mood today...
-- Bert
Bert Gunter
Nonclinical Statistics
7-7374
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Henrik Bengtsson
Sent: Tuesday, November 20
Maura:
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Maura E Monville
Sent: Wednesday, November 21, 2007 12:46 PM
To: [EMAIL PROTECTED]
Subject: [R] How can I save a plot ?
I recently installed R 2.6 on Linux/SuSE
When I was running the previous
style, to give a much wider thesaurus -- to quote another
pertinent comment -- to help in searching? It certainly works pretty well
for Google (and for my research this morning on slow cookers).
-- Bert
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business
Better yet! Thanks Gabor.
-- Bert
-Original Message-
From: Gabor Grothendieck [mailto:[EMAIL PROTECTED]
Sent: Thursday, November 29, 2007 12:10 PM
To: Bert Gunter
Cc: r-help@r-project.org
Subject: Re: [R] MS Excel Data
On Nov 29, 2007 3:01 PM, Bert Gunter [EMAIL PROTECTED] wrote
Read ?panel.superpose carefully and pay attention to the comments about the
group.number argument passed to panel.groups. You can test the value of this
argument and adjust the plot accordingly.
Bert Gunter
Genentech Nonclinical Statistics
7-7374
-Original Message-
From: [EMAIL
... but the best option is not to do this kind of technicolor extravaganza
at all!
See ?dotplot (in lattice) and ?dotchart for better alternatives.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of hadley
, but that is
inevitable (if R was told it should use a cyrillic character set for
Russian, it won't automatically switch to French when the flash drive is
stuck into a French computer).
If I am wrong about any of this, **PLEASE CORRECT**
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From
of tests make sense for each.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of T.K.
Sent: Tuesday, December 04, 2007 2:46 PM
To: Zembower, Kevin
Cc: r-help@r-project.org
Subject: Re: [R] Learning to do
Try This:
plot(A,main=paste(This is my really long title and,\n,it's so long that
I can see
just about half of it., sep = ))
-- Bert Gunter
Genentech
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
RSiteSearch(RExcel)
Please do at least start by using already available R information resources.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of [EMAIL PROTECTED]
Sent: Wednesday, December 12, 2007 1:12 PM
Are there lots of missing values in the data? If so, my guess would be that
na.omit(h1) leaves you with no data. Have you checked this?
Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of SNN
Sent: Wednesday, December 26, 2007 2:49 PM
Please read the posting guide and provide a simple reproducible example as
it asks you to. Very likely a loop is not even needed.
-- Bert Gunter
Genentech, Inc.
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Mike Jones
Sent: Thursday, December 27, 2007
The semantic gymnastics below aren't necessary (nor are necessarily they a
bad idea). This is basically FAQ 7.21. So, e.g.
z - c(x1,x2)
mod - lm(y ~ get(z[1]) + get(z[2]))
Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Charilaos
] - leftlab
axlab[length(axlab)] - rightlab
axis(1, at=1:3, labels=axlab, line=2.25, lty=0)
Cheers,
Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of John Maindonald
Sent: Sunday, December 30, 2007 3:45 PM
To: [EMAIL
(mydata,ncol=2))
for(i in seq(mydf))class(mydf[[i]]) - class(mydata)
Cheers,
Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Thomas Pujol
Sent: Thursday, January 03, 2008 9:02 AM
To: [EMAIL PROTECTED]
Subject: [R] retaining formatting
: evaluation nested too deeply: infinite recursion /
options(expressions=)?
## I don't have a clue. Any words of wisdom would be appreciated.
Bert Gunter
Genentech
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
: Thursday, January 03, 2008 10:57 AM
To: Bert Gunter
Cc: r-help@r-project.org
Subject: Re: [R] Infinite recursion error that I don't understand
The error is in printing ...
foo - data.frame(I(z))
works. It is not hard to create a data frame you cannot print, and this
is one. It could be considered
()
[[5]][[3]]
list()
Cheers,
Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Charilaos Skiadas
Sent: Thursday, January 03, 2008 1:46 PM
To: Peter Waltman
Cc: r-help@r-project.org
Subject: Re: [R] Suggestion on how to make permanent
See FAQ 7.12 and the R Language Definition manual on scoping.
R uses lexical scoping which means here that s does not exist in the
enclosing environment of inner(), which is the global environment.
Bert Gunter
Genentech Nonclinical Statistics
7-7374
-Original Message-
From: [EMAIL
Why do you want to do this? X[,cname] gives the column named cname of
matrix X (as a vector, unless drop=FALSE). The $ operator on data frames is
essentially equivalent to this, anyway (see ? Extract).
-- Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL
You may also wish to consider runmed(), which is a more robust running
median, and is automatically available as part of the stats package.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Jared O'Connell
are identical and unfilled.
sessionInfo()
R version 2.6.1 Patched (2007-12-03 r43574)
i386-pc-mingw32
Bert Gunter
Genentech Nonclinical Statistics
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
?unlist
mean(unlist(test))
Bert Gunter
Genentch Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of B. Bogart
Sent: Friday, January 11, 2008 11:29 AM
To: r-help@r-project.org
Subject: [R] How to calculate the mean of all values
What it probably means is that you can't make a silk purse out of a sow's
ear: all those missings in your your data do not allow you to estimate the
model(s) that you've specified.
If you don't know what this means, I suggest you seek the help of a local
statistician.
--
Bert Gunter
Genentech
Thanks Hadley. Witty! Profound! Concise! I think this is definitely a
Fortunes candidate.
-- Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of hadley wickham
Sent: Thursday, January 17, 2008 9:56 AM
To: Wittner
. Since a matrix is just a vector with a dim
attribute, you're done; but if you want to show it as a vector, just use,
e.g.
as.vector(outer(z,z,/))
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Milton Cezar
Try this:
z1 - c(1,-1,0,99)[Data$A]
z2 - c(-1,1,0,99)[Data$A]
Data$new - ifelse(Data$B == 1, z1,z2)
(This does not generalize, however, and assumes that the values you gave are
exact. Gabor's approach would work more generally)
-- Bert Gunter
Genentech Nonclinical Statistics
-Original
and outs of computer languages). But I do appreciate
that the folks who develop and maintain the core R language do such a
careful job of documenting things that even dolts like me can figure it out
if we make the effort.
Cheers,
Bert Gunter
-Original Message-
From: [EMAIL PROTECTED
to the plot.default method that will be called by plot.formula to do the
plot.
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Schools Statistics Poster Competition
Sent: Thursday, January 24, 2008 5:58 PM
To: R-help mailing list
Subject: [R] Using
Or you can use lattice's xyplot(..., type = c(p,g))
?panel.xyplot provides details
-- Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of John Kane
Sent: Friday, January 25, 2008 9:49 AM
To: John Lande; r-help@r-project.org
Subject
estimation or
perhaps talking with your local statistician?
-- Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Mike Lawrence
Sent: Thursday, March 12, 2009 5:40 PM
To: g...@ucalgary.ca
Cc
Edna:
Please make at least a minimal effort to answer such questions before
posting.
is.vector(function(x)x) ## FALSE
as.vector(function(x)x) ## try it
or even
is.vector(plot)
-- Bert Gunter
Genentech
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r
Not surprising at all -- expected!
is.function(vector) ##TRUE
-- Bert
-Original Message-
From: Wacek Kusnierczyk [mailto:waclaw.marcin.kusnierc...@idi.ntnu.no]
Sent: Thursday, March 19, 2009 1:03 PM
To: Bert Gunter
Cc: Edna Bell; R help
Subject: Re: [R] function question
Bert Gunter
Inline Below.
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Wacek Kusnierczyk
Sent: Sunday, March 22, 2009 2:16 AM
To: rkevinbur...@charter.net
Cc: r-help@r
a tradeoff between
computational accuracy (especially for ill-conditioned matrices) and speed
which your remarks seem to neglect.
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org
.
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Greg Snow
Sent: Thursday, March 26, 2009 9:25 AM
To: Florin Maican; r-help@r-project.org
Subject: Re: [R] programming
Good point, Greg. So to handle i as a vector, you'd probably want something
like
f - function(x,i)outer(x,i,+)
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf
You are misusing source.
Write a function to do what you want. An Introduction to R documents how.
Have you read it?
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org
such a
thing? category.
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Patrick Giraudoux
Sent: Friday, March 27, 2009 12:39 PM
To: r-h...@stat.math.ethz.ch
Cc: Francis
:
A data frame is a list of variables of the same length with unique row
names, given class 'data.frame'.
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Wacek Kusnierczyk
Sent
Folks:
I do not wish to agree or disagree with the criticisms of either the speed
or possible design flaws of [. But let's at least see what the docs say
about the issues, using the simple example you provided:
m = matrix(1:9, 3, 3)
md = data.frame(m)
md[1]
# the first column
?cut
?quantile (perhaps, to define the breaks)
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Dan Dube
Sent: Monday, April 06, 2009 12:45 PM
To: r-help@r-project.org
not provide all of the details, either. As a result, I still occasionally
get bitten by these indexing subtleties (which, of course, may just be due
to my dummheit).
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r
I **think** the explanation is:
When sending code to R through TINN-R, TINN-R **sources** the code into R.
Sourced code is **not** autoprinted. Indeed, ?source explicitly tells you
this.
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r
, of course.
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Ben Bolker
Sent: Wednesday, April 08, 2009 12:47 PM
To: r-help@r-project.org
Subject: Re: [R] MLE
Gabe:
Don't be silly. lme is not appropriate here --you have only one stratum!.
Use lm and you'll see what's going on. It also looks like you should do some
reading up on linear models. VR's MASS or Peter Dalgaard's INTRO to
Statistics with R might be places to start.
Incidentallyy, Version
plot(..., axes=FALSE)
axis(...)
?plot; ?axis
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Gundala Viswanath
Sent: Friday, April 10, 2009 8:26 AM
To: r-h
complex functionality is difficult. Computer scientists, web
designers, etc. have whole curricula devoted to this.
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf
Below.
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of jimm-pa...@gmx.de
Sent: Tuesday, April 14, 2009 10:08 AM
To: r-help@r-project.org
Subject: [R] Forcing
How about:
xx[,-match(x2,names(xx))]
or
xx[,names(xx) != x2]
etc.
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Erin Hodgess
Sent: Tuesday, April 14, 2009 10:39
Everything in R is an object! ?str
z - t.test(x,y,...)
str(z)
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Jun Shen
Sent: Wednesday, April 15, 2009 1:30 PM
To: r
Folks:
Not to be picky, but depending on exactly what's meant by attributes, I
think it's impossible:
x - 1:6
attr(x,length) - length(x)
xx - rep(x,length=10)
xx[7:10] - NA
attr(xx,length)
NULL
attr(x,length)
[1] 6
So attributes aren't preserved. The whole point of object orientation and
Perhaps a good R WIKI topic, though there's no assurance that info there
would be consulted, of course...
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of baptiste
Is this homework? If so, you need to read the text and/or class notes more
carefully.
-- Bert Gunter
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Vemuri, Aparna
Sent: Monday, April 20, 2009 4:26 PM
To: r-help@r-project.org
dangerous to generalize from an example or two (as I
recently experienced to my own chagrin).
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Jason Liao
Sent
... or ?polygon in base graphics
-- Bert Gunter
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Dieter Menne
Sent: Wednesday, April 22, 2009 8:23 AM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] plot confidence intervals as shaded
?do.call
... as in do.call(rbind,X)
(A very useful and powerful feature os the S language: Computing on the
language. See VR's S Programming for an informative discussion (there may
well be others, of course).
-- Bert
-Original Message-
From: r-help-boun...@r-project.org
?traceback
options(error) ( ?options)
?debug
?try
?tryCatch
R has facilities to help you with such problems. Please use them -- and then
repost with more specific info if you still cannot solve it.
-- Bert Gunter
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help
, of course, but
perhaps one worth being reminded of in such situations.
As always, thanks for **your** knowledgeable summary of exactly these
matters, Frank.
-- Bert Gunter
Genetech
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf
Don't know, but try checking the Graphical Models Task View on CRAN.
-- Bert Gunter
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Talita Perciano
Sent: Monday, April 27, 2009 1:38 PM
To: r-help@r-project.org; R help
Subject
...)
Cheers to all,
Bert Gunter
DEFINITELY MY OWN OPINIONS HERE!
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of David Winsemius
Sent: Saturday, September 27, 2008 5:34 PM
To: Darin Brooks
Cc: [EMAIL PROTECTED]; [EMAIL PROTECTED]
Subject: Re: [R] FW: logistic
imho better: ?xyplot using the groups argument.
-- Bert
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Kingsford Jones
Sent: Thursday, October 02, 2008 11:17 AM
To: Michel PETITJEAN
Cc: r-help@r-project.org
Subject: Re: [R] 3 curves / 1 plot
I am sure you will get helpful answers. I am almost as sure that you
shouldn't be doing this. I suggest you consult with your local statistician.
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Gang Chen
Sent: Thursday, October 02, 2008
Seek and ye shall find ... Check the RGUI's link on the other web page on
CRAN.
If you are on Windows, there is some simple built-in GUI functionality.
?winMenuAdd, ?select.list and the links therein will get you started there.
Cheers,
Bert Gunter
-Original Message-
From: [EMAIL
, it's values are those of
trellis.par.get(superpose.symbol) I believe. Also, group.number appears to
be undefined in your code.
HTH
Cheers,
Bert Gunter
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
is, as appears to be the case here, that even the most technical
aspects of the discipline can be made manifest to anyone with half a brain
and a stat 101 course under their belt.
I think we owe Doug Bates a little more respect than that!
Cheers,
Bert Gunter
__
R
Please use R's search tools before posting:
?RsiteSearch
RSiteSearch(mcfadden,restr=func)
RSiteSearch(discrete choice model,restr=func)
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Richard Palmer
Sent: Tuesday, October 14, 2008 2:01
of a
difference to be worth the effort. Random number generation is so efficient
in R that avoiding loops rarely matters.
Also see ?replicate for a way to perhaps write cleaner code (but still using
hidden interpreted loops).
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL
to exclude each row where Net=0. Could this be done in a simple way?
Yes. help([)
Also see An Introduction to R. Do your homework before posting, please.
-- Bert Gunter
Thanks in advance
Jim
A
Business.Unit Event1 Net Date
1General Fraud 170.000 2006-01-01
2
CART is a commercial package and not part of R.
R has several packages that do various kinds of regression and
classification trees. Try:
RSiteSearch(Classification Tree,restr=func)
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of [EMAIL
and informative manner[1], I would be delighted. Thanks!
-- I didn't consider your questions naive or impertiment. I leave it to you
to judge whether I have been friendly or helpful.
Cheers,
Bert Gunter
[1] Or, translated from craven, please don't ignite a flame war on my
no-account account
Regards
David:
IMHO, this is beginning to be commercial abuse of this list. Please no more
advertising.
(Others may have different opinions).
-- Bert
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of [EMAIL PROTECTED]
Sent: Tuesday, October 28, 2008 9:33 AM
To:
Why don't you ask your someone at work to help, as you appear to have
local resources handy?
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Laura Bonnett
Sent: Tuesday, October 28, 2008 9:33 AM
To: r-help@r-project.org
Subject: [R] Dose
David et. al:
My sentiments are the same, but my apology to you (and I'm glad to learn
that it's not you!).
-- Bert
_
From: David Smith [mailto:[EMAIL PROTECTED]
Sent: Tuesday, October 28, 2008 11:09 AM
To: Bert Gunter
Subject: Re: [R] Recommended R books by XLSolutions
RSiteSearch(ewma,rest=func)
would have told you there already exists an ewma implementation in the qcc
package.
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of [EMAIL PROTECTED]
Sent: Friday, October 31, 2008 10:04 AM
To: [EMAIL
... and also note the drop argument. Adding drop = TRUE to the call is
usually what is wanted, e.g.
df$factorC - with(df, interaction(factorA, factorB), drop = TRUE )
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Marc Schwartz
Sent
Take a look at the AlgDesign package.
Of course, lm(), or whatever be your favorite version of multiple
regression, will fit the results.
Obviously (I think) you'll have to write your own plotting code or perhaps
make use of anything available in the other packages Ben mentioned.
-- Bert Gunter
])] - foo[,3]
See ?[ for details.
-- Bert Gunter
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of stephen sefick
Sent: Wednesday, November 05, 2008 2:52 PM
To: Philipp Pagel
Cc: r-help@r-project.org
Subject: Re: [R] Efficient way to fill a matrix
#reshape
23
svn rev46779
language R
version.string R version 2.8.0 Patched (2008-10-23 r46779)
Cheers,
Bert Gunter
__
R-help@r
Ummm... as today is still Armistice day (in my time zone, anyway), maybe we
should call a truce and end this flame war...
Cheers,
Bert
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Berwin A Turlach
Sent: Tuesday, November 11, 2008 9:31 AM
To: Duncan
automatically; other arguments in lapply are passed in as **additional**
arguments to fn, of which there are none here.
Compare:
fn - function(i,a)i^2+a
lapply(1:4,fn) ## error
lapply(1:4,fn,5)
Cheers,
Bert
Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED
into R. Totally innocuous. But just a guess
without the info as above.
-- Bert Gunter
Genentech
__
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
sample of 73? binom.test(12,72) gives us .088
to .27 for an exact 2 sided interval (and a P value of 2.2e-16 for the null
= 0).
Seems rather convincing -- and simple -- to me!
-- Bert Gunter
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf
a prior guess if
you like for a Null. Which is exactly the point that Marc Schwartz made --
that is, that the data are probably completely useless to answer the
question of interest because the researcher messed up the design.
-- Bert Gunter
-Original Message-
From: markle
boxplot is an (S3) generic function: ?UseMethod
1. Uses the default method.
2. Uses the formula method, since biomass~clipping is a formula: ?~ ;
?formula; ?boxplot.
See also an Introduction to R where these matters are also explained.
-- Bert Gunter, Genentech
-Original Message
At the risk of beating the decaying skeletal remains of the stone dead
horse, a one-liner:
t(apply(mat,1,function(x)c(unique(x),rep(NA,sum(duplicated(x))
(probably more efficient as a 2-liner that calls duplicated/unique only
once, though)
-- Bert Gunter, Genentech
-Original Message
See the warn.conflicts arguments of require() and library(). People
thought about this a long time ago.
-- Bert Gunter
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Carl Witthoft
Sent: Friday, February 13, 2009 4:05 PM
To: r-help
of the intimation that there is an objective
statistical meaning to this term, which there is not) for each.
Statistical algorithms cannot replace careful thinking. Sorry about that.
-- Bert Gunter
Genentech
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org
...
I suppose the clean way to do this would be to define a cartesian product of
two factors with the induced lexicographic order (is there a standard
function for doing this?):
Of course. ?interaction.
-- Bert Gunter
Genentch
__
R-help@r
You did not read the tree graph correctly. Mortality is **not** positively
related to incidence. You're reading the tree backwards. Read the output
of summary() on your rpart fit object for clarity.
-- Bert Gunter
Genentech
-Original Message-
From: r-help-boun...@r-project.org
Nothing wrong with prior suggestions, but strictly speaking, (fully) sorting
the vector is unnecessary.
y[y quantile(y, 1- p/length(y))]
will do it without the (complete) sort. (But sorting is so efficient anyway,
I don't think you could notice any difference).
-- Bert Gunter
Genentech
into the psychology of
surprise that I don't have).
-- Bert
-Original Message-
From: macra...@gmail.com [mailto:macra...@gmail.com] On Behalf Of Stavros
Macrakis
Sent: Tuesday, February 24, 2009 2:12 PM
To: Peterko
Cc: Bert Gunter; r-help@r-project.org
Subject: Re: [R] bigest part of vector
On Tue, Feb 24
-- or perhaps even the section on Additive Models in the latest edition of
VR's MASS (still a useul book to have in one's library, IMHO, although,
like me, it's getting grayer)
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r
process.
This is much much more than you intimate.
Cheers to all,
Bert Gunter
Genentech Nonclinical Biostatistics
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R
Below. Brief summary is: You **need** to consult a statistician. You know
far too little statistics to do statistical analysis on your own.
-- Bert
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun
... but in general
?strip.default
will explain these matters. In particular, note the style argument.
Cheers,
Bert Gunter
Genentech Nonclinical Biostatistics
650-467-7374
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Elena
FWIW
Indeed! And IMHO such a nice example of the power and beauty of the R
language.
-- Bert
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Greg Snow
Sent: Tuesday, November 25, 2008 10:00 AM
To: [EMAIL PROTECTED]; r-help@r-project.org
Subject: Re: [R]
401 - 500 of 5096 matches
Mail list logo