similar to: binary representation of an integer?

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

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

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

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

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

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

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 all, Basically, I have data in the format of (up to 1 gig in size) text files containing stuff like: F34060F81000F28055F8A000F2E05EF8F000F34 (...) The data is basically strings denoting hex values (9 = 9, A = 10, B = 11, ...) organised in fixed, small blocks. What I want to do is to read in a specified segment of the string, break it up into blocks, and convert it into a vector of integers

> 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

I installed rw0641 from the binary rw0641b.zip but not the help files (not rw0641h, not rw0641w ). I do not like the MDI so I edited the ./etc/Rconsole to read MDI = no and I commented out #toolbar = yes #statusbar = yes and I even changed pagerstyle = singlewindows Now when I start Rgui, and do > ?ls causes R to crash, instead of reporting "help page not found". Is there

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)

Hi, all, I want to use asterisk SIP as a proxy, and two softphone (Ubiquity SIP Phone) as user agent, like below: Softphone1<-------------->Asterisk SIP<------------>Softphone2 (User Agent) (Proxy) (User Agent) 155.69.xx.xx 155.69.yy.yy 155.69.zz.zz zhou mysipproxy.com