search for: wittkowski

...Est(0:nP,nP:0, 0:nX,nX:0, 0:nQ,nQ:0)^2) return(1-sum(tb[tb[,2]<c(Est(pP,qP,pX,qX,pQ,qQ)^2),1])) } #-------------------------------------------------------------------------- # Note: The numbers in parentheses in inline comments refer to the equation # numbers in # Wittkowski KM, Liu X # A Statistically Valid Alternative to the TDT # Hum Hered 2002, 54:157-164; 2004, 58:60-61 (discussion) #-------------------------------------------------------------------------- Knut M. Wittkowski, PhD,DSc ------------------------------------------ The Rockefeller Univers...

...,] > dim(SNP2) [1] 3178 93 > dim(pmin(SNP1,SNP2)) [1] 3178 93 > GT[1,,] <- pmin(SNP1,SNP2) > dim(GT) NULL # why?????????????????????????????????????? > GT[2,,] <- pmax(SNP1,SNP2) Error in GT[2, , ] <- pmax(SNP1, SNP2) : incorrect number of subscripts Knut M. Wittkowski, PhD,DSc ---------------------------------------------------- Sr Research Associate, The Rockefeller University Biometrician, The Rockefeller University Hospital Director, Biostatistics/Epidemiology/Research Design Center for Clinical and Translational Science (CCTS) apl. Professor, Eberhard...

...ot;) Wilcoxon rank sum test: p-value = 0.05556 R> wilcox.exact(x1,x2,alternative="two.sided") Exact Wilcoxon rank sum test: p-value = 0.05556 S> wilcox.test(x1,x2,alternative="two.sided") Exact Wilcoxon rank sum test: p-value = 0.0556 Knut M. Wittkowski, PhD,DSc ------------------------------------------ The Rockefeller University, GCRC 1230 York Ave #121B, Box 322, NY,NY 10021 +1(212)327-7175, +1(212)327-8450 (Fax) kmw at rockefeller.edu http://www.rucares.org/statist/

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

When all else fails, read the help page... ?fivenum says to look at ?boxplot.stats, and the "Details" section of ?boxplot.stats has, well, details. Tukey had reasons to call those hinges rather than quartiles. Andy > -----Original Message----- > From: Knut M. Wittkowski [mailto:kmw at rockefeller.edu] > Sent: Thursday, May 29, 2003 5:58 PM > To: Matthias Kirschner > Cc: R-help at stat.math.ethz.ch > Subject: Re: [R] R summary (and quantiles) > > > Matthias > > The function "fivenum", defines quantiles by assuming that >...

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

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

David A. Paul wrote: > I don't have the reference, but a biologist friend of mine once > showed me a refereed journal article that purported to demonstrate > numerical errors made by MSExcel. This would have been Excel97 or > Excel2000... In any case, the journal's scope was biological in > nature and the article was of interest since Excel is heavily used in > that

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,

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

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