Re: [R] lmer coefficient distributions and p values

2007-08-15 Thread Steven McKinney
ethz.ch/mailman/listinfo/r-sig-mixed-models Also check the latest documentation for lme4 and the lmer() and lmer2() functions at http://cran.r-project.org/ in the Packages ... lme4 pages. Hope this helps Steven McKinney Statistician Molecular Oncology and Breast Cancer Program British Columbia

Re: [R] Systematically biased count data regression model

2007-08-09 Thread Steven McKinney
omewhat related to your situation, and it has plenty of good references Application of Shrinkage Techniques in Logistic Regression Analysis: A Case Study E. W. Steyerberg Statistica Neerlandica, 2001, vol. 55, issue 1, pages 76-88 Steven McKinney Statistician Molecular Oncology and Breast

Re: [R] FW: Selecting undefined column of a data frame (was[BioC]read.phenoData vs read.AnnotatedDataFrame)

2007-08-04 Thread Steven McKinney
") : undefined columns selected > sessionInfo() R version 2.5.1 (2007-06-27) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base

Re: [R] FW: Selecting undefined column of a data frame (was [BioC]read.phenoData vs read.AnnotatedDataFrame)

2007-08-03 Thread Steven McKinney
Hi Bert, > -Original Message- > From: Bert Gunter [mailto:[EMAIL PROTECTED] > Sent: Fri 8/3/2007 3:19 PM > To: Steven McKinney; r-help@stat.math.ethz.ch > Subject: RE: [R] FW: Selecting undefined column of a data frame (was > [BioC]read.phenoData vs read.Annotated

Re: [R] FW: Selecting undefined column of a data frame (was [BioC] read.phenoData vs read.AnnotatedDataFrame)

2007-08-03 Thread Steven McKinney
oding/debugging implications, and as you point out, "[.data.frame is one of the most complex functions in R" so please bear with me. This change in behaviour has taken away a side-effect debugging tool, discussed below. > > > On Fri, 3 Aug 2007, Steven McKinney wrote:

Re: [R] FW: Selecting undefined column of a data frame (was [BioC] read.phenoData vs read.AnnotatedDataFrame)

2007-08-03 Thread Steven McKinney
> -Original Message- > From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] > Sent: Fri 8/3/2007 1:05 PM > To: Steven McKinney > Cc: r-help@stat.math.ethz.ch > Subject: Re: [R] FW: Selecting undefined column of a data frame (was [BioC] > read.phenoData vs read.Annotated

Re: [R] FW: Selecting undefined column of a data frame (was [BioC] read.phenoData vs read.AnnotatedDataFrame)

2007-08-03 Thread Steven McKinney
er than after? -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Fri 8/3/2007 12:25 PM To: Steven McKinney Cc: r-help@stat.math.ethz.ch Subject: Re: [R] FW: Selecting undefined column of a data frame (was [BioC] read.phenoData vs read.AnnotatedDataFrame) You a

Re: [R] FW: Selecting undefined column of a data frame (was [BioC]read.phenoData vs read.AnnotatedDataFrame)

2007-08-03 Thread Steven McKinney
erpc-apple-darwin8.9.1 locale: en_CA.UTF-8/en_CA.UTF-8/en_CA.UTF-8/C/en_CA.UTF-8/en_CA.UTF-8 attached base packages: [1] "stats" "graphics" "grDevices" "utils" "datasets" "methods" "base" other attached packages:

[R] FW: Selecting undefined column of a data frame (was [BioC] read.phenoData vs read.AnnotatedDataFrame)

2007-08-03 Thread Steven McKinney
UTF-8/en_CA.UTF-8/en_CA.UTF-8/C/en_CA.UTF-8/en_CA.UTF-8 attached base packages: [1] "stats" "graphics" "grDevices" "utils" "datasets" "methods" "base" other attached packages: plotrix lme4 Matr

Re: [R] Extracting a website text content using R

2007-08-01 Thread Steven McKinney
>-Original Message- >From: [EMAIL PROTECTED] on behalf of Am Stat >Sent: Wed 8/1/2007 2:19 PM >To: r-help@stat.math.ethz.ch >Subject: [R] Extracting a website text content using R >Dear useR, >Just wandering whether it is possible that there is any function in R could >let me get the t

[R] Design matrix question

2007-05-16 Thread Steven McKinney
del > model.matrix(~ A * B - B, df) (Intercept) Aa Ax:Bb Aa:Bb 1 1 1 0 1 2 1 1 0 0 3 1 0 0 0 4 1 0 0 0 attr(,"assign") [1] 0 1 2 2 attr(,"contrasts") attr(,"contrasts")$A [1]

Re: [R] lmer function

2007-05-14 Thread Steven McKinney
powerpc-apple-darwin8.9.1 locale: C attached base packages: [1] "stats" "graphics" "grDevices" "utils" "datasets" "methods" [7] "base" other attached packages: lme4 Matrix lattice "0.997

Re: [R] reodering factor

2007-05-03 Thread Steven McKinney
B 20 10.6319798 BB > lm(y ~ grp2, data = mydata) Call: lm(formula = y ~ grp2, data = mydata) Coefficients: (Intercept)grp2A 10.66 -10.17 > Steven McKinney Statistician Molecular Oncology and Breast Cancer Program British Columbia Cancer Research Centre

Re: [R] Computing an ordering on subsets of a data frame

2007-04-19 Thread Steven McKinney
} > Y A B C 3 1 4 3 7 3 2 2 9 3 3 3 6 2 1 2 5 2 1 1 1 1 2 1 2 1 3 2 8 3 1 1 4 2 3 3 > oY <- order(Y$A) > Y[oY,] A B C 3 1 4 3 1 1 2 1 2 1 3 2 6 2 1 2 5 2 1 1 4 2 3 3 7 3 2 2 9 3 3 3 8 3 1 1 > HTH Steven McKinney Statistician Molecular Oncology and Breast Cancer Program Briti

Re: [R] extracting intercept from ppr fit

2007-04-18 Thread Steven McKinney
rms =1) > ## how to extract a = 10 from data.ppr? > data.ppr$yb [1] 9.973964 > a <- 210 > data$z <- evalq(a + atan (x + y) + rnorm (n), data) > data.ppr <- ppr(z ~ x + y, data=data, nterms =1) > data.ppr$yb [1] 209.9773 > HTH Steven McKinney Statistician Molecular

Re: [R] Extracting approximate Wald test (Chisq) fromcoxph(..frailty)

2007-04-17 Thread Steven McKinney
score cat("\nIterations:", object$iter[1], "outer,", object$iter[2], "Newton-Raphson\n") if (length(print2)) { for (i in 1:length(print2)) cat("", print2[i], "\n") } if (is.null(object$df)) df <- s

Re: [R] problem with RCurl install on Unix

2007-03-21 Thread Steven McKinney
quot;grDevices" "utils" "datasets" "methods" [7] "base" other attached packages: DBI "0.1-12" > Steven McKinney Statistician Molecular Oncology and Breast Cancer Program British Columbia Cancer Research Centre email: [EMAIL

Re: [R] replacing all NA's in a dataframe with zeros...

2007-03-14 Thread Steven McKinney
E 5 TRUE FALSE TRUE TRUE FALSE 6 FALSE TRUE TRUE FALSE FALSE > mydata.df[is.na(mydata.df)] <- 0 > mydata.df V1 V2 V3 V4 V5 1 1 0 1 1 1 2 1 0 0 0 1 3 0 0 1 0 0 4 0 0 0 0 1 5 0 1 0 0 1 6 1 0 0 1 1 > Steven McKinney Statistician Molecular Oncology and Br

Re: [R] dendrogram - got it , just need to label :)

2007-03-09 Thread Steven McKinney
e able to identify nodes deeper in the tree. This is aided by having access to internal node labels/names and being able to extract internal nodes by those labels/names. Best Steven McKinney Statistician Molecular Oncology and Breast Cancer Program British Columbia Cancer Research Centr

Re: [R] Double-banger function names: preferences and suggestions

2007-02-26 Thread Steven McKinney
of the graphical codes, left-arrow mutated to the underscore of ASCII-1967. It may have had earlier, or other, meanings, but for some early programming languages it was "assignment", eg. c ? b + a "C is assigned the sum of B and A". Steven McKinney Statistician Mo

Re: [R] How to plot two graphs on one single plot?

2007-02-23 Thread Steven McKinney
ean=0, sd=1) > > x2=rnorm(25, mean=0, sd=1) > y2=dnorm(x2, mean=0, sd=1) > plot(x1, y1, type='p', xlim=range(x1,x2), ylim=range(y1, y2), xlab='x', > ylab='y') > points(x2, y2, type='p', col="red", xlab='x', ylab=&#

Re: [R] Google Custom Search Engine for R

2007-02-23 Thread Steven McKinney
Whereas "R" is very generic, "CRAN" is much less so. I've had very good luck adding CRAN to my search terms, e.g. try to Google cran 3d scatterplot This produces all R-related hits on the first Google page. Hope this helps Steven McKinney Statistician Molecular O

[R] Bug/problem reporting: Possible to modify posting guide FAQ?

2006-08-28 Thread Steven McKinney
If users post a bug or problem issue to an R-based news group (R-devel, R-help, BioC - though BioC is far more forgiving) they get yelled at for not reading the posting guide and FAQ. "Please *_do_* read the FAQ, the posting guide, ..." the yellers do say. So I read the BioC FAQ and it says...

[R] How do I modify an exported function in a locked environment?

2006-07-20 Thread Steven McKinney
iting, I see the original copy. But if I reinvoke the editor via fixInNamespace(), I do see my modification. Where is my copy residing? How do I push it out to replace the exported copy? Is this the proper way to modify a package function? Are there other ways? I've searched webpages, R news

[R] scripts to process array CGH data files from NimbleGen

2006-06-12 Thread Steven McKinney
and use to create an aCGH object. Any info appreciated. Best Steve McKinney Steven McKinney Statistician Molecular Oncology and Breast Cancer Program British Columbia Cancer Research Centre [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list