Re: [R] Weighted Mann-Whitney-Wilcoxon-Test

nder complex sampling. Biometrika, 100, 831-842. -thomas -- Thomas Lumley Professor of Biostatistics University of Auckland __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] problem when I Call C subfunction in void function

gt; $Px > [1] 4 10 10 8 > > $tailleP > [1] 4 > > $res > [1] 4 0 0 0 > > I haven't problem in "essai" function but I can't correctly return "Px" > vector. > I d'ont understand why I get only the first number (number 4 in my exemp

Re: [R] Can you get the DEFT from svyratio?

No, you can't (at the moment), though it shouldn't be too hard to extend. I can't run your example, though. I get: Error in eval(expr, envir, enclos) : object 'M' not found -thomas Thomas Lumley Professor of Biostatistics University of Auckland ___

Re: [R] Warning when using survey:::svyglm

cative of a problem? Thanks, Ben __ 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 and provide commented, minimal, self-contained, re

Re: [R] Attempting to make a custom color spectrum to use in heatmap.2

rifying before you can give me what you feel is an appropriate response. Thanks again for your time. Regards, Daniel __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-p

Re: [R] produce variable on the fly

s, and then give them a value, for example: vect = c(10:20) for (i in 1:10) { cat(paste("VAR",i,sep="")) = vect[i] } This is almost a candidate for a FAQ: see It *is* FAQ 7.21. -thomas Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECT

Re: [R] ignoring zeros or converting to NA

t.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 and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED]

Re: [R] ignoring zeros or converting to NA

-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __ R-help@r-project.org mailing l

Re: [R] Quantile regression with complex survey data

nimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do r

Re: [R] Survey Design / Rake questions

od reference for such things as "PSUs", "cluster sampling", and so on. -- Stas Kolenikov, also found at http://stas.kolenikov.name Small print: I use this email account for mailing lists only. __ R-help@r-project.org mailing list

Re: [R] Help Regarding 'integrate'

ailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle

Re: [R] Large data sets with R (binding to hadoop available?)

sure of the current status of the HDF5 support in Bioconductor, though. -thomas Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __ R-help@r-project.org mailing list

Re: [R] filtering out data

he posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle _

Re: [R] Survey Design / Rake questions

Woodman str(EBSurvey$NumStn) Factor w/ 12 levels "1","2","3","4",..: 10 12 4 12 8 1 8 8 12 4 ... EBSurvey$NumStn[1:5] [1] 10 12 4 12 8 Levels: 1 2 3 4 5 6 7 8 9 10 11 12 str(ByEBNum$StnTraveld) Factor w/ 12 levels "1","10","11",..: 1 5 6 7 8 9 10 11 12 2 ... ByEBNum$StnTraveld[1:5] [1] 1 2 3 4 5 Levels: 1 10 11 12 2 3 4 5 6 7 8 9 ***

Re: [R] SQL Primer for R

der results of a query as you read them: SELECT foo, bar FROM table_name ORDER BY baz -thomas Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __ R-help@r-project.

Re: [R] How to do a meta-analysis plot

On Mon, 25 Aug 2008, Jorge Ivan Velez wrote: My problem is that I don't have the raw data as rmeta _requires_ and, even when I have my data set in the _same_ (?) format that summary(a), when I tried plot(mydata) it doesn't work. No, it doesn't _require_ that. If you just want a forest plot the

Re: [R] SQL Primer for R

. dbClearResult -thomas Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] svymeans question

r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __

Re: [R] Survey Design / Rake questions

;,"2","3","4",..: 10 12 4 12 8 1 8 8 12 4 ... EBSurvey$NumStn[1:5] [1] 10 12 4 12 8 Levels: 1 2 3 4 5 6 7 8 9 10 11 12 str(ByEBNum$StnTraveld) int [1:12] 1 2 3 4 5 6 7 8 9 10 ... ByEBNum$StnTraveld[1:5] [1] 1 2 3 4 5 *

Re: [R] Survey Design / Rake questions

quot;Balboa","De Soto",..: 11 2 5 8 6 1 12 7 10 13 ... all(levels(EBSurvey$lineon)==StnName) [1] TRUE # str(EBDesign$NumStn) NULL str(EBSurvey$NumStn) Factor w/ 12 levels "1","2","3","4",..: 10 12 4 12 8 1 8 8 12 4 ... str(ByEBNum$StnTraveld) Factor w/ 12 levels "1","2","3","4",..: 1 2 3 4 5 6 7 8

Re: [R] How is the binding for a super assignment made visible?

p level in your package. The superassignments would then modify that variable. -thomas Thomas Lumley Assoc. Professor, Biostatistics [EMAIL PROTECTED] University of Washington, Seattle __ R-help@r-project.org mailing

Re: [R] Several simple but hard tasks to do with R

ction and sends it to the screen. The R sink("analysis.log",split=TRUE) is almost identical to the Stata log using analysis.log with sink() to close the log and flush to disk at the end of the session. -thomas Thomas Lumley Assoc. Professor, Biostatis

Re: [R] R Memory Usage Concerns

l use more than the object size for the data set, because it makes temporary copies of things. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-pr

Re: [R] Factor Analysis function source code required

ost texts and some other programs would lead one to believe." If SPSS and R are fitting the same factor analysis model you can compare the fitted likelihood to see which package has found a better solution. -thomas Thomas Lumley Assoc. Professor, Biostatistics t

Re: [R] identical(length(x), 1) returns FALSE, but print(length(x)) is 1, length(x)==1 is TRUE, and is.integer(lenght(x)) is TRUE????

al, just equal, so identical(length(x), 1) is FALSE and length(x)==1 is TRUE. identical(length(x), 1L) should be TRUE, since 1L is the way to specify an integer value of 1. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniver

Re: [R] wilcox.test p-value = 0

do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle

Re: [R] Teasing out logrank differences *between* groups using survdiff or something else?

ead the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle _

Re: [R] What is the time complexity of 'match()'?

etimes used for quite large m and n. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

Re: [R] Problems with the commands FUNCTION and DERIV to build a polynomial

ou really don't need deriv() to differentiate a polynomial, and so you can use ordinary lexical scope rather than substitution, for a much tidier answer derivative <- function(coef){ function(x) coef[2]+x*(2*coef[3]+x*3*coef[4]) } -thomas Thomas Lumley A

Re: [R] Fastest Way to Divide Elements of Row With Its RowSum

ame and 1 to the number of entries if dat is a matrix, not from 1 to the number of rows. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-proj

Re: [R] WMF conversion...

://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-proj

Re: [R] cluster/distance large matrix (fwd)

On Thu, 11 Feb 2010, Christian Hennig wrote: >It is well know that hierarchical methods are problematic with too large >dissimilarity matrices; even if you resolve the memory problem, the number of >operations required is enormous. There is at least one exception to this. Single-linkage hierar

Re: [R] R ANOVA gives diferent results than SPSS

. Either you specified a different model or different tests for the same model in the two systems, or you are interpreting the output incorrectly, or the results are different. Without more detail it is hard to be sure, but the first two possibilities seem more likely. -thomas Thomas Lumley

Re: [R] difftime result for days not an integer?

riods, due to daylight saving time: 0.958 is 23/24. Would you just round off? Yes, or use as.Date() if you only want to consider whole days R> as.Date('2004-08-05')-as.Date('2001-01-03') Time difference of 1310 days -thomas Thomas Lumley

Re: [R] expression() help

". Thanks, quote("Average PM"[10]) or quote(Average~PM[10]) seem to work. There is an example of the first style in ?plotmath. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle

Re: [R] R 64-bit memory.limit() is max of ?

algorithms that scale slower than linearly in the sample size will get very painful. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.o

Re: [R] survey package question

their SEs be computed for subgroups? You can also use subset(), which is what svyby() does internally svyquantile(~api00, quantile=c(0.25,0.5,0.75), subset(dclus1, stype=="E")) 0.25 0.5 0.75 api00 553 652 729 -thomas Thomas Lumley Assoc.

Re: [R] How to parse the arguments from a function call and evaluate them in a dataframe?

to create a model frame that will contain the variables in the formula, and columns `(ftime)` and `(fstatus)` for the other arguments. If you use formulas for ftime and fstatus you would have to call model.frame() multiple times, which is a bit more work. You would also need to use na.action=na.pa

[R] chairs for stat computing sessions at JSM

. The rewards are also pretty minimal: the gratitude of the ASA and possibly of the speakers. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-hel

Re: [R] How to take out the content of character string

s you a string with code in it, the only way to parse it is with parse(). The problem with eval(parse()) is when people use it to do macro processing or other things that are much better done with expressions. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum

Re: [R] How to take out the content of character string

On Wed, 10 Mar 2010, Henrique Dallazuanna wrote: On Wed, Mar 10, 2010 at 7:50 PM, Thomas Lumley wrote: On Wed, 10 Mar 2010, baptiste auguie wrote: Hi, it's generally considered a bad practice but try this, eval(parse(text=AA)) library(fortunes) fortune(106) HTH, baptiste On 10

Re: [R] mmap error-12, macbook pro

www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-projec

Re: [R] constrOptim - error: initial value not feasible

istinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle

Re: [R] Does S inherit the enhancements in R language?

at a function is desired. I think this has also been true in S-PLUS for a long time, but I'm not completely sure. There are just a few places where problems can occur, such as do.call() with an unquoted function name, where a user-defined variable can get in the way. -thomas Thomas

Re: [R] using a condition given as string in subset function - how?

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 and provide commented, minimal, self-contained, reproducible code. [[alternative HTML version deleted]]

Re: [R] Software lifecycle for R releases (aka practical limits of support for older versions)

Based just on R itself a longer update delay might be ok, but CRAN doesn't supply binaries of new or updated packages for old versions of R. Many packages will become seriously outdated much faster than base R. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum

Re: [R] questions about meta-analysis

read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle

Re: [R] multiple regression w/ no intercept; strange results

that your definition loses in the no-intercept case. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.

Re: [R] Large Stata file Import in R

worse in terms of memory than using read.dta, although it does provide the possibility of reading the data in chunks. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Se

Re: [R] Stata file and R Interaction :File Size Problem in Import

h of these least squares fits to a data set that was too big to fit in memory. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-proje

Re: [R] Large Stata file Import in R

in memory. Stata's -compress- option just chooses smaller data types, eg, byte instead of integer. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-he

Re: [R] system function question

= TRUE the value is 256 times the error code returned by the command, and if wait = FALSE it is 0 (the conventional success value). -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __

Re: [R] clogit comparison between Stata and R

On Wed, 8 Jul 2009, David Hugh-Jones wrote: Hello all I'm moving back and forth between stata and R at the moment - of course, using R whenever possible :-) I'm running conditional logits on some panel data and I get slightly different results and different N in the two programs. That's prob

Re: [R] (-8)^(1/3) == NaN?

solve? __ 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 and provide commented, minimal, self-contained, reproducible code. Thomas Lumley

Re: [R] what is meaning of the bubbles in boxplots?

www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-proje

[R] suggestions for JSM 2010 invited sessions

ed posters, and other aspects of the program are also welcome, though less urgent. For suggestions on topics other than statistical computing, you can find the rest of the Program Committee at http://www.amstat.org/meetings/jsm/2010/index.cfm?fuseaction=program -thomas

Re: [R] Random # generator accuracy

f ways to draw PPS samples without replacement. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

Re: [R] return.replicates in survey pkg

bute of the result. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

Re: [R] large dataset

the nearly ten years I've been saying it. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

Re: [R] Design of a survey using the "survey" package

"a guide to analysis using R", and you get what it says on the tin. You might look at the 'sampling' package, which has techniques for taking samples and which is used in some simulations in the book. -thomas Thomas Lumley Assoc. Professor, Biostatist

Re: [R] large dataset

On Mon, 29 Mar 2010, Gabor Grothendieck wrote: On Mon, Mar 29, 2010 at 4:12 PM, Thomas Lumley wrote: On Sun, 28 Mar 2010, kMan wrote: This was *very* useful for me when I dealt with a 1.5Gb text file http://www.csc.fi/sivut/atcsc/arkisto/atcsc3_2007/ohjelmistot_html/R_and_la rge_data

Re: [R] large dataset

f working with raw text files is something more than mere prescription. Sincerely, KeithC. -Original Message- From: Thomas Lumley [mailto:tlum...@u.washington.edu] Sent: Monday, March 29, 2010 2:56 PM To: Gabor Grothendieck Cc: kMan; r-help; n.via...@libero.it Subject: Re: [R] large dataset On

Re: [R] predicted time length differs from survfit.coxph:

sion deleted]] __ 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 and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, B

Re: [R] A questionb about the Wilcoxon signed rank test

85, 13.77, 295.17, 191.54, 126.44, 84.83) ? __ Make your browsing faster, safer, and easier with the new Internet Explorer[[elided Yahoo spam]] com/ca/internetexplorer/ [[alternative HTML version deleted]] Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washin

Re: [R] Help transfrom R to C

m. __ 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 and provide commented, minimal, self-contained, reproducible code. Thomas Lumley

Re: [R] SAS and R on multiple operating systems

would have thought Linux virtualization would be the most straightforward approach to using R. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-p

Re: [R] logistic regression in an incomplete dataset

enetics Institute d.campb...@ucl.ac.uk Tel. ext. 020 31084006, int. 54006 __ 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 and provide commented, minimal, self-containe

Re: [R] median of two groups

g/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing l

Re: [R] help in function in R akin to macro in SAS

s reserved use in the language. -- Stuart Luppescu University of Chicago __ 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 and provide commented, minimal, self-contained, reproducibl

Re: [R] Odd results with %% and conserving memory

as to be off by at least 1, no matter how high the numerical precision. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list

Re: [R] median of two groups

sometime people like me a bit confused! Regards, Cheba 2010/4/6 Thomas Lumley None of them.  - mood.test() looks promising until you read the help page and see that it does not do Mood's test for equality of quantiles, it does Mood's test for equality of scale para

Re: [R] Struggeling with svydesign()

be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher Indeed. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.or

Re: [R] Struggeling with svydesign()

to sampling probability). Also, how did you do the sampling? It's quite hard to do unequal probability sampling without replacement (the R sample() function doesn't actually do it, though the sampling package does). -thomas Thomas Lumley Assoc. Professor, Biost

Re: [R] Error in leaps.setup

been stripped off by the mailing list. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/l

Re: [R] erasing an area of a graph

ovide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/li

Re: [R] How to run Shapiro-Wilk test for each grouped

from a unique normal distribution, and only one normality test on the residuals is necessary (no need for a loop). Or, as some of us would say, *no* normality test on the residuals is necessary, which is even simpler. -thomas Thomas Lumley Assoc. Professor, Biostatistics

Re: [R] glmer with non integer weights

quasilikelihood chapter in McCullagh and Nelder, where the observations are the proportion of damage on a set of leaves. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle

Re: [R] Any chance R will ever get beyond the 2^31-1 vector size limit?

lock at a time.  Etc. Duncan Murdoch Best, Matt -- Matthew C Keller Asst. Professor of Psychology University of Colorado at Boulder www.matthewckeller.com __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://

Re: [R] Weights in binomial glm

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-contained, reproducible code. Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity

Re: [R] plotting RR, 95% CI as table and figure in same plot

MMP) 325 9th Avenue, 2HH-15 Box 359911 Seattle, WA 98104? 206-897-4210 http://www.chammp.org (Thurs) __ 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/pos

Re: [R] nls for piecewise linear regression not converging to least square

ith your graph. If I do points(x1,fitted(nlmod),col="red") I get points that are on the horizontal line segment, but then go through the data nicely on the right. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUnive

Re: [R] Results from clogit out of range?

p://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ 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/p

Re: [R] Generating a model fitness when score using svyglm?

Now, the deviance returned by svyglm() is scaled to the sample size, so if the survey design isn't informative it should be more or less in the same ballpark as a deviance from an independent sample, and the usual BIC calculation might give somewhat helpful results. -thomas Tho

Re: [R] Generating a model fitness score when using svyglm?

mpled data, then applying the same formula to svyglm() output will give a reasonable approximation to the same thing, but I wouldn't put much weight in small differences. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversit

[R] [R-pkgs] survey package (3.18)

PPS designs without replacement There is also experimental support for parallel processing using the 'multicore' package. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washingto

Re: [R] Fortran vs R

cluding in R) written in Fortran. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/

Re: [R] Newbie Question About Setting Plot Axis Limits

beyond 15. Is there some way to have the plot command >respect the requested ylim and xlim settings? par(xaxs="i",yaxs="i") plot(c(1,1), ylim=c(0, 15), xlim=c(0, 13)) -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.e

Re: [R] scaled Schoenfeld residuals

alues. Because their purpose in life is to be smoothed against time to get an estimate of the parameter as a function of time (plot.cox.zph). -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Se

Re: [R] xtable - print - suppress output

, otherwise it is the captured output. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Function to find prime numbers

create a vector of small primes and use my_array %in% smallprimes. For example, the first 1000 primes are at http://primes.utm.edu/lists/small/1000.txt and the first 1 are on the same site. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum

Re: [R] svy / weighted regression

, self-contained, reproducible code. David Winsemius, MD Heritage Laboratories West Hartford, CT __ 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/postin

Re: [R] post-hoc test with kruskal.test()

___ 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 and provide commented, minimal, self-contained, reproducible code. Thomas Lumley Assoc. Pro

Re: [R] Generating plots with two lines, each line plotted with a different color

(Surv(Age,Died) ~ Rx, data = GVHDdata),lty=c(1,2),lwd=2) ,col=c("tomato","goldenrod") -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle ___

Re: [R] Merging two data frames with 3 common variables makes duplicated rows

45F clas 03:46:15 645F free 02:20:15 745H clas 02:30:07 845H free 01:59:36 938F clas 02:43:17 10 38F free 02:24:46 Lines 9 and 1 appear to be the same in meil2, as do 2 and 10. If the 16 rows consist of two repeats of 8 rows that would expla

Re: [R] plot(survfit(fitCox)) graph shows one line - should show two

On Sun, 10 May 2009, John Sorkin wrote: R 2.8.1 Windows XP I am trying to plot the results of a coxph using plot(survfit()). The plot should, I believe, show two lines one for survival in each of two treatment (Drug) groups, however my plot shows only one line. What am I doing wrong? Expect

Re: [R] How does one plot the -log(log(survival))

I would use cox.zph() and plot.cox.zph(), which give valid formal tests and reasonable estimates of how the log hazard ratio varies with time. Still, if you want to do log-log plots, the survival package is up to the task. -thomas Thomas Lumley Assoc. Professor, Biosta

Re: [R] (no subject)

and binom.test() gives an exact confidence interval. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.

Re: [R] Error with regsubset in leaps package - vcov and all.best option (plus calculating VIFs for subsets)

iance matrices for any of the best models. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum...@u.washington.eduUniversity of Washington, Seattle __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-h

Re: [R] Sample size calculation proportions with EpiR: Discrepancy to other calculators

;). epi.studysize(,method="cohort") doesn't seem exactly appropriate, since judging from the example on the help page the inputs are supposed to be cumulative incidence rather than probabilities. -thomas Thomas Lumley Assoc. Professor, Biostatistics tlum..

Re: [R] (OT) Does pearson correlation assume bivariate normality of the data?

alienetworks.com/srtest.cfm Do you know how to write? http://garbl.home.comcast.net/~garbl/stylemanual/e.htm#e-mail __ 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/posti

  1   2   3   4   5   6   >