similar to: New Ranking Scheme - Help

Hi all, I am a biostatistician and I have developed my own ranking system for clinical data. I would like to test the efficiency of it w.r.t. to other ranking systems. I would like to simulate the data and after assigning ranks to my observed scores(after neglecting dropouts), observe the type I error. If I want to do a Kruuskal Wallis type of test, what test statistic should I use to test for a

Hi all, I am a biostatistician and I have developed my own ranking system for clinical data. I would like to test the efficiency of it w.r.t. to other ranking systems. I would like to simulate the data and after assigning ranks to my observed scores(after neglecting dropouts), observe the type I error. If I want to do a Kruuskal Wallis type of test, what test statistic should I use to test for a

Hello, I need help in performing a Van_der_Waerden normal scores test in R. I have two arrays of scores(final on therapy scores from drug and placebo) and want to use the normal scores procdeure to test for significance. (observations are unequal in number - due to dropouts). Could you please help me out with the coding or let me know if there is a package that can be used (for example,

Hello, I need some help in performing a Van_der_Waerden normal scores test in R. I have two arrays of scores(final on therapy scores from drug and placebo) and want to use the normal scores procdeure to test for significance. (observations are unequal in number - due to dropouts). Thanks in advance Regards ~S

Hello All, I want "Ranks' Table in R as like in SPSS ouput in the given link. http://www.statisticssolutions.com/methods-chapter/statistical-tests/kruskal-wallis-test/ Is the code is already available? Please let me know. Thanks, Lawrence

Hello all, Here goes one of my first functions. I want to make a nonparametric multiple sample comparison with unequal sample sizes (see Zar?s Biostatistical Analysis, 3rd. Ed., pg. 201 Example 10.11, pg. 288 Example 11.10). In the real world, I want to compare samples of fish length captured with different fishing gears. After using the Kruskal-Wallis test I want to check the differences

how do I check the rank of a matrix ? say A= 1 0 0 0 1 0 then rank(A)=2 what is this function? thanks I did try help.search("rank"), but all the returned help information seem irrelevant to what I want. I would like to know how people search for help information like this. rank(base) Sample Ranks SignRank(stats) Distribution of the

Hello, I have some patient data for my masters thesis with three groups (n=16, 19 & 20) I have completed compiling the results of 7 tests, for which one of these tests the variances are unequal. I wish to perform an ANOVA between the three groups but for the one test with unequal variance (<0.001 by both bartlett and levene's test) I am not sure what to do. I thought i would run

Well fellow R users, I throw myself on your mercy. Help me, the unworthy, satisfy my employer, the ungrateful. My feeble ramblings follow... I've searched R-Help, the R Website and done a GOOGLE without success for a one way ANOVA procedure to analyse data that are both non-normal in nature and which exhibit unequal variances and unequal sample sizes across the 4 treatment levels. My

Dear fellow R users, Keywords: Kruskal-Wallis, Post-Hoc, pair-wise comparisons, Nemenyi-Damico-Wolfe-Dunn test, coin package, oneway_test I am using the "oneway_test" function in the R package "coin" and I am obtaining results which I cannot believe are accurate. I do not wish to waste anyone's time and so if the following problem is rather trivial, I apologize, however I

Hi, Is anyone aware of a rank-based, non-paired test such as the Krustal-Wallis, that can accommodate weights? Alternatively, would it make sense to simulate a dataset by duplicating observations in proportion to their weight, and then using the Krustal-Wallis test? thanks! Dylan

Suppose I have a data frame with two columns ``district'' and ``score'' --- score is numeric; district may be considered categorical. I wish to append to this data frame a third column whose entries are the ranks of ``score'' ***within*** district. I've tried fiddling about with tapply() and by() but the result is a list whose i-th component consists of the ranks of

Hi there, I am doing multiple comparisons for data that is not normally distributed. For this purpose I tried both functions kruskal{agricolae} and kruskalmc{pgirmess}. It confuses me that these functions do not yield the same results although they are doing the same thing, don't they? Can anyone tell my why this happens and which function I can trust? kruskalmc() tells me that there are no

I would like to make a ranking operation more efficient if possible. The goal is to rank a set of points representing objective function values such that points which are "dominated" by no others have rank 1, those which are dominated by one other point have rank 2, etc. In the example with two dimensions below, objective functions 1 and 2 are to be minimized. Points a-e are

Hello List, I'm trying to prepare some lecture notes on non parametric methods, and I can't manually reproduce the results of the wilcox.test function for ordinal data. The data I'm using are from David Howell's website, available here http://www.uvm.edu/~dhowell/StatPages/More_Stuff/OrdinalChisq/OrdinalChiSq.html If I run the wilcox.test function on the data I get a p-value of

Hi, I'm playing with moving some of my lattice graphics into ggplot2, and I'd like to ask how to achieve a couple of things, both of which are fully illustrated in self-contained code (and mostly minimal, although that left quite a bit) following this written description. 1. I quite often like to use a 'ghosted' reference across facets - for example, in my example program below,

In http://tolstoy.newcastle.edu.au/R/e2/help/06/10/3064.html a method was given for converting a frequency table to an expanded data frame representing each observation as a set of factors. A slightly modified version was later included in the NCStats package, only on http://rforge.net/ (and it has too many dependencies to be useful). I've tried to make it more general, allowing an input

I think I've gotten my Emacs/Sweave/R system set up correctly, thanks to Vincent and Jim, but I haven't been successful getting my first document produced. I'm trying to use one of Friedrich Leisch's examples, http://www.ci.tuwien.ac.at/~leisch/Sweave/example-1.Snw. I cut and pasted the text into a document sweaveexample.Rnw in Emacs. It seemed to be processed successfully with R:

Dear All, I met problems when doing contrast and now really need some help in the model below: Fit=gam(y~treat+SEQUENCE+PERIOD+SEX+s(x),data=dat, random=list(SUBJID=~1),correlation=corAR1(form=~1|SUBJID)) And error message keeps coming out when I want to compare the differences between treatments: Diff=contrast(Fit, list(treat=treatment[-placebo.pos]),list(treat="Placebo"),

Dear all, The incredibly useful Hmisc package provides a method to generate summary tables that can be typeset in latex. The Alzola and Harrell book "An introduction to S and the Hmisc and Design libraries" provides an example that generates mean and quartiles for continuous variables, and numbers and percentages for count variables: summary() with method = 'reverse'. I