similar to: Non-parametric test for repeated measures and post-hoc single comparisons in R?

Dear List, after updating the exactRanksumTests package I receive a warning that the package is not developed any further and that one should consider the coin package. I don't find the signed rank test in the coin package, only the Wilcoxon Mann Whitney U-Test. I only found a signed rank test in the stats package (wilcox.test) which is able to calculate the exact pvalues but unfortunately

Dear R list members, I am seeking for an elegant way for matching NA elements of a matrix. Everthing I tried, the result was a vector. Here ist an example with vectors that works: a and b are vectors of same lenght representing pairs of data (a[i] b[i]) with some NA elements: ## R Code a <- (1:5) a[2] <- NA b <- (6:10) b[3] <- NA a 1 NA 3 4 5 b 6 7 NA 9 10 ## /R

Dear Everyone, I am using R 2.6.2 on my Redhat AS4. I installed the RPM offered on the cran website. I would like to use R from my bash console with the following command: /path/to/R --vanilla --slave --args < /data/myscript.R this script invokes a range of functions and tools, among others to process a bunch of raster images in a loop using GDAL. The thing is that releasing this command in

dd if=/dev/zero of=sl1 bs=512 count=256000 dd if=/dev/zero of=sl2 bs=512 count=256000 dd if=/dev/zero of=sl3 bs=512 count=256000 dd if=/dev/zero of=sl4 bs=512 count=256000 zpool create -m /export/test1 test1 raidz /export/sl1 /export/sl2 /export/sl3 zpool add -f test1 /export/sl4 dd if=/dev/zero of=sl4 bs=512 count=256000 zpool scrub test1 panic. and message like on image. This message posted

Dear list, has anyone written a package/function in R for computing a point- biserial resp. biserial correlation? Thanks in advance Bernd

IIRC, uncorrectable bitrot even in a nonessential file detected by ZFS used to cause a kernel panic. Bug ID 4924238 was closed with the claim that bitrot-induced panics is not a bug, but the description did mention an open bug ID 4879357, which suggests that it''s considered a bug after all. Can somebody clarify the intended behavior? For example, if I''m running Solaris in a VM,

Hallo! I have data of the following design: NSubj were measured at Baseline (visit 1) and at 3 following time points (visit 2, visit 3, visit 4). There is or is not a treatment. Most interesting is the question if there is a difference in treatment between the results of visit 4 and baseline. (The other time points are also of interest.) The level of significance is alpha=0.0179 (because of an

Greetings, I have two questions that I could not answer from the documentation. A - ecdf and confidence intervals : Is there a (simple) way to generate confidence intervals (95%) for a ecdf? B - cross-validation of rpart trees : a colleague is using S to generate decision tree and mentioned to me the use of cross-validation. Is this function enabled in R ? if so, how should one proceed to

I need to use a red-black tree, which package provides that data structure? -- Best regards Rune Schjellerup Philosof Ph.d-stipendiat, Forskningsenheden for Biostatistik Telefon: 6550 3607 E-mail: rphilosof at health.sdu.dk Adresse: J.B. Winsl?wsvej 9, 5000 Odense C SYDDANSK UNIVERSITET _______________________________________________________________ * Campusvej 55 * 5230 * Odense M * 6550

How come the k is 3 in all of this output? I expected it to be equal to r. > tmp3 <- lapply(1:3, function(k) function(r) print(paste(r, "<- r | k ->", k))) > for (i in 1:3) { tmp3[[i]](i) } [1] "1 <- r | k -> 3" [1] "2 <- r | k -> 3" [1] "3 <- r | k -> 3" -- Med venlig hilsen Rune Schjellerup Philosof Ph.d-stipendiat,

I often find myself making lists similar to this list(var1=var1, var2=var2) It doesn't seem list has an option, to make it use the name of the variable as name in the list. Is there another function that does this? -- Med venlig hilsen Rune Schjellerup Philosof Ph.d-stipendiat, Forskningsenheden for Biostatistik Telefon: 6550 3607 E-mail: rphilosof at health.sdu.dk Adresse: J.B.

Hi dear R-users I try to reproduce the steps included in a LDA. Concerning the eigenvectors there is a difference to SPSS. In my textbook (Bortz) it says, that the matrix with the eigenvectors V usually are not normalized to the length of 1, but in the way that the following holds (SPSS does the same thing): t(Vstar)%*%Derror%*%Vstar = I where Vstar are the normalized eigenvectors. Derror

Hi all, I have been trying to reproduce an analysis from Douglas Montgomery?s book on design and analysis of experiments. Table 6.10 of example 6.2 on page 246, gives a table as follows: > NPK <- expand.grid(A=mp,B=mp,C=mp,D=mp) > Rate <- c(45,71,48,65,68,60,80,65,43,100,45,104,75,86,70,96) > filtration <- cbind(NPK,Rate) > filtration A B C D Rate 1 - - - - 45 2

I am trying to use the friedman.test() on a data frame, d, but I am receiving the following error message: > d AW HS IAC WA 1 6 8 3 5 2 2 2 3 6 3 7 7 8 3 4 8 5 4 5 .... 20 2 5 2 7 21 7 7 6 7 22 7 8 6 8 23 6 8 4 5 24 5 7 5 2 > friedman.test(d) Error in any(is.na(groups)) : argument "groups" is missing, with no default I think I

Hi, I would need to extend the Friedman test to a replicated design. Currently the function: friedman.test(y, ...) only works for unreplicated designs. I found in Conover 1999 "Practical Nonparamteric statistics" an extension of the formula to my case. Nevertheless, other sources, like Sheskin 2000 "Parametric and Nonparametric statistical Procedures" and Daniel 1990

Hi, I was trying to draw an empirical distribution function with uniform confidence bands. So I tried to find a way to calculate values of the Kolmogorov-Smirnov Distribution but failed. I guess it must be hidden somewhere (since the ks-test is implemented), but I was unable to find it. Is there any way to do this? Thanks Leif Boysen

Dear List, I am looking for a function `Pseudo standard error' (PSE), which is available in S+ DOX (design of experiemnt) module - Is there a similar function available in R? Reference for PSE function is in the paper: 'Quick and easy analysis of unreplicated factorials' by Russell V. Lenth, Technometrics, 1989, 31, 4, 469-473. Thanks. -Nitin

Hi, I'm currently working on porting some SAS scripts to R, and hence need to do the same calculation (and get the same results) as SAS in order to make the transition easier for users of the script. In the script, I'm dealing with a two-factorial repeated-measures anova. I'll try to give you a short overview of the setup: - two between-cell factors: facBetweenROI (numbering

Hallo, is there any extension of the pairwise Wilcoxon test to a dependent samples layout with replicates (or, in other terms, a one-way layout with blocking and replicates)? The Wilcoxon method with matched pairs works for the case of dependent samples with one observation per block, while the Mann-Whitney test works for independent samples, thus one single block and replicated observations. Is

Hola a todos, Acudo a vosotros después de darle muchas vueltas al tema en Internet y no localizar la solución para algo que creo trivial. Necesito aplicar el test de Quade a un conjunto de datos en el que existe más de una covariable. El ejemplo de dicho test en R es el siguiente: ## Conover (1999, p. 375f): ## Numbers of five brands of a new hand lotion sold in seven stores ## during one week. y