similar to: Rounding problem R vs Excel.

Hi again Sam, Sorry for not replying sooner. I have been ignoring everything for a while to work on the AMMP analysis tool and the continuing saga of migrating the legacy AMMP data. Our consultant has written some programs that take a VERY long time to run :( Anyway, I now know Java and a certain amount of JSP. I hope to be able to get around to this soon. It is all a bit hard now with our

For numbers ending in 5 Excel always rounds up and I want to reproduce this in R rather than use the round function. I have tried: x<-as.integer(x*100+0.5)/100 and x<-floor(x*100+0.5)/100 However, some values of x cause problems e.g x<-floor(4.145*100+0.5)/100 4.14 I have tried breaking it down into steps and force the results to 10 significant figures in the following function,

Your surest bet is to look into S+, not R, since the former does 'integrate' with Microsoft applications. To what extent, I don' know. This is a desirable feature (IMO) and a reason to use S+ rather then R (already discussed in different thread). Try the S+ newsgroup. > -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at

Dear all i use R only a few days and don't understand the difference between fivenum(x) und summary(x). > x [1] 20.77 22.56 22.71 22.99 26.39 27.08 27.32 27.33 27.57 27.81 28.69 29.36 [13] 30.25 31.89 32.88 33.23 33.28 33.40 33.52 33.83 33.95 34.82 > fivenum(x) [1] 20.770 27.080 29.025 33.280 34.820 > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 20.77 27.14

Hi, I'm wondering if anyone can point me to a function that will allow me to do a ranking that treats ties differently than rank() provides for? I'd like a method that will assign to the elements of each tie group the largest rank. An example: For the vector 'v', I'd like the method to return 'rv' v: 1 2 3 3 3 4 5 5 6 7 rv: 1 2 5 5 5 6 8 8 9 10 Thanks,

Prof Brian Ripley wrote: > > Yes, we need one. I am having problems with Erich's latest add-in: > it works about every other time, I think because is has some memory of the > previous one. I don't think Excel 2000 removes add-ins cleanly. > brian, we had a stress test of the addin yesterday election ight forecast, the systen was running about 8 hours without any problems.

The function 'wilcox.test' in R and S gives (almost) identical results (see below). 'qwilcox' however, does not: > qwilcox(p,5,5) p: 0.025 0.975 -------------------- R> 3 22 S> 18 37 I originally wanted to ask a questions, but then I found the answer. Given the confusion I run into, I wonder if this experience is worth reporting. The

Hello, I am analysing several samples whose sizes are from 9 to 110. I would like to test their distribution with R, whether they are normal or not. I wonder which test for normality from R should I use . Thank you for help. Samuel Samuel BERTRAND Doctorant Laboratoire de Biomecanique LBM - ENSAM - CNRS UMR 8005 151, bd de l'Hopital 75013 PARIS Tel. +33 (0) 1 44 24 64 53 Fax +33 (0) 1

Hi, I would like to do a MannWitney test. Can anyone help me with the propper command? Thanks, Margarida

Hello, I've tried to analyze some data with a CMH test. My 3 dimensional contingency tables are 2x2xN where N is usually between 10 and 100. The problem is that there may be 2 strata with opposite counts (the 2x2 contigency table for these are reversed), producing opposite odds ratios that cancle out in the overall statistics. These opposite counts are very important for my analysis, since

Hi, I use ks.test or lillie.test to verify a normal distribution. It's performed for a group My users use SigmaStat software and a One Way ANOVA on several groups In the result page there is a probability value to determine if Normality test is failed or passed So, how can i retrieve this probability value on several groups? Is there another function in R to verify normality on several

Hi. I'm going to introduce the R-package for a group of medical doctors later this week and is a little confused about there use of a test named "willcoxon-pratt" for testing if the clinical and biochemical markers has decreased significantly after the use of some drugs for a group of patients. Looking into the R-functions I would in R recommand using a matched-pairs Wilcoxon

I am trying to `cut' a continuous variable into contiguous classes containing approximately an equal number of observations. I thought quantile() was the appropriate function to use in order to find the breakpoints, but I end up with classes of different sizes - see example below. Does anybody have an explanation for that? And what is the `recommended' way of computing what I am looking

Hi, I wonder if anybody could help me in converting this easy SAS program into R. (I'm still trying to do that!) PROC IMPORT OUT= WORK.CHLA_italian DATAFILE= "C:\Documents and Settings\carleal\My Documents\REBECCA\stat\sas\All&nutrients.xls" DBMS=EXCEL2000 REPLACE; GETNAMES=YES; RUN; data chla_italian; set chla_italian;

I remember in one slide of Prof. Ripley's presentation overhead, he said the most popular data analysis software is excel. So is there any resource or tutorial on this topic? Thank you so much!

Last Friday, Gregory Chaitin (http://www.umcs.maine.edu/~chaitin/lm.html) mentioned that there can be no proof that a given code is the shortest for a problem, even within a language. Still, the script below, a replacement of the "TDT", one of the most frequently used tests in genetics (http://mustat.rockefeller.edu under "downloads") may get close. It contains a few

G'day, I'm?calculating LSM for the following model, and am finding that R and SAS give different answers. Whilst the error is at the second or third decimal, the percentage error can be quite large. I'm using the effects library (Version: 2.0-10) on R?2.11.1 in the following manner: options(contrasts=c("contr.helmert","contr.poly"))

Hello all... I am really new to statistics and I am trying to figure out a way to apply Chauvenet's criterion using the t-distribution on a set of numbers in perl. I was unable to find a TDIST and TINV function for perl. I am getting these functions from Excel. So, I figured that I would install R and call it from perl, overkill for what I need I know. I am having a hard time figuring

Why does the assignment of a 3178x93 object to another 3178x93 object remove the dimension attribute? > GT <- array(dim = c(6,nrow(InData),ncol(InSNPs))) > dim(GT) [1] 6 3178 93 > SNP1 <- InSNPs[InData[,"C1"],] > dim(SNP1) [1] 3178 93 > SNP2 <- InSNPs[InData[,"C2"],] > dim(SNP2) [1] 3178 93 > dim(pmin(SNP1,SNP2)) [1] 3178 93

Hi, I just try to draw multiple plots in one page using R, I used par command. For example I have 7 plots, but instead of arranging them in the default way plot1 plot2 plot3 plot4 plot5 plot6 plot7 I want them in this order plot1 plot2 plot3 plot4 plot5 plot6 plot7 Could somebody tell me how to do this, please? Thanks so many. Suyan