Re: [R] pairing columns based on a value

2014-12-17 Thread Bert Gunter
y be more inclined than I am to help, you would certainly increase the likelihood by first doing some homework and showing us code that you tried. Although, by that time, you probably will have figured it out for yourself. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374

Re: [R] Help: The Difference Between Workspace, Script, and History

2014-12-17 Thread Bert Gunter
I suggest that you have a look at "An Introduction to R" for a better description of how things work than your current tutorial apparently gives. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledg

Re: [R] how to make this get command work?

2014-12-17 Thread Bert Gunter
?"[" Read the docs! Go thru an R tutorial. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Wed, Dec 17, 2014 at 7:24 PM, ce wrote:

Re: [R] list of lists, is this element empty

2014-12-20 Thread Bert Gunter
h the length of object [1] ## outer gets completely flummoxed, as it should! > expand.grid(z,z) Var1 Var2 1 NULL NULL 2 NULL NULL 3 NA NULL 4 NULL NULL 5 NULL NULL 6 NA NULL 7 NULL NA 8 NULL NA 9 NA NA ## and expand.grid gets confused, as it probably should. :-) Cheers, Bert Be

Re: [R] Automating regression

2014-12-22 Thread Bert Gunter
lse will groc what you mean -- or you could define yourself more clearly. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Dec 22, 2014 a

Re: [R] Automating regression

2014-12-22 Thread Bert Gunter
.. that should have been either myProc <- function(FUN, ...) do.call(FUN,list(...)) or myProc <- function(FUN, ...) FUN(...) My other comments still apply. Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. An

Re: [R] Help re. Fitting non-linear dynamic model to time series data

2014-12-22 Thread Bert Gunter
with statistics than R coding. Cheers, Bert Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Dec 22, 2014 at 12:14 PM, Ciaran Roberts w

Re: [R] ave(x, y, FUN=length) produces character output when x is character

2014-12-24 Thread Bert Gunter
actly what you claim is wanted without trying to invent square wheels. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Wed, Dec 24, 2014 at

Re: [R] ave(x, y, FUN=length) produces character output when x is character

2014-12-24 Thread Bert Gunter
y due to historical quirks -- evolution is not necessarily orderly. But that's the way it is; that's the way it's documented; and tutorials will tell you about this (that's how I learned). So please stop guessing and intuiting and read the docs to understand how things w

Re: [R] ave(x, y, FUN=length) produces character output when x is character

2014-12-25 Thread Bert Gunter
On Thu, Dec 25, 2014 at 12:41 PM, Duncan Murdoch wrote: > Would you say a cube contains a polygon, or a cube is a polygon? Neither, actually. I'd say a cube is a polyhedron or a square is a polygon. :-) But point taken, of course. Cheers, Bert Bert Gunter Genentech Non

Re: [R] ave(x, y, FUN=length) produces character output when x is character

2014-12-25 Thread Bert Gunter
[1] "a" "b" "c" $class [1] "factor" > is.vector(f) [1] FALSE > attributes(f) <- NULL > f [1] 1 2 3 > is.vector(f) [1] TRUE Don't you think it's time to call a halt to this? Cheers, Bert Bert Gunter Genentech Nonclinical Biostat

Re: [R] ANOVA test

2014-12-26 Thread Bert Gunter
This is a statistical question primarily and, as such, is off topic here. Either consult a local statistical expert or post to a statistical site like stats.stackexchange.com . Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Informati

Re: [R] CRF Usage example in R

2015-01-09 Thread Bert Gunter
ternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and > more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contai

Re: [R] CRF Usage example in R

2015-01-09 Thread Bert Gunter
LEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowle

Re: [R] R vs. RStudio?

2015-01-10 Thread Bert Gunter
SE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge

Re: [R] Reading large data matrix in R

2015-01-12 Thread Bert Gunter
software. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Jan 12, 2015 at 3:17 AM, Madhavan BL wrote: > Hello, > > While I

Re: [R] overlapping coefficient bidimensional distribution

2015-01-12 Thread Bert Gunter
?intersect or, more generally, ?match Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Jan 12, 2015 at 5:09 AM, Meli Massimili

Re: [R] nonmonotonic glm?

2015-01-12 Thread Bert Gunter
... but do realize that after you have looked at the data to determine the appropriate modeling approach, no statistical inference (significance tests, confidence intervals, etc.) should be done on the model used. Or more precisely, any that is done is wrong. Cheers, Bert Bert Gunter Genentech

Re: [R] any r package can handle factor levels not in the test set

2015-01-13 Thread Bert Gunter
ing a search engine) to see whether/how such methodology is implemented in R. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Jan 13,

Re: [R] Mac Yosemite R.app console typing / autorepeat performance

2015-01-13 Thread Bert Gunter
This clearly should go to the r-sig-mac list, not r-help. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Jan 13, 2015 at 10:44 A

Re: [R] Probably dumb question about lists

2015-01-14 Thread Bert Gunter
erving the class attribute of tl. Ergo you are getting the basic plot method and not the plot.tile.list method. HTH. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom.&

Re: [R] Extracting elements out of list in list in list

2015-01-16 Thread Bert Gunter
) to remove them > unlist(extr(x,"A")) f1.x1.A f1.x2.A f2.x3.A f2.x4.A 11 12 13 14 I would welcome any possibly "slicker" versions of the above. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. In

Re: [R] Proportion of equal entries in dist()?

2015-01-20 Thread Bert Gunter
large" matrices. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Jan 20, 2015 at 10:10 AM, Adams, Jean wrote: > Jorge, > > I

Re: [R] Re-order levels of a categorical (factor) variable

2015-01-21 Thread Bert Gunter
ep(c("pre", "day2","day10"))) > xyplot(y~x,data=z) ## x axis order is day 10, day2, pre > levels(z$x) [1] "day10" "day2" "pre" > z$x <- factor(as.character(z$x),levels=c(levels(z$x)[3:1])) ## explicitly > defines lev

Re: [R] How to split a data table into multiple tables based on column (unique value) and save them separately

2015-01-22 Thread Bert Gunter
nt, so you might start there. Hadley Wickham's widely used plyr package provides what many consider to be much more usable interfaces and flexibility for this sort of thing, so you may wish to learn and use this instead. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-737

Re: [R] ggplot courses?

2015-01-22 Thread Bert Gunter
Check: http://www.rstudio.com/ (where Hadley works) -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Thu, Jan 22, 2015 at 10:14 AM, Erin Hodg

Re: [R] sum of grouped elements of vector

2015-01-23 Thread Bert Gunter
. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Fri, Jan 23, 2015 at 12:45 AM, Kathryn Lord wrote: > Dear R users, > > I have a quick q

Re: [R] Testing for significant differences between groups in multiple linear regression

2015-01-23 Thread Bert Gunter
. ... oh, and while you're at it, please read the posting guide for this list (see link below) and, in particular, DO NOT POST IN HTML, which, as you can see here, often becomes a mess on this **plain text** mailing list. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467

Re: [R] GLM for proportional data

2015-01-24 Thread Bert Gunter
wrong list. Post on stats.stackexchange.com -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sat, Jan 24, 2015 at 1:53 AM, Mauricio Go

Re: [R] How to add a line of "mean" to the bwplot (lattice package)

2015-01-25 Thread Bert Gunter
ers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Jan 25, 2015 at 7:06 AM, Jun Shen wrote: > Dear list, > > The bwplot generates box plo

Re: [R] Sum function and missing values --- need to mimic SAS sum function

2015-01-27 Thread Bert Gunter
Huh?? > ifelse(TRUE, a <- 2L, a <- 3L) [1] 2 > a [1] 2 Please clarify. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, J

Re: [R] panel.xyplot and panel.loess using two different groupings

2015-01-27 Thread Bert Gunter
you might want to add checks in the panel.groups function for this. HTH. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Jan 27, 2015 at

Re: [R] panel.xyplot and panel.loess using two different groupings

2015-01-27 Thread Bert Gunter
ra package, as it provides a ggplot-like syntax ("+") and facilities to handle duplicated parameters in situations like this. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly no

Re: [R] How to add error bars to a line xyplot (lattice package)

2015-01-27 Thread Bert Gunter
hey should be model based, obtained,e.g. by the relevant predict() method. But this is OT for here, so I won't rant further. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisd

Re: [R] Paste every two columns together

2015-01-28 Thread Bert Gunter
ith a " " (space): paste0(ID1[i],ID1[i+1],collapse= " ") Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Wed, Jan 28, 201

Re: [R] repeated measures ANOVA

2015-01-30 Thread Bert Gunter
Mixed effects models. Packages nlme and lme4 among others. Ask further on r-sig-mixed-models list. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford St

Re: [R] save program results and graphs to one file

2015-02-01 Thread Bert Gunter
But in addition to what Jeff noted, see ?save and ?save.image (noting that that the resulting .Rdata file can only be read by R). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainl

Re: [R] (Cross-)validating GLM including factors

2015-02-02 Thread Bert Gunter
This is not an r-help issue. Post on a statistics list like stats.stackexchange.com instead for a variety of opinions.* Cheers, Bert *Which will probably be useless given your paucity of data. But that's just *my* useless opinion. Bert Gunter Genentech Nonclinical Biostatistics (650) 467

Re: [R] Still trying to avoid loops

2015-02-04 Thread Bert Gunter
tapply() (of which by() is essentially a wrapper) **is** a (disguised) loop (at the R level, of course). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clif

Re: [R] Computing RLE and NUSE Scores

2015-02-06 Thread Bert Gunter
Ask on the BioConductor list, not here. Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Fri, Feb 6, 2015 at 9:27 PM, Beverly Nguyen wrote: > I am

Re: [R] prediction intervals for robust regression

2015-02-11 Thread Bert Gunter
;re looking for (and it's not already available). Further questions should go to a statistics site like stats.stackexchange.com, as statistical questions are off topic here. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Informat

Re: [R] How to get the error and error variance after HB using bayesm

2015-02-13 Thread Bert Gunter
approaches. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Fri, Feb 13, 2015 at 11:46 AM, ARNAB KR MAITY wrote: > Hello Michael, > > I h

Re: [R] Comparing gam models fitted using ti()

2015-02-16 Thread Bert Gunter
I wonder if this might be better asked on a statistical list -- e.g. stats.stackexchange.com -- as this seems to involve complex statistical model comparison issues, which are normally OT here. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not inform

Re: [R] split dataframe to several dataframes in R

2015-02-17 Thread Bert Gunter
Inline. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Feb 17, 2015 at 1:33 AM, Chris Campbell wrote: > # Assuming you want to cre

[R] Fwd: How to analyse nonlinear response to categorical and quantitative explanatory variables?

2015-02-19 Thread Bert Gunter
Statistical methodology questions are generally off topic here. This list is about R programming. I suggest you post to a statistics site like stats.stackexchange.com instead. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Informati

Re: [R] Density plots for Chi-square and F-distribution on my data

2015-02-19 Thread Bert Gunter
I think you might do well to consult a local statistical expert, as you appear to be out of your statistical depth here. Also, is this homework? If so, this is not a homework site (although posters sometimes get help anyway). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650

Re: [R] multiple parameter optimization with optim()

2015-02-20 Thread Bert Gunter
This is not the proper venue for a discussion of the mathematics of optimization, no matter that it is interesting. Please take it off list. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowled

Re: [R] Subsetting a list of lists using lapply

2015-02-20 Thread Bert Gunter
How can you expect a solution if you cannot specify the problem? -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Fri, Feb 20, 2015 at 6:1

Re: [R] Split a dataframe by rownames and/or colnames

2015-02-20 Thread Bert Gunter
I think ?tapply and friends: ?by ?aggregate ?ave is what you want. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Fri, Feb 20, 2015 at 9:

Re: [R] glm.fit when family argument is not a "family" object

2015-02-22 Thread Bert Gunter
... But glm(1~1, family = gaussian) ## works. As does family = "gaussian" (R 3.1.2 ) Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stol

Re: [R] glm.fit when family argument is not a "family" object

2015-02-22 Thread Bert Gunter
Exactly. Just a matter of fixing the glm.fit Help page. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Feb 22, 2015 at 3:11 PM, David

Re: [R] how to standardize raster image (GeoTiff) in R?

2015-02-23 Thread Bert Gunter
much much faster. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Feb 23, 2015 at 3:18 PM, Angela Smith wrote: > Hi R user, > Would

Re: [R] Best Mac for R

2015-02-25 Thread Bert Gunter
What does this have to do with R? Does the answer not depend on what you intend to do with your laptop, e.g the sorts of data you deal with, of which we have no idea? Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information i

Re: [R] Dummy variable in ARIMA

2015-02-26 Thread Bert Gunter
Inline. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Thu, Feb 26, 2015 at 8:02 AM, Mikael Olai Milhøj wrote: > Hi all > > I

Re: [R] Dummy variable in ARIMA

2015-02-26 Thread Bert Gunter
Dig deeper. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Thu, Feb 26, 2015 at 8:29 AM, Mikael Olai Milhøj wrote: > Hi. > > First of

Re: [R] (no subject)

2015-03-01 Thread Bert Gunter
1. Do not multiple post -- it's rude. 2. Use plain text emails, not HTML (although it makes no difference here). 3. Homework? -- This list has a no homework policy. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Informat

Re: [R] R CMD wants (unnecessary) packages ?

2015-03-01 Thread Bert Gunter
stall these core packages, as they may be necessary for key functionality that you use, directly or indirectly. And apologies if I have misinterpreted... Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And kn

Re: [R] How to decide " weight" in WLS model in R ?

2015-03-02 Thread Bert Gunter
Angela: These are statistical, not R, issues I believe, and you appear to be out of your depth statistically here. I suggest you talk to a local statistical resource or, if you can't find such help, post on a statistical site like stats.stackexchange.com. Cheers, Bert Bert Gunter Gene

Re: [R] R 3.1.2 using a custom function in aggregate() function on Windows 7 OS 64bit

2015-03-04 Thread Bert Gunter
What do you think dat$a is? I recommend that you spend some time with an R tutorial if you plan to use R. Your code is pretty bad. Examples: use of the ifelse construction instead of if ... else; return() Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data i

Re: [R] R 3.1.2 using a custom function in aggregate() function on Windows 7 OS 64bit

2015-03-05 Thread Bert Gunter
Sorry, Jeff. aggregate() is generic. >From ?aggregate: "## S3 method for class 'data.frame' aggregate(x, by, FUN, ..., simplify = TRUE)" Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not

Re: [R] R 3.1.2 using a custom function in aggregate() function on Windows 7 OS 64bit

2015-03-05 Thread Bert Gunter
components of by, and FUN is applied to each such subset with further arguments in ... passed to it." As I read this, the argument of FUN is a data frame that is a subset of the original frame, defined by the by variable values. No? -- Bert Bert Gunter Genentech Nonclinical Biostatistics

Re: [R] R 3.1.2 using a custom function in aggregate() function on Windows 7 OS 64bit

2015-03-05 Thread Bert Gunter
Well, I obviously don't use it either, as I'm just quoting the docs. I either use by(), or tapply(). -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Cliff

Re: [R] Special case of data classification

2015-03-07 Thread Bert Gunter
Have you done any homework for yourself? CRAN task views would be an obvious first place to look: http://cran.r-project.org/web/views/ Both the environmental and spatio-temporal views seem like possibilities. If that doesn't suit, try searching on suitable keywords Cheers, Bert Bert G

Re: [R] Date extract Year

2015-03-08 Thread Bert Gunter
Please read ?as.Date carefully. Your argument appears to be numeric (??) and you do not seem to have specified an origin: "as.Date will accept numeric data (the number of days since an epoch), but only if origin is supplied." Cheers, Bert Bert Gunter Genentech Nonclinical Biostati

Re: [R] Cluster mapping data

2015-03-08 Thread Bert Gunter
Have you looked at the "Cluster" task view on CRAN? http://cran.r-project.org/web/views/ -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll

Re: [R] Date extract Year

2015-03-08 Thread Bert Gunter
... and that is the source of your difficulties, as I and others have told you. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Mar 8, 201

Re: [R] Returning to parent function

2015-03-17 Thread Bert Gunter
gt; >> You need to explain what you want to accomplish rather than posting >> failed code. >> >> > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > P

Re: [R] Joining two datasets - recursive procedure?

2015-03-21 Thread Bert Gunter
1. Still not sure what you mean, but maybe look at ?ave and ?tapply, for which ave() is a wrapper. 2. You still need to heed the rest of Jeff's advice. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not know

Re: [R] Joining two datasets - recursive procedure?

2015-03-21 Thread Bert Gunter
z <- rnorm(nrow(f1)) ## or anything you want z1 <- f1$v4 + z - with(f1,ave(z,v1,v2,FUN=mean)) aggregate(v4~v1,f1,sum) aggregate(z1~v1,f1,sum) aggregate(v4~v2,f1,sum) aggregate(z1~v2,f1,sum) aggregate(v4~v3,f1,sum) aggregate(z1~v3,f1,sum) Cheers, Bert Bert Gunter Genentech Noncl

Re: [R] Joining two datasets - recursive procedure?

2015-03-21 Thread Bert Gunter
... or cleaner: z1 <- with(f1,v4 + z -ave(z,v1,v2,FUN=mean)) Just for curiosity, was this homework? (in which case I should probably have not provided you an answer -- that is, assuming that I HAVE provided an answer). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650)

Re: [R] Joining two datasets - recursive procedure?

2015-03-22 Thread Bert Gunter
se that it relies on the implementation of factor() rather than its API -- is: y <- f1$v3 ## to simplify the notation; could be done using with() z <- (c(29,2.567)/table(y))[c(y)] Then proceed to z1 as I previously described -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650)

Re: [R] Joining two datasets - recursive procedure?

2015-03-22 Thread Bert Gunter
what you are trying to do. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Mar 22, 2015 at 7:55 AM, Bert Gunter wrote: > I would have thou

Re: [R] Why can't I access this type?

2015-03-22 Thread Bert Gunter
ames and the links therein for more info. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Mar 22, 2015 at 7:39 AM, Yves S. Garret wrote: &g

Re: [R] My R stopped working

2015-03-22 Thread Bert Gunter
restart R. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Sun, Mar 22, 2015 at 2:11 AM, adeela uaf wrote: > Hi, > I was downloading som

Re: [R] Fwd: Joining two datasets - recursive procedure?

2015-03-22 Thread Bert Gunter
for v1, which is not supposed to change according to my understanding of your specification. So I'm done. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford

Re: [R] Conversion of Matlab code to an R code

2015-03-23 Thread Bert Gunter
asks. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Mar 23, 2015 at 8:10 AM, Abhinaba Roy wrote: > Hi, > > Can a Matlab code be co

Re: [R] multiple break in univariate series

2015-03-28 Thread Bert Gunter
Once you have looked at the data and chosen change points to test based on the data, the tests for change points are invalid (unless you make appropriate adjustments for post hoc tests). And no, I am not making this up. Consult any competent statistician. Cheers, Bert Bert Gunter Genentech

Re: [R] changing column labels for data frames inside a list

2015-03-30 Thread Bert Gunter
Sarah's statement is correct. So is yours. They are not contradictory, and I believe Sarah's point was that the OP needed to learn the appropriate syntax. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not

[R] Fwd: changing column labels for data frames inside a list

2015-03-30 Thread Bert Gunter
Sorry, I failed to cc the list. -- Bert -- Forwarded message -- From: Bert Gunter Date: Mon, Mar 30, 2015 at 8:36 AM Subject: Re: [R] changing column labels for data frames inside a list To: Vikram Chhatre You really really need to spend (more?) time with a good R

Re: [R] Calculating Kendall's tau

2015-03-31 Thread Bert Gunter
This sounds like homework. Homework is discouraged on this list (but you might get lucky). Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll O

Re: [R] Calculating Kendall's tau

2015-03-31 Thread Bert Gunter
chance of getting useful help. -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Mar 31, 2015 at 9:24 AM, Desta Yoseph wrote: > Dear Bert

Re: [R] about model.matrix

2015-04-02 Thread Bert Gunter
27987551&sr=1-5&keywords=linear+models+in+statistics ... but there are dozens. Better yet, consult a local statistical expert for help. Trying to educate yourself is laudable, but may be unrealistic. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is

Re: [R] Aggregating daily rainfall raster data to bimontly data

2015-04-05 Thread Bert Gunter
-monthly, depending on your meaning of "bi") without knowing the months. Aggregating by every n days would be easy, but that's probably not what you want. Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not know

Re: [R] Verify that a grid is uniform

2015-04-06 Thread Bert Gunter
of ==, e.g. all(diff(z) - z[2] + z[1] < tol) Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Apr 6, 2015 at 10:11 AM, Marc La

Re: [R] Verify that a grid is uniform

2015-04-06 Thread Bert Gunter
... correction: you need to use absolute value for the comparison, of course. all(abs(diff(z) - z[2] + z[1]) < tol) -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom.&q

Re: [R] sort adjacency matrix

2015-04-06 Thread Bert Gunter
Not quite, David. If I understand the OP's query, he wants the ties to be broken by the "lexicographic" order (with apologies if I have misused this term) of the 1's within the rows. Makes things a bit more interesting. Have at it! Cheers, Bert Bert Gunter Genentech Noncl

Re: [R] Verify that a grid is uniform

2015-04-06 Thread Bert Gunter
Does not min(abs(diff(z))) give you the scaling you need to set a tolerance? -- Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Mon, Apr 6, 2015

Re: [R] nested for loops too slow

2015-04-12 Thread Bert Gunter
ecessarily denies the wisdom of the approach you've suggested, however. It may indeed be considerably faster, but timing will have to tell. I am just trying to correct ​(again) ​ the ​ ​ widely held misperception ​that yo​ u ​ seem to​ express ​.​ Cheers, Bert Bert Gunter Genentech Nonclinica

Re: [R] : automated levene test and other tests for variable datasets

2015-04-14 Thread Bert Gunter
searching. Choose one that suits your tastes/learning style and have at it! Cheers, Bert Bert Gunter Genentech Nonclinical Biostatistics (650) 467-7374 "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." Clifford Stoll On Tue, Apr 14, 2015

Re: [R] Formatting summary() output to a file

2018-07-27 Thread Bert Gunter
Not quite sure what you mean here. R is open source, so > print.summaryDefault ## at the prompt. It's in base R, so no package:: prefix needed will give you the code used for formatting. You can then do the same. Cheers, Bert Bert Gunter "The trouble with having an open

Re: [R] (no subject)

2018-07-30 Thread Bert Gunter
How can one possibly answer this without knowing the structure of your dataset? -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On

Re: [R] New post for Rhelp

2018-08-01 Thread Bert Gunter
Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Aug 1, 2018 at 8:40 AM, Edoardo Silvestri < silvestri.cas...@gmail.com> wrote:

Re: [R] Philip Morris International - Windows10 migration assessment

2018-08-02 Thread Bert Gunter
R is free and open source. Your queries are inappropriate for this list, which is about help for programming in R. Please go here and follow the relevant links to answer your questions: https://www.r-project.org/ Cheers, Bert Bert Gunter "The trouble with having an open mind is that p

Re: [R] Combinations of true/false values where one pair is mutually exclusive

2018-08-02 Thread Bert Gunter
Logic: !(E == "fail" & F == "fail) <==> (E == "pass" | F == "pass") -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "

Re: [R] kSamples ad.test question

2018-08-02 Thread Bert Gunter
, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Aug 2, 2018 at 12:00 PM, Andras Farkas via R-help < r-help@r-project.org> wr

Re: [R] Lag function

2018-08-07 Thread Bert Gunter
(= hours) and -24 moves it forward 24 hours. This means that the day 2 observations are those from day 1, etc., which is usually what is wanted for lag. But you decide. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things

Re: [R] Submit your own R package - @examples

2018-08-08 Thread Bert Gunter
This should be posted on the r-package-devel list rather than here: https://stat.ethz.ch/mailman/listinfo/r-package-devel Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in

[R] Fwd: Help with finalfit and knitr

2018-08-08 Thread Bert Gunter
; Fortune nomination! Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Aug 8, 2018 at 8:06 AM, Jeff Newmiller wrote: > R syn

Re: [R] Warnings using SuperLearner?

2018-08-09 Thread Bert Gunter
, though there are certainy wiser folks than I out there. Also, this is a plain text maiing list. Please post in plain text, not HTML. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breat

Re: [R] Trying to Generalize a Function in R

2018-08-09 Thread Bert Gunter
" I am thinking that there should be an R command to tell me about the structure of series" ?str ## perhaps also/instead ?summary Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka

<    4   5   6   7   8   9   10   11   12   13   >