similar to: Wilcoxon paired signed rank test and continuity correction

Hi Running R v2.0 on SuSe linux 8.2. I'm trying to build a package (which built perfectly on Windows...) on Linux, and I ran: R CMD check mypackage I got: * checking mypackage-maual.tex ... ERROR LaTeX errors when creating DVI version This typically indicates Rd problems OK, there are no problems with my Rd - I got 3 warnings but they were all expected, and the whole package builds find

Hi First a simple question to do with plot(). How do I change the x-axis labels on a plot? For example, I am plotting each row of a matrix, and I want each row to be a line on my graph. Simple really. Eg: plot(sg[1,], type="l") When I do this, the x-axis is labelled 1:38, as I have 38 columns in my matrix. When I do: plot(sg[1,order(sg[1,])], type="l") Ideally I would

Replacing the NA?s with eg. 1 : > a 1 2 3 4 1 20 50 10 80 2 NA 19 NA 49 3 NA 32 NA 61 4 45 101 44 190 a[try(is.na(a)) == TRUE] <- 1 > a 1 2 3 4 1 20 50 10 80 2 1 19 1 49 3 1 32 1 61 4 45 101 44 190 I this can help you, Matthias > -----Urspr?ngliche Nachricht----- > Von: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at

Hi I'm seeing some "odd" behaviour with cbind(). My code is: > cat <- read.table("cogs_category.txt", sep="\t", header=TRUE, quote=NULL, colClasses="character") > colnames(cat) [1] "Code" "Description" > is.factor(cat$Code) [1] FALSE > is.factor(cat$Description) [1] FALSE > is.factor(rainbow(nrow(cat))) [1]

OK, dumb question, and it is probably in the docs somewhere, but after 12 months working with R and quite a while looking at the docs, I still don't know (or have forgotten) how to replace all NA values in a matrix at once with some other value. I can do it column by column using is.na(), but I can't figure out how to do it for the whole matrix. My apologies, I am ashamed ;-) Michael

Hi I have recently installed and implemented Xvfb (X virtual frame buffer) so that I can create jpegs using R over CGI (SUSE Linux 8.1 and Apache 1.3). I have noticed that in order to do this, a file (Rplots.ps) is created in my cgi-bin directory everytime a cgi script is run. This could cause problems though as I have a multi-user system where it is possible that two different users will run

Hi You may remember a while ago I was having problems with creating jpegs using R over CGI. I solved this on my LAPTOP (*grins*) by - executing "xhost +localhost" from the command line - putting "SetEnv DISPLAY=:0.0" in my httpd.conf (apache) And all worked well. Now I'm trying to move the whole system onto my shiny new server - and encountering yet more problems

Hi When I use plot(hclust(dist..)...)...) etc to create a dendrogram of a hierarchial cluster analysis, I end up with a vertical tree. What do I need to do to get a horizontal tree? Also, my users are used to seeing trees who's leaves all "end" at the same place (eg. Like in minitab). Is this possible in R? Thanks Mick Michael Watson Head of Informatics Institute for Animal

Hi I have seen a few posts to this list regarding problems accessing the x11() device over cgi - namely, when trying to create a graphic using the jpeg() function, everything is fine from the command line but it won't work over cgi, producing the error: "Unable to open connection to X11 display" Has anyone actually solved this particular problem satisfactorily? Please reply

Hi I am running an R script that creates 100s of graphs, and I need to use the greek CAPITAL letter delta in the mtext() function. I got as far as expression(delta) but this gives me the lowercase version. Can anyone help? I'm using R 2.7 on Windows XP Mick Head of Informatics Institute for Animal Health Compton Berks RG20 7NN 01635 578411

Following on from my previous mail! plot(1:10,1:10) mtext("title", side=3, adj=0, font=3, cex=1.5) This works as expected and puts the font in italics. tag <- "A)" suffix <- "genea::" plot(1:10,1:10) mtext(bquote(.(tag) ~ Delta * .(suffix)), side=3, adj=0, font=3, cex=1.5) Here, the font isn't in italics, it's normal. I presume this is some

Hi all, I have a distribution, and take a sample of it. Then I compare that sample with the mean of the population like here in "Wilcoxon signed rank test with continuity correction": > wilcox.test(Sample,mu=mean(All), alt="two.sided") Wilcoxon signed rank test with continuity correction data: AlphaNoteOnsetDists V = 63855, p-value = 0.0002093 alternative hypothesis:

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, I'm using R software to evaluate Wilcoxon Rank Sum Test and I' getting one Warning message as this: > C1dea_com [1] 1.000 0.345 0.200 0.208 0.508 0.480 0.545 0.563 0.451 0.683 0.380 0.913 1.000 0.506 > C1dea_sem [1] 1.000 0.665 0.284 0.394 0.509 0.721 0.545 0.898 0.744 0.683 0.382 0.913 1.000 0.970 > wilcox.test(C1dea_sem,C1dea_com, paired = TRUE, alternative =

Hi, I am using wilcox.test function to test the difference between the means of two samples. The data points are paired, so I am using a paired test. There is one strange case. Sample A has a higher mean than a sample B. However, wilcox.test function says that sample B has a significantly higher "mean rank" than sample A. How is it possible? Here is the code (data file is attached):

Dear all, I am trying to run Wilcoxon Rank-Sum Test with Bonferroni's correction. I have two lists: l0, l1: mapply(function(x,y)wilcox.test(x,y)$p.value, l0, l1) How do I run Bonferroni's correction on mapply? Any help is much apperciated. Thanks, -Raj

Hi there, I'm a total newbie to R. I'd like to use a Wilcoxon Rank Sum test to compare two populations of values. Further, I'd like to do this simultaneously for 114 sets of values. The two populations are C and N. The different sets of values have arbitrary names (I'll call them a, b, c etc). The set-up is as follows: a b c d .... C 2 C 3 C 5 C 9 C 4

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 have a question concerning the Wilcoxon signed-rank test, and specifically, which R subroutine I should use for my particular dataset. There are three different commands in R (that I'm aware of) that calculate the Wilcoxon signed-rank test; wilcox.test, wilcox.exact, and wilcoxsign_test. When I run the three commands on the same dataset, I get different p-values. I'm hoping that

Hi guys,   I have two data sets of prices: endprice0, endprice1   I use the Wilcox test:   wilcox.test(endprice0, endprice1, paired = TRUE, alternative = "two.sided",  conf.int = T, conf.level = 0.9)   The result is with V = 1819, p-value = 0.8812.   Then I calculated the z-value of the test: z-value = -2.661263. The corresponding p-value is: p-value = 0.003892, which is different from