similar to: Biobass, SAGx, and Jonckheere-Terpstra test

Is it possible to implement the Jonckheere-Terpstra test for ordered alternatives using the coin package: Conditional Inference Procedures in a Permutation Test Framework? I found jonckheere.test{clinfun}, but it uses a normal approximation when ties are present in the data. To make this concrete, I've include a small dataset. Thanks. --Dale Hollander and Wolfe, 1999 Table 6.6, pg 205

Hello, can anybody here explain what a Jonckheere-Terpstra test is and whether it is implemented in R? I just know it's a non-parametric test, otherwise I've no clue about it ;-( . Are there alternatives to this test? thanks for help, Arne

Dear All: I am a newbie to R and as such i am posting this request for help. I am trying to use R to compute the "Calinski Harabasz (CH) Index ". The CH Index is available in the library SAGx. The version of R i am running is 2.2.1. I have my data in a CSV format which i read into R using the read.table() command. After the data has been read i am loading in the "SAGx"

Full_Name: Alexandra R Mendes de Almeida Version: 2.4.1 OS: Linux Submission from: (NULL) (201.51.225.154) R-Masters The packahe SGAx is not on the repository. See the message returned by R when i tried to install it: > install.packages("SAGx") --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done Warning in download.packages(unique(pkgs),

Hi! I have large matrices, one column per variable and I have a vector of factors / grouping symbols. Then I am computing subtotals for the groups but it takes pretty long and thus I wanted to ask if there is a better way to do it or if this is already the best way: subTotals <- function(x, groupvec) do.call("rbind",lapply(split(x,groupvec),colSums,na.rm=T)) Thanks reading my

Has anyone implemented Tibrishani's gap statistic in R or S plus? If so I would greatly appreciate the relevant script file. Any help will be much appreciated Kind regards Dr Graham Leask Economics and Strategy Group Aston Business School Aston University Aston Triangle Birmingham B4 7ET Tel: Direct line 0121 204 3150 email g.leask@aston.ac.uk [[alternative HTML version deleted]]

I wrote an R wrapper function (phylpro) around a C function (Rphylpro). The first time I'm running my function, it runs with no errors. The second time I'm trying to run it, I get an error message with the first argument to .Call garbled. Set up: > dyn.load("Phylpro.so") > source("phylpro.R") > WinHalfWidth<-30 > permReps<-10 > breaks<-c(548,

Hi, I am doing kernel density plots, and am trying to make the lines thicker. I comparing three groups, in sm.density.compare. I tried changing lwd to make the line sthicker right on the density compare call, but was not able to do it. There is not an option in sm.options to specify line thickness, as well as cex.ylab or cex.xlab- I tried it and it does not change the thickness of the lines.

Hi, I am using this subset statement and it works outside a function. LIS[[i]]<- lapply(LI, subset, select=cov[[i]]) However, wrapped inside a function this statement produces the same values for every LIS[[1]] which is only the first subset of LI. Does anyone know why is not working correctly inside a function? ff = factor(covariate) nLev <- nlevels(ff) cov <-

Hello R useRs, I have a function which returns a list of functions : freq1 <- function(x) { lev <- unique(x[!is.na(x)]) nlev <- length(lev) args <- alist(x=) if (nlev == 1) { body <- c("{", "sum(!is.na(x))", "}") f <- function() {} formals(f) <- as.pairlist(args) body(f) <- parse(text = body) namef <-

I apologize if this is not appropriate for this mailing list. In R, there is already functionality to drop unused factor levels. However, I am proposing the code below that I wrote. In some occasions, it was faster than applying function 'factor'. In any case, there is no restriction for anyone to use the code below. droplevels2 <- function(x) { if (is.null(levels(x))) stop("no

Michael, Got any solutions to this. We are confronted with the same problem only SUSE 9.x ,Samba , W200 workstations. Only few of them have this problem. (exactly as how you describe it). Mark

Hi all, Small problem in generating the line charts. Question: Legend for the first graph is coming wrong., for second graph correctly. Please fix the legend postion at the down of graph. Plesae give me the solution. Thank you Devarayalu Orange1 <- structure(list(REFID = c(7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9), ARM = c(1, 1, 1, 1, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 2,

Full_Name: Robert Kushler Version: 2.4.1 OS: Windows XP Submission from: (NULL) (69.245.71.40) In the constant leverage case, plot #5 is not correctly produced. The labels on the x-axis are sorted correctly by magnitude of the fitted value, but the data are plotted in the original factor order. I changed facval[ord] <- facval xx <- facval

I want to overlay 50 denisty plots on a single plot. For each plot there are 10,000 data points and i want the empirical density of the data. I have not been able to find an easy way to achieve this (I have scoured the manula and website so sorry if I missed it!), does anyone have any suggestions? Thankyou.

Hello again, First, thanks for the help that got the latest plotrix package finished. I had been planning to write something about packages since Scott Waichler offered the gantt.chart function. Then Ben Bolker (who helped me to write the axis.break function) asked if I would be willing to include some of his plotting functions and almost immediately after that Sander Oom kindly donated the

Hello again, First, thanks for the help that got the latest plotrix package finished. I had been planning to write something about packages since Scott Waichler offered the gantt.chart function. Then Ben Bolker (who helped me to write the axis.break function) asked if I would be willing to include some of his plotting functions and almost immediately after that Sander Oom kindly donated the

Dear list, I'm having second thoughts after solving a very trivial problem: I want to extend the relevel() function to reorder an arbitrary number of levels of a factor in one go. I could not find a trivial way of using the code obtained by getS3method("relevel","factor"). Instead, I thought of solving the problem in a recursive manner (possibly after reading

I use the following function to select a subgroup from a data vector (usually to calculate descriptive stats). select.range<-function (groupvec, min, max, data) { if (nargs() > 3) { min.cond <- groupvec >= min max.cond <- groupvec < max cond <- min.cond & max.cond selected <- na.remove(ifelse(cond, data, NA))

Hi folks, I am tearing my hair out on this one. I am using an example from Pinheiro and Bates. ### this works data(Orange) mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange ) ### This works mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange, fixed = Asymp + xmid + scal ~ 1, start =