search for: echip

...[2]<-"us" > citizen [1] uk us no au uk us us Levels: as au no uk us ** [STICKY "as"] ** > citizen[2]<-"gs" Warning message: invalid factor level, NAs generated in: [<-.factor(*tmp*, 2, value = "gs") -- Bob Wheeler --- (Reply to: bwheeler@echip.com) ECHIP, Inc. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !)...

...running Windows NT4. Rcmd does not seen to read the command line. Rterm and other things work well. For example "Rcmd check --help" fails with the command interpreter trying to find a program or file called "check." What have I missed? -- Bob Wheeler --- (Reply to: bwheeler at echip.com) ECHIP, Inc. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To...

...ation has little effect the variability may be large, while if it has a dramatic effect (kills everything!), there may be no variability to speak of. John Maindonald. On 22 Feb 2005, at 10:06 PM, r-help-request at stat.math.ethz.ch wrote: > To: "'Bob Wheeler'" <bwheeler at echip.com> > Cc: r-help at stat.math.ethz.ch > Subject: RE: [R] power.anova.test for interaction effects > Reply-To: akniss at uwyo.edu > > > It's a rather complex model. A 37*4 factorial (37 cultivars[var]; 4 > herbicide treatments[trt]) with three replications[rep] was carr...

Dear R-helpers, Would anybody have a compiled binary version (for Windows) of the Dopt library that is in the development site?. best regards, IOsu -------------------------------------------------------------- Jes?s Mar?a Fr?as Celayeta Dublin Institute of Technology School of Food Sci. and Env. Health. Cathal Bruha St., Dublin 1. Ireland Phone: +353 1 4024459 Fax: +353 1 4024495

Hello everyone. If I do f <- function(n){max(rnorm(n))} plot(sapply(rep(5000,4000),f)) #[this takes my PC about 30 seconds] then I get something quite unexpected: gaps in the distribution. For me, the most noticable one is at about 3.6. Do others get this? Is it an optical illusion? It can't be right, can it? Or maybe I just don't understand the good ol' Gaussian very

Dear R-people Does someone have a routine to ngenerate inverse-gaussian random numbers. I am thinking of something similar to rinvgauss, pinvgauss etc. in S-plus. best regards Helgi -- Helgi Tomasson FAX: 354-552-6806 University of Iceland PHONE:354-525-4571 Faculty of Economics and Business Administration

Hi folks: Is there a lookup function that returns the variate given the cumulative probability for an object returned by the density(...) function? > mydata _ as.vector(mymatrix) > mydata.density _ density(mydata) > mydata.p80 _ lookup(mydata.density, p=0.8) # is there any function to accomplish this task? Thanks. Rajiv. -------- Rajiv Prasad, Postdoctoral Research Associate,

Dear People, Consider the following fragment of R code choose <- function(n,r) { if( is.integer(n) && is.integer(r) && n > 0 && r >= 0 ) { .C("choose",as.double(n),as.double(r),comb = double(1))$comb } else stop("n must be a positive integer and r a non-negative integer.") } This is a practice function (n choose r), which I

...or this purpose. I'd normally submit this sort of thing to beta test, but I guess the R users are the beta testers, so have at it. It would probably be best to send the bug reports directly to me rather then clutter up the newsgroups. -- Bob Wheeler --- http://www.bobwheeler.com/ ECHIP, Inc. --- Randomness comes in bunches. _______________________________________________ R-packages mailing list R-packages at stat.math.ethz.ch https://www.stat.math.ethz.ch/mailman/listinfo/r-packages

...et.al (2005). A comment on the implementation of the Ziggurat Method. Jour. Stat. Sci. 12-7, 1-4. The fix takes care of problems for very large sets of random values. I tested it a bit on small samples and it seems ok. It will be on CRAN shortly. -- Bob Wheeler --- http://www.bobwheeler.com/ ECHIP, Inc. --- Randomness comes in bunches. _______________________________________________ R-packages mailing list R-packages at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-packages

...u> The t-distribution family may help. These have zero skewness and kurtosis ranging from the normal to infinite. He doen't know offhand the relationship between degrees of freedom and kurtosis, but it should be easy to look up (or evaluate by simulation). 6. Bob Wheeler <bwheeler at echip.com> We can use the Johnson system. See the SuppDists package. Abdul Kudus ===================== Dept. of Statistics Bandung Islamic University Indonesia ==================== -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://w...

...u> The t-distribution family may help. These have zero skewness and kurtosis ranging from the normal to infinite. He doen't know offhand the relationship between degrees of freedom and kurtosis, but it should be easy to look up (or evaluate by simulation). 6. Bob Wheeler <bwheeler at echip.com> We can use the Johnson system. See the SuppDists package. Abdul Kudus ===================== Dept. of Statistics Bandung Islamic University Indonesia ==================== -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://w...

...or this purpose. I'd normally submit this sort of thing to beta test, but I guess the R users are the beta testers, so have at it. It would probably be best to send the bug reports directly to me rather then clutter up the newsgroups. -- Bob Wheeler --- http://www.bobwheeler.com/ ECHIP, Inc. --- Randomness comes in bunches. _______________________________________________ R-packages mailing list R-packages at stat.math.ethz.ch https://www.stat.math.ethz.ch/mailman/listinfo/r-packages

cor.test and SuppDists give me different P-values for the same Spearman's rho. Which is correct, or am I doing something wrong? > x <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1) > y <- c( 2.6, 3.1, 2.5, 5.0, 3.6, 4.0, 5.2, 2.8, 3.8) > cor.test(x,y,method="spearman") Spearman's rank correlation rho data: x and y S = 48, p-value =

Hello I am searching ghyper package (generalized hypergeometric distributions). Does anyone can send it to me? Regards from Mexico Lizbeth Román [[alternative HTML version deleted]]

Hello, My name is Lisa and I'm a statistician at Princess Margare Hospital. I wonder if there is any function in R that calculate the Normal rankits based on Blom's approximation? Thank you very much Lisa Wang Msc. Princess Margaret hospital Toronto, Ca

...s a single folder "doc," but apparently any folder causes the problem. If inst is empty, the project checks OK. This was not a problem before 1.9. I've checked the documentation, but don't see a change. What am I missing. -- Bob Wheeler --- http://www.bobwheeler.com/ ECHIP, Inc. --- Randomness comes in bunches.

Hi list, I am building a regression model with categorical predictor variable coded by treatment contrasts. The summary of the regression model shows that some levels are significant while others are not. The significant ones show that they are statistically significant from the basis factor (at level 0). How do we generally deal with those insignificant levels? If they are used for

Hello, Could somebody point me, is there any function in R which returns expected values of order statistics for normal distribution? I have been looking and couldn't find it. Thanks! Anna

...NEW_NUMERIC(1)); NUMERIC_POINTER(rvector)[0]=0.5; SET_ELEMENT(alist,2,rvector); UNPROTECT(1); PROTECT(kvector=NEW_INTEGER(1)); INTEGER_POINTER(kvector)[0]=k; SET_ELEMENT(alist,3,kvector); UNPROTECT(1); UNPROTECT(1); return alist; } -- Bob Wheeler --- http://www.bobwheeler.com/ ECHIP, Inc. --- Randomness comes in bunches.