similar to: survival5

I downloaded binary of survival5 for windows95 (also the splines, date). The function survsum shows up in content, has a help page, but the function itself is missing. Can we cut and paste the survsum from survival4 ? Mai Zhou P.S. I am not picking on R, not even complaining. I think R is great and I am try to contribute my (small) part to make it even better in way of bug report (or

Dear r-help, I ran output <- lm(formu, data=mydata) I want to look at the condition number (to see if the matrix is near singular). How? Also, I use the function stepAIC from MASS to select models, how can I see the condition numbers in each step? While I am at it. I find a minor bug: R 1.7.1 on window XP, when you copy and paste by click one button, then the input stopped, until you

Dear R community, I wrote and someone else helped me type a 5-page R-info in Chinese, covering the most basic information on how to get and install and start using R on Wintel PC. This could serve those not fluent in/intimidated by English but want to start using R. (I met a lot of those in China) <http://www.ms.uky.edu/~mai/ZhongWen.htm> also a MSword version of the same (may print

OK, I am forwarding this to the r-help. In the original message I simply hit reply and I thought I send to ALL r-help list. (apparent not) also I recieved another piece of mail.....about the RNG in Splus 4.0 or later indicating it was modified and getting worse. I only have Splus 3.4 for HP. So I cannot try it. This makes the idea of using the pre-generated RN on CDROM more appealing.. Mai

Full_Name: Mai Zhou Version: 1.9.1 OS: Win XP Professional Submission from: (NULL) (12.222.227.93) > power.t.test(n=25, delta=0.1, sig.level=1.1, strict=TRUE, type="one.sample") One-sample t test power calculation n = 25 delta = 0.1 sd = 1 sig.level = 1.1 power = 1.088311 alternative = two.sided ### power can never be

Dear helper, I am using linux (redhat) R 1.5.0, in the batch mode. The device png() do not seem to be available. A look at the help says png() is available if I use --gui option. Is there a way to get around it? (to get png("Rplot.png") in batch mode?) The reason I need batch mode is that I am working to get the Rcgi to run. Thanks, Mai Zhou mai at ms.uky.edu

survexp in survival5 doesn't seem to work for me. see below: > library(survival5) Attaching Package "package:survival5": The following object(s) are masked from package:base : sort.list > library(chron) > data(ratetables) > survexp(~ratetable(year=julian(6,1,1991), + sex=1,age=35*365.24),times=(0:30)/6*365.24) Error in as.character(as.date(c(min(R[, 3]),

Hi, I was thinking to ask the University computing Labs here to install R and I was told they need a statement that the software (R) is Y2K compliant. Which I cannot find in CRAN. (in addition to the licensing info). Is it a big issue to put such a statement in R? or this has legal implications that R do not want as a "free" software? If possible I would like to see the statement on

Dear R-Helpers: What other references are there on the capabilities of the survival5 package other than the help files and the chapter on survival analysis in every edition of Modern Applied Statistics with S? I'm thinking of something like "An Introduction to Survival Analysis in R" with worked examples that might complement or extend the chapter in MASS. Thanks, Spencer

R 0.64 on windows NT 4.0 Sometimes I got an error message by doing this > .Random.seed <- c(1, 1:2) > .Random.seed [1] 1 1 2 > runif(5) Warning: Wrong length .Random.seed; forgot initial RNGkind? set to Wichmann-Hill[1] 0.02253721 0.84832584 ........ Sometimes I do not get error message: > .Random.seed <- c(1, 1:2) > .Random.seed [1] 1 1 2 > runif(1) [1] 0.5641106 >

Dear list, Did anybody have something like multiroot(), similar to uniroot() only for several varibles/function? I know the difficulty several varibles bring. But for nice functions, it can be useful. We are thinking of writing one....I thought to ask first. Mai Zhou -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello all. I've been attempting to utilize the "survival" pkg ( http://cran.r-project.org/web/packages/survival/index.html), while reading through this guide (http://www.ms.uky.edu/~mai/Rsurv.pdf). I figured working through the guide would be the best way to go, before attempting my own data. I tried to utilize the Kaplain-Meier estimator as shown in the guide:

Anyone using the penalised partial likelihood routines in survival5 should update their version. A bug has been fixed in the S package: in coxph() models with penalised likelihood and strata it was possible in some circumstances to get an infinite loop or perhaps an incorrect answer. The new version (2.3) is on cran.r-project.org and will percolate through CRAN in the next few days. -thomas

Anyone using the penalised partial likelihood routines in survival5 should update their version. A bug has been fixed in the S package: in coxph() models with penalised likelihood and strata it was possible in some circumstances to get an infinite loop or perhaps an incorrect answer. The new version (2.3) is on cran.r-project.org and will percolate through CRAN in the next few days. -thomas

Dear R-help, > prop.test(9, 137, p=0.066) > prop.test(9, 137, p=0.05) give two different 95% confidence intervals. I thought the confidence interval calculation should be independent of testing calculations (and thus the null hypothesis)? Splus 2000 has similar problems but give slightly different answer. Using R1.3.0 on windows. Mai Zhou

Windows 9*/NT version, on an NT box R 0.63.0 Beta (Nov 13, 1998) > library(survival4) > test <- function(x,d) { survfit(Surv(x,d) ~1) } > test( rexp(50), sample(c(0,1), size=50, replace=T) ) Error: Object "x" not found > Other functions like in the ctest library works out OK, the x and d got passed to the inner function. And here it did not got passed to

R 0630 for windows > library(ctest) > binom.test(7,10,p=0.3, alternative="two.sided") returns a p-value of =< 2.2e-016 and a warning In Splus 3.4 > binom.test(7,10,p=0.3, alternative="two.sided") returns a p-value of 0.0106 I think it is the max(v[v<=(1+eps)*PVAL]) causing the problem... max() of an empty vector....... Mai Z

Hi, everyone, I was downloading the rm0641b.zip etc. (thanks to Guido Masarotto and Brian Ripley) I noticed rw0641sp.zip there, which do not seem mentioned in the News file. My guess is that it is not a binary file? Is is a source patch file(?) for compiling on Wintel PC? Mai Zhou -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Is there a function or other ways in R that can find the binary representation of an integer? (or list all those representations for integers 0:63, for example) OR do I have to built such a table? i.e. something like binary(3)= 11 binary(4)= 100 etc. Thanks. Mai Zhou -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list

I have installed R version 1 for windows and have installed the survival5 package in the library folder, but the base package only seems to partially acknowledge survival5's presence. How have I screwed up? The diagnostics are as follows... >library(survival5) ** no warning given.. seems OK ** >help(package=survival5) **this gives the contents of the package correctly,