similar to: R-alpha: another ctest question

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:

Hi R users! I have the following problem: how appropriate is my aov model under the violation of anova assumptions? Example: a<-c(1,1,1,1,1,1,1,1,1,1,2,2,2,3,3,3,3,3,3,3) b<-c(101,1010,200,300,400, 202, 121, 234, 55,555,66,76,88,34,239, 30, 40, 50,50,60) z<-data.frame(a, b) fligner.test(z$b, factor(z$a)) aov(z$b~factor(z$a))->ll TukeyHSD(ll) Now from the aov i found that my model

Hello, Is there like the kruskal wallis test in relation to ANOVA (no restrictions on normallity and variance homogenity) something (in R) for MANOVA? thanks -- Dr.Nicolaas Busscher Universit?t GH Kassel Nordbahnhofstrasse: 1a, D-37213 Witzenhausen Phone: 0049-(0)5542-98-1715, Fax: 0049-(0)5542-98-1713

I keep getting the 'memory.size' error message when I run a program I have been writing. It always it cannot allocate a vector of a certain size. I believe the error comes in the code fragement below where I have multiple arrays that could be taking up space. Does anyone know a good way around this? w1 <- outer(xk$xk1, data[,x1], function(y,z) abs(z-y)) w2 <- outer(xk$xk2,

Consider a data frame (x) with 2 variables, x1 and x2, having equal values. It looks like: x1 x2 1 1 2 2 3 3 Now, consider a second data frame (xk): xk1 xk2 0.5 0.5 1.0 0.5 1.5 0.5 2.0 0.5 0.5 1 1.0 1 1.5 1 2.0 1 0.5 1.5 1.0 1.5 1.5 1.5 2.0 1.5 0.5 2 1.0 2 1.5 2 2.0 2 I have written code to calculate some differences between these

Hi All, I've reproduced the example from Prof. Friedrich Leisch's webpage. When I write sweave("Example-1.Snw") OR sweave("Example-1.Rnw"), (yes, I renamed them). I get the following error: Writing to file example-1.tex Processing code chunks ... 1 : echo term verbatim Error: chunk 1 Error in library(ctest) : there is no package called 'ctest' Also while

useR's, I am trying trying to find out if there is a faster way to do a certain computation. I have successfully used FOR loops and the apply function to do this, but it can take some time to fully compute, but I was wondering if anyone may know of a different function or way to do this: > x [1] 1 2 3 4 5 > xk [1] 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 I want to do:

Here is a drop-in replacement for the R incomplete beta function. src/math/pbeta.c It is a slightly modified version of the cephes library one from Netlib. In the few cases I tried it seems to give at least 14 digit agreement with the one in S-PLUS (its hard to get more). I'm not sure what performance is like. I'd like to know if it helps with some of the problems which have been

I am using R 2.0.2 on a WinXP I am trying to get the code of the Kruskal-Wallis test but > kruskal.test function (x, ...) UseMethod("kruskal.test") <environment: namespace:stats> > ls(3) [1] "acf" "acf2AR" "add.scope" .............. [181] "kruskal.test" "ks.test"

Hi all, I am performing hundreds of kruskal wallis tests and trying to figure out how to create a file of the p.values I obtain. This is the code I use for the tests: A2<-kruskal.test(X2~treatment) A3<-kruskal.test(X3~treatment) A4<-kruskal.test(X4~treatment) A5<-kruskal.test(X5~treatment) A6<-kruskal.test(X6~treatment) A7<-kruskal.test(X7~treatment)

Hello, I think I am right in saying that a 2 sample wilcox.test is equal to a 2 sample kruskal.test and a 2 sample t.test is equal to a 2 sample anova. This is also stated in the ?kruskal.test man page: The Wilcoxon rank sum test (wilcox.test) as the special case for two samples; lm together with anova for performing one-way location analysis under normality assumptions; with Student's t

I am contemplating improving my ctest package as follows: * Add exact p,q,r,s functions for the Wilcoxon distribution, and change the test accordingly (make `exact' work). * Make Fisher's test work for tables larger than 2 by 2. * Perhaps add an exact unconditional test for 2 by 2 tables? * Perhaps add something on estimating/testing relative risk and odds? As clearly I'd like to

Recently, I ran a series of Kruskal-Wallace tests [kruskal.test()] using by() to group by site Output is a list: >Herb.KW Herb.df$ID: 10-1 Kruskal-Wallis rank sum test data: Indicator_Rating by Year Kruskal-Wallis chi-squared = 15.24, df = 7, p-value = 0.03302 ----------------------------------------------------------------------------------------------------- Herb.df$ID: 18-1

Hi, I am trying to do a kruskal wallis test on two lists, fVisited and cSN: fVisited[[1]] [1] 0.17097623 0.30376141 0.17577266 0.14951855 0.03959753 0.08096217 0.05744888 0.02196257.... .... cSN[[1]] [1] 0.08557303 0.36477791 0.19601252 0.12981040 0.05351320 0.10385542 0.03539577 0.03106175.... So if I just want to do a test on just one of the entries this is simple enough:

Hi, ?kruskal.test() ?a<- c(2,4,5,2,7) ?b<- c(2,2,6) ?c<- c(3,7,9,3) ?kruskal.test(list(a,b,c)) # ?#? Kruskal-Wallis rank sum test # #data:? list(a, b, c) #Kruskal-Wallis chi-squared = 2.003, df = 2, p-value = 0.3673 A.K. Hi I need an expression in R to apply a kruskal.test to this data (for example). a ? a ? a ? a ? ?a ? ?b ? b ? ?b ? ?c ? ?c ? c ? ?c 2 ?4 ? ?5 ? 2 ? ?7 ? ?2 ? 2 ?

Hi, I have got two groups of samples; and for every row, I wish to calculate Kruskal-Wallis' p-value. In the example below, and the stars (****) show where I am struggling to design and put things together. Any help would be appreciated. myFile <- data.frame(Sample_1a = 1:10, Sample_1b = 2:11, Sample_1c = 3:12, Sample_2a=4:13, Sample_2b=7:16, row.names=LETTERS[1:10]) groups <- rep

Hi all My data table (g) contains a continues data column (plant.height) and other columns (columns 8 to 57), each with number of levels of different factors. ANOVA test was done and the p-values were extracted as follos: a <- function(x) anova(lm(plant.height ~ x))$"Pr(>F)"[1] r<- apply(g[,8:57],2,a) If I try to do a Kruskal-Wallis test : kw <- function(x)

typically, the apply family wants you to use vectors to run functions on. However, I have a function, kruskal.test, that requires 2 arguments. kruskal.test(Herb.df$Score,Herb.df$Year) This easily computes the KW ANOVA statistic for any difference across years.... However, my data has multiple sites on which KW needs to be run... here's the data: Herb.df<-

Hi there, I'm rying to 'turn' an Schunk in an .Rnw file(Xemacs-21.4.13, ESS-5.2.8, R-2.1, miktex-2.4.1705). Has anyone got the isorot package to work with Sweave? JC example test.Rnw: \documentclass[a4paper]{article} \usepackage{Sweave} \usepackage{isorot} \rotdriver{dvips} \clockwise \title{Sweave Example 1} \author{apologies to Friedrich Leisch } \begin{document} \maketitle

I submit a couple options for addressing bug 16719: kruskal.test documentation for formula. https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=16719 disallow-character.diff changes the documentation and error message to indicate that factors are accepted. allow-character.diff changes the kruskal.test functions to convert character vectors to factors; documentation is updated accordingly. I