similar to: bug in power.t.test( ) (PR#7245)

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

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

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

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

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

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

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

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

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 >

When legend() is used with the angle argument as follows, not only the boxes beside the legend text, but also the whole legend box is filled with shading lines. I think this is not intended: plot(1:10) legend(8, 4, c("A", "B"), angle=c(10, 80), fill=NULL, density=20) I'd suggest as a fix (legend.R of R-1.5.0): 25c25 < rect2 <- function(left, top, dx, dy,

Problem: neither fg or bg nor color work properly in dotchart. version: R-1.7.1 for windows code which shows the errors: x <- matrix(rnorm(16),ncol=2,dimnames=list(paste("a",1:8,sep=""),c("before","afte r"))) dotchart(x,fg="blue",bg="lightgrey") dotchart(x,color=c("red","blue")) Dr Ian J Wilson Lecturer in

The package ipred is uploaded to CRAN. The main focus of the package is the calculation of improved predictors in classification tasks. Misclassification error can be improved by bootstrap aggregated classification trees and/or the framework of indirect classification. Furthermore, a unified interface for the estimation of misclassification error completes the features of ipred. We try to make

The package ipred is uploaded to CRAN. The main focus of the package is the calculation of improved predictors in classification tasks. Misclassification error can be improved by bootstrap aggregated classification trees and/or the framework of indirect classification. Furthermore, a unified interface for the estimation of misclassification error completes the features of ipred. We try to make

Code like df <- data.frame(x=1:10) y <- 20:29 eval(quote(x+y), env=df) does what you might expect: it looks for x and y in the data.frame, and when it doesn't find y there, it looks in the parent environment. However, sometimes I'd like to construct a single environment out of df, so that I can pass it to nested functions and get the same behaviour. Right now, I get the wrong

Code like df <- data.frame(x=1:10) y <- 20:29 eval(quote(x+y), env=df) does what you might expect: it looks for x and y in the data.frame, and when it doesn't find y there, it looks in the parent environment. However, sometimes I'd like to construct a single environment out of df, so that I can pass it to nested functions and get the same behaviour. Right now, I get the wrong

Full_Name: Tomas Larsson Version: 2.6.0 OS: Windows XP Submission from: (NULL) (198.208.251.24) This is not a bug, it is a performance issue but I think it should have an easy fix. I have a large matrix (about 2,000,000 by 20), when I type colnames(x) it takes a long time to get the result. However, if I select just the first couple of rows of the matrix I don't have to wait for the

Full_Name: Guan-Hua Huang Version: 2.0.1 OS: Linux Submission from: (NULL) (140.113.114.123) I install the package cluster by using install.packages("cluster"). After install it, it runs fine for function clara, but it does not work for function fanny. I did the following things: library(cluster) set.seed(21) x <- rbind(cbind(rnorm(10, 0, 0.5), rnorm(10, 0, 0.5)),