Re: [R] Normal plot and xyplot side by side in one figure?

2009-06-18 Thread Bert Gunter
= 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

Re: [R] Using plotmath expressions in lattice key text

2007-11-16 Thread Bert Gunter
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

Re: [R] display basename

2007-11-20 Thread Bert Gunter
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

Re: [R] display basename

2007-11-20 Thread Bert Gunter
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

Re: [R] How can I save a plot ?

2007-11-21 Thread Bert Gunter
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

Re: [R] Packages - a great resource, but hard to find the right one

2007-11-25 Thread Bert Gunter
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

Re: [R] MS Excel Data

2007-11-29 Thread Bert Gunter
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

Re: [R] Set panel background color in a multi-group trellis plot

2007-11-30 Thread Bert Gunter
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

Re: [R] using color coded colorbars for bar plots

2007-12-03 Thread Bert Gunter
... 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

Re: [R] Is R portable?

2007-12-04 Thread Bert Gunter
, 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

Re: [R] Learning to do randomized block design analysis

2007-12-04 Thread Bert Gunter
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

Re: [R] Fitting large titles in a plot

2007-12-06 Thread Bert Gunter
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

Re: [R] R and Excel Interface

2007-12-12 Thread Bert Gunter
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

Re: [R] Principal Components Analysis

2007-12-27 Thread Bert Gunter
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

Re: [R] Conditionally incrementing a loop counter

2007-12-27 Thread Bert Gunter
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

Re: [R] refering to variable names in lm where the variable name isin another variable

2007-12-31 Thread Bert Gunter
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

Re: [R] Symbolic substitution in parallel; use infinity symbol?

2007-12-31 Thread Bert Gunter
] - 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

Re: [R] retaining formatting when converting a vector to amatrix/data.frame?

2008-01-03 Thread Bert Gunter
(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

[R] Infinite recursion error that I don't understand

2008-01-03 Thread Bert Gunter
: 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

Re: [R] Infinite recursion error that I don't understand

2008-01-03 Thread Bert Gunter
: 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

Re: [R] Suggestion on how to make permanent changes to a singleobject in a list?

2008-01-03 Thread Bert Gunter
() [[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

Re: [R] subsetting

2008-01-04 Thread Bert Gunter
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

Re: [R] slow access to matrix dimnames

2008-01-04 Thread Bert Gunter
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

Re: [R] moving or running average

2008-01-07 Thread Bert Gunter
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

Re: [R] lattice color problem with symbols: bug?

2008-01-11 Thread Bert Gunter
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

Re: [R] How to calculate the mean of all values in a list or dataframe

2008-01-11 Thread Bert Gunter
?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

Re: [R] Error in solve.default message using lme

2008-01-14 Thread Bert Gunter
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

Re: [R] things that are difficult/impossible to do in SAS orSPSSbut simple in R

2008-01-17 Thread Bert Gunter
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

Re: [R] Res: vector generation

2008-01-17 Thread Bert Gunter
. 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

Re: [R] recoding one variable into another - but differently fordifferent cases

2008-01-22 Thread Bert Gunter
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

Re: [R] Using bquote: question

2008-01-24 Thread Bert Gunter
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

Re: [R] Using bquote: question

2008-01-24 Thread Bert Gunter
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

Re: [R] plotting gridlines

2008-01-25 Thread Bert Gunter
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

Re: [R] How to find maximum values on the density function of arandom variable

2009-03-12 Thread Bert Gunter
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

Re: [R] function question

2009-03-19 Thread Bert Gunter
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

Re: [R] function question

2009-03-19 Thread Bert Gunter
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

Re: [R] variance/mean

2009-03-23 Thread 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

Re: [R] very fast OLS regression?

2009-03-25 Thread Bert Gunter
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

Re: [R] programming creating different functions in a loop

2009-03-26 Thread Bert Gunter
. 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

Re: [R] programming creating different functions in a loop

2009-03-26 Thread Bert Gunter
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

Re: [R] Passing parameters to R script in Rgui

2009-03-27 Thread Bert Gunter
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

Re: [R] nls, convergence and starting values

2009-03-27 Thread Bert Gunter
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

Re: [R] Matrix max by row

2009-03-29 Thread Bert Gunter
: 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

Re: [R] [.data.frame and lapply

2009-03-29 Thread Bert Gunter
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

Re: [R] bucketing observations

2009-04-06 Thread Bert Gunter
?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

Re: [R] Sequences

2009-04-07 Thread Bert Gunter
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

Re: [R] Tinn-R pdf()

2009-04-08 Thread Bert Gunter
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

Re: [R] MLE for bimodal distribution

2009-04-08 Thread Bert Gunter
, 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

Re: [R] Anova interaction not tested

2009-04-09 Thread Bert Gunter
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

Re: [R] Coloring X and Y axis

2009-04-10 Thread Bert Gunter
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

Re: [R] Building GUI for custom R application

2009-04-14 Thread Bert Gunter
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

Re: [R] Forcing the extrapolation of loess through the origin

2009-04-14 Thread Bert Gunter
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

Re: [R] excluding a column from a data frame

2009-04-15 Thread Bert Gunter
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

Re: [R] How to extract p-value into a table from a groups of t test

2009-04-15 Thread Bert Gunter
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

Re: [R] Extending a vector to length n

2009-04-16 Thread Bert Gunter
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

Re: [R] The assign(paste(...,i),...) idiom

2009-04-20 Thread Bert Gunter
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

Re: [R] Fitting linear models

2009-04-20 Thread Bert Gunter
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

Re: [R] My surprising experience in trying out REvolution's R

2009-04-21 Thread Bert Gunter
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

Re: [R] plot confidence intervals as shaded band

2009-04-22 Thread Bert Gunter
... 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

Re: [R] rbind data frames stored in a list

2009-04-22 Thread Bert Gunter
?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

Re: [R] Floating simulation error

2009-04-23 Thread Bert Gunter
?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

Re: [R] Multiple Imputation in mice/norm

2009-04-27 Thread Bert Gunter
, 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

Re: [R] Generate ramified structures

2009-04-27 Thread Bert Gunter
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

Re: [R] FW: logistic regression

2008-09-28 Thread Bert Gunter
...) 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

Re: [R] 3 curves / 1 plot

2008-10-02 Thread Bert Gunter
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

Re: [R] t.test() on a list

2008-10-02 Thread Bert Gunter
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

Re: [R] Creating GUIs for R

2008-10-07 Thread Bert Gunter
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

Re: [R] panel.groups: use group.number to define colors

2008-10-07 Thread Bert Gunter
, 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

Re: [R] lme and lmer df's and F-statistics again

2008-10-07 Thread Bert Gunter
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

Re: [R] does any r package work with mcfadden type discrete choicemodels?...

2008-10-14 Thread Bert Gunter
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

Re: [R] Loop avoidance in simulating a vector

2008-10-16 Thread Bert Gunter
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

Re: [R] Exclude rows in table

2008-10-21 Thread Bert Gunter
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

Re: [R] Problem in Cart

2008-10-22 Thread Bert Gunter
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

Re: [R] Think I'm sure, but confirm: lme4 vs. nlme

2008-10-22 Thread Bert Gunter
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

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread Bert Gunter
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:

Re: [R] Dose Profile

2008-10-28 Thread Bert Gunter
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

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread Bert Gunter
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

Re: [R] smooth function advice

2008-10-31 Thread Bert Gunter
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

Re: [R] How to generate a new factor variable by two other factor variables

2008-11-04 Thread Bert Gunter
... 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

Re: [R] response surfaces

2008-11-04 Thread Bert Gunter
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

Re: [R] Efficient way to fill a matrix

2008-11-05 Thread 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

[R] esoteric inconsistency -- intended or not?

2008-11-07 Thread Bert Gunter
23 svn rev46779 language R version.string R version 2.8.0 Patched (2008-10-23 r46779) Cheers, Bert Gunter __ R-help@r

Re: [R] Variable passed to function not used in function in select=... in subset

2008-11-11 Thread Bert Gunter
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

Re: [R] Error with lapply

2008-11-20 Thread Bert Gunter
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

Re: [R] meaning of warning messages

2009-02-09 Thread Bert Gunter
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

Re: [R] OT: A test with dependent samples.

2009-02-10 Thread Bert Gunter
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

Re: [R] OT: A test with dependent samples.

2009-02-10 Thread Bert Gunter
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

Re: [R] Tilde in boxplots

2009-02-12 Thread Bert Gunter
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

Re: [R] Extending each element in a list, or rbind()-ing arrays of different length without recycling

2009-02-12 Thread Bert Gunter
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

Re: [R] question about tool collisions

2009-02-13 Thread Bert Gunter
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

Re: [R] Outlier Detection for timeseries

2009-02-16 Thread Bert Gunter
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

Re: [R] Applying functions to partitions

2009-02-16 Thread Bert Gunter
... 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

Re: [R] why results from regression tree (rpart) are totallyinconsistent with ordinary regression

2009-02-23 Thread Bert Gunter
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

Re: [R] bigest part of vector

2009-02-24 Thread Bert Gunter
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

Re: [R] bigest part of vector

2009-02-24 Thread Bert Gunter
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

Re: [R] User manual or tutorial for GAM package

2009-02-26 Thread Bert Gunter
-- 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

Re: [R] Inference for R Spam

2009-03-04 Thread Bert Gunter
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

Re: [R] Is there any method to identify the distribution of a given dataset?

2009-03-09 Thread Bert Gunter
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

Re: [R] 2 Simple Lattice Plot Questions

2009-03-10 Thread Bert Gunter
... 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

Re: [R] Line color based on data values?

2008-11-25 Thread Bert Gunter
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]

<    1   2   3   4   5   6   7   8   9   10   >