similar to: Y2K compliance of R

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

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 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

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

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:

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

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

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 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

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

Hi, I'd like to execute simple commands and functions in R through a website, is there any service like this somewhere? I only found http://www.osvisions.com/r-online/ but it does not work (last update 2003) and the links to releated websites only give errors (if I calculate 7+3). Thanks for help & hints, Thomas

Hi, dear R-help, I want to fit a Buckley-James censored regression without intercept. The function bj() inside the Design library have to have an intercept, I try Surv(y,d)~ 0 + x , or -1+x, or x-1, none works. Any suggestions? Is there another BJ() code out there that do not have to have an intercept? Or may be it is easier to modify the bj()? I need the estimator only, no need of var estimate

I have commited fixes for 0.64.1 for (1) From: Mai Zhou <mai@ms.uky.edu> > .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 ........ Here the length of the seed was being tested before the kind was picked out, so the length of the previous type was used. (2)

Gregory Hosler <greg@camelot.eno.ericsson.se> wrote: > Everything was running smoothly, and then I applied the solaris "y2k" patches. > Now, most accesses to samba will hang for a minute or 2, before completing. In > the event log I see messages like: > > The redirector has times out a request to <server name> > > This did not used to happen.

Dear R-help, Does anybody know where can I find R/Splus code for computing PREdiction Sum of Squares (PRESS) in a linear regression model? I have a large regression model (~100 predictors, ~30,000 cases) and good prediction is the main goal. I remember there is a faster way to compute it (no need to repeatedly fit models), but if there is existing code...... Thanks. Mai Z

I came across this: shouldn't the last value be a more exact zero? It did not do that with 1 - sum( rep(0.1, 10) ) > 1 - cumsum( rep(0.1, 10) ) [1] 9.000000e-01 8.000000e-01 7.000000e-01 6.000000e-01 5.000000e-01 4.000000e-01 3.000000e-01 2.000000e-01 [9] 1.000000e-01 1.110223e-16 > version _ platform

> data(package=survival5) Warning: no data listing found Warning: package `survival5' contains datasets but no index The data loads OK. Mai Zhou -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

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

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